(* 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[     30917,        826]
NotebookOptionsPosition[     30159,        798]
NotebookOutlinePosition[     30582,        814]
CellTagsIndexPosition[     30539,        811]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"O", "-", 
     RowBox[{"Atom", " ", "and", " ", "homologues", " ", "to", " ", 
      RowBox[{"Ar", "[", 
       RowBox[{"+", "10"}], "]"}], " ", "1", "s", 
      RowBox[{"(", "2", ")"}], "2", "s", 
      RowBox[{"(", "2", ")"}], " ", "2", "p", 
      RowBox[{"(", "4", ")"}]}]}], ",", " ", 
    RowBox[{"p", 
     RowBox[{"(", "1", ")"}]}], ",", " ", 
    RowBox[{"p", 
     RowBox[{"(", "2", ")"}], " ", "different", " ", "size"}], ",", 
    RowBox[{
     RowBox[{
     "with", " ", "exchange", " ", "correction", " ", "page", " ", "40"}], 
     ";", " ", 
     RowBox[{"11.12", 
      RowBox[{".2011", "/", "28.04"}], ".2013"}]}]}], " ", "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"Clear", "[", 
     RowBox[{"z", ",", "sig"}], "]"}], ";", " ", 
    RowBox[{"sig", "=", "0.29"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"T", "=", 
     RowBox[{
      RowBox[{"0.976", "*", 
       RowBox[{"9", "/", 
        RowBox[{"(", 
         RowBox[{"4", "*", 
          RowBox[{"P", "^", "2"}]}], ")"}]}]}], "+", 
      RowBox[{"9", "*", 
       RowBox[{"4", "/", 
        RowBox[{"(", 
         RowBox[{"4", "*", 
          RowBox[{"Q", "^", "2"}]}], ")"}]}]}], "+", 
      RowBox[{"4", "/", 
       RowBox[{"R", "^", "2"}]}], "+", 
      RowBox[{"4", "/", 
       RowBox[{"S", "^", "2"}]}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Vne", "=", 
     RowBox[{
      RowBox[{
       RowBox[{"-", "3"}], "*", 
       RowBox[{"z", "/", "P"}]}], "-", 
      RowBox[{"3", "*", 
       RowBox[{"z", "/", "Q"}]}], "-", 
      RowBox[{"2", "*", 
       RowBox[{"z", "/", "R"}]}], "-", 
      RowBox[{"2", "*", 
       RowBox[{"z", "/", "S"}]}]}]}], ";"}], "\[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[{"1.5", "*", 
       RowBox[{"sig", "/", "S"}]}], "+", 
      RowBox[{"4", "*", 
       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", "*", 
       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[{"4", "*", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"1", "-", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{"3", "/", "20"}], ")"}], "*", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{"P", "/", "S"}], ")"}], "^", "2"}]}], "+", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{"1", "/", "32"}], ")"}], "*", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{"P", "/", "S"}], ")"}], "^", "3"}]}]}], ")"}], "/", 
        "S"}]}], "+", 
      RowBox[{"4", "*", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"1", "-", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{"3", "/", "20"}], ")"}], "*", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{"Q", "/", "S"}], ")"}], "^", "2"}]}], "+", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{"1", "/", "32"}], ")"}], "*", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{"Q", "/", "S"}], ")"}], "^", "3"}]}]}], ")"}], "/", 
        "S"}]}], "+", 
      RowBox[{"1.517767", "/", "R"}], "+", 
      RowBox[{"1.517767", "/", "S"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"func", "=", 
     RowBox[{"T", "+", "Vne", "+", "Vee"}]}], ";", " ", 
    RowBox[{"Vir", " ", "=", " ", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"Vne", "+", "Vee"}], ")"}], "/", "T"}]}], ";"}]}]}]], "Input",
 CellChangeTimes->{{3.5761312847167473`*^9, 3.576131295589967*^9}, {
   3.5761313718429003`*^9, 3.576131380095315*^9}, {3.5761314680326695`*^9, 
   3.576131473414679*^9}, {3.576131505082735*^9, 3.576131508857941*^9}, {
   3.576131550291614*^9, 3.57613157653086*^9}, {3.576131668602222*^9, 
   3.576131691456262*^9}, {3.5761317303939304`*^9, 3.57613173585394*^9}, {
   3.576131887221006*^9, 3.576131897142623*^9}, {3.576131931337883*^9, 
   3.576131942975504*^9}, {3.5761319922559905`*^9, 3.5761320081836185`*^9}, {
   3.576132261559263*^9, 3.576132328686181*^9}, 3.576132365736246*^9, {
   3.5761324620820155`*^9, 3.576132462924417*^9}, {3.5761324940464716`*^9, 
   3.5761324980712786`*^9}, {3.576132621748296*^9, 3.5761326404839287`*^9}, {
   3.5761341272913404`*^9, 3.5761341524073844`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"t", " ", "=", 
  RowBox[{"Table", "[", 
   RowBox[{
    RowBox[{"FindMinimum", "[", 
     RowBox[{"func", ",", 
      RowBox[{"{", 
       RowBox[{"P", ",", "0.1", ",", "0.2"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"Q", ",", "0.4", ",", "2.0"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"R", ",", "0.2", ",", "0.8"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"S", ",", "0.2", ",", "0.84"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"Method", "\[Rule]", "Automatic"}], ",", 
        RowBox[{"MaxIterations", "\[Rule]", "800"}]}], "}"}]}], "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"z", ",", "8", ",", "18"}], "}"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.5761312847167473`*^9, 3.576131295589967*^9}, {
   3.5761313718429003`*^9, 3.576131380095315*^9}, {3.5761314680326695`*^9, 
   3.576131473414679*^9}, {3.576131505082735*^9, 3.576131508857941*^9}, {
   3.576131550291614*^9, 3.57613157653086*^9}, {3.576131668602222*^9, 
   3.576131691456262*^9}, {3.5761317303939304`*^9, 3.57613173585394*^9}, {
   3.576131887221006*^9, 3.576131897142623*^9}, {3.576131931337883*^9, 
   3.576131942975504*^9}, {3.5761319922559905`*^9, 3.5761320081836185`*^9}, {
   3.576132261559263*^9, 3.576132328686181*^9}, {3.576132365736246*^9, 
   3.5761323988395042`*^9}, 3.5761324339863663`*^9, {3.57613253287494*^9, 
   3.576132579503422*^9}, {3.5761326607327642`*^9, 3.5761327003724337`*^9}, {
   3.5761327906341925`*^9, 3.5761328172010393`*^9}, {3.5761328674331274`*^9, 
   3.5761330469114428`*^9}, {3.576133718196222*^9, 3.576133726261436*^9}, {
   3.5761337715015154`*^9, 3.576133793513154*^9}, {3.576133849735653*^9, 
   3.5761338744616966`*^9}, {3.5761339511982307`*^9, 3.576133954536637*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "74.97919781546207`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.19016103313136445`"}], ",", 
       RowBox[{"Q", "\[Rule]", "2.00613980390547`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.7991308946707711`"}], ",", 
       RowBox[{"S", "\[Rule]", "0.8340156998994251`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "98.94804940780332`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.16837331852380474`"}], ",", 
       RowBox[{"Q", "\[Rule]", "1.1392028956639269`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.7210223669522138`"}], ",", 
       RowBox[{"S", "\[Rule]", "0.7637319360584117`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "126.56856801396133`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.15105639582380012`"}], ",", 
       RowBox[{"Q", "\[Rule]", "0.9479339217745721`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.6136974331817198`"}], ",", 
       RowBox[{"S", "\[Rule]", "0.6424140485901305`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "157.74184884346596`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.13696446588749034`"}], ",", 
       RowBox[{"Q", "\[Rule]", "0.8164429820582584`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.5325572343612969`"}], ",", 
       RowBox[{"S", "\[Rule]", "0.5531087429559473`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "192.46694807302367`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.12527464775687008`"}], ",", 
       RowBox[{"Q", "\[Rule]", "0.7176633703394801`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.470135435458793`"}], ",", 
       RowBox[{"S", "\[Rule]", "0.4855453221796489`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "230.74327835280695`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.11542100141690369`"}], ",", 
       RowBox[{"Q", "\[Rule]", "0.6404282593247813`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.4207340118584999`"}], ",", 
       RowBox[{"S", "\[Rule]", "0.4327267621998089`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "272.570419387621`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.10700243225641431`"}], ",", 
       RowBox[{"Q", "\[Rule]", "0.5782831713533064`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.38069328627820354`"}], ",", 
       RowBox[{"S", "\[Rule]", "0.3902878194986044`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "317.94805897007274`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.09972707193611059`"}], ",", 
       RowBox[{"Q", "\[Rule]", "0.5271957960864525`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.34760331335404854`"}], ",", 
       RowBox[{"S", "\[Rule]", "0.3554611120253503`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "366.8759595370991`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.09337752959951504`"}], ",", 
       RowBox[{"Q", "\[Rule]", "0.4844186272782812`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.31979629899225614`"}], ",", 
       RowBox[{"S", "\[Rule]", "0.32632666697388263`"}]}], "}"}]}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "419.35393789179125`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.08778632759816841`"}], ",", 
       RowBox[{"Q", "\[Rule]", "0.4481117913167129`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.29611015926291645`"}], ",", 
       RowBox[{"S", "\[Rule]", "0.3016501513979271`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "475.3818488937672`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.08282655232738126`"}], ",", 
       RowBox[{"Q", "\[Rule]", "0.416866054429945`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.2756977039977251`"}], ",", 
       RowBox[{"S", "\[Rule]", "0.2804375478878805`"}]}], "}"}]}], "}"}]}], 
  "}"}]], "Output",
 CellChangeTimes->{
  3.576132593917847*^9, {3.576132647893942*^9, 3.5761326770815935`*^9}, 
   3.5761327075172462`*^9, 3.5761328243146515`*^9, {3.5761328715359344`*^9, 
   3.5761330515602508`*^9}, 3.5761337322518463`*^9, 3.576133801063567*^9, 
   3.576133879750106*^9, 3.576133960854648*^9, 3.576134301418846*^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", "-", "7"}], "}"}], ",", "2"}], "]"}], "]"}]}], ",", 
       "10"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"z", " ", ",", "8", ",", "18"}], "}"}]}], "]"}], 
   "\[IndentingNewLine]", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"N", "[", 
      RowBox[{
       RowBox[{"Vne", " ", "/.", " ", 
        RowBox[{"t", "[", 
         RowBox[{"[", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{"z", "-", "7"}], "}"}], ",", "2"}], "]"}], "]"}]}], ",", 
       "10"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"z", " ", ",", "8", ",", "18"}], "}"}]}], "]"}], 
   "\[IndentingNewLine]", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"N", "[", 
      RowBox[{
       RowBox[{"Vee", " ", "/.", " ", 
        RowBox[{"t", "[", 
         RowBox[{"[", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{"z", "-", "7"}], "}"}], ",", "2"}], "]"}], "]"}]}], ",", 
       "10"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"z", " ", ",", "8", ",", "18"}], "}"}]}], "]"}], 
   "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{"Etot", " ", "in", " ", "eV"}], " ", "*)"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"N", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"t", "[", 
          RowBox[{"[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{"z", "-", "7"}], "}"}], ",", "1"}], "]"}], "]"}], "*", 
         "27.21138"}], ",", "10"}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"z", ",", "8", ",", "18"}], "}"}]}], "]"}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"PA", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"t", "[", 
        RowBox[{"[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{"z", "-", "7"}], "}"}], ",", "2", ",", "1"}], "]"}], "]"}],
        ",", 
       RowBox[{"{", 
        RowBox[{"z", ",", "8", ",", "18"}], "}"}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"QA", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"t", "[", 
        RowBox[{"[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{"z", "-", "7"}], "}"}], ",", "2", ",", "2"}], "]"}], "]"}],
        ",", 
       RowBox[{"{", 
        RowBox[{"z", ",", "8", ",", "18"}], "}"}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"RA", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"t", "[", 
        RowBox[{"[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{"z", "-", "7"}], "}"}], ",", "2", ",", "3"}], "]"}], "]"}],
        ",", 
       RowBox[{"{", 
        RowBox[{"z", ",", "8", ",", "18"}], "}"}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"SA", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"t", "[", 
        RowBox[{"[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{"z", "-", "7"}], "}"}], ",", "2", ",", "4"}], "]"}], "]"}],
        ",", 
       RowBox[{"{", 
        RowBox[{"z", ",", "8", ",", "18"}], "}"}]}], "]"}]}], 
    ";"}]}]}]], "Input",
 CellChangeTimes->{{3.5761312847167473`*^9, 3.576131295589967*^9}, {
   3.5761313718429003`*^9, 3.576131380095315*^9}, {3.5761314680326695`*^9, 
   3.576131473414679*^9}, {3.576131505082735*^9, 3.576131508857941*^9}, {
   3.576131550291614*^9, 3.57613157653086*^9}, {3.576131668602222*^9, 
   3.576131691456262*^9}, {3.5761317303939304`*^9, 3.57613173585394*^9}, {
   3.576131887221006*^9, 3.576131897142623*^9}, {3.576131931337883*^9, 
   3.576131942975504*^9}, {3.5761319922559905`*^9, 3.5761320081836185`*^9}, {
   3.576132261559263*^9, 3.576132328686181*^9}, 3.576132365736246*^9, {
   3.5761342884708233`*^9, 3.576134294882435*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"-", "2.0000096890553416`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1.999997875994375`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "2.0000004828377635`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1.9999967435699681`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "2.000000228722352`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "2.0000076117523142`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1.9999974101664466`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "2.000001498399218`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1.9999857074852636`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1.999996028701908`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1.9999987114177247`"}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.5761343014344463`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"-", "177.37813841843325`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "232.59176166198563`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "293.97101692472427`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "362.4429336442146`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "438.00947194947986`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "520.6710648668254`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "610.4351694148108`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "707.2915131936256`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "811.2550929731952`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "922.3024430683325`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1040.4514813247192`"}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.5761343014344463`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", "27.4204692580674`", "}"}], ",", 
   RowBox[{"{", "34.69545267971906`", "}"}], ",", 
   RowBox[{"{", "40.833942008856496`", "}"}], ",", 
   RowBox[{"{", "46.95872228031603`", "}"}], ",", 
   RowBox[{"{", "53.075619824915506`", "}"}], ",", 
   RowBox[{"{", "59.18626450852584`", "}"}], ",", 
   RowBox[{"{", "65.29362472572296`", "}"}], ",", 
   RowBox[{"{", "71.39587166588933`", "}"}], ",", 
   RowBox[{"{", "77.49793024399403`", "}"}], ",", 
   RowBox[{"{", "83.59290189864282`", "}"}], ",", 
   RowBox[{"{", "89.68717096777105`", "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.5761343014344463`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"Table", "[", 
   RowBox[{
    RowBox[{"Show", "[", 
     RowBox[{
      RowBox[{"Graphics", "[", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"Circle", "[", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{"0", ",", "0"}], "}"}], ",", 
           RowBox[{"P", "/.", " ", 
            RowBox[{"PA", "[", 
             RowBox[{"[", "i", "]"}], "]"}]}]}], "]"}], ",", 
         RowBox[{"Circle", "[", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{"0", ",", "0"}], "}"}], ",", 
           RowBox[{"Q", " ", "/.", " ", 
            RowBox[{"QA", "[", 
             RowBox[{"[", "i", "]"}], "]"}]}]}], "]"}], ",", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"Dashing", "[", 
            RowBox[{"{", 
             RowBox[{"0.01", ",", "0.01"}], "}"}], "]"}], ",", 
           RowBox[{"Circle", "[", 
            RowBox[{
             RowBox[{"{", 
              RowBox[{"0", ",", 
               RowBox[{
                RowBox[{"-", "R"}], "/.", " ", 
                RowBox[{"RA", "[", 
                 RowBox[{"[", "i", "]"}], "]"}]}]}], "}"}], ",", 
             RowBox[{"R", " ", "/.", " ", 
              RowBox[{"RA", "[", 
               RowBox[{"[", "i", "]"}], "]"}]}]}], "]"}], ",", 
           RowBox[{"Circle", "[", 
            RowBox[{
             RowBox[{"{", 
              RowBox[{
               RowBox[{
                RowBox[{"-", "R"}], "/.", " ", 
                RowBox[{"RA", "[", 
                 RowBox[{"[", "i", "]"}], "]"}]}], ",", "0"}], "}"}], ",", 
             RowBox[{"R", " ", "/.", " ", 
              RowBox[{"RA", "[", 
               RowBox[{"[", "i", "]"}], "]"}]}]}], "]"}], ",", 
           RowBox[{"Circle", "[", 
            RowBox[{
             RowBox[{"{", 
              RowBox[{
               RowBox[{"R", "/.", " ", 
                RowBox[{"RA", "[", 
                 RowBox[{"[", "i", "]"}], "]"}]}], ",", "0"}], "}"}], ",", 
             RowBox[{"R", " ", "/.", " ", 
              RowBox[{"RA", "[", 
               RowBox[{"[", "i", "]"}], "]"}]}]}], "]"}], ",", 
           RowBox[{"Circle", "[", 
            RowBox[{
             RowBox[{"{", 
              RowBox[{"0", ",", 
               RowBox[{"R", "/.", " ", 
                RowBox[{"RA", "[", 
                 RowBox[{"[", "i", "]"}], "]"}]}]}], "}"}], ",", 
             RowBox[{"R", " ", "/.", " ", 
              RowBox[{"RA", "[", 
               RowBox[{"[", "i", "]"}], "]"}]}]}], "]"}], ",", 
           RowBox[{"Circle", "[", 
            RowBox[{
             RowBox[{"{", 
              RowBox[{"0", ",", "0"}], "}"}], ",", 
             RowBox[{"S", " ", "/.", " ", 
              RowBox[{"SA", "[", 
               RowBox[{"[", "i", "]"}], "]"}]}]}], "]"}]}], "}"}]}], "}"}], 
       "]"}], ",", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"AspectRatio", " ", "\[Rule]", " ", "Automatic"}], ",", 
        RowBox[{"Axes", " ", "->", " ", "True"}], ",", 
        RowBox[{"GridLines", " ", "->", " ", "Automatic"}], ",", " ", 
        RowBox[{"PlotRange", " ", "\[Rule]", " ", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{
             RowBox[{"-", "2"}], ",", "2"}], "}"}], ",", 
           RowBox[{"{", 
            RowBox[{
             RowBox[{"-", "2"}], ",", "2"}], "}"}]}], "}"}]}], ",", " ", 
        RowBox[{"Frame", " ", "->", " ", "True"}]}], "}"}]}], "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"i", ",", "1", ",", "11"}], "}"}]}], "]"}], 
  "\[IndentingNewLine]"}]], "Input",
 CellChangeTimes->{{3.5761312847167473`*^9, 3.576131295589967*^9}, {
   3.5761313718429003`*^9, 3.576131380095315*^9}, {3.5761314680326695`*^9, 
   3.576131473414679*^9}, {3.576131505082735*^9, 3.576131508857941*^9}, {
   3.576131550291614*^9, 3.57613157653086*^9}, {3.576131668602222*^9, 
   3.576131691456262*^9}, {3.5761317303939304`*^9, 3.57613173585394*^9}, {
   3.576131887221006*^9, 3.576131897142623*^9}, {3.576131931337883*^9, 
   3.576131942975504*^9}, {3.5761319922559905`*^9, 3.5761320081836185`*^9}, {
   3.576132261559263*^9, 3.576132328686181*^9}, 3.576132365736246*^9, 
   3.5761342884708233`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   GraphicsBox[{CircleBox[{0, 0}, 0.19016103313136445`], 
     CircleBox[{0, 0}, 2.00613980390547], 
     {Dashing[{0.01, 0.01}], 
      CircleBox[{0, -0.7991308946707711}, 0.7991308946707711], 
      CircleBox[{-0.7991308946707711, 0}, 0.7991308946707711], 
      CircleBox[{0.7991308946707711, 0}, 0.7991308946707711], 
      CircleBox[{0, 0.7991308946707711}, 0.7991308946707711], 
      CircleBox[{0, 0}, 0.8340156998994251]}},
    AspectRatio->Automatic,
    Axes->True,
    Frame->True,
    GridLines->Automatic,
    PlotRange->{{-2, 2}, {-2, 2}}], ",", 
   GraphicsBox[{CircleBox[{0, 0}, 0.16837331852380474`], 
     CircleBox[{0, 0}, 1.1392028956639269`], 
     {Dashing[{0.01, 0.01}], 
      CircleBox[{0, -0.7210223669522138}, 0.7210223669522138], 
      CircleBox[{-0.7210223669522138, 0}, 0.7210223669522138], 
      CircleBox[{0.7210223669522138, 0}, 0.7210223669522138], 
      CircleBox[{0, 0.7210223669522138}, 0.7210223669522138], 
      CircleBox[{0, 0}, 0.7637319360584117]}},
    AspectRatio->Automatic,
    Axes->True,
    Frame->True,
    GridLines->Automatic,
    PlotRange->{{-2, 2}, {-2, 2}}], ",", 
   GraphicsBox[{CircleBox[{0, 0}, 0.15105639582380012`], 
     CircleBox[{0, 0}, 0.9479339217745721], 
     {Dashing[{0.01, 0.01}], 
      CircleBox[{0, -0.6136974331817198}, 0.6136974331817198], 
      CircleBox[{-0.6136974331817198, 0}, 0.6136974331817198], 
      CircleBox[{0.6136974331817198, 0}, 0.6136974331817198], 
      CircleBox[{0, 0.6136974331817198}, 0.6136974331817198], 
      CircleBox[{0, 0}, 0.6424140485901305]}},
    AspectRatio->Automatic,
    Axes->True,
    Frame->True,
    GridLines->Automatic,
    PlotRange->{{-2, 2}, {-2, 2}}], ",", 
   GraphicsBox[{CircleBox[{0, 0}, 0.13696446588749034`], 
     CircleBox[{0, 0}, 0.8164429820582584], 
     {Dashing[{0.01, 0.01}], 
      CircleBox[{0, -0.5325572343612969}, 0.5325572343612969], 
      CircleBox[{-0.5325572343612969, 0}, 0.5325572343612969], 
      CircleBox[{0.5325572343612969, 0}, 0.5325572343612969], 
      CircleBox[{0, 0.5325572343612969}, 0.5325572343612969], 
      CircleBox[{0, 0}, 0.5531087429559473]}},
    AspectRatio->Automatic,
    Axes->True,
    Frame->True,
    GridLines->Automatic,
    PlotRange->{{-2, 2}, {-2, 2}}], ",", 
   GraphicsBox[{CircleBox[{0, 0}, 0.12527464775687008`], 
     CircleBox[{0, 0}, 0.7176633703394801], 
     {Dashing[{0.01, 0.01}], 
      CircleBox[{0, -0.470135435458793}, 0.470135435458793], 
      CircleBox[{-0.470135435458793, 0}, 0.470135435458793], 
      CircleBox[{0.470135435458793, 0}, 0.470135435458793], 
      CircleBox[{0, 0.470135435458793}, 0.470135435458793], 
      CircleBox[{0, 0}, 0.4855453221796489]}},
    AspectRatio->Automatic,
    Axes->True,
    Frame->True,
    GridLines->Automatic,
    PlotRange->{{-2, 2}, {-2, 2}}], ",", 
   GraphicsBox[{CircleBox[{0, 0}, 0.11542100141690369`], 
     CircleBox[{0, 0}, 0.6404282593247813], 
     {Dashing[{0.01, 0.01}], 
      CircleBox[{0, -0.4207340118584999}, 0.4207340118584999], 
      CircleBox[{-0.4207340118584999, 0}, 0.4207340118584999], 
      CircleBox[{0.4207340118584999, 0}, 0.4207340118584999], 
      CircleBox[{0, 0.4207340118584999}, 0.4207340118584999], 
      CircleBox[{0, 0}, 0.4327267621998089]}},
    AspectRatio->Automatic,
    Axes->True,
    Frame->True,
    GridLines->Automatic,
    PlotRange->{{-2, 2}, {-2, 2}}], ",", 
   GraphicsBox[{CircleBox[{0, 0}, 0.10700243225641431`], 
     CircleBox[{0, 0}, 0.5782831713533064], 
     {Dashing[{0.01, 0.01}], 
      CircleBox[{0, -0.38069328627820354`}, 0.38069328627820354`], 
      CircleBox[{-0.38069328627820354`, 0}, 0.38069328627820354`], 
      CircleBox[{0.38069328627820354`, 0}, 0.38069328627820354`], 
      CircleBox[{0, 0.38069328627820354`}, 0.38069328627820354`], 
      CircleBox[{0, 0}, 0.3902878194986044]}},
    AspectRatio->Automatic,
    Axes->True,
    Frame->True,
    GridLines->Automatic,
    PlotRange->{{-2, 2}, {-2, 2}}], ",", 
   GraphicsBox[{CircleBox[{0, 0}, 0.09972707193611059], 
     CircleBox[{0, 0}, 0.5271957960864525], 
     {Dashing[{0.01, 0.01}], 
      CircleBox[{0, -0.34760331335404854`}, 0.34760331335404854`], 
      CircleBox[{-0.34760331335404854`, 0}, 0.34760331335404854`], 
      CircleBox[{0.34760331335404854`, 0}, 0.34760331335404854`], 
      CircleBox[{0, 0.34760331335404854`}, 0.34760331335404854`], 
      CircleBox[{0, 0}, 0.3554611120253503]}},
    AspectRatio->Automatic,
    Axes->True,
    Frame->True,
    GridLines->Automatic,
    PlotRange->{{-2, 2}, {-2, 2}}], ",", 
   GraphicsBox[{CircleBox[{0, 0}, 0.09337752959951504], 
     CircleBox[{0, 0}, 0.4844186272782812], 
     {Dashing[{0.01, 0.01}], 
      CircleBox[{0, -0.31979629899225614`}, 0.31979629899225614`], 
      CircleBox[{-0.31979629899225614`, 0}, 0.31979629899225614`], 
      CircleBox[{0.31979629899225614`, 0}, 0.31979629899225614`], 
      CircleBox[{0, 0.31979629899225614`}, 0.31979629899225614`], 
      CircleBox[{0, 0}, 0.32632666697388263`]}},
    AspectRatio->Automatic,
    Axes->True,
    Frame->True,
    GridLines->Automatic,
    PlotRange->{{-2, 2}, {-2, 2}}], ",", 
   GraphicsBox[{CircleBox[{0, 0}, 0.08778632759816841], 
     CircleBox[{0, 0}, 0.4481117913167129], 
     {Dashing[{0.01, 0.01}], 
      CircleBox[{0, -0.29611015926291645`}, 0.29611015926291645`], 
      CircleBox[{-0.29611015926291645`, 0}, 0.29611015926291645`], 
      CircleBox[{0.29611015926291645`, 0}, 0.29611015926291645`], 
      CircleBox[{0, 0.29611015926291645`}, 0.29611015926291645`], 
      CircleBox[{0, 0}, 0.3016501513979271]}},
    AspectRatio->Automatic,
    Axes->True,
    Frame->True,
    GridLines->Automatic,
    PlotRange->{{-2, 2}, {-2, 2}}], ",", 
   GraphicsBox[{CircleBox[{0, 0}, 0.08282655232738126], 
     CircleBox[{0, 0}, 0.416866054429945], 
     {Dashing[{0.01, 0.01}], 
      CircleBox[{0, -0.2756977039977251}, 0.2756977039977251], 
      CircleBox[{-0.2756977039977251, 0}, 0.2756977039977251], 
      CircleBox[{0.2756977039977251, 0}, 0.2756977039977251], 
      CircleBox[{0, 0.2756977039977251}, 0.2756977039977251], 
      CircleBox[{0, 0}, 0.2804375478878805]}},
    AspectRatio->Automatic,
    Axes->True,
    Frame->True,
    GridLines->Automatic,
    PlotRange->{{-2, 2}, {-2, 2}}]}], "}"}]], "Output",
 CellChangeTimes->{
  3.576132593933447*^9, {3.576132647909542*^9, 3.5761326770971932`*^9}, 
   3.5761327075328465`*^9, 3.5761328243302517`*^9, {3.5761328715515347`*^9, 
   3.5761330515758505`*^9}, 3.5761337325014467`*^9, 3.5761338010791674`*^9, 
   3.5761338797657056`*^9, 3.576133960870248*^9, 3.5761343014656463`*^9}]
}, Open  ]]
},
WindowSize->{808, 815},
WindowMargins->{{Automatic, 277}, {Automatic, 13}},
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[557, 20, 5823, 167, 252, "Input"],
Cell[CellGroupData[{
Cell[6405, 191, 1755, 33, 52, "Input"],
Cell[8163, 226, 4550, 108, 232, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[12750, 339, 4117, 125, 212, "Input"],
Cell[16870, 466, 885, 25, 52, "Output"],
Cell[17758, 493, 883, 25, 52, "Output"],
Cell[18644, 520, 659, 14, 52, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[19340, 539, 4181, 104, 172, "Input"],
Cell[23524, 645, 6619, 150, 881, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)