(* 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[     20803,        590]
NotebookOptionsPosition[     19859,        557]
NotebookOutlinePosition[     20282,        573]
CellTagsIndexPosition[     20239,        570]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell["\<\
Ne-Atom and homologues 1s(2)2s(2) 2p(6) , all p equal size, with exchange \
correction page 40; 24.12.2012\
\>", "Text",
 CellChangeTimes->{{3.565285708771845*^9, 3.565285709801447*^9}, {
  3.5652857568979297`*^9, 3.5652857982848024`*^9}, {3.565285856628905*^9, 
  3.5652858574401064`*^9}, {3.565285894817772*^9, 3.56528589956018*^9}, {
  3.5653268341468997`*^9, 3.5653268721017666`*^9}, {3.565326947169098*^9, 
  3.5653269759043484`*^9}, {3.565332760327036*^9, 3.5653327614658375`*^9}, {
  3.5653331650385466`*^9, 3.5653331986254053`*^9}, {3.565333237984275*^9, 
  3.5653333056259937`*^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[{"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"}], "}"}]}], 
   "]"}]}], "\[IndentingNewLine]"}], "Input",
 CellChangeTimes->{{3.565285708771845*^9, 3.565285709801447*^9}, {
  3.5652857568979297`*^9, 3.5652857982848024`*^9}, {3.565285856628905*^9, 
  3.5652858574401064`*^9}, {3.565285894817772*^9, 3.56528589956018*^9}, {
  3.5653268341468997`*^9, 3.5653268721017666`*^9}, {3.565326947169098*^9, 
  3.5653269759043484`*^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.5653345488401775`*^9, 3.5653345786362295`*^9}}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"(*", " ", 
  RowBox[{"Virial", " ", "ratio"}], " ", "*)"}]], "Input",
 CellChangeTimes->{{3.565285708771845*^9, 3.565285709801447*^9}, {
  3.5652857568979297`*^9, 3.5652857982848024`*^9}, {3.565285856628905*^9, 
  3.5652858574401064`*^9}, {3.565285894817772*^9, 3.56528589956018*^9}, {
  3.5653268341468997`*^9, 3.5653268721017666`*^9}, {3.565326947169098*^9, 
  3.5653269849523644`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"Table", "[", 
  RowBox[{
   RowBox[{"N", "[", 
    RowBox[{
     RowBox[{"Vir", " ", "/.", " ", 
      RowBox[{"t", "[", 
       RowBox[{"[", 
        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[{"{", 
          RowBox[{"z", "-", "9"}], "}"}], ",", "2"}], "]"}], "]"}]}], ",", 
     "10"}], "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"z", " ", ",", "10", ",", "18"}], "}"}]}], 
  "]"}], "\[IndentingNewLine]", 
 RowBox[{"Table", "[", 
  RowBox[{
   RowBox[{"N", "[", 
    RowBox[{
     RowBox[{"Vee", " ", "/.", " ", 
      RowBox[{"t", "[", 
       RowBox[{"[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{"z", "-", "9"}], "}"}], ",", "2"}], "]"}], "]"}]}], ",", 
     "10"}], "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"z", " ", ",", "10", ",", "18"}], "}"}]}], 
  "]"}], "\[IndentingNewLine]"}], "Input",
 CellChangeTimes->{{3.565285708771845*^9, 3.565285709801447*^9}, {
  3.5652857568979297`*^9, 3.5652857982848024`*^9}, {3.565285856628905*^9, 
  3.5652858574401064`*^9}, {3.565285894817772*^9, 3.56528589956018*^9}, {
  3.5653268341468997`*^9, 3.5653268721017666`*^9}, {3.565326947169098*^9, 
  3.5653269849523644`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"-", "1.9999999999999964`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1.9999999999630884`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1.9999999999999936`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1.9999999954166932`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1.9999999999806777`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1.9999999999904157`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1.9999999999998836`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1.9999999999964242`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1.9999999999988225`"}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{{3.5653345488557773`*^9, 3.5653345786518297`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"-", "310.83962654506786`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "385.6667237935207`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "469.17780773429934`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "561.9234344188847`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "661.3820931128294`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "768.7124464586709`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "884.0137357751692`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1007.303537477298`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1138.5874621119622`"}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{{3.5653345488557773`*^9, 3.5653345786518297`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", "51.02863534087648`", "}"}], ",", 
   RowBox[{"{", "59.71714466347031`", "}"}], ",", 
   RowBox[{"{", "69.08642149012044`", "}"}], ",", 
   RowBox[{"{", "79.48454657456786`", "}"}], ",", 
   RowBox[{"{", "88.4728987692344`", "}"}], ",", 
   RowBox[{"{", "97.31349415602482`", "}"}], ",", 
   RowBox[{"{", "106.11617929625945`", "}"}], ",", 
   RowBox[{"{", "114.90184869317312`", "}"}], ",", 
   RowBox[{"{", "123.67769075767262`", "}"}]}], "}"}]], "Output",
 CellChangeTimes->{{3.5653345488557773`*^9, 3.5653345786518297`*^9}}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"(*", " ", 
  RowBox[{"Etot", " ", "in", " ", "eV"}], " ", "*)"}]], "Input",
 CellChangeTimes->{{3.565285708771845*^9, 3.565285709801447*^9}, {
  3.5652857568979297`*^9, 3.5652857982848024`*^9}, {3.565285856628905*^9, 
  3.5652858574401064`*^9}, {3.565285894817772*^9, 3.56528589956018*^9}, {
  3.5653268341468997`*^9, 3.5653268721017666`*^9}, {3.565326947169098*^9, 
  3.565326993033179*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"Table", "[", 
   RowBox[{
    RowBox[{"N", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"t", "[", 
        RowBox[{"[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{"z", "-", "9"}], "}"}], ",", "1"}], "]"}], "]"}], "*", 
       "27.21138"}], ",", "10"}], "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"z", ",", "10", ",", "18"}], "}"}]}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"PA", "=", 
  RowBox[{"Table", "[", 
   RowBox[{
    RowBox[{"P", " ", "/.", " ", 
     RowBox[{"t", "[", 
      RowBox[{"[", 
       RowBox[{"i", ",", "2"}], "]"}], "]"}]}], ",", 
    RowBox[{"{", 
     RowBox[{"i", ",", "1", ",", "9"}], "}"}]}], 
   "]"}]}], "\[IndentingNewLine]", 
 RowBox[{"QA", "=", 
  RowBox[{"Table", "[", 
   RowBox[{
    RowBox[{"Q", "  ", "/.", " ", 
     RowBox[{"t", "[", 
      RowBox[{"[", 
       RowBox[{"i", ",", "2"}], "]"}], "]"}]}], ",", 
    RowBox[{"{", 
     RowBox[{"i", ",", "1", ",", "9"}], "}"}]}], 
   "]"}]}], "\[IndentingNewLine]", 
 RowBox[{"RA", "=", 
  RowBox[{"Table", "[", 
   RowBox[{
    RowBox[{"R", " ", "/.", " ", 
     RowBox[{"t", "[", 
      RowBox[{"[", 
       RowBox[{"i", ",", "2"}], "]"}], "]"}]}], ",", 
    RowBox[{"{", 
     RowBox[{"i", ",", "1", ",", "9"}], "}"}]}], "]"}]}]}], "Input",
 CellChangeTimes->{{3.565285708771845*^9, 3.565285709801447*^9}, {
   3.5652857568979297`*^9, 3.5652857982848024`*^9}, {3.565285856628905*^9, 
   3.5652858574401064`*^9}, {3.565285894817772*^9, 3.56528589956018*^9}, {
   3.5653268341468997`*^9, 3.5653268721017666`*^9}, {3.565326947169098*^9, 
   3.565326993033179*^9}, {3.5653270417364645`*^9, 3.565327146755849*^9}, {
   3.565327396200287*^9, 3.5653274410503654`*^9}, {3.5653280284537973`*^9, 
   3.565328048905433*^9}, {3.565328087905502*^9, 3.565328089200304*^9}, {
   3.56532818927448*^9, 3.5653282015673018`*^9}, {3.56532824045817*^9, 
   3.5653283594083786`*^9}, {3.5653288706680765`*^9, 3.5653289157677555`*^9}, 
   3.565332684635703*^9, {3.5653327310301847`*^9, 3.565332731560585*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "0.15398421744914842`", ",", "0.13969192675004327`", ",", 
   "0.12782468311414308`", ",", "0.11781805218407913`", ",", 
   "0.1092617149883094`", ",", "0.10186175056008497`", ",", 
   "0.0953990268490955`", ",", "0.0897062911586054`", ",", 
   "0.08465380331289858`"}], "}"}]], "Output",
 CellChangeTimes->{{3.565334548886977*^9, 3.56533457866743*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "1.5885139078383728`", ",", "1.2913825954113067`", ",", "1.02443385856661`",
    ",", "0.74977754076882`", ",", "0.6466370326233891`", ",", 
   "0.5772654286442558`", ",", "0.5235313645874632`", ",", 
   "0.47979765093440596`", ",", "0.44320969672752136`"}], "}"}]], "Output",
 CellChangeTimes->{{3.565334548886977*^9, 3.56533457866743*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "0.6177358807009392`", ",", "0.5327797060592275`", ",", 
   "0.4724387728144514`", ",", "0.4360242734454343`", ",", 
   "0.39616563882769434`", ",", "0.3614697967847427`", ",", 
   "0.33197443067159177`", ",", "0.3067791621916184`", ",", 
   "0.28506814451379575`"}], "}"}]], "Output",
 CellChangeTimes->{{3.565334548886977*^9, 3.56533457866743*^9}}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"Do", "[", 
  RowBox[{
   RowBox[{"Show", "[", 
    RowBox[{
     RowBox[{"Graphics", "[", 
      RowBox[{"{", 
       RowBox[{
        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[{"Dashing", "[", 
           RowBox[{"{", 
            RowBox[{"0.02", ",", " ", "0.02"}], "}"}], "]"}], ",", " ", 
          RowBox[{"Circle", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{"0", ",", " ", 
              RowBox[{"-", 
               RowBox[{"RA", "[", 
                RowBox[{"[", "i", "]"}], "]"}]}]}], "}"}], ",", " ", 
            RowBox[{"RA", "[", 
             RowBox[{"[", "i", "]"}], "]"}]}], "]"}], ",", " ", "\n", 
          "      ", 
          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[{"Circle", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{"0", ",", " ", 
              RowBox[{"RA", "[", 
               RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", " ", 
            RowBox[{"RA", "[", 
             RowBox[{"[", "i", "]"}], "]"}]}], "]"}], ",", " ", "\n", 
          "      ", 
          RowBox[{"Circle", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{"0", ",", " ", "0"}], "}"}], ",", " ", 
            RowBox[{"RA", "[", 
             RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}], "}"}]}], "}"}], "]"}],
      ",", " ", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"AspectRatio", " ", "->", " ", "Automatic"}], ",", " ", "\n", 
       "    ", 
       RowBox[{"Axes", " ", "->", " ", "True"}], ",", " ", 
       RowBox[{"GridLines", " ", "->", " ", "Automatic"}], ",", " ", "\n", 
       "    ", 
       RowBox[{"PlotRange", " ", "->", " ", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"-", "2"}], ",", " ", "2"}], "}"}], ",", " ", 
          RowBox[{"{", 
           RowBox[{
            RowBox[{"-", "2"}], ",", " ", "2"}], "}"}]}], "}"}]}], ",", " ", 
       RowBox[{"Frame", " ", "->", " ", "True"}]}], "}"}]}], "]"}], ",", " ", 
   RowBox[{"{", 
    RowBox[{"i", ",", " ", "1", ",", " ", "9"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.565285708771845*^9, 3.565285709801447*^9}, {
   3.5652857568979297`*^9, 3.5652857982848024`*^9}, {3.565285856628905*^9, 
   3.5652858574401064`*^9}, 3.565285894817772*^9, {3.5653271780807037`*^9, 
   3.565327253366436*^9}, 3.565327613087468*^9, {3.565327667344363*^9, 
   3.565327878256734*^9}, {3.565328125236367*^9, 3.56532813834039*^9}, {
   3.565328384742823*^9, 3.565328389516432*^9}, {3.5653284227912903`*^9, 
   3.565328433602109*^9}, {3.565328951975419*^9, 3.565329081408847*^9}, 
   3.5653326246691976`*^9},
 FormatType->"InputForm"]
},
WindowSize->{663, 815},
WindowMargins->{{Automatic, 431}, {Automatic, 5}},
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, 601, 10, 49, "Text"],
Cell[CellGroupData[{
Cell[1183, 34, 4044, 123, 252, "Input"],
Cell[5230, 159, 2988, 78, 192, "Output"]
}, Open  ]],
Cell[8233, 240, 412, 7, 31, "Input"],
Cell[CellGroupData[{
Cell[8670, 251, 1546, 47, 92, "Input"],
Cell[10219, 300, 771, 21, 31, "Output"],
Cell[10993, 323, 765, 21, 52, "Output"],
Cell[11761, 346, 584, 12, 52, "Output"]
}, Open  ]],
Cell[12360, 361, 416, 7, 31, "Input"],
Cell[CellGroupData[{
Cell[12801, 372, 2046, 54, 92, "Input"],
Cell[14850, 428, 395, 8, 52, "Output"],
Cell[15248, 438, 383, 7, 52, "Output"],
Cell[15634, 447, 392, 8, 52, "Output"]
}, Open  ]],
Cell[16041, 458, 3814, 97, 134, "Input"]
}
]
*)

(* End of internal cache information *)