(* 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[     33441,        861]
NotebookOptionsPosition[     32643,        832]
NotebookOutlinePosition[     33066,        848]
CellTagsIndexPosition[     33023,        845]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"Ne", "-", 
     RowBox[{"Atom", " ", "and", " ", "homologues", " ", "1", "s", 
      RowBox[{"(", "2", ")"}], "2", "s", 
      RowBox[{"(", "2", ")"}], " ", "2", "p", 
      RowBox[{"(", "6", ")"}]}]}], " ", ",", " ", 
    RowBox[{"all", " ", "p", " ", "equal", " ", "size", " ", 
     RowBox[{"(", "?", ")"}]}], ",", 
    RowBox[{
     RowBox[{
     "with", " ", "exchange", " ", "correction", " ", "page", " ", "40"}], 
     ";", " ", 
     RowBox[{"8.12", 
      RowBox[{".2011", "/", "10.04"}], ".2013"}]}]}], " ", "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"Clear", "[", 
     RowBox[{
     "T", ",", "Vne", ",", " ", "Vee", ",", "func", ",", "Vir", ",", "P", ",",
       "Q", ",", "R", ",", "t", ",", "z", ",", "sig", ",", "max"}], "]"}], 
    ";", " ", 
    RowBox[{"sig", "=", "0.285"}], ";", " ", 
    RowBox[{"max", "=", "12"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"T", "=", 
     RowBox[{
      RowBox[{"0.992", "*", 
       RowBox[{"9", "/", 
        RowBox[{"(", 
         RowBox[{"4", "*", 
          RowBox[{"P", "^", "2"}]}], ")"}]}]}], "+", 
      RowBox[{"0.999", "*", "9", "*", 
       RowBox[{"4", "/", 
        RowBox[{"(", 
         RowBox[{"4", "*", 
          RowBox[{"Q", "^", "2"}]}], ")"}]}]}], "+", 
      RowBox[{"12", "/", 
       RowBox[{"R", "^", "2"}]}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Vne", "=", 
     RowBox[{
      RowBox[{
       RowBox[{"-", "3"}], "*", 
       RowBox[{"z", "/", "P"}]}], "-", 
      RowBox[{"3", "*", 
       RowBox[{"z", "/", "Q"}]}], "-", 
      RowBox[{"6", "*", 
       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[{"12", "*", 
       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[{"4.5", "*", 
       RowBox[{"sig", "/", "R"}]}], "+", 
      RowBox[{"12", "*", 
       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[{"6.582336", "/", "R"}]}]}], ";"}], " ", 
   RowBox[{"(*", " ", 
    RowBox[{"8.48528", "/", "R"}], " ", "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"func", "=", 
     RowBox[{"T", "+", "Vne", "+", "Vee"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    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.07", ",", "0.16"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"Q", ",", "0.3", ",", "1.6"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"R", ",", "0.23", ",", "0.63"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"MaxIterations", "\[Rule]", "800"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"Method", "\[Rule]", "Automatic"}], "}"}]}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"z", ",", "10", ",", "21"}], "}"}]}], "]"}]}]}]}]], "Input",
 CellChangeTimes->{{3.5650905841503596`*^9, 3.5650906167700167`*^9}, {
   3.5650906713233128`*^9, 3.565090682914133*^9}, {3.5650908040951457`*^9, 
   3.5650908556844363`*^9}, {3.5650909164933434`*^9, 3.5650909189737473`*^9}, 
   3.565091060341196*^9, {3.5650911731137943`*^9, 3.5650911987602386`*^9}, {
   3.5650977150648775`*^9, 3.5650977777457876`*^9}, {3.56509786437274*^9, 
   3.5650978692087483`*^9}, {3.5650979228416424`*^9, 3.565098002058582*^9}, {
   3.5650981295420055`*^9, 3.5650981774028893`*^9}, {3.565098221301367*^9, 
   3.565098223937772*^9}, {3.5650982722978563`*^9, 3.565098279489469*^9}, 
   3.5650983268823524`*^9, {3.5650984202017164`*^9, 3.5650984208569174`*^9}, {
   3.5746057867575455`*^9, 3.574605836069232*^9}, {3.574605926580591*^9, 
   3.5746059374538107`*^9}, {3.5746092646432543`*^9, 
   3.5746092733168697`*^9}, {3.5746094926688547`*^9, 3.574609496210061*^9}, {
   3.5746095739137974`*^9, 3.574609582353412*^9}, {3.574609677295179*^9, 
   3.5746096775291796`*^9}, {3.574609866835512*^9, 3.5746099292512217`*^9}, {
   3.5746099721200967`*^9, 3.574609972806498*^9}, {3.574610100804723*^9, 
   3.574610263154208*^9}, {3.574610310952692*^9, 3.5746103551943693`*^9}, {
   3.574610387829627*^9, 3.5746103886720285`*^9}, {3.5746104792769876`*^9, 
   3.5746104969050183`*^9}, {3.57461053177108*^9, 3.57461054889991*^9}, {
   3.574610587619178*^9, 3.57461064028487*^9}, {3.574610709533392*^9, 
   3.5746107138233995`*^9}, {3.5746108731464796`*^9, 
   3.5746108762820845`*^9}, {3.5746109163273554`*^9, 3.574610923612568*^9}, {
   3.5746113636113405`*^9, 3.574611369913752*^9}, {3.576143373886182*^9, 
   3.576143483070773*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "129.38259834303952`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.15345587139622607`"}], ",", 
       RowBox[{"Q", "\[Rule]", "1.6010380759400116`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.6212428915437588`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "162.45693750348283`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.139153091770794`"}], ",", 
       RowBox[{"Q", "\[Rule]", "1.3015090469763464`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.5352792600072026`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "199.5471494650992`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.1272870456796772`"}], ",", 
       RowBox[{"Q", "\[Rule]", "1.0350319503299297`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.47407331201124425`"}]}], "}"}]}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "240.7550989351488`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.11728966965060303`"}], ",", 
       RowBox[{"Q", "\[Rule]", "0.7531404879903423`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.43764118177483596`"}]}], "}"}]}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "286.04558693253426`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.10875653174508455`"}], ",", 
       RowBox[{"Q", "\[Rule]", "0.6484484857689474`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.39762982861344637`"}]}], "}"}]}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "335.36277807567984`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.10135787081911282`"}], ",", 
       RowBox[{"Q", "\[Rule]", "0.5784548729068761`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.3626738319176281`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "388.7011148808075`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.09490890714492488`"}], ",", 
       RowBox[{"Q", "\[Rule]", "0.5244047743191271`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.33299720864285964`"}]}], "}"}]}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "446.05889780758264`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.0892303055238827`"}], ",", 
       RowBox[{"Q", "\[Rule]", "0.4804698927161771`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.3076536503937662`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "507.4353216713465`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.08419213442898728`"}], ",", 
       RowBox[{"Q", "\[Rule]", "0.4437279503153145`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.2858235861926292`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "572.8299175364587`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.07969111416683325`"}], ",", 
       RowBox[{"Q", "\[Rule]", "0.412457029720728`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.2668482887343572`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "642.2423770176188`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.07564733674124068`"}], ",", 
       RowBox[{"Q", "\[Rule]", "0.3853729210326859`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.2502189080072669`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "715.6724815800172`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.07199288077147124`"}], ",", 
       RowBox[{"Q", "\[Rule]", "0.36176030161490674`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.23552246818552905`"}]}], "}"}]}], "}"}]}], 
  "}"}]], "Output",
 CellChangeTimes->{3.5746109380425935`*^9, 3.574610997041897*^9, 
  3.574611184881827*^9, 3.57461123230591*^9, 3.5761434438523045`*^9, 
  3.5761434922123895`*^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", "-", "9"}], "}"}], ",", "2"}], "]"}], "]"}]}], ",", 
       "10"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"z", " ", ",", "10", ",", 
       RowBox[{"9", "+", "max"}]}], "}"}]}], "]"}], "\[IndentingNewLine]", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"N", "[", 
      RowBox[{
       RowBox[{"Vne", " ", "/.", " ", 
        RowBox[{"t", "[", 
         RowBox[{"[", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{"z", "-", "9"}], "}"}], ",", "2"}], "]"}], "]"}]}], ",", 
       "10"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"z", " ", ",", "10", ",", 
       RowBox[{"9", "+", "max"}]}], "}"}]}], "]"}], "\[IndentingNewLine]", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"N", "[", 
      RowBox[{
       RowBox[{"Vee", " ", "/.", " ", 
        RowBox[{"t", "[", 
         RowBox[{"[", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{"z", "-", "9"}], "}"}], ",", "2"}], "]"}], "]"}]}], ",", 
       "10"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"z", " ", ",", "10", ",", 
       RowBox[{"9", "+", "max"}]}], "}"}]}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.5650905841503596`*^9, 3.5650906167700167`*^9}, {
   3.5650906713233128`*^9, 3.565090682914133*^9}, {3.5650908040951457`*^9, 
   3.5650908556844363`*^9}, {3.5650909164933434`*^9, 
   3.5650909189737473`*^9}, {3.565091060341196*^9, 3.565091078874028*^9}, {
   3.565091206123452*^9, 3.5650912083386555`*^9}, {3.5746092956405087`*^9, 
   3.574609335779379*^9}, {3.5746098161354227`*^9, 3.574609841407467*^9}, {
   3.5746109807866683`*^9, 3.5746109859190774`*^9}, 3.57461106735122*^9, {
   3.57461115246497*^9, 3.574611170607802*^9}, {3.574611213461077*^9, 
   3.5746112265339003`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"-", "2.000000475832554`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1.999999386457429`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1.9999995402954787`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "2.0000184423243166`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "2.0002145032418968`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "2.000001870731068`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "2.0000012576231976`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1.9999999813321274`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1.999998568564546`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "2.0000001608278573`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1.999999607874518`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "2.0000005165121735`"}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.5746109380893936`*^9, 3.5746109970886974`*^9, 
  3.574611184928627*^9, 3.5746112323371105`*^9, 3.5761434442891054`*^9, 
  3.5761434922435894`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"-", "310.81436601546875`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "385.8041895627786`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "469.4821036374029`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "562.521469115292`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "662.205445762744`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "769.9217527898759`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "885.5711579802546`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1009.2424369974066`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1140.9417205626896`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1280.6669269167019`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1428.4275738033832`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1584.2157951641143`"}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.5746109380893936`*^9, 3.5746109970886974`*^9, 
  3.574611184928627*^9, 3.5746112323371105`*^9, 3.5761434442891054`*^9, 
  3.5761434922435894`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", "52.04923089381258`", "}"}], ",", 
   RowBox[{"{", "60.89021488150466`", "}"}], ",", 
   RowBox[{"{", "70.38771297443549`", "}"}], ",", 
   RowBox[{"{", "81.01571124672597`", "}"}], ",", 
   RowBox[{"{", "90.17561644479864`", "}"}], ",", 
   RowBox[{"{", "99.19682401091057`", "}"}], ",", 
   RowBox[{"{", "108.1694170575638`", "}"}], ",", 
   RowBox[{"{", "117.12463305527044`", "}"}], ",", 
   RowBox[{"{", "126.0703508580465`", "}"}], ",", 
   RowBox[{"{", "135.0071839707781`", "}"}], ",", 
   RowBox[{"{", "143.94256792844533`", "}"}], ",", 
   RowBox[{"{", "152.87120165743792`", "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.5746109380893936`*^9, 3.5746109970886974`*^9, 
  3.574611184928627*^9, 3.5746112323371105`*^9, 3.5761434442891054`*^9, 
  3.5761434922435894`*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"(*", " ", 
  RowBox[{"Etot", " ", "in", " ", "eV"}], " ", "*)"}]], "Input",
 CellChangeTimes->{{3.5650905841503596`*^9, 3.5650906167700167`*^9}, {
  3.5650906713233128`*^9, 3.565090682914133*^9}, {3.5650908040951457`*^9, 
  3.5650908556844363`*^9}, {3.5650909164933434`*^9, 3.5650909189737473`*^9}, {
  3.565091060341196*^9, 3.565091067579609*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"PA", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"P", " ", "/.", " ", 
      RowBox[{"t", "[", 
       RowBox[{"[", 
        RowBox[{"i", ",", "2"}], "]"}], "]"}]}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "max"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"QA", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"Q", " ", "/.", " ", 
      RowBox[{"t", "[", 
       RowBox[{"[", 
        RowBox[{"i", ",", "2"}], "]"}], "]"}]}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "max"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"RA", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"R", " ", "/.", " ", 
      RowBox[{"t", "[", 
       RowBox[{"[", 
        RowBox[{"i", ",", "2"}], "]"}], "]"}]}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "max"}], "}"}]}], "]"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.5650905841503596`*^9, 3.5650906167700167`*^9}, {
   3.5650906713233128`*^9, 3.565090682914133*^9}, {3.5650908040951457`*^9, 
   3.5650908556844363`*^9}, {3.5650909164933434`*^9, 
   3.5650909189737473`*^9}, {3.565091060341196*^9, 3.5650910731176186`*^9}, {
   3.565097584804649*^9, 3.5650975899994583`*^9}, {3.5650976564711747`*^9, 
   3.5650976693879976`*^9}, {3.5650987312506623`*^9, 
   3.5650987434654837`*^9}, {3.565098960711466*^9, 3.56509896351947*^9}, {
   3.5650992007490873`*^9, 3.565099203135891*^9}, {3.565099728076813*^9, 
   3.5650997446440425`*^9}, {3.5650997851105137`*^9, 3.565099863469451*^9}, {
   3.565099907024728*^9, 3.565099907211928*^9}, {3.565099947522399*^9, 
   3.5650999484116*^9}, 3.5651002560285406`*^9, 3.5651004699673166`*^9, {
   3.5651142606546755`*^9, 3.565114270139492*^9}, {3.5746053843080387`*^9, 
   3.574605384464039*^9}, 3.574607264376541*^9, {3.574607420298815*^9, 
   3.574607421546817*^9}, {3.5746074897657366`*^9, 3.574607556643054*^9}, {
   3.5746076010407324`*^9, 3.574607602335534*^9}, {3.5746076352515926`*^9, 
   3.574607644440008*^9}, 3.5746076911932907`*^9, 3.574607728992157*^9, {
   3.5746077986462793`*^9, 3.574607823231922*^9}, {3.574608622889327*^9, 
   3.574608665649002*^9}, 3.5746087344919233`*^9, {3.574608797734434*^9, 
   3.5746088048012466`*^9}, {3.5746093421909904`*^9, 
   3.5746093582434187`*^9}, {3.574610808718366*^9, 3.5746108452692304`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"Table", "[", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Show", "[", 
     RowBox[{
      RowBox[{"Graphics", "[", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"Thickness", "[", "0.015", "]"}], ",", 
           RowBox[{"Circle", "[", 
            RowBox[{
             RowBox[{"{", 
              RowBox[{"0", ",", " ", "0"}], "}"}], ",", " ", 
             RowBox[{"PA", "[", 
              RowBox[{"[", "i", "]"}], "]"}]}], "]"}], ",", " ", 
           RowBox[{"Circle", "[", 
            RowBox[{
             RowBox[{"{", 
              RowBox[{"0", ",", " ", "0"}], "}"}], ",", " ", 
             RowBox[{"QA", "[", 
              RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}], "}"}], ",", 
         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", "]"}], "]"}]}], "]"}]}], "}"}]}], "}"}]}], 
        "}"}], "]"}], ",", " ", "\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[{"-", "1.6"}], ",", " ", "1.6"}], "}"}], ",", " ", 
           RowBox[{"{", 
            RowBox[{
             RowBox[{"-", "1.6"}], ",", " ", "1.6"}], "}"}]}], "}"}]}], ",", 
        " ", 
        RowBox[{"Frame", " ", "->", " ", "True"}]}], "}"}]}], "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"i", ",", "1", ",", "max"}], "}"}]}], "]"}], 
  "\[IndentingNewLine]"}]], "Input",
 CellChangeTimes->{{3.5650905841503596`*^9, 3.5650906167700167`*^9}, {
   3.5650906713233128`*^9, 3.565090682914133*^9}, {3.5650908040951457`*^9, 
   3.5650908556844363`*^9}, {3.5650909164933434`*^9, 
   3.5650909189737473`*^9}, {3.565091060341196*^9, 3.5650910970792603`*^9}, {
   3.565091263016752*^9, 3.5650912652007556`*^9}, {3.565091306478428*^9, 
   3.565091306915229*^9}, {3.5650987968175774`*^9, 3.5650988590772867`*^9}, {
   3.565099041613208*^9, 3.5650990568388343`*^9}, {3.5746068435190015`*^9, 
   3.574606848651411*^9}, {3.574606879586265*^9, 3.5746069350599623`*^9}, {
   3.5746071115118723`*^9, 3.574607140325123*^9}, {3.5746080941731987`*^9, 
   3.574608110802828*^9}, {3.5746081704417324`*^9, 3.5746081717833347`*^9}, 
   3.5746082427010593`*^9, {3.5746082974727554`*^9, 3.5746083203267956`*^9}, 
   3.574608402944541*^9, 3.5746093666362333`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   GraphicsBox[{
     {Thickness[0.015], CircleBox[{0, 0}, 0.15345587139622607`], 
      CircleBox[{0, 0}, 1.6010380759400116`]}, 
     {Thickness[0.01], 
      {RGBColor[1, 0, 0], 
       CircleBox[{0, -0.6212428915437588}, 0.6212428915437588], 
       CircleBox[{0, 0.6212428915437588}, 0.6212428915437588]}, 
      {RGBColor[1, 0, 1], 
       CircleBox[{-0.6212428915437588, 0}, 0.6212428915437588], 
       CircleBox[{0.6212428915437588, 0}, 0.6212428915437588]}, 
      {RGBColor[0, 1, 1], CircleBox[{0, 0}, 0.6212428915437588]}}},
    AspectRatio->Automatic,
    Axes->True,
    ColorOutput->Automatic,
    Frame->True,
    GridLines->Automatic,
    PlotRange->{{-1.6, 1.6}, {-1.6, 1.6}}], ",", 
   GraphicsBox[{
     {Thickness[0.015], CircleBox[{0, 0}, 0.139153091770794], 
      CircleBox[{0, 0}, 1.3015090469763464`]}, 
     {Thickness[0.01], 
      {RGBColor[1, 0, 0], 
       CircleBox[{0, -0.5352792600072026}, 0.5352792600072026], 
       CircleBox[{0, 0.5352792600072026}, 0.5352792600072026]}, 
      {RGBColor[1, 0, 1], 
       CircleBox[{-0.5352792600072026, 0}, 0.5352792600072026], 
       CircleBox[{0.5352792600072026, 0}, 0.5352792600072026]}, 
      {RGBColor[0, 1, 1], CircleBox[{0, 0}, 0.5352792600072026]}}},
    AspectRatio->Automatic,
    Axes->True,
    ColorOutput->Automatic,
    Frame->True,
    GridLines->Automatic,
    PlotRange->{{-1.6, 1.6}, {-1.6, 1.6}}], ",", 
   GraphicsBox[{
     {Thickness[0.015], CircleBox[{0, 0}, 0.1272870456796772], 
      CircleBox[{0, 0}, 1.0350319503299297`]}, 
     {Thickness[0.01], 
      {RGBColor[1, 0, 0], 
       CircleBox[{0, -0.47407331201124425`}, 0.47407331201124425`], 
       CircleBox[{0, 0.47407331201124425`}, 0.47407331201124425`]}, 
      {RGBColor[1, 0, 1], 
       CircleBox[{-0.47407331201124425`, 0}, 0.47407331201124425`], 
       CircleBox[{0.47407331201124425`, 0}, 0.47407331201124425`]}, 
      {RGBColor[0, 1, 1], CircleBox[{0, 0}, 0.47407331201124425`]}}},
    AspectRatio->Automatic,
    Axes->True,
    ColorOutput->Automatic,
    Frame->True,
    GridLines->Automatic,
    PlotRange->{{-1.6, 1.6}, {-1.6, 1.6}}], ",", 
   GraphicsBox[{
     {Thickness[0.015], CircleBox[{0, 0}, 0.11728966965060303`], 
      CircleBox[{0, 0}, 0.7531404879903423]}, 
     {Thickness[0.01], 
      {RGBColor[1, 0, 0], 
       CircleBox[{0, -0.43764118177483596`}, 0.43764118177483596`], 
       CircleBox[{0, 0.43764118177483596`}, 0.43764118177483596`]}, 
      {RGBColor[1, 0, 1], 
       CircleBox[{-0.43764118177483596`, 0}, 0.43764118177483596`], 
       CircleBox[{0.43764118177483596`, 0}, 0.43764118177483596`]}, 
      {RGBColor[0, 1, 1], CircleBox[{0, 0}, 0.43764118177483596`]}}},
    AspectRatio->Automatic,
    Axes->True,
    ColorOutput->Automatic,
    Frame->True,
    GridLines->Automatic,
    PlotRange->{{-1.6, 1.6}, {-1.6, 1.6}}], ",", 
   GraphicsBox[{
     {Thickness[0.015], CircleBox[{0, 0}, 0.10875653174508455`], 
      CircleBox[{0, 0}, 0.6484484857689474]}, 
     {Thickness[0.01], 
      {RGBColor[1, 0, 0], 
       CircleBox[{0, -0.39762982861344637`}, 0.39762982861344637`], 
       CircleBox[{0, 0.39762982861344637`}, 0.39762982861344637`]}, 
      {RGBColor[1, 0, 1], 
       CircleBox[{-0.39762982861344637`, 0}, 0.39762982861344637`], 
       CircleBox[{0.39762982861344637`, 0}, 0.39762982861344637`]}, 
      {RGBColor[0, 1, 1], CircleBox[{0, 0}, 0.39762982861344637`]}}},
    AspectRatio->Automatic,
    Axes->True,
    ColorOutput->Automatic,
    Frame->True,
    GridLines->Automatic,
    PlotRange->{{-1.6, 1.6}, {-1.6, 1.6}}], ",", 
   GraphicsBox[{
     {Thickness[0.015], CircleBox[{0, 0}, 0.10135787081911282`], 
      CircleBox[{0, 0}, 0.5784548729068761]}, 
     {Thickness[0.01], 
      {RGBColor[1, 0, 0], 
       CircleBox[{0, -0.3626738319176281}, 0.3626738319176281], 
       CircleBox[{0, 0.3626738319176281}, 0.3626738319176281]}, 
      {RGBColor[1, 0, 1], 
       CircleBox[{-0.3626738319176281, 0}, 0.3626738319176281], 
       CircleBox[{0.3626738319176281, 0}, 0.3626738319176281]}, 
      {RGBColor[0, 1, 1], CircleBox[{0, 0}, 0.3626738319176281]}}},
    AspectRatio->Automatic,
    Axes->True,
    ColorOutput->Automatic,
    Frame->True,
    GridLines->Automatic,
    PlotRange->{{-1.6, 1.6}, {-1.6, 1.6}}], ",", 
   GraphicsBox[{
     {Thickness[0.015], CircleBox[{0, 0}, 0.09490890714492488], 
      CircleBox[{0, 0}, 0.5244047743191271]}, 
     {Thickness[0.01], 
      {RGBColor[1, 0, 0], 
       CircleBox[{0, -0.33299720864285964`}, 0.33299720864285964`], 
       CircleBox[{0, 0.33299720864285964`}, 0.33299720864285964`]}, 
      {RGBColor[1, 0, 1], 
       CircleBox[{-0.33299720864285964`, 0}, 0.33299720864285964`], 
       CircleBox[{0.33299720864285964`, 0}, 0.33299720864285964`]}, 
      {RGBColor[0, 1, 1], CircleBox[{0, 0}, 0.33299720864285964`]}}},
    AspectRatio->Automatic,
    Axes->True,
    ColorOutput->Automatic,
    Frame->True,
    GridLines->Automatic,
    PlotRange->{{-1.6, 1.6}, {-1.6, 1.6}}], ",", 
   GraphicsBox[{
     {Thickness[0.015], CircleBox[{0, 0}, 0.0892303055238827], 
      CircleBox[{0, 0}, 0.4804698927161771]}, 
     {Thickness[0.01], 
      {RGBColor[1, 0, 0], 
       CircleBox[{0, -0.3076536503937662}, 0.3076536503937662], 
       CircleBox[{0, 0.3076536503937662}, 0.3076536503937662]}, 
      {RGBColor[1, 0, 1], 
       CircleBox[{-0.3076536503937662, 0}, 0.3076536503937662], 
       CircleBox[{0.3076536503937662, 0}, 0.3076536503937662]}, 
      {RGBColor[0, 1, 1], CircleBox[{0, 0}, 0.3076536503937662]}}},
    AspectRatio->Automatic,
    Axes->True,
    ColorOutput->Automatic,
    Frame->True,
    GridLines->Automatic,
    PlotRange->{{-1.6, 1.6}, {-1.6, 1.6}}], ",", 
   GraphicsBox[{
     {Thickness[0.015], CircleBox[{0, 0}, 0.08419213442898728], 
      CircleBox[{0, 0}, 0.4437279503153145]}, 
     {Thickness[0.01], 
      {RGBColor[1, 0, 0], 
       CircleBox[{0, -0.2858235861926292}, 0.2858235861926292], 
       CircleBox[{0, 0.2858235861926292}, 0.2858235861926292]}, 
      {RGBColor[1, 0, 1], 
       CircleBox[{-0.2858235861926292, 0}, 0.2858235861926292], 
       CircleBox[{0.2858235861926292, 0}, 0.2858235861926292]}, 
      {RGBColor[0, 1, 1], CircleBox[{0, 0}, 0.2858235861926292]}}},
    AspectRatio->Automatic,
    Axes->True,
    ColorOutput->Automatic,
    Frame->True,
    GridLines->Automatic,
    PlotRange->{{-1.6, 1.6}, {-1.6, 1.6}}], ",", 
   GraphicsBox[{
     {Thickness[0.015], CircleBox[{0, 0}, 0.07969111416683325], 
      CircleBox[{0, 0}, 0.412457029720728]}, 
     {Thickness[0.01], 
      {RGBColor[1, 0, 0], 
       CircleBox[{0, -0.2668482887343572}, 0.2668482887343572], 
       CircleBox[{0, 0.2668482887343572}, 0.2668482887343572]}, 
      {RGBColor[1, 0, 1], 
       CircleBox[{-0.2668482887343572, 0}, 0.2668482887343572], 
       CircleBox[{0.2668482887343572, 0}, 0.2668482887343572]}, 
      {RGBColor[0, 1, 1], CircleBox[{0, 0}, 0.2668482887343572]}}},
    AspectRatio->Automatic,
    Axes->True,
    ColorOutput->Automatic,
    Frame->True,
    GridLines->Automatic,
    PlotRange->{{-1.6, 1.6}, {-1.6, 1.6}}], ",", 
   GraphicsBox[{
     {Thickness[0.015], CircleBox[{0, 0}, 0.07564733674124068], 
      CircleBox[{0, 0}, 0.3853729210326859]}, 
     {Thickness[0.01], 
      {RGBColor[1, 0, 0], 
       CircleBox[{0, -0.2502189080072669}, 0.2502189080072669], 
       CircleBox[{0, 0.2502189080072669}, 0.2502189080072669]}, 
      {RGBColor[1, 0, 1], 
       CircleBox[{-0.2502189080072669, 0}, 0.2502189080072669], 
       CircleBox[{0.2502189080072669, 0}, 0.2502189080072669]}, 
      {RGBColor[0, 1, 1], CircleBox[{0, 0}, 0.2502189080072669]}}},
    AspectRatio->Automatic,
    Axes->True,
    ColorOutput->Automatic,
    Frame->True,
    GridLines->Automatic,
    PlotRange->{{-1.6, 1.6}, {-1.6, 1.6}}], ",", 
   GraphicsBox[{
     {Thickness[0.015], CircleBox[{0, 0}, 0.07199288077147124], 
      CircleBox[{0, 0}, 0.36176030161490674`]}, 
     {Thickness[0.01], 
      {RGBColor[1, 0, 0], 
       CircleBox[{0, -0.23552246818552905`}, 0.23552246818552905`], 
       CircleBox[{0, 0.23552246818552905`}, 0.23552246818552905`]}, 
      {RGBColor[1, 0, 1], 
       CircleBox[{-0.23552246818552905`, 0}, 0.23552246818552905`], 
       CircleBox[{0.23552246818552905`, 0}, 0.23552246818552905`]}, 
      {RGBColor[0, 1, 1], CircleBox[{0, 0}, 0.23552246818552905`]}}},
    AspectRatio->Automatic,
    Axes->True,
    ColorOutput->Automatic,
    Frame->True,
    GridLines->Automatic,
    PlotRange->{{-1.6, 1.6}, {-1.6, 1.6}}]}], "}"}]], "Output",
 CellChangeTimes->{3.5746109381205935`*^9, 3.574610997119897*^9, 
  3.574611184975427*^9, 3.57461123238391*^9, 3.5761434463483086`*^9, 
  3.5761434922903895`*^9}]
}, Open  ]]
},
WindowSize->{827, 815},
WindowMargins->{{484, Automatic}, {Automatic, 26}},
PrivateNotebookOptions->{"VersionedStylesheet"->{"Default.nb"[8.] -> False}},
FrontEndVersion->"9.0 for Microsoft Windows (64-bit) (January 25, 2013)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[579, 22, 6418, 164, 272, "Input"],
Cell[7000, 188, 4045, 106, 252, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[11082, 299, 2083, 55, 92, "Input"],
Cell[13168, 356, 1080, 29, 52, "Output"],
Cell[14251, 387, 1077, 29, 52, "Output"],
Cell[15331, 418, 836, 17, 52, "Output"]
}, Open  ]],
Cell[16182, 438, 372, 6, 31, "Input"],
Cell[16557, 446, 2396, 52, 72, "Input"],
Cell[CellGroupData[{
Cell[18978, 502, 4936, 116, 212, "Input"],
Cell[23917, 620, 8710, 209, 833, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)