(* Content-type: application/vnd.wolfram.mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 9.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 157, 7] NotebookDataLength[ 31982, 864] NotebookOptionsPosition[ 31182, 835] NotebookOutlinePosition[ 31607, 851] CellTagsIndexPosition[ 31564, 848] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[TextData[StyleBox["(*N-Atom and homologues 1s(2)2s(2) 2p(3) concentric \ spherical,with exchange correction page 40; 6.12.2011 *)", FontFamily->"Courier New"]], "Input", CellChangeTimes->{{3.5662329377492924`*^9, 3.566232983706973*^9}, { 3.566233035452264*^9, 3.566233065435517*^9}, {3.5662332314354086`*^9, 3.566233239828223*^9}, {3.566233339808799*^9, 3.5662333604944353`*^9}, 3.5662334004149055`*^9, 3.566233505231489*^9}, FormatType->"TextForm"], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"Clear", "[", RowBox[{"z", ",", "sig"}], "]"}], ";", " ", RowBox[{"sig", "=", "0.205"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"T", "=", RowBox[{ RowBox[{"9", "/", RowBox[{"(", RowBox[{"4", "*", RowBox[{"P", "^", "2"}]}], ")"}]}], "+", RowBox[{"9", "*", RowBox[{"4", "/", RowBox[{"(", RowBox[{"4", "*", RowBox[{"Q", "^", "2"}]}], ")"}]}]}], "+", RowBox[{"6", "/", RowBox[{"R", "^", "2"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vne", "=", RowBox[{ RowBox[{ RowBox[{"-", "3"}], "*", RowBox[{"z", "/", "P"}]}], "-", RowBox[{"3", "*", RowBox[{"z", "/", "Q"}]}], "-", RowBox[{"3", "*", RowBox[{"z", "/", "R"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vee", "=", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"2", "/", "Q"}], ")"}], "*", RowBox[{"(", RowBox[{"3", "-", RowBox[{ RowBox[{"(", RowBox[{"6", "/", "10"}], ")"}], "*", RowBox[{ RowBox[{"(", RowBox[{"P", "/", "Q"}], ")"}], "^", "2"}]}]}], ")"}]}], "-", RowBox[{"6", "*", "sig", "*", RowBox[{ RowBox[{"P", "^", "2"}], "/", RowBox[{"Q", "^", "3"}]}]}], " ", "+", " ", RowBox[{"3", "*", RowBox[{"sig", "/", "P"}]}], "+", RowBox[{"3", "*", RowBox[{"sig", "/", "Q"}]}], "+", RowBox[{"6", "*", RowBox[{ RowBox[{"(", RowBox[{"1", "-", RowBox[{ RowBox[{"(", RowBox[{"3", "/", "20"}], ")"}], "*", RowBox[{ RowBox[{"(", RowBox[{"P", "/", "R"}], ")"}], "^", "2"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"1", "/", "32"}], ")"}], "*", RowBox[{ RowBox[{"(", RowBox[{"P", "/", "R"}], ")"}], "^", "3"}]}]}], ")"}], "/", "R"}]}], "+", RowBox[{"6", "*", RowBox[{ RowBox[{"(", RowBox[{"1", "-", RowBox[{ RowBox[{"(", RowBox[{"3", "/", "20"}], ")"}], "*", RowBox[{ RowBox[{"(", RowBox[{"Q", "/", "R"}], ")"}], "^", "2"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"1", "/", "32"}], ")"}], "*", RowBox[{ RowBox[{"(", RowBox[{"Q", "/", "R"}], ")"}], "^", "3"}]}]}], ")"}], "/", "R"}]}], "+", RowBox[{"1.645584", "/", "R"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"func", "=", RowBox[{"T", "+", "Vne", "+", "Vee"}]}], ";", " ", RowBox[{"Vir", " ", "=", " ", RowBox[{ RowBox[{"(", RowBox[{"Vne", "+", "Vee"}], ")"}], "/", "T"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{"Etot", ",", " ", "P", ",", " ", "Q"}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{"t", " ", "=", " ", RowBox[{"Table", "[", RowBox[{ RowBox[{"FindMinimum", "[", RowBox[{"func", ",", RowBox[{"{", RowBox[{"P", ",", "0.22"}], "}"}], ",", RowBox[{"{", RowBox[{"Q", ",", "1.41"}], "}"}], ",", RowBox[{"{", RowBox[{"R", ",", "1.20"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"MaxIterations", "\[Rule]", "800"}], ",", " ", RowBox[{"Method", " ", "\[Rule]", "\"\\""}]}], "}"}]}], "]"}], ",", RowBox[{"{", RowBox[{"z", ",", "7", ",", "18"}], "}"}]}], "]"}]}]}], "Input", CellChangeTimes->{{3.5662329377492924`*^9, 3.566232983706973*^9}, { 3.566233035452264*^9, 3.566233065435517*^9}, {3.5662332314354086`*^9, 3.566233239828223*^9}, {3.566233339808799*^9, 3.5662333604944353`*^9}, 3.5662334004149055`*^9, 3.5662334336897635`*^9, {3.5662334712234297`*^9, 3.5662334720658307`*^9}, {3.566233741400304*^9, 3.5662337423831053`*^9}, { 3.566233776235165*^9, 3.5662339244198256`*^9}, {3.5662339597538877`*^9, 3.5662339608302894`*^9}, {3.566234379331424*^9, 3.5662343864918365`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "54.658142264226285`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.2210649259732874`"}], ",", RowBox[{"Q", "\[Rule]", "1.4139278082118352`"}], ",", RowBox[{"R", "\[Rule]", "1.2074327021227957`"}]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "74.4785045269774`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.19273932823620674`"}], ",", RowBox[{"Q", "\[Rule]", "1.2039252853802211`"}], ",", RowBox[{"R", "\[Rule]", "0.8826552542228082`"}]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "97.56313787960495`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.1708486158800319`"}], ",", RowBox[{"Q", "\[Rule]", "1.016528139283902`"}], ",", RowBox[{"R", "\[Rule]", "0.713974515356104`"}]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "123.9023982086206`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.15341956034790566`"}], ",", RowBox[{"Q", "\[Rule]", "0.8729411962866754`"}], ",", RowBox[{"R", "\[Rule]", "0.6030246576350509`"}]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "153.4943850491857`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.13921389588243155`"}], ",", RowBox[{"Q", "\[Rule]", "0.7633157280569021`"}], ",", RowBox[{"R", "\[Rule]", "0.5228069425569312`"}]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "186.3383072169185`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.1274132944241037`"}], ",", RowBox[{"Q", "\[Rule]", "0.6776520079939099`"}], ",", RowBox[{"R", "\[Rule]", "0.4617342683443797`"}]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "222.4336966904474`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.11745490408870395`"}], ",", RowBox[{"Q", "\[Rule]", "0.609084616175056`"}], ",", RowBox[{"R", "\[Rule]", "0.413572499219011`"}]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "261.780236965286`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.10893879524989715`"}], ",", RowBox[{"Q", "\[Rule]", "0.553035628118999`"}], ",", RowBox[{"R", "\[Rule]", "0.3745771973920501`"}]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "304.37770022464434`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.10157297068870219`"}], ",", RowBox[{"Q", "\[Rule]", "0.5063941876492405`"}], ",", RowBox[{"R", "\[Rule]", "0.34234034141791986`"}]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "350.2259164035074`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.09513925227682629`"}], ",", RowBox[{"Q", "\[Rule]", "0.4669888918398166`"}], ",", RowBox[{"R", "\[Rule]", "0.3152358329258193`"}]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "399.3247553149391`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.08947134065170427`"}], ",", RowBox[{"Q", "\[Rule]", "0.4332636446740114`"}], ",", RowBox[{"R", "\[Rule]", "0.2921233619462681`"}]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "451.6741153560028`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.08444026225663692`"}], ",", RowBox[{"Q", "\[Rule]", "0.404076281771941`"}], ",", RowBox[{"R", "\[Rule]", "0.2721785497620972`"}]}], "}"}]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{ 3.566233444313382*^9, 3.566233476886239*^9, 3.566233748186316*^9, { 3.566233782381576*^9, 3.566233929786235*^9}, 3.5662339660874987`*^9, 3.5662341994795084`*^9, 3.566234318241717*^9, 3.5662343908130445`*^9, 3.5662348042405715`*^9, 3.566234994170905*^9, 3.5662354527769113`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{"Virial", " ", "ratio"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"N", "[", RowBox[{ RowBox[{"Vir", " ", "/.", " ", RowBox[{"t", "[", RowBox[{"[", RowBox[{ RowBox[{"{", RowBox[{"z", "-", "6"}], "}"}], ",", "2"}], "]"}], "]"}]}], ",", "10"}], "]"}], ",", RowBox[{"{", RowBox[{"z", " ", ",", "7", ",", "18"}], "}"}]}], "]"}], "\[IndentingNewLine]", RowBox[{"Table", "[", RowBox[{ RowBox[{"N", "[", RowBox[{ RowBox[{"Vne", " ", "/.", " ", RowBox[{"t", "[", RowBox[{"[", RowBox[{ RowBox[{"{", RowBox[{"z", "-", "6"}], "}"}], ",", "2"}], "]"}], "]"}]}], ",", "10"}], "]"}], ",", RowBox[{"{", RowBox[{"z", " ", ",", "7", ",", "18"}], "}"}]}], "]"}], "\[IndentingNewLine]", RowBox[{"Table", "[", RowBox[{ RowBox[{"N", "[", RowBox[{ RowBox[{"Vee", " ", "/.", " ", RowBox[{"t", "[", RowBox[{"[", RowBox[{ RowBox[{"{", RowBox[{"z", "-", "6"}], "}"}], ",", "2"}], "]"}], "]"}]}], ",", "10"}], "]"}], ",", RowBox[{"{", RowBox[{"z", " ", ",", "7", ",", "18"}], "}"}]}], "]"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{"Etot", " ", "in", " ", "eV"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"Table", "[", RowBox[{ RowBox[{"N", "[", RowBox[{ RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{ RowBox[{"{", RowBox[{"z", "-", "6"}], "}"}], ",", "1"}], "]"}], "]"}], "*", "27.21138"}], ",", "10"}], "]"}], ",", RowBox[{"{", RowBox[{"z", ",", "7", ",", "18"}], "}"}]}], "]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"PA", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"P", " ", "/.", " ", RowBox[{"t", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "12"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"QA", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"Q", " ", "/.", " ", RowBox[{"t", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "12"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"RA", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"R", " ", "/.", " ", RowBox[{"t", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "12"}], "}"}]}], "]"}]}], ";"}]}]}]], "Input", CellChangeTimes->{{3.5662329377492924`*^9, 3.566232983706973*^9}, { 3.566233035452264*^9, 3.566233065435517*^9}, {3.5662332314354086`*^9, 3.566233239828223*^9}, {3.566233339808799*^9, 3.5662333604944353`*^9}, { 3.5662334004149055`*^9, 3.5662334369657693`*^9}, {3.566235340035513*^9, 3.5662354402344894`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"-", "1.999999999999996`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "1.9999999940576103`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "1.9999999999999953`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "2.0000000000000004`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "1.9999999964308768`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "1.9999999999997582`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "1.999999997954952`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "1.9999999999999951`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "1.9999999999814375`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "1.9999999999999987`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "1.999999999999964`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "1.9999999920174028`"}], "}"}]}], "}"}]], "Output", CellChangeTimes->{ 3.5662334443289824`*^9, 3.5662334769018393`*^9, 3.5662337482019157`*^9, { 3.566233782397176*^9, 3.566233929801835*^9}, 3.566233966103099*^9, 3.5662341994951086`*^9, 3.566234318241717*^9, 3.5662343908286448`*^9, 3.5662348042405715`*^9, 3.5662349942021055`*^9, 3.566235452792511*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"-", "127.23923373695581`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "171.64598862442463`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "222.41211487825376`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "279.65799622782794`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "343.39855462695016`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "413.6366188792659`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "490.3731160984445`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "573.6084784888344`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "663.3429642307925`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "759.576744968487`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "862.3099456667937`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "971.5426655858216`"}], "}"}]}], "}"}]], "Output", CellChangeTimes->{ 3.5662334443289824`*^9, 3.5662334769018393`*^9, 3.5662337482019157`*^9, { 3.566233782397176*^9, 3.566233929801835*^9}, 3.566233966103099*^9, 3.5662341994951086`*^9, 3.566234318241717*^9, 3.5662343908286448`*^9, 3.5662348042405715`*^9, 3.5662349942021055`*^9, 3.566235452792511*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", "17.922949208503045`", "}"}], ",", RowBox[{"{", "22.68897912788953`", "}"}], ",", RowBox[{"{", "27.285839119043427`", "}"}], ",", RowBox[{"{", "31.85319981058681`", "}"}], ",", RowBox[{"{", "36.40978398073842`", "}"}], ",", RowBox[{"{", "40.960004445383774`", "}"}], ",", RowBox[{"{", "45.505722262662125`", "}"}], ",", RowBox[{"{", "50.04800455826108`", "}"}], ",", RowBox[{"{", "54.58756377585393`", "}"}], ",", RowBox[{"{", "59.124912161471784`", "}"}], ",", RowBox[{"{", "63.660435036901106`", "}"}], ",", RowBox[{"{", "68.19443126828355`", "}"}]}], "}"}]], "Output", CellChangeTimes->{ 3.5662334443289824`*^9, 3.5662334769018393`*^9, 3.5662337482019157`*^9, { 3.566233782397176*^9, 3.566233929801835*^9}, 3.566233966103099*^9, 3.5662341994951086`*^9, 3.566234318241717*^9, 3.5662343908286448`*^9, 3.5662348042405715`*^9, 3.5662349942021055`*^9, 3.566235452808111*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"-", "1487.3234792459218`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "2026.6628885153023`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "2654.8276188343243`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "3371.555240566094`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "4176.79403943971`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "5070.522486236312`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "6052.727845448506`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "7123.401504552444`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "8282.537264338882`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "9530.130497104074`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "10866.177660281828`"}], "}"}], ",", RowBox[{"{", RowBox[{"-", "12290.675989116027`"}], "}"}]}], "}"}]], "Output", CellChangeTimes->{ 3.5662334443289824`*^9, 3.5662334769018393`*^9, 3.5662337482019157`*^9, { 3.566233782397176*^9, 3.566233929801835*^9}, 3.566233966103099*^9, 3.5662341994951086`*^9, 3.566234318241717*^9, 3.5662343908286448`*^9, 3.5662348042405715`*^9, 3.5662349942021055`*^9, 3.566235452808111*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"Table", "[", "\[IndentingNewLine]", RowBox[{ RowBox[{"Show", "[", RowBox[{ RowBox[{"Graphics", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"Thickness", "[", "0.01", "]"}], ",", " ", "\n", " ", RowBox[{"{", RowBox[{ RowBox[{"RGBColor", "[", RowBox[{"1", ",", "0", ",", "0"}], "]"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", " ", RowBox[{"-", RowBox[{"RA", "[", RowBox[{"[", "i", "]"}], "]"}]}]}], "}"}], ",", " ", RowBox[{"RA", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}], ",", " ", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", " ", RowBox[{"RA", "[", RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", " ", RowBox[{"RA", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}], "}"}], ",", " ", "\n", " ", RowBox[{"{", RowBox[{ RowBox[{"RGBColor", "[", RowBox[{"1", ",", "0", ",", "1"}], "]"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", RowBox[{"RA", "[", RowBox[{"[", "i", "]"}], "]"}]}], ",", " ", "0"}], "}"}], ",", " ", RowBox[{"RA", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}], ",", " ", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"RA", "[", RowBox[{"[", "i", "]"}], "]"}], ",", " ", "0"}], "}"}], ",", " ", "\n", " ", RowBox[{"RA", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"RGBColor", "[", RowBox[{"0", ",", "1", ",", "1"}], "]"}], ",", " ", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", " ", "0"}], "}"}], ",", " ", RowBox[{"RA", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}], "}"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Thickness", "[", "0.015", "]"}], ",", RowBox[{"{", RowBox[{ RowBox[{"RGBColor", "[", RowBox[{"0", ",", "0", ",", "0"}], "]"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", " ", "0"}], "}"}], ",", " ", RowBox[{"PA", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}], ",", " ", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", " ", "0"}], "}"}], ",", " ", RowBox[{"QA", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}], "}"}]}], "}"}]}], "}"}], "]"}], ",", " ", "\n", " ", RowBox[{"{", RowBox[{ RowBox[{"AspectRatio", " ", "->", " ", "Automatic"}], ",", " ", "\n", " ", RowBox[{"Axes", " ", "->", " ", "True"}], ",", " ", RowBox[{"GridLines", " ", "->", " ", "Automatic"}], ",", RowBox[{"ColorOutput", "\[Rule]", " ", "Automatic"}], ",", " ", "\n", " ", RowBox[{"PlotRange", " ", "->", " ", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "2.4"}], ",", "2.4"}], "}"}], ",", " ", RowBox[{"{", RowBox[{ RowBox[{"-", "2.4"}], ",", " ", "2.4"}], "}"}]}], "}"}]}], ",", " ", RowBox[{"Frame", " ", "->", " ", "True"}]}], "}"}]}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "12"}], "}"}]}], "]"}], "\[IndentingNewLine]"}]], "Input", CellChangeTimes->{{3.5662329377492924`*^9, 3.566232983706973*^9}, { 3.566233035452264*^9, 3.566233065435517*^9}, {3.5662332314354086`*^9, 3.566233239828223*^9}, {3.566233339808799*^9, 3.5662333604944353`*^9}, { 3.5662334004149055`*^9, 3.566233411490925*^9}, {3.5662341880134883`*^9, 3.5662341936294985`*^9}, {3.5662343131093082`*^9, 3.5662343132965083`*^9}, 3.566234759936494*^9, {3.5662347927745514`*^9, 3.5662347929461517`*^9}, { 3.566234953376834*^9, 3.566234975872073*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ GraphicsBox[{ {Thickness[0.01], {RGBColor[1, 0, 0], CircleBox[{0, -1.2074327021227957`}, 1.2074327021227957`], CircleBox[{0, 1.2074327021227957`}, 1.2074327021227957`]}, {RGBColor[1, 0, 1], CircleBox[{-1.2074327021227957`, 0}, 1.2074327021227957`], CircleBox[{1.2074327021227957`, 0}, 1.2074327021227957`]}, {RGBColor[0, 1, 1], CircleBox[{0, 0}, 1.2074327021227957`]}}, {Thickness[0.015], {RGBColor[0, 0, 0], CircleBox[{0, 0}, 0.2210649259732874], CircleBox[{0, 0}, 1.4139278082118352`]}}}, AspectRatio->Automatic, Axes->True, ColorOutput->Automatic, Frame->True, GridLines->Automatic, PlotRange->{{-2.4, 2.4}, {-2.4, 2.4}}], ",", GraphicsBox[{ {Thickness[0.01], {RGBColor[1, 0, 0], CircleBox[{0, -0.8826552542228082}, 0.8826552542228082], CircleBox[{0, 0.8826552542228082}, 0.8826552542228082]}, {RGBColor[1, 0, 1], CircleBox[{-0.8826552542228082, 0}, 0.8826552542228082], CircleBox[{0.8826552542228082, 0}, 0.8826552542228082]}, {RGBColor[0, 1, 1], CircleBox[{0, 0}, 0.8826552542228082]}}, {Thickness[0.015], {RGBColor[0, 0, 0], CircleBox[{0, 0}, 0.19273932823620674`], CircleBox[{0, 0}, 1.2039252853802211`]}}}, AspectRatio->Automatic, Axes->True, ColorOutput->Automatic, Frame->True, GridLines->Automatic, PlotRange->{{-2.4, 2.4}, {-2.4, 2.4}}], ",", GraphicsBox[{ {Thickness[0.01], {RGBColor[1, 0, 0], CircleBox[{0, -0.713974515356104}, 0.713974515356104], CircleBox[{0, 0.713974515356104}, 0.713974515356104]}, {RGBColor[1, 0, 1], CircleBox[{-0.713974515356104, 0}, 0.713974515356104], CircleBox[{0.713974515356104, 0}, 0.713974515356104]}, {RGBColor[0, 1, 1], CircleBox[{0, 0}, 0.713974515356104]}}, {Thickness[0.015], {RGBColor[0, 0, 0], CircleBox[{0, 0}, 0.1708486158800319], CircleBox[{0, 0}, 1.016528139283902]}}}, AspectRatio->Automatic, Axes->True, ColorOutput->Automatic, Frame->True, GridLines->Automatic, PlotRange->{{-2.4, 2.4}, {-2.4, 2.4}}], ",", GraphicsBox[{ {Thickness[0.01], {RGBColor[1, 0, 0], CircleBox[{0, -0.6030246576350509}, 0.6030246576350509], CircleBox[{0, 0.6030246576350509}, 0.6030246576350509]}, {RGBColor[1, 0, 1], CircleBox[{-0.6030246576350509, 0}, 0.6030246576350509], CircleBox[{0.6030246576350509, 0}, 0.6030246576350509]}, {RGBColor[0, 1, 1], CircleBox[{0, 0}, 0.6030246576350509]}}, {Thickness[0.015], {RGBColor[0, 0, 0], CircleBox[{0, 0}, 0.15341956034790566`], CircleBox[{0, 0}, 0.8729411962866754]}}}, AspectRatio->Automatic, Axes->True, ColorOutput->Automatic, Frame->True, GridLines->Automatic, PlotRange->{{-2.4, 2.4}, {-2.4, 2.4}}], ",", GraphicsBox[{ {Thickness[0.01], {RGBColor[1, 0, 0], CircleBox[{0, -0.5228069425569312}, 0.5228069425569312], CircleBox[{0, 0.5228069425569312}, 0.5228069425569312]}, {RGBColor[1, 0, 1], CircleBox[{-0.5228069425569312, 0}, 0.5228069425569312], CircleBox[{0.5228069425569312, 0}, 0.5228069425569312]}, {RGBColor[0, 1, 1], CircleBox[{0, 0}, 0.5228069425569312]}}, {Thickness[0.015], {RGBColor[0, 0, 0], CircleBox[{0, 0}, 0.13921389588243155`], CircleBox[{0, 0}, 0.7633157280569021]}}}, AspectRatio->Automatic, Axes->True, ColorOutput->Automatic, Frame->True, GridLines->Automatic, PlotRange->{{-2.4, 2.4}, {-2.4, 2.4}}], ",", GraphicsBox[{ {Thickness[0.01], {RGBColor[1, 0, 0], CircleBox[{0, -0.4617342683443797}, 0.4617342683443797], CircleBox[{0, 0.4617342683443797}, 0.4617342683443797]}, {RGBColor[1, 0, 1], CircleBox[{-0.4617342683443797, 0}, 0.4617342683443797], CircleBox[{0.4617342683443797, 0}, 0.4617342683443797]}, {RGBColor[0, 1, 1], CircleBox[{0, 0}, 0.4617342683443797]}}, {Thickness[0.015], {RGBColor[0, 0, 0], CircleBox[{0, 0}, 0.1274132944241037], CircleBox[{0, 0}, 0.6776520079939099]}}}, AspectRatio->Automatic, Axes->True, ColorOutput->Automatic, Frame->True, GridLines->Automatic, PlotRange->{{-2.4, 2.4}, {-2.4, 2.4}}], ",", GraphicsBox[{ {Thickness[0.01], {RGBColor[1, 0, 0], CircleBox[{0, -0.413572499219011}, 0.413572499219011], CircleBox[{0, 0.413572499219011}, 0.413572499219011]}, {RGBColor[1, 0, 1], CircleBox[{-0.413572499219011, 0}, 0.413572499219011], CircleBox[{0.413572499219011, 0}, 0.413572499219011]}, {RGBColor[0, 1, 1], CircleBox[{0, 0}, 0.413572499219011]}}, {Thickness[0.015], {RGBColor[0, 0, 0], CircleBox[{0, 0}, 0.11745490408870395`], CircleBox[{0, 0}, 0.609084616175056]}}}, AspectRatio->Automatic, Axes->True, ColorOutput->Automatic, Frame->True, GridLines->Automatic, PlotRange->{{-2.4, 2.4}, {-2.4, 2.4}}], ",", GraphicsBox[{ {Thickness[0.01], {RGBColor[1, 0, 0], CircleBox[{0, -0.3745771973920501}, 0.3745771973920501], CircleBox[{0, 0.3745771973920501}, 0.3745771973920501]}, {RGBColor[1, 0, 1], CircleBox[{-0.3745771973920501, 0}, 0.3745771973920501], CircleBox[{0.3745771973920501, 0}, 0.3745771973920501]}, {RGBColor[0, 1, 1], CircleBox[{0, 0}, 0.3745771973920501]}}, {Thickness[0.015], {RGBColor[0, 0, 0], CircleBox[{0, 0}, 0.10893879524989715`], CircleBox[{0, 0}, 0.553035628118999]}}}, AspectRatio->Automatic, Axes->True, ColorOutput->Automatic, Frame->True, GridLines->Automatic, PlotRange->{{-2.4, 2.4}, {-2.4, 2.4}}], ",", GraphicsBox[{ {Thickness[0.01], {RGBColor[1, 0, 0], CircleBox[{0, -0.34234034141791986`}, 0.34234034141791986`], CircleBox[{0, 0.34234034141791986`}, 0.34234034141791986`]}, {RGBColor[1, 0, 1], CircleBox[{-0.34234034141791986`, 0}, 0.34234034141791986`], CircleBox[{0.34234034141791986`, 0}, 0.34234034141791986`]}, {RGBColor[0, 1, 1], CircleBox[{0, 0}, 0.34234034141791986`]}}, {Thickness[0.015], {RGBColor[0, 0, 0], CircleBox[{0, 0}, 0.10157297068870219`], CircleBox[{0, 0}, 0.5063941876492405]}}}, AspectRatio->Automatic, Axes->True, ColorOutput->Automatic, Frame->True, GridLines->Automatic, PlotRange->{{-2.4, 2.4}, {-2.4, 2.4}}], ",", GraphicsBox[{ {Thickness[0.01], {RGBColor[1, 0, 0], CircleBox[{0, -0.3152358329258193}, 0.3152358329258193], CircleBox[{0, 0.3152358329258193}, 0.3152358329258193]}, {RGBColor[1, 0, 1], CircleBox[{-0.3152358329258193, 0}, 0.3152358329258193], CircleBox[{0.3152358329258193, 0}, 0.3152358329258193]}, {RGBColor[0, 1, 1], CircleBox[{0, 0}, 0.3152358329258193]}}, {Thickness[0.015], {RGBColor[0, 0, 0], CircleBox[{0, 0}, 0.09513925227682629], CircleBox[{0, 0}, 0.4669888918398166]}}}, AspectRatio->Automatic, Axes->True, ColorOutput->Automatic, Frame->True, GridLines->Automatic, PlotRange->{{-2.4, 2.4}, {-2.4, 2.4}}], ",", GraphicsBox[{ {Thickness[0.01], {RGBColor[1, 0, 0], CircleBox[{0, -0.2921233619462681}, 0.2921233619462681], CircleBox[{0, 0.2921233619462681}, 0.2921233619462681]}, {RGBColor[1, 0, 1], CircleBox[{-0.2921233619462681, 0}, 0.2921233619462681], CircleBox[{0.2921233619462681, 0}, 0.2921233619462681]}, {RGBColor[0, 1, 1], CircleBox[{0, 0}, 0.2921233619462681]}}, {Thickness[0.015], {RGBColor[0, 0, 0], CircleBox[{0, 0}, 0.08947134065170427], CircleBox[{0, 0}, 0.4332636446740114]}}}, AspectRatio->Automatic, Axes->True, ColorOutput->Automatic, Frame->True, GridLines->Automatic, PlotRange->{{-2.4, 2.4}, {-2.4, 2.4}}], ",", GraphicsBox[{ {Thickness[0.01], {RGBColor[1, 0, 0], CircleBox[{0, -0.2721785497620972}, 0.2721785497620972], CircleBox[{0, 0.2721785497620972}, 0.2721785497620972]}, {RGBColor[1, 0, 1], CircleBox[{-0.2721785497620972, 0}, 0.2721785497620972], CircleBox[{0.2721785497620972, 0}, 0.2721785497620972]}, {RGBColor[0, 1, 1], CircleBox[{0, 0}, 0.2721785497620972]}}, {Thickness[0.015], {RGBColor[0, 0, 0], CircleBox[{0, 0}, 0.08444026225663692], CircleBox[{0, 0}, 0.404076281771941]}}}, AspectRatio->Automatic, Axes->True, ColorOutput->Automatic, Frame->True, GridLines->Automatic, PlotRange->{{-2.4, 2.4}, {-2.4, 2.4}}]}], "}"}]], "Output", CellChangeTimes->{ 3.5662334443757825`*^9, 3.5662334769330397`*^9, 3.566233748233116*^9, { 3.566233782428376*^9, 3.566233929833035*^9}, 3.566233966134299*^9, 3.566234199526309*^9, 3.566234318272917*^9, 3.5662343908754444`*^9, 3.5662348042873716`*^9, 3.5662349942333055`*^9, 3.5662354528393116`*^9}] }, Open ]] }, WindowSize->{774, 815}, WindowMargins->{{562, Automatic}, {-39, Automatic}}, PrivateNotebookOptions->{"VersionedStylesheet"->{"Default.nb"[8.] -> False}}, FrontEndVersion->"9.0 for Microsoft Windows (64-bit) (November 20, 2012)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[557, 20, 468, 7, 48, "Input"], Cell[CellGroupData[{ Cell[1050, 31, 4044, 122, 212, "Input"], Cell[5097, 155, 4167, 105, 252, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[9301, 265, 3270, 101, 192, "Input"], Cell[12574, 368, 1229, 31, 31, "Output"], Cell[13806, 401, 1224, 31, 52, "Output"], Cell[15033, 434, 986, 19, 52, "Output"], Cell[16022, 455, 1224, 31, 52, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[17283, 491, 4756, 116, 252, "Input"], Cell[22042, 609, 9124, 223, 833, "Output"] }, Open ]] } ] *) (* End of internal cache information *)