(* 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[     24716,        694]
NotebookOptionsPosition[     23695,        658]
NotebookOutlinePosition[     24275,        679]
CellTagsIndexPosition[     24232,        676]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[TextData[StyleBox["(*Ne-Atom and homologues 1s(2)2s(2) 2p(6), all p \
equal size, with exchange correction page 40; 25.12.2012*)",
 FontFamily->"Courier New"]], "Text",
 CellChangeTimes->{{3.5654359755976987`*^9, 3.565435976689701*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"Clear", "[", 
   RowBox[{"z", ",", "sig"}], "]"}], ";", " ", 
  RowBox[{"sig", "=", "0.24"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"T", "=", 
   RowBox[{
    RowBox[{"9", "/", 
     RowBox[{"(", 
      RowBox[{"4", "*", 
       RowBox[{"P", "^", "2"}]}], ")"}]}], "+", 
    RowBox[{"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[{
   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[{
   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.15"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"Q", ",", "1.6"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"R", ",", "0.6"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"MaxIterations", "\[Rule]", "800"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"Method", "\[Rule]", "\"\<Newton\>\""}], "}"}]}], "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"z", ",", "10", ",", "18"}], "}"}]}], "]"}]}]}], "Input",
 CellChangeTimes->{3.565347278648552*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "129.90549560209547`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.15398421744914842`"}], ",", 
       RowBox[{"Q", "\[Rule]", "1.5885139078383728`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.6177358807009392`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "162.97478956201732`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.13969192675004327`"}], ",", 
       RowBox[{"Q", "\[Rule]", "1.2913825954113067`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.5327797060592275`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "200.04569312208884`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.12782468311414308`"}], ",", 
       RowBox[{"Q", "\[Rule]", "1.02443385856661`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.4724387728144514`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "241.2194433693671`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.11781805218407913`"}], ",", 
       RowBox[{"Q", "\[Rule]", "0.74977754076882`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.4360242734454343`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "286.45459716903`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.1092617149883094`"}], ",", 
       RowBox[{"Q", "\[Rule]", "0.6466370326233891`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.39616563882769434`"}]}], "}"}]}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "335.6994761497142`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.10186175056008497`"}], ",", 
       RowBox[{"Q", "\[Rule]", "0.5772654286442558`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.3614697967847427`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "388.9487782394322`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.0953990268490955`"}], ",", 
       RowBox[{"Q", "\[Rule]", "0.5235313645874632`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.33197443067159177`"}]}], "}"}]}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "446.20084439126464`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.0897062911586054`"}], ",", 
       RowBox[{"Q", "\[Rule]", "0.47979765093440596`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.3067791621916184`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "507.45488567684595`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.08465380331289858`"}], ",", 
       RowBox[{"Q", "\[Rule]", "0.44320969672752136`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.28506814451379575`"}]}], "}"}]}], "}"}]}], 
  "}"}]], "Output",
 CellChangeTimes->{
  3.5654322794695463`*^9, {3.5654323231964226`*^9, 3.565432348764868*^9}, {
   3.56543240690617*^9, 3.5654324480122423`*^9}, 3.5654351182513933`*^9, 
   3.5654351887011166`*^9, {3.5654352375604024`*^9, 3.565435249853224*^9}, {
   3.5654354295811396`*^9, 3.5654354355715504`*^9}, 3.5654355591861677`*^9, 
   3.565435704640823*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"(*", " ", 
  RowBox[{"Virial", " ", "ratio"}], " ", "*)"}]], "Input"],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"Table", "[", 
  RowBox[{
   RowBox[{"N", "[", 
    RowBox[{
     RowBox[{"Vir", " ", "/.", " ", 
      RowBox[{"t", "[", 
       RowBox[{"[", 
        RowBox[{
         RowBox[{"z", "-", "9"}], ",", "2"}], "]"}], "]"}]}], ",", "10"}], 
    "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"z", ",", "10", ",", "18"}], "}"}]}], 
  "]"}], "\[IndentingNewLine]", 
 RowBox[{"Table", "[", 
  RowBox[{
   RowBox[{"N", "[", 
    RowBox[{
     RowBox[{"Vne", " ", "/.", " ", 
      RowBox[{"t", "[", 
       RowBox[{"[", 
        RowBox[{
         RowBox[{"z", "-", "9"}], ",", "2"}], "]"}], "]"}]}], ",", "10"}], 
    "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"z", ",", "10", ",", "18"}], "}"}]}], 
  "]"}], "\[IndentingNewLine]", 
 RowBox[{"Table", "[", 
  RowBox[{
   RowBox[{"N", "[", 
    RowBox[{
     RowBox[{"Vee", " ", "/.", " ", 
      RowBox[{"t", "[", 
       RowBox[{"[", 
        RowBox[{
         RowBox[{"z", "-", "9"}], ",", "2"}], "]"}], "]"}]}], ",", "10"}], 
    "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"z", ",", "10", ",", "18"}], "}"}]}], "]"}]}], "Input",
 CellChangeTimes->{
  3.5653472854657636`*^9, {3.5653474173639956`*^9, 3.56534748266571*^9}, {
   3.5653476185107484`*^9, 3.565347687197669*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", "1.9999999999999964`"}], ",", 
   RowBox[{"-", "1.9999999999630884`"}], ",", 
   RowBox[{"-", "1.9999999999999936`"}], ",", 
   RowBox[{"-", "1.9999999954166932`"}], ",", 
   RowBox[{"-", "1.9999999999806777`"}], ",", 
   RowBox[{"-", "1.9999999999904157`"}], ",", 
   RowBox[{"-", "1.9999999999998836`"}], ",", 
   RowBox[{"-", "1.9999999999964242`"}], ",", 
   RowBox[{"-", "1.9999999999988225`"}]}], "}"}]], "Output",
 CellChangeTimes->{
  3.565432279500746*^9, {3.565432323212023*^9, 3.5654323487804675`*^9}, {
   3.5654324069217696`*^9, 3.565432448043442*^9}, 3.565435118282593*^9, 
   3.565435188716717*^9, {3.5654352375916023`*^9, 3.565435249868824*^9}, {
   3.56543542959674*^9, 3.5654354355871506`*^9}, 3.565435559201768*^9, 
   3.565435704687623*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", "310.83962654506786`"}], ",", 
   RowBox[{"-", "385.6667237935207`"}], ",", 
   RowBox[{"-", "469.17780773429934`"}], ",", 
   RowBox[{"-", "561.9234344188847`"}], ",", 
   RowBox[{"-", "661.3820931128294`"}], ",", 
   RowBox[{"-", "768.7124464586709`"}], ",", 
   RowBox[{"-", "884.0137357751692`"}], ",", 
   RowBox[{"-", "1007.303537477298`"}], ",", 
   RowBox[{"-", "1138.5874621119622`"}]}], "}"}]], "Output",
 CellChangeTimes->{
  3.565432279500746*^9, {3.565432323212023*^9, 3.5654323487804675`*^9}, {
   3.5654324069217696`*^9, 3.565432448043442*^9}, 3.565435118282593*^9, 
   3.565435188716717*^9, {3.5654352375916023`*^9, 3.565435249868824*^9}, {
   3.56543542959674*^9, 3.5654354355871506`*^9}, 3.565435559201768*^9, 
   3.565435704687623*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "51.02863534087648`", ",", "59.71714466347031`", ",", "69.08642149012044`", 
   ",", "79.48454657456786`", ",", "88.4728987692344`", ",", 
   "97.31349415602482`", ",", "106.11617929625945`", ",", 
   "114.90184869317312`", ",", "123.67769075767262`"}], "}"}]], "Output",
 CellChangeTimes->{
  3.565432279500746*^9, {3.565432323212023*^9, 3.5654323487804675`*^9}, {
   3.5654324069217696`*^9, 3.565432448043442*^9}, 3.565435118282593*^9, 
   3.565435188716717*^9, {3.5654352375916023`*^9, 3.565435249868824*^9}, {
   3.56543542959674*^9, 3.5654354355871506`*^9}, 3.565435559201768*^9, 
   3.565435704687623*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"(*", " ", 
  RowBox[{"Etot", " ", "in", " ", "eV"}], " ", "*)"}]], "Input"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"Table", "[", 
   RowBox[{
    RowBox[{"N", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"t", "[", 
        RowBox[{"[", 
         RowBox[{"i", ",", "1"}], "]"}], "]"}], "*", "27.21138"}], ",", 
      "10"}], "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"i", ",", "1", ",", "9"}], "}"}]}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"PA", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"P", " ", "/.", " ", 
      RowBox[{"t", "[", 
       RowBox[{"[", 
        RowBox[{"i", ",", "2"}], "]"}], "]"}]}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "9"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"QA", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"Q", "  ", "/.", " ", 
      RowBox[{"t", "[", 
       RowBox[{"[", 
        RowBox[{"i", ",", "2"}], "]"}], "]"}]}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "9"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"RA", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"R", " ", "/.", " ", 
      RowBox[{"t", "[", 
       RowBox[{"[", 
        RowBox[{"i", ",", "2"}], "]"}], "]"}]}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "9"}], "}"}]}], "]"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.565347493398529*^9, 3.565347525316185*^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", ",", "9"}], "}"}]}], "]"}], 
  "\[IndentingNewLine]"}]], "Input",
 CellChangeTimes->{{3.5654355069104757`*^9, 3.565435545723344*^9}, 
   3.5654356945320053`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   GraphicsBox[{
     {Thickness[0.015], CircleBox[{0, 0}, 0.15398421744914842`], 
      CircleBox[{0, 0}, 1.5885139078383728`]}, 
     {Thickness[0.01], 
      {RGBColor[1, 0, 0], 
       CircleBox[{0, -0.6177358807009392}, 0.6177358807009392], 
       CircleBox[{0, 0.6177358807009392}, 0.6177358807009392]}, 
      {RGBColor[1, 0, 1], 
       CircleBox[{-0.6177358807009392, 0}, 0.6177358807009392], 
       CircleBox[{0.6177358807009392, 0}, 0.6177358807009392]}, 
      {RGBColor[0, 1, 1], CircleBox[{0, 0}, 0.6177358807009392]}}},
    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.13969192675004327`], 
      CircleBox[{0, 0}, 1.2913825954113067`]}, 
     {Thickness[0.01], 
      {RGBColor[1, 0, 0], 
       CircleBox[{0, -0.5327797060592275}, 0.5327797060592275], 
       CircleBox[{0, 0.5327797060592275}, 0.5327797060592275]}, 
      {RGBColor[1, 0, 1], 
       CircleBox[{-0.5327797060592275, 0}, 0.5327797060592275], 
       CircleBox[{0.5327797060592275, 0}, 0.5327797060592275]}, 
      {RGBColor[0, 1, 1], CircleBox[{0, 0}, 0.5327797060592275]}}},
    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.12782468311414308`], 
      CircleBox[{0, 0}, 1.02443385856661]}, 
     {Thickness[0.01], 
      {RGBColor[1, 0, 0], 
       CircleBox[{0, -0.4724387728144514}, 0.4724387728144514], 
       CircleBox[{0, 0.4724387728144514}, 0.4724387728144514]}, 
      {RGBColor[1, 0, 1], 
       CircleBox[{-0.4724387728144514, 0}, 0.4724387728144514], 
       CircleBox[{0.4724387728144514, 0}, 0.4724387728144514]}, 
      {RGBColor[0, 1, 1], CircleBox[{0, 0}, 0.4724387728144514]}}},
    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.11781805218407913`], 
      CircleBox[{0, 0}, 0.74977754076882]}, 
     {Thickness[0.01], 
      {RGBColor[1, 0, 0], 
       CircleBox[{0, -0.4360242734454343}, 0.4360242734454343], 
       CircleBox[{0, 0.4360242734454343}, 0.4360242734454343]}, 
      {RGBColor[1, 0, 1], 
       CircleBox[{-0.4360242734454343, 0}, 0.4360242734454343], 
       CircleBox[{0.4360242734454343, 0}, 0.4360242734454343]}, 
      {RGBColor[0, 1, 1], CircleBox[{0, 0}, 0.4360242734454343]}}},
    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.1092617149883094], 
      CircleBox[{0, 0}, 0.6466370326233891]}, 
     {Thickness[0.01], 
      {RGBColor[1, 0, 0], 
       CircleBox[{0, -0.39616563882769434`}, 0.39616563882769434`], 
       CircleBox[{0, 0.39616563882769434`}, 0.39616563882769434`]}, 
      {RGBColor[1, 0, 1], 
       CircleBox[{-0.39616563882769434`, 0}, 0.39616563882769434`], 
       CircleBox[{0.39616563882769434`, 0}, 0.39616563882769434`]}, 
      {RGBColor[0, 1, 1], CircleBox[{0, 0}, 0.39616563882769434`]}}},
    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.10186175056008497`], 
      CircleBox[{0, 0}, 0.5772654286442558]}, 
     {Thickness[0.01], 
      {RGBColor[1, 0, 0], 
       CircleBox[{0, -0.3614697967847427}, 0.3614697967847427], 
       CircleBox[{0, 0.3614697967847427}, 0.3614697967847427]}, 
      {RGBColor[1, 0, 1], 
       CircleBox[{-0.3614697967847427, 0}, 0.3614697967847427], 
       CircleBox[{0.3614697967847427, 0}, 0.3614697967847427]}, 
      {RGBColor[0, 1, 1], CircleBox[{0, 0}, 0.3614697967847427]}}},
    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.0953990268490955], 
      CircleBox[{0, 0}, 0.5235313645874632]}, 
     {Thickness[0.01], 
      {RGBColor[1, 0, 0], 
       CircleBox[{0, -0.33197443067159177`}, 0.33197443067159177`], 
       CircleBox[{0, 0.33197443067159177`}, 0.33197443067159177`]}, 
      {RGBColor[1, 0, 1], 
       CircleBox[{-0.33197443067159177`, 0}, 0.33197443067159177`], 
       CircleBox[{0.33197443067159177`, 0}, 0.33197443067159177`]}, 
      {RGBColor[0, 1, 1], CircleBox[{0, 0}, 0.33197443067159177`]}}},
    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.0897062911586054], 
      CircleBox[{0, 0}, 0.47979765093440596`]}, 
     {Thickness[0.01], 
      {RGBColor[1, 0, 0], 
       CircleBox[{0, -0.3067791621916184}, 0.3067791621916184], 
       CircleBox[{0, 0.3067791621916184}, 0.3067791621916184]}, 
      {RGBColor[1, 0, 1], 
       CircleBox[{-0.3067791621916184, 0}, 0.3067791621916184], 
       CircleBox[{0.3067791621916184, 0}, 0.3067791621916184]}, 
      {RGBColor[0, 1, 1], CircleBox[{0, 0}, 0.3067791621916184]}}},
    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.08465380331289858], 
      CircleBox[{0, 0}, 0.44320969672752136`]}, 
     {Thickness[0.01], 
      {RGBColor[1, 0, 0], 
       CircleBox[{0, -0.28506814451379575`}, 0.28506814451379575`], 
       CircleBox[{0, 0.28506814451379575`}, 0.28506814451379575`]}, 
      {RGBColor[1, 0, 1], 
       CircleBox[{-0.28506814451379575`, 0}, 0.28506814451379575`], 
       CircleBox[{0.28506814451379575`, 0}, 0.28506814451379575`]}, 
      {RGBColor[0, 1, 1], CircleBox[{0, 0}, 0.28506814451379575`]}}},
    AspectRatio->Automatic,
    Axes->True,
    ColorOutput->Automatic,
    Frame->True,
    GridLines->Automatic,
    PlotRange->{{-1.6, 1.6}, {-1.6, 1.6}}]}], "}"}]], "Output",
 CellChangeTimes->{3.5654355592485676`*^9, 3.565435704734423*^9}]
}, Open  ]]
},
WindowSize->{761, 815},
WindowMargins->{{Automatic, 314}, {Automatic, -47}},
PrintingCopies->1,
PrintingPageRange->{32000, 32000},
PrintingOptions->{"Magnification"->1.,
"PaperOrientation"->"Portrait",
"PaperSize"->{595.2, 841.8}},
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, 241, 3, 51, "Text"],
Cell[CellGroupData[{
Cell[823, 27, 3763, 120, 212, "Input"],
Cell[4589, 149, 3268, 83, 192, "Output"]
}, Open  ]],
Cell[7872, 235, 93, 2, 31, "Input"],
Cell[CellGroupData[{
Cell[7990, 241, 1236, 41, 72, "Input"],
Cell[9229, 284, 816, 17, 31, "Output"],
Cell[10048, 303, 810, 17, 52, "Output"],
Cell[10861, 322, 653, 12, 31, "Output"]
}, Open  ]],
Cell[11529, 337, 99, 2, 31, "Input"],
Cell[11631, 341, 1377, 46, 92, "Input"],
Cell[CellGroupData[{
Cell[13033, 391, 4166, 106, 252, "Input"],
Cell[17202, 499, 6477, 156, 597, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)