(* 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[     46590,       1074]
NotebookOptionsPosition[     45229,       1028]
NotebookOutlinePosition[     45741,       1047]
CellTagsIndexPosition[     45698,       1044]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{"O2", " ", "with", " ", 
    RowBox[{"Linnet", "'"}], "s", " ", "clouds", "  ", "13.05", ".2013"}], 
   " ", "*)"}], "\n", 
  RowBox[{
   RowBox[{
    RowBox[{"Clear", "[", 
     RowBox[{
     "a", ",", "z", ",", "sig1", ",", "sig2", ",", "sig3", ",", "sig4", ",", 
      "k1", ",", "k2", ",", "k3", ",", "k4", ",", "k5", ",", "nc", ",", "vee",
       ",", "vne", ",", "vnn", ",", "R1", ",", "R2", ",", "R3", ",", "R4", 
      ",", "R5", ",", "\n", "xc", ",", "yc", ",", "zc", ",", "xn", ",", "yn", 
      ",", "zn", ",", "ow", ",", "cn", ",", "rr", ",", "d1", ",", "d2", ",", 
      "d3", ",", "d4", ",", "d5", ",", "i", ",", "j", ",", "d", ",", "t", ",",
       "sq", ",", "Ekin", ",", "Epot"}], "]"}], ";"}], "\n", 
   RowBox[{
    RowBox[{"z", "=", "8.0"}], ";"}], "\n", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{"pi", "=", "0.0"}], ";"}], "*)"}], " ", 
   RowBox[{"(*", " ", 
    RowBox[{
     RowBox[{
      RowBox[{"distance", " ", "of", " ", "bananas", " ", "from", " ", "C"}], 
      "-", 
      RowBox[{"C", " ", "axis"}]}], ",", " ", 
     RowBox[{"pi", " ", "=", 
      RowBox[{"0", " ", "is", " ", "regular", " ", "tetraedric"}]}]}], " ", 
    "*)"}], "\n", 
   RowBox[{
    RowBox[{"nc", "=", "13"}], ";"}], "  ", 
   RowBox[{"(*", " ", 
    RowBox[{"number", " ", "of", " ", "clouds"}], " ", "*)"}], "\n", 
   RowBox[{
    RowBox[{"sig1", "=", "0.3"}], ";", " ", 
    RowBox[{"sig2", "=", "0.3"}], ";", " ", 
    RowBox[{"sig3", "=", "sig2"}], ";"}], " ", 
   RowBox[{"(*", " ", 
    RowBox[{
     RowBox[{"screening", " ", 
      RowBox[{"const", ".", " ", "from"}], " ", "e"}], "-", 
     RowBox[{"e", " ", "interaction", " ", "in", " ", "doubly", " ", 
      RowBox[{"occ", ".", " ", "clouds"}]}]}], " ", "*)"}], "\n", 
   RowBox[{
    RowBox[{"k1", "=", "0.99"}], ";", " ", 
    RowBox[{"k2", "=", "0.9"}], ";", " ", 
    RowBox[{"k3", "=", "0.9"}], ";", " ", 
    RowBox[{"k4", "=", "1.0"}], ";", " ", 
    RowBox[{"k5", "=", "1.3"}], ";"}], 
   RowBox[{"(*", " ", 
    RowBox[{
     RowBox[{
     "parameters", " ", "for", " ", "kinetic", " ", "energy", " ", "of", " ", 
      "clouds"}], ";", " ", 
     RowBox[{"k", "=", 
      RowBox[{"1.0", " ", 
       RowBox[{"Kimball", "'"}], "s", " ", "lowest", " ", "value"}]}]}], " ", 
    "*)"}], "\n", 
   RowBox[{
    RowBox[{"bohr", "=", "0.529177"}], ";"}], " ", "\n", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{"R2", "=", "1.1"}], ";", " ", 
     RowBox[{"R3", "=", "0.98"}], ";"}], "*)"}], 
   RowBox[{
    RowBox[{"R4", "=", "1.3"}], ";", " ", 
    RowBox[{"R5", "=", "1.15"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{"R1", "=", "0.18"}], ";", " ", 
     RowBox[{"R2", "=", "1.0"}], ";", " ", 
     RowBox[{"R3", "=", "0.96"}], ";", " ", 
     RowBox[{"R4", "=", "1.0"}], ";", " ", 
     RowBox[{"R5", "=", "1.0"}], ";"}], "*)"}]}]}]], "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5773415047046976`*^9, 3.5773415137059135`*^9}, {
   3.5773416513449545`*^9, 3.577341669752987*^9}, {3.5773417239318824`*^9, 
   3.577341893956581*^9}, 3.5773419407566633`*^9, {3.5773420175555983`*^9, 
   3.577342056165666*^9}, {3.5773420929973307`*^9, 3.5773421056489525`*^9}, {
   3.5773421396414127`*^9, 3.5773423470593767`*^9}, {3.577342382081438*^9, 
   3.577342441158742*^9}, {3.5773425028880506`*^9, 3.577342523682887*^9}, {
   3.5773425980794177`*^9, 3.5773426626791315`*^9}, {3.5773426933487854`*^9, 
   3.5773426941131864`*^9}, {3.577342738354864*^9, 3.577342889035529*^9}, {
   3.5773429412644205`*^9, 3.57734295234044*^9}, {3.577342991824109*^9, 
   3.5773430231333647`*^9}, {3.5773430561274223`*^9, 3.577343143737176*^9}, {
   3.5773431993668737`*^9, 3.577343258334977*^9}, 3.5773434426249013`*^9, {
   3.577343480042567*^9, 3.5773435125686235`*^9}, {3.5773436063903885`*^9, 
   3.5773436272008247`*^9}, {3.5773436617392855`*^9, 
   3.5773437808610945`*^9}, {3.577343868193248*^9, 3.577343879659268*^9}, {
   3.5773439214829416`*^9, 3.577344118838888*^9}, {3.5773441503509436`*^9, 
   3.577344170740179*^9}, {3.5773446920694942`*^9, 3.577344693286297*^9}, 
   3.5773526894956493`*^9, {3.5773528388503113`*^9, 3.577352838897112*^9}, {
   3.577353457328998*^9, 3.577353503364679*^9}, {3.5773556679320807`*^9, 
   3.5773557165573664`*^9}, 3.5773607112650003`*^9, {3.577363145384075*^9, 
   3.5773631641821084`*^9}, 3.5773713130222435`*^9, {3.5773735851350346`*^9, 
   3.577373695474028*^9}, {3.577373964730501*^9, 3.5773739677725067`*^9}, 
   3.5773740596254673`*^9, {3.577374092245125*^9, 3.5773741045067463`*^9}, {
   3.5773741853616886`*^9, 3.577374230414568*^9}, 3.5773742637986264`*^9, {
   3.5773746139412413`*^9, 3.5773746339092765`*^9}, {3.5773750495252066`*^9, 
   3.5773750633000307`*^9}, 3.5773754561555204`*^9, {3.5773757271123962`*^9, 
   3.577375728063998*^9}, 3.577375963078411*^9, {3.5773761205918875`*^9, 
   3.5773762048008356`*^9}, 3.5774291695117903`*^9, {3.5774424964722877`*^9, 
   3.5774424977046895`*^9}, {3.5774428197828555`*^9, 3.577442825336465*^9}, {
   3.577442927267044*^9, 3.5774429281406455`*^9}, {3.5774429766411304`*^9, 
   3.577442992833959*^9}, {3.5774430583072743`*^9, 3.577443066715689*^9}, {
   3.5774436977055974`*^9, 3.5774437122604227`*^9}, 3.5774438419902506`*^9, {
   3.577446912512444*^9, 3.577446913651246*^9}, {3.577447169460495*^9, 
   3.577447179070112*^9}, {3.577449454744052*^9, 3.57744946514927*^9}, {
   3.577449615377534*^9, 3.5774496616940155`*^9}, {3.577450578242425*^9, 
   3.5774506245901065`*^9}, {3.577451618688603*^9, 3.5774516195310044`*^9}, {
   3.5774516668614874`*^9, 3.5774516695602922`*^9}, 3.577455410676982*^9, {
   3.5774557532847834`*^9, 3.5774558029396706`*^9}, {3.5774565586985984`*^9, 
   3.5774565589793987`*^9}, 3.577457075496306*^9, {3.577457117273179*^9, 
   3.577457148941235*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{"C", " ", "He"}], "-", "shells"}], " ", "*)"}], "\n", 
  RowBox[{
   RowBox[{
    RowBox[{"ow", "=", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"-", "2"}], ",", 
       RowBox[{"-", "2"}], ",", 
       RowBox[{"-", "2"}], ",", 
       RowBox[{"-", "2"}], ",", 
       RowBox[{"-", "2"}], ",", 
       RowBox[{"-", "1"}], ",", 
       RowBox[{"-", "1"}], ",", 
       RowBox[{"-", "1"}], ",", 
       RowBox[{"-", "1"}], ",", 
       RowBox[{"-", "0.5"}], ",", 
       RowBox[{"-", "0.5"}], ",", 
       RowBox[{"-", "0.5"}], ",", 
       RowBox[{"-", "0.5"}]}], "}"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"cn", "=", 
     RowBox[{"{", " ", 
      RowBox[{"z", ",", "z"}], "}"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Ekin", " ", "=", " ", 
     RowBox[{"2", "*", 
      RowBox[{"(", 
       RowBox[{"2.25", "*", 
        RowBox[{"k1", "/", 
         RowBox[{"R1", "^", "2"}]}]}], ")"}]}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"vee", "=", 
     RowBox[{"2", "*", 
      RowBox[{"(", 
       RowBox[{"3.0", "*", 
        RowBox[{"sig1", "/", "R1"}]}], ")"}]}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"vne", "=", 
     RowBox[{
      RowBox[{"-", "2"}], "*", 
      RowBox[{"(", 
       RowBox[{"3.0", "*", 
        RowBox[{"z", "/", "R1"}]}], ")"}]}]}], ";"}], "\n", 
   RowBox[{"(*", " ", 
    RowBox[{"bonding", " ", "pairs"}], " ", "*)"}], "\n", 
   RowBox[{
    RowBox[{"Ekin", " ", "=", " ", 
     RowBox[{"Ekin", " ", "+", " ", 
      RowBox[{"2.25", "*", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"2", "*", 
          RowBox[{"k2", "/", 
           RowBox[{"R2", "^", "2"}]}]}], "+", 
         RowBox[{"k3", "/", 
          RowBox[{"R3", "^", "2"}]}], "+", 
         RowBox[{"2", "*", 
          RowBox[{"k4", "/", 
           RowBox[{"R4", "^", "2"}]}]}]}], ")"}]}]}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"vee", "=", 
     RowBox[{"vee", "+", 
      RowBox[{"3.0", "*", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"2", "*", 
          RowBox[{"sig2", "/", "R2"}]}], "+", 
         RowBox[{"sig3", "/", "R3"}]}], ")"}]}]}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"d1", "=", 
     RowBox[{"R1", "+", "R3"}]}], ";"}], "  ", "\n", 
   RowBox[{
    RowBox[{"d2", "=", 
     RowBox[{"d1", "+", "R1", "+", "R2"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"xn", "=", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"-", "d1"}], ",", "d1", ",", "0", ",", "0", ",", "0", ",", "0",
        ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", 
       "0"}], "}"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"yn", "=", 
     RowBox[{"{", 
      RowBox[{
      "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", 
       "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0"}], "}"}]}], 
    ";"}], "\n", 
   RowBox[{
    RowBox[{"zn", "=", 
     RowBox[{"{", 
      RowBox[{
      "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", 
       "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0"}], "}"}]}], 
    ";"}], "\n", 
   RowBox[{"(*", " ", 
    RowBox[{
    "cloud", " ", "coordinates", " ", "in", " ", "terms", " ", "of", " ", 
     "radii"}], " ", "*)"}], "        ", "\n", 
   RowBox[{
    RowBox[{"xc", "=", 
     RowBox[{"{", 
      RowBox[{"d1", ",", 
       RowBox[{"-", "d1"}], ",", "d2", ",", 
       RowBox[{"-", "d2"}], ",", "0", ",", "d1", ",", " ", "d1", ",", 
       RowBox[{"-", "d1"}], ",", 
       RowBox[{"-", "d1"}], ",", "d1", ",", "d1", ",", 
       RowBox[{"-", "d1"}], ",", 
       RowBox[{"-", "d1"}]}], "}"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"yc", "=", 
     RowBox[{"{", 
      RowBox[{
      "0", ",", "0", ",", "0", ",", "0", ",", "                ", "0", ",", 
       "R4", ",", 
       RowBox[{"-", "R4"}], ",", " ", "R4", ",", 
       RowBox[{"-", "R4"}], ",", " ", "0", ",", "      ", "0", ",", "  ", "0",
        ",", "      ", "0"}], "}"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"zc", "=", 
     RowBox[{"{", 
      RowBox[{
      "0", ",", "0", ",", "0", ",", "0", ",", "                ", "0", ",", 
       "0", ",", "      ", "0", ",", "    ", "0", ",", "      ", "0", ",", 
       " ", "R5", ",", 
       RowBox[{"-", "R5"}], ",", "     ", "R5", ",", 
       RowBox[{"-", "R5"}]}], "}"}]}], ";"}], "\n", 
   RowBox[{"(*", " ", 
    RowBox[{"cc", ":", " ", 
     RowBox[{
      RowBox[{"sum", " ", "of", " ", "cloud"}], "-", 
      RowBox[{"cloud", " ", "potential", " ", "energies"}]}]}], " ", "*)"}], 
   "\n", 
   RowBox[{"For", "[", 
    RowBox[{
     RowBox[{"i", " ", "=", " ", "1"}], ",", " ", 
     RowBox[{"i", " ", "<", " ", "nc"}], ",", " ", 
     RowBox[{"i", "++"}], ",", "\n", "  ", 
     RowBox[{"For", "[", 
      RowBox[{
       RowBox[{"j", " ", "=", " ", 
        RowBox[{"i", "+", "1"}]}], ",", " ", 
       RowBox[{"j", " ", "\[LessEqual]", " ", "nc"}], ",", " ", 
       RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", 
       RowBox[{"vee", " ", "=", " ", 
        RowBox[{"vee", " ", "+", " ", 
         RowBox[{
          RowBox[{"ow", "[", 
           RowBox[{"[", "i", "]"}], "]"}], "*", 
          RowBox[{
           RowBox[{"ow", "[", 
            RowBox[{"[", "j", "]"}], "]"}], "/", 
           RowBox[{"Sqrt", "[", 
            RowBox[{
             RowBox[{
              RowBox[{"(", 
               RowBox[{
                RowBox[{"xc", "[", 
                 RowBox[{"[", "i", "]"}], "]"}], "-", 
                RowBox[{"xc", "[", 
                 RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}], "+", 
             RowBox[{
              RowBox[{"(", 
               RowBox[{
                RowBox[{"yc", "[", 
                 RowBox[{"[", "i", "]"}], "]"}], "-", 
                RowBox[{"yc", "[", 
                 RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}], "+", 
             RowBox[{
              RowBox[{"(", 
               RowBox[{
                RowBox[{"zc", "[", 
                 RowBox[{"[", "i", "]"}], "]"}], "-", 
                RowBox[{"zc", "[", 
                 RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}]}], 
            "]"}]}]}]}]}]}], "]"}]}], "]"}], "\n", 
   RowBox[{
    RowBox[{"vee", "=", 
     RowBox[{"vee", "-", 
      RowBox[{"2", "*", 
       RowBox[{"0.25", "/", 
        RowBox[{"(", 
         RowBox[{"2", "*", "R5"}], ")"}]}]}]}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{"nn", ":", " ", 
     RowBox[{
      RowBox[{"sum", " ", "of", " ", "nuclei"}], "-", 
      RowBox[{"nuclei", " ", "potential", " ", "energies"}]}]}], " ", "*)"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"vnn", " ", "=", 
     RowBox[{
      RowBox[{"z", "^", "2"}], "/", 
      RowBox[{"(", 
       RowBox[{"2", "*", "d1"}], ")"}]}]}], ";"}], "\n", 
   RowBox[{"(*", " ", 
    RowBox[{"cn", ":", " ", 
     RowBox[{
      RowBox[{"sum", " ", "of", " ", "cloud"}], "-", 
      RowBox[{"nuclei", " ", "potential", " ", "energies"}]}]}], " ", "*)"}], 
   "\n", 
   RowBox[{"For", "[", 
    RowBox[{
     RowBox[{"i", " ", "=", " ", "2"}], ",", " ", 
     RowBox[{"i", " ", "\[LessEqual]", " ", "nc"}], ",", " ", 
     RowBox[{"i", "++"}], ",", "\n", "  ", 
     RowBox[{
      RowBox[{"sq", " ", "=", " ", 
       RowBox[{"Sqrt", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"xc", "[", 
             RowBox[{"[", "i", "]"}], "]"}], "-", 
            RowBox[{"xn", "[", 
             RowBox[{"[", "1", "]"}], "]"}]}], ")"}], "^", "2"}], "+", 
         RowBox[{
          RowBox[{"yc", "[", 
           RowBox[{"[", "i", "]"}], "]"}], "^", "2"}], "+", 
         RowBox[{
          RowBox[{"zc", "[", 
           RowBox[{"[", "i", "]"}], "]"}], "^", "2"}]}], "]"}]}], ";", "\n", 
      "    ", 
      RowBox[{"If", "[", " ", 
       RowBox[{
        RowBox[{"sq", " ", "\[NotEqual]", " ", "0"}], ",", 
        RowBox[{"vne", " ", "=", " ", 
         RowBox[{"vne", " ", "+", " ", 
          RowBox[{"2", "*", "z", "*", 
           RowBox[{
            RowBox[{"ow", "[", 
             RowBox[{"[", "i", "]"}], "]"}], "/", "sq"}]}]}]}]}], "]"}]}]}], 
    "]"}], "\n", 
   RowBox[{
    RowBox[{"Epot", "=", 
     RowBox[{"vne", "+", "vee", "+", "vnn"}]}], " ", ";"}], "\n", 
   RowBox[{
    RowBox[{"func", "=", 
     RowBox[{"Ekin", "+", "Epot"}]}], ";"}], "\n", 
   RowBox[{"t", " ", "=", " ", 
    RowBox[{"FindMinimum", "[", 
     RowBox[{"func", ",", 
      RowBox[{"{", 
       RowBox[{"R1", ",", "0.2"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"R2", ",", "1.1"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"R3", ",", "0.98"}], "}"}], ",", 
      RowBox[{"(*", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"R4", ",", "1.3"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"R5", ",", "1.15"}], "}"}], ","}], "*)"}], 
      RowBox[{"Method", " ", "\[Rule]", "\"\<Newton\>\""}], ",", " ", 
      RowBox[{"MaxIterations", " ", "->", " ", "500"}]}], 
     "]"}]}]}]}]], "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5773415047046976`*^9, 3.5773415137059135`*^9}, {
   3.5773416513449545`*^9, 3.577341669752987*^9}, {3.5773417239318824`*^9, 
   3.577341893956581*^9}, 3.5773419407566633`*^9, {3.5773420175555983`*^9, 
   3.577342056165666*^9}, {3.5773420929973307`*^9, 3.5773421056489525`*^9}, {
   3.5773421396414127`*^9, 3.5773423470593767`*^9}, {3.577342382081438*^9, 
   3.577342441158742*^9}, {3.5773425028880506`*^9, 3.577342523682887*^9}, {
   3.5773425980794177`*^9, 3.5773426626791315`*^9}, {3.5773426933487854`*^9, 
   3.5773426941131864`*^9}, {3.577342738354864*^9, 3.577342889035529*^9}, {
   3.5773429412644205`*^9, 3.57734295234044*^9}, {3.577342991824109*^9, 
   3.5773430231333647`*^9}, {3.5773430561274223`*^9, 3.577343143737176*^9}, {
   3.5773431993668737`*^9, 3.577343258334977*^9}, 3.5773434426249013`*^9, {
   3.577343480042567*^9, 3.5773435125686235`*^9}, {3.5773436063903885`*^9, 
   3.5773436272008247`*^9}, {3.5773436617392855`*^9, 
   3.5773437808610945`*^9}, {3.577343868193248*^9, 3.577343879659268*^9}, {
   3.5773439214829416`*^9, 3.577344118838888*^9}, {3.5773441503509436`*^9, 
   3.577344170740179*^9}, {3.5773446920694942`*^9, 3.577344693286297*^9}, 
   3.5773526894956493`*^9, {3.577352838897112*^9, 3.577352838928312*^9}, 
   3.577362343573867*^9, {3.577362462570876*^9, 3.577362524440585*^9}, 
   3.5773625762794757`*^9, {3.5773626257315626`*^9, 3.577362806021079*^9}, {
   3.5773631418428693`*^9, 3.577363185554146*^9}, 3.5773632603562775`*^9, {
   3.577371455731294*^9, 3.5773714994893713`*^9}, {3.5773724898259106`*^9, 
   3.577372494162718*^9}, {3.5773730043772144`*^9, 3.5773730238460484`*^9}, {
   3.5773731837619295`*^9, 3.577373238128025*^9}, {3.577373271434083*^9, 
   3.577373393987899*^9}, {3.577373430913163*^9, 3.5773734337835684`*^9}, {
   3.577373480833251*^9, 3.5773734946236753`*^9}, {3.5773739993625617`*^9, 
   3.5773740373018284`*^9}, {3.5773743739348197`*^9, 
   3.5773743842152376`*^9}, {3.5773744685645857`*^9, 3.577374501527444*^9}, {
   3.5773746732993455`*^9, 3.577374707073405*^9}, {3.577375258362773*^9, 
   3.577375352243738*^9}, 3.577375468307942*^9, {3.577375746784031*^9, 
   3.577375845345004*^9}, {3.577375895483492*^9, 3.5773759185091324`*^9}, 
   3.577376045072155*^9, 3.5773761022150555`*^9, {3.577376336184266*^9, 
   3.577376368523123*^9}, {3.577416169561963*^9, 3.5774161709503655`*^9}, {
   3.5774290920108542`*^9, 3.577429100450469*^9}, {3.577429155814966*^9, 
   3.577429156813368*^9}, {3.5774425034142995`*^9, 3.5774426597733746`*^9}, {
   3.5774426930950327`*^9, 3.5774427236242867`*^9}, {3.577442792748008*^9, 
   3.577442798395218*^9}, 3.577442841279693*^9, {3.5774430051891813`*^9, 
   3.577443009635189*^9}, {3.5774430800225124`*^9, 3.577443093625736*^9}, {
   3.577443239392392*^9, 3.5774432451488028`*^9}, {3.5774469888589783`*^9, 
   3.5774470085618124`*^9}, {3.577447132940831*^9, 3.5774471595388775`*^9}, {
   3.577449541277404*^9, 3.5774495651610456`*^9}, {3.577449860999965*^9, 
   3.5774498666939754`*^9}, {3.5774500984479823`*^9, 
   3.5774501088844004`*^9}, {3.577451894887088*^9, 3.577451972668825*^9}, {
   3.5774553353600492`*^9, 3.5774553614276953`*^9}, {3.577456222674008*^9, 
   3.577456231300823*^9}, {3.577456262485278*^9, 3.5774563142149687`*^9}, {
   3.5774564314179745`*^9, 3.5774564367843843`*^9}, {3.577456616169099*^9, 
   3.57745662209711*^9}, {3.577456952084489*^9, 3.577457050863863*^9}}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"FindMinimum", "::", "lstol"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"The line search decreased the step size to within the \
tolerance specified by AccuracyGoal and PrecisionGoal but was unable to find \
a sufficient decrease in the function. You may need more than \
\[NoBreak]\\!\\(MachinePrecision\\)\[NoBreak] digits of working precision to \
meet these tolerances. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \
ButtonStyle->\\\"Link\\\", ButtonFrame->None, \
ButtonData:>\\\"paclet:ref/message/FindMinimum/lstol\\\", ButtonNote -> \
\\\"FindMinimum::lstol\\\"]\\)\"\>"}]], "Message", "MSG",
 CellChangeTimes->{3.577457157006449*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", "142.9375353891195`"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"R1", "\[Rule]", "0.19285714279201807`"}], ",", 
     RowBox[{"R2", "\[Rule]", "1.2802758809663549`*^7"}], ",", 
     RowBox[{"R3", "\[Rule]", "1.8364951733505525`*^7"}]}], "}"}]}], 
  "}"}]], "Output",
 CellChangeTimes->{3.5774504994154863`*^9, 3.577450631890919*^9, 
  3.5774507316062946`*^9, 3.577451627580619*^9, 3.5774516740063*^9, 
  3.5774520186109056`*^9, 3.5774551196116705`*^9, 3.577455371240113*^9, 
  3.577455500065139*^9, 3.577455581466082*^9, 3.577455683287461*^9, 
  3.5774560517133083`*^9, 3.5774563553522415`*^9, 3.5774564425719943`*^9, 
  3.5774565674502134`*^9, 3.5774566327675285`*^9, 3.5774566931552343`*^9, 
  3.577457157006449*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", "func", " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"Ekin", " ", "/.", " ", 
    RowBox[{"t", "[", 
     RowBox[{"[", "2", "]"}], "]"}]}], "\n", 
   RowBox[{"vne", "  ", "/.", " ", 
    RowBox[{"t", "[", 
     RowBox[{"[", "2", "]"}], "]"}]}], " ", "\n", 
   RowBox[{"vee", "  ", "/.", " ", 
    RowBox[{"t", "[", 
     RowBox[{"[", "2", "]"}], "]"}]}], " ", "\n", 
   RowBox[{"vnn", "  ", "/.", " ", 
    RowBox[{"t", "[", 
     RowBox[{"[", "2", "]"}], "]"}]}], "\n", 
   RowBox[{"Vir", " ", "=", " ", 
    RowBox[{
     RowBox[{
      RowBox[{"-", "Epot"}], "/", "Ekin"}], "  ", "/.", " ", 
     RowBox[{"t", "[", 
      RowBox[{"[", "2", "]"}], "]"}]}]}], "\n", 
   RowBox[{
    RowBox[{"2", "*", "d1", "*", "bohr"}], "  ", "/.", " ", 
    RowBox[{"t", "[", 
     RowBox[{"[", "2", "]"}], "]"}]}]}]}]], "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5773415047046976`*^9, 3.5773415137059135`*^9}, {
   3.5773416513449545`*^9, 3.577341669752987*^9}, {3.5773417239318824`*^9, 
   3.577341893956581*^9}, 3.5773419407566633`*^9, {3.5773420175555983`*^9, 
   3.577342056165666*^9}, {3.5773420929973307`*^9, 3.5773421056489525`*^9}, {
   3.5773421396414127`*^9, 3.5773423470593767`*^9}, {3.577342382081438*^9, 
   3.577342441158742*^9}, {3.5773425028880506`*^9, 3.577342523682887*^9}, {
   3.5773425980794177`*^9, 3.5773426626791315`*^9}, {3.5773426933487854`*^9, 
   3.5773426941131864`*^9}, {3.577342738354864*^9, 3.577342889035529*^9}, {
   3.5773429412644205`*^9, 3.57734295234044*^9}, {3.577342991824109*^9, 
   3.5773430231333647`*^9}, {3.5773430561274223`*^9, 3.577343143737176*^9}, {
   3.5773431993668737`*^9, 3.577343258334977*^9}, 3.5773434426249013`*^9, {
   3.577343480042567*^9, 3.5773435125686235`*^9}, {3.5773436063903885`*^9, 
   3.5773436272008247`*^9}, {3.5773436617392855`*^9, 
   3.5773437808610945`*^9}, {3.577343868193248*^9, 3.577343879659268*^9}, {
   3.5773439214829416`*^9, 3.577344086609232*^9}, {3.5773444013153844`*^9, 
   3.5773444098953996`*^9}, {3.5773527214601054`*^9, 3.577352754235763*^9}, {
   3.577352839084312*^9, 3.577352839115512*^9}, {3.5773562861299667`*^9, 
   3.577356291496376*^9}, {3.5773717712106485`*^9, 3.5773718737028284`*^9}, {
   3.5774560203260527`*^9, 3.577456044802496*^9}}],

Cell[BoxData["122.44049975216302`"], "Output",
 CellChangeTimes->{3.5774504994154863`*^9, 3.577450631906519*^9, 
  3.5774507316218944`*^9, 3.5774516275962186`*^9, 3.5774516740219*^9, 
  3.5774520186265054`*^9, 3.5774551196272707`*^9, 3.5774553712557125`*^9, 
  3.577455500080739*^9, 3.577455581466082*^9, 3.5774556833030605`*^9, 
  3.577456051728908*^9, 3.5774563553678412`*^9, 3.5774564426967945`*^9, 
  3.577456567559414*^9, 3.5774566327987285`*^9, 3.5774566931708345`*^9, 
  3.577457157022049*^9}],

Cell[BoxData[
 RowBox[{"-", "287.4173170665802`"}]], "Output",
 CellChangeTimes->{3.5774504994154863`*^9, 3.577450631906519*^9, 
  3.5774507316218944`*^9, 3.5774516275962186`*^9, 3.5774516740219*^9, 
  3.5774520186265054`*^9, 3.5774551196272707`*^9, 3.5774553712557125`*^9, 
  3.577455500080739*^9, 3.577455581466082*^9, 3.5774556833030605`*^9, 
  3.577456051728908*^9, 3.5774563553678412`*^9, 3.5774564426967945`*^9, 
  3.577456567559414*^9, 3.5774566327987285`*^9, 3.5774566931708345`*^9, 
  3.577457157022049*^9}],

Cell[BoxData["22.039280182848312`"], "Output",
 CellChangeTimes->{3.5774504994154863`*^9, 3.577450631906519*^9, 
  3.5774507316218944`*^9, 3.5774516275962186`*^9, 3.5774516740219*^9, 
  3.5774520186265054`*^9, 3.5774551196272707`*^9, 3.5774553712557125`*^9, 
  3.577455500080739*^9, 3.577455581466082*^9, 3.5774556833030605`*^9, 
  3.577456051728908*^9, 3.5774563553678412`*^9, 3.5774564426967945`*^9, 
  3.577456567559414*^9, 3.5774566327987285`*^9, 3.5774566931708345`*^9, 
  3.577457157022049*^9}],

Cell[BoxData["1.7424494291250709`*^-6"], "Output",
 CellChangeTimes->{3.5774504994154863`*^9, 3.577450631906519*^9, 
  3.5774507316218944`*^9, 3.5774516275962186`*^9, 3.5774516740219*^9, 
  3.5774520186265054`*^9, 3.5774551196272707`*^9, 3.5774553712557125`*^9, 
  3.577455500080739*^9, 3.577455581466082*^9, 3.5774556833030605`*^9, 
  3.577456051728908*^9, 3.5774563553678412`*^9, 3.5774564426967945`*^9, 
  3.577456567559414*^9, 3.5774566327987285`*^9, 3.5774566931708345`*^9, 
  3.577457157022049*^9}],

Cell[BoxData["2.1674040507711534`"], "Output",
 CellChangeTimes->{3.5774504994154863`*^9, 3.577450631906519*^9, 
  3.5774507316218944`*^9, 3.5774516275962186`*^9, 3.5774516740219*^9, 
  3.5774520186265054`*^9, 3.5774551196272707`*^9, 3.5774553712557125`*^9, 
  3.577455500080739*^9, 3.577455581466082*^9, 3.5774556833030605`*^9, 
  3.577456051728908*^9, 3.5774563553678412`*^9, 3.5774564426967945`*^9, 
  3.577456567559414*^9, 3.5774566327987285`*^9, 3.5774566931708345`*^9, 
  3.577457157022049*^9}],

Cell[BoxData["1.9436620331073634`*^7"], "Output",
 CellChangeTimes->{3.5774504994154863`*^9, 3.577450631906519*^9, 
  3.5774507316218944`*^9, 3.5774516275962186`*^9, 3.5774516740219*^9, 
  3.5774520186265054`*^9, 3.5774551196272707`*^9, 3.5774553712557125`*^9, 
  3.577455500080739*^9, 3.577455581466082*^9, 3.5774556833030605`*^9, 
  3.577456051728908*^9, 3.5774563553678412`*^9, 3.5774564426967945`*^9, 
  3.577456567559414*^9, 3.5774566327987285`*^9, 3.5774566931708345`*^9, 
  3.577457157022049*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"rr", "=", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
     "R1", ",", "R1", ",", "R2", ",", "R2", ",", "R3", ",", "R4", ",", "R4", 
      ",", "R4", ",", "R4", ",", "R5", ",", "R5", ",", "R5", ",", "R5"}], 
     "}"}], " ", "/.", " ", 
    RowBox[{"t", "[", 
     RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}]], "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5773415047046976`*^9, 3.5773415137059135`*^9}, {
   3.5773416513449545`*^9, 3.577341669752987*^9}, {3.5773417239318824`*^9, 
   3.577341893956581*^9}, 3.5773419407566633`*^9, {3.5773420175555983`*^9, 
   3.577342056165666*^9}, {3.5773420929973307`*^9, 3.5773421056489525`*^9}, {
   3.5773421396414127`*^9, 3.5773423470593767`*^9}, {3.577342382081438*^9, 
   3.577342441158742*^9}, {3.5773425028880506`*^9, 3.577342523682887*^9}, {
   3.5773425980794177`*^9, 3.5773426626791315`*^9}, {3.5773426933487854`*^9, 
   3.5773426941131864`*^9}, {3.577342738354864*^9, 3.577342889035529*^9}, {
   3.5773429412644205`*^9, 3.57734295234044*^9}, {3.577342991824109*^9, 
   3.5773430231333647`*^9}, {3.5773430561274223`*^9, 3.577343143737176*^9}, {
   3.5773431993668737`*^9, 3.577343258334977*^9}, 3.5773434426249013`*^9, {
   3.577343480042567*^9, 3.5773435125686235`*^9}, {3.5773436063903885`*^9, 
   3.5773436272008247`*^9}, {3.5773436617392855`*^9, 
   3.5773437808610945`*^9}, {3.577343868193248*^9, 3.577343879659268*^9}, {
   3.5773439214829416`*^9, 3.577344086609232*^9}, {3.5773444013153844`*^9, 
   3.5773444098953996`*^9}, {3.5773527214601054`*^9, 3.577352754235763*^9}, {
   3.577352839084312*^9, 3.577352839115512*^9}, {3.5773562861299667`*^9, 
   3.577356291496376*^9}, {3.5773717712106485`*^9, 3.577371904793683*^9}, {
   3.577372601444106*^9, 3.577372628525754*^9}, 3.577375878432662*^9, 
   3.577442909545413*^9, {3.5774471930165367`*^9, 3.5774472038429556`*^9}, {
   3.5774495179865627`*^9, 3.577449527440179*^9}, {3.577455105602846*^9, 
   3.577455114822462*^9}, {3.577455475619896*^9, 3.577455493403927*^9}}],

Cell[BoxData[
 RowBox[{"(*", " ", 
  RowBox[{
   RowBox[{"projection", " ", "on", " ", "xy"}], "-", 
   RowBox[{"plane", " ", "of", " ", "molecule"}]}], " ", "*)"}]], "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5773413317003937`*^9, 3.5773413362088013`*^9}, {
   3.577352839115512*^9, 3.577352839131112*^9}, {3.5773715990175457`*^9, 
   3.5773716034635534`*^9}, {3.5773716711988726`*^9, 
   3.5773716785932856`*^9}, {3.5773719148557005`*^9, 
   3.5773720320119066`*^9}, {3.5773720625879602`*^9, 3.577372089326407*^9}, {
   3.5773721384508934`*^9, 3.577372247229884*^9}, {3.5773723618276854`*^9, 
   3.577372407098965*^9}, {3.577372640428575*^9, 3.57737267169103*^9}, {
   3.5773727054338894`*^9, 3.5773727803452206`*^9}, {3.577372855490553*^9, 
   3.5773728739765854`*^9}, {3.5773729151762576`*^9, 3.577372935331493*^9}, 
   3.577373068493327*^9}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"plot0", "=", 
   RowBox[{
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"Graphics3D", "[", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"Darker", "[", 
          RowBox[{"Red", ",", "1"}], "]"}], ",", 
         RowBox[{"Sphere", "[", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{
             RowBox[{"xc", "[", 
              RowBox[{"[", "i", "]"}], "]"}], ",", 
             RowBox[{"yc", "[", 
              RowBox[{"[", "i", "]"}], "]"}], ",", 
             RowBox[{"zc", "[", 
              RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", 
           RowBox[{"rr", "[", 
            RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}], "}"}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "2"}], "}"}]}], "]"}], " ", "/.", " ", 
    RowBox[{"t", "[", 
     RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"plot1", "=", 
   RowBox[{
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"Graphics3D", "[", 
       RowBox[{"{", 
        RowBox[{"Sphere", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"xc", "[", 
             RowBox[{"[", "i", "]"}], "]"}], ",", 
            RowBox[{"yc", "[", 
             RowBox[{"[", "i", "]"}], "]"}], ",", 
            RowBox[{"zc", "[", 
             RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", 
          RowBox[{"rr", "[", 
           RowBox[{"[", "i", "]"}], "]"}]}], "]"}], "}"}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "3", ",", "5"}], "}"}]}], "]"}], " ", "/.", " ", 
    RowBox[{"t", "[", 
     RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"plot2", "=", 
   RowBox[{
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"Graphics3D", "[", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"Opacity", "[", "0.5", "]"}], ",", 
         RowBox[{"Lighter", "[", 
          RowBox[{"Green", ",", "0.4"}], "]"}], ",", 
         RowBox[{"Sphere", "[", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{
             RowBox[{"xc", "[", 
              RowBox[{"[", "i", "]"}], "]"}], ",", 
             RowBox[{"yc", "[", 
              RowBox[{"[", "i", "]"}], "]"}], ",", 
             RowBox[{"zc", "[", 
              RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", 
           RowBox[{"rr", "[", 
            RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}], "}"}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "6", ",", "9"}], "}"}]}], "]"}], " ", "/.", " ", 
    RowBox[{"t", "[", 
     RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"plot2a", "=", 
    RowBox[{
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"Graphics3D", "[", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"Opacity", "[", "0.5", "]"}], ",", 
          RowBox[{"Lighter", "[", 
           RowBox[{"Cyan", ",", "0.3"}], "]"}], ",", 
          RowBox[{"Sphere", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{"xc", "[", 
               RowBox[{"[", "i", "]"}], "]"}], ",", 
              RowBox[{"yc", "[", 
               RowBox[{"[", "i", "]"}], "]"}], ",", 
              RowBox[{"zc", "[", 
               RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", 
            RowBox[{"rr", "[", 
             RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}], "}"}], "]"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "10", ",", "nc"}], "}"}]}], "]"}], " ", "/.", " ", 
     RowBox[{"t", "[", 
      RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], 
  "\n"}], "\[IndentingNewLine]", 
 RowBox[{"Show", "[", 
  RowBox[{"plot0", ",", "plot1", ",", "plot2", ",", "plot2a", ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"AspectRatio", " ", "\[Rule]", " ", "Automatic"}], ",", 
     RowBox[{"Axes", " ", "->", " ", "True"}], ",", 
     RowBox[{"PlotRange", " ", "\[Rule]", " ", "All"}], ",", " ", 
     RowBox[{"AxesLabel", " ", "->", " ", 
      RowBox[{"{", 
       RowBox[{"xl", ",", " ", "yl", ",", " ", "zl"}], "}"}]}], ",", 
     RowBox[{"SphericalRegion", "\[Rule]", "True"}]}], "}"}]}], 
  "]"}], "\n"}], "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5773413317003937`*^9, 3.5773413362088013`*^9}, {
   3.577352839115512*^9, 3.577352839131112*^9}, {3.5773715990175457`*^9, 
   3.5773716034635534`*^9}, {3.5773716711988726`*^9, 
   3.5773716785932856`*^9}, {3.5773719148557005`*^9, 
   3.5773720320119066`*^9}, {3.5773720625879602`*^9, 3.577372089326407*^9}, {
   3.5773721384508934`*^9, 3.577372247229884*^9}, {3.5773723618276854`*^9, 
   3.577372407098965*^9}, {3.577372640428575*^9, 3.57737267169103*^9}, {
   3.5773727054338894`*^9, 3.5773727803452206`*^9}, {3.577372855490553*^9, 
   3.5773728739765854`*^9}, {3.5773729151762576`*^9, 3.577372935331493*^9}, 
   3.577373068493327*^9, 3.577375863940237*^9, {3.577442873244149*^9, 
   3.5774428788913593`*^9}, {3.5774433556593966`*^9, 
   3.5774433846442475`*^9}, {3.5774434296971264`*^9, 
   3.5774434305863285`*^9}, {3.5774434802880154`*^9, 3.5774435257464952`*^9}, 
   3.577443563108561*^9, 3.5774436149318523`*^9, {3.577443975932086*^9, 
   3.5774439769772882`*^9}, {3.5774443030334606`*^9, 
   3.5774444499077187`*^9}, {3.577444584567155*^9, 3.5774445963139753`*^9}, {
   3.577444660835689*^9, 3.577444676014516*^9}, {3.577444738118225*^9, 
   3.577444777929495*^9}, 3.5774449044145164`*^9, {3.5774460837297883`*^9, 
   3.577446116255845*^9}, {3.577446165957533*^9, 3.5774461885775723`*^9}, {
   3.577446237358858*^9, 3.5774463455450478`*^9}, {3.577446401939147*^9, 
   3.5774464300347967`*^9}, 3.5774465521986113`*^9, {3.577446588484275*^9, 
   3.5774466383887625`*^9}, 3.5774507241494813`*^9, {3.577451828867772*^9, 
   3.5774518404585924`*^9}, {3.5774520011856747`*^9, 3.577452007347686*^9}, {
   3.5774555692980604`*^9, 3.5774555760684724`*^9}, {3.5774556601526203`*^9, 
   3.5774556771410503`*^9}, {3.5774566838108177`*^9, 3.577456688662427*^9}}],

Cell[BoxData[
 Graphics3DBox[{
   {RGBColor[0, 0, 0], 
    SphereBox[{1.8364951926362667`*^7, 0, 0}, 0.19285714279201807`]}, 
   {RGBColor[0, 0, 0], 
    SphereBox[{-1.8364951926362667`*^7, 0, 0}, 0.19285714279201807`]}, 
   SphereBox[{3.116771092888336*^7, 0, 0}, 1.2802758809663549`*^7], 
   SphereBox[{-3.116771092888336*^7, 0, 0}, 1.2802758809663549`*^7], 
   SphereBox[{0, 0, 0}, 1.8364951733505525`*^7], 
   {RGBColor[0.4, 1., 0.4], Opacity[0.5], 
    SphereBox[{1.8364951926362667`*^7, 1.3, 0}, 1.3]}, 
   {RGBColor[0.4, 1., 0.4], Opacity[0.5], 
    SphereBox[{1.8364951926362667`*^7, -1.3, 0}, 1.3]}, 
   {RGBColor[0.4, 1., 0.4], Opacity[0.5], 
    SphereBox[{-1.8364951926362667`*^7, 1.3, 0}, 1.3]}, 
   {RGBColor[0.4, 1., 0.4], Opacity[0.5], 
    SphereBox[{-1.8364951926362667`*^7, -1.3, 0}, 1.3]}, 
   {RGBColor[0.3, 1., 1.], Opacity[0.5], 
    SphereBox[{1.8364951926362667`*^7, 0, 1.15}, 1.15]}, 
   {RGBColor[0.3, 1., 1.], Opacity[0.5], 
    SphereBox[{1.8364951926362667`*^7, 0, -1.15}, 1.15]}, 
   {RGBColor[0.3, 1., 1.], Opacity[0.5], 
    SphereBox[{-1.8364951926362667`*^7, 0, 1.15}, 1.15]}, 
   {RGBColor[0.3, 1., 1.], Opacity[0.5], 
    SphereBox[{-1.8364951926362667`*^7, 0, -1.15}, 1.15]}},
  AspectRatio->Automatic,
  AutomaticImageSize->True,
  Axes->True,
  AxesLabel->{
    FormBox["xl", TraditionalForm], 
    FormBox["yl", TraditionalForm], 
    FormBox["zl", TraditionalForm]},
  ImageSize->{360., 349.2921617131196},
  PlotRange->All,
  SphericalRegion->True,
  ViewPoint->{1.2310063697582567`, 2.8487895411714437`, -1.3487110393730777`},
  ViewVertical->{0.09844965406807721, 
   0.2726635767396057, -1.5199635220184475`}]], "Output",
 CellChangeTimes->{3.5774504994466867`*^9, 3.5774506319533195`*^9, 
  3.5774507316530943`*^9, 3.5774516279082193`*^9, 3.5774516740531*^9, 
  3.5774520186577053`*^9, 3.577455119658471*^9, 3.577455371302513*^9, 
  3.577455500111939*^9, 3.577455581497282*^9, 3.577455683334261*^9, 
  3.577456051760108*^9, 3.577456355399041*^9, 3.577456445083599*^9, 
  3.5774565699618177`*^9, 3.5774566328455286`*^9, 3.577456693202035*^9, 
  3.5774571570532494`*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"(*", " ", 
  RowBox[{
   RowBox[{"projection", " ", "on", " ", "xy"}], "-", "plane"}], " ", 
  "*)"}]], "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5773413317003937`*^9, 3.5773413362088013`*^9}, {
  3.577352839115512*^9, 3.577352839131112*^9}, {3.5773715990175457`*^9, 
  3.5773716034635534`*^9}, {3.5773716711988726`*^9, 3.5773716785932856`*^9}, {
  3.5773719148557005`*^9, 3.5773720320119066`*^9}, {3.5773720625879602`*^9, 
  3.577372089326407*^9}, {3.5773721384508934`*^9, 3.577372247229884*^9}, {
  3.5773723618276854`*^9, 3.577372407098965*^9}, {3.577372640428575*^9, 
  3.57737267169103*^9}, {3.5773727054338894`*^9, 3.5773727803452206`*^9}, {
  3.577372855490553*^9, 3.5773728739765854`*^9}, {3.5773729151762576`*^9, 
  3.577372935331493*^9}, {3.577373068493327*^9, 3.577373077494543*^9}, {
  3.5774504640814247`*^9, 3.5774504745490427`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"plot3", "=", 
    RowBox[{
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"Graphics", "[", 
        RowBox[{"{", 
         RowBox[{"Circle", "[", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{
             RowBox[{"xc", "[", 
              RowBox[{"[", "i", "]"}], "]"}], ",", 
             RowBox[{"yc", "[", 
              RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", 
           RowBox[{"rr", "[", 
            RowBox[{"[", "i", "]"}], "]"}]}], "]"}], "}"}], "]"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "nc"}], "}"}]}], "]"}], " ", "/.", " ", 
     RowBox[{"t", "[", 
      RowBox[{"[", "2", "]"}], "]"}]}]}], ";", 
   RowBox[{"Show", "[", 
    RowBox[{"plot3", ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"AspectRatio", " ", "\[Rule]", " ", "Automatic"}], ",", 
       RowBox[{"Axes", " ", "->", " ", "True"}], ",", 
       RowBox[{"GridLines", " ", "->", " ", "Automatic"}], ",", " ", 
       RowBox[{"PlotRange", " ", "\[Rule]", " ", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"-", "4"}], ",", "4"}], "}"}], ",", 
          RowBox[{"{", 
           RowBox[{
            RowBox[{"-", "3"}], ",", "3"}], "}"}]}], "}"}]}], ",", " ", 
       RowBox[{"Frame", " ", "->", " ", "True"}]}], "}"}]}], "]"}]}], 
  "\[IndentingNewLine]"}]], "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5773413317003937`*^9, 3.5773413362088013`*^9}, {
   3.577352839115512*^9, 3.577352839131112*^9}, {3.5773715990175457`*^9, 
   3.5773716034635534`*^9}, {3.5773716711988726`*^9, 
   3.5773716785932856`*^9}, {3.5773719148557005`*^9, 
   3.5773720320119066`*^9}, {3.5773720625879602`*^9, 3.577372089326407*^9}, {
   3.5773721384508934`*^9, 3.577372247229884*^9}, {3.5773723618276854`*^9, 
   3.577372407098965*^9}, {3.577372640428575*^9, 3.57737267169103*^9}, {
   3.5773727054338894`*^9, 3.5773727803452206`*^9}, {3.577372855490553*^9, 
   3.5773728739765854`*^9}, {3.5773729151762576`*^9, 3.577372935331493*^9}, {
   3.577373068493327*^9, 3.577373077494543*^9}, 3.577375870741849*^9, {
   3.5774428830877666`*^9, 3.577442884476169*^9}, {3.5774502155573883`*^9, 
   3.5774503400456066`*^9}, {3.5774503829300823`*^9, 3.5774503991853104`*^9}}],

Cell[BoxData[
 GraphicsBox[{CircleBox[{1.8364951926362667`*^7, 0}, 0.19285714279201807`], 
   CircleBox[{-1.8364951926362667`*^7, 0}, 0.19285714279201807`], 
   CircleBox[{3.116771092888336*^7, 0}, 1.2802758809663549`*^7], 
   CircleBox[{-3.116771092888336*^7, 0}, 1.2802758809663549`*^7], 
   CircleBox[{0, 0}, 1.8364951733505525`*^7], 
   CircleBox[{1.8364951926362667`*^7, 1.3}, 1.3], 
   CircleBox[{1.8364951926362667`*^7, -1.3}, 1.3], 
   CircleBox[{-1.8364951926362667`*^7, 1.3}, 1.3], 
   CircleBox[{-1.8364951926362667`*^7, -1.3}, 1.3], 
   CircleBox[{1.8364951926362667`*^7, 0}, 1.15], 
   CircleBox[{1.8364951926362667`*^7, 0}, 1.15], 
   CircleBox[{-1.8364951926362667`*^7, 0}, 1.15], 
   CircleBox[{-1.8364951926362667`*^7, 0}, 1.15]},
  AspectRatio->Automatic,
  Axes->True,
  Frame->True,
  GridLines->Automatic,
  PlotRange->{{-4, 4}, {-3, 3}}]], "Output",
 CellChangeTimes->{3.577450499540287*^9, 3.5774506319689193`*^9, 
  3.5774507317310944`*^9, 3.5774516279082193`*^9, 3.5774516740687003`*^9, 
  3.5774520187513056`*^9, 3.5774551197364707`*^9, 3.577455371318113*^9, 
  3.577455500205539*^9, 3.5774555816064825`*^9, 3.577455683443461*^9, 
  3.577456051869308*^9, 3.5774563555082417`*^9, 3.5774564455671997`*^9, 
  3.577456570429819*^9, 3.577456632954729*^9, 3.577456693280035*^9, 
  3.577457157068849*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{"projection", " ", "on", " ", "xz"}], "-", "plane"}], " ", "*)"}],
   "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"plot4", "=", 
     RowBox[{
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"Graphics", "[", 
         RowBox[{"{", 
          RowBox[{"Circle", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{"xc", "[", 
               RowBox[{"[", "i", "]"}], "]"}], ",", 
              RowBox[{"zc", "[", 
               RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", 
            RowBox[{"rr", "[", 
             RowBox[{"[", "i", "]"}], "]"}]}], "]"}], "}"}], "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "1", ",", "nc"}], "}"}]}], "]"}], " ", "/.", " ", 
      RowBox[{"t", "[", 
       RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], "\n", 
   RowBox[{"Show", "[", 
    RowBox[{"plot4", ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"AspectRatio", " ", "\[Rule]", " ", "Automatic"}], ",", 
       RowBox[{"Axes", " ", "->", " ", "True"}], ",", 
       RowBox[{"GridLines", " ", "->", " ", "Automatic"}], ",", " ", 
       RowBox[{"PlotRange", " ", "\[Rule]", " ", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"-", "4"}], ",", "4"}], "}"}], ",", 
          RowBox[{"{", 
           RowBox[{
            RowBox[{"-", "3"}], ",", "3"}], "}"}]}], "}"}]}], ",", " ", 
       RowBox[{"Frame", " ", "->", " ", "True"}]}], "}"}]}], 
    "]"}]}]}]], "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5773413317003937`*^9, 3.5773413362088013`*^9}, {
   3.577352839115512*^9, 3.577352839131112*^9}, {3.5773715990175457`*^9, 
   3.5773716034635534`*^9}, {3.5773716711988726`*^9, 
   3.5773716785932856`*^9}, {3.5773719148557005`*^9, 
   3.5773720320119066`*^9}, {3.5773720625879602`*^9, 3.577372089326407*^9}, {
   3.5773721384508934`*^9, 3.577372247229884*^9}, {3.5773723618276854`*^9, 
   3.577372407098965*^9}, {3.577372640428575*^9, 3.57737267169103*^9}, {
   3.5773727054338894`*^9, 3.5773727803452206`*^9}, {3.577372855490553*^9, 
   3.5773728739765854`*^9}, {3.5773729151762576`*^9, 3.577372935331493*^9}, {
   3.577373068493327*^9, 3.577373077494543*^9}, 3.577375870741849*^9, {
   3.5774428830877666`*^9, 3.577442884476169*^9}, {3.5774502155573883`*^9, 
   3.5774503400456066`*^9}, {3.5774503829300823`*^9, 3.5774503991853104`*^9}, 
   3.577450442132186*^9, {3.577450485905863*^9, 3.577450491911873*^9}}],

Cell[BoxData[
 GraphicsBox[{CircleBox[{1.8364951926362667`*^7, 0}, 0.19285714279201807`], 
   CircleBox[{-1.8364951926362667`*^7, 0}, 0.19285714279201807`], 
   CircleBox[{3.116771092888336*^7, 0}, 1.2802758809663549`*^7], 
   CircleBox[{-3.116771092888336*^7, 0}, 1.2802758809663549`*^7], 
   CircleBox[{0, 0}, 1.8364951733505525`*^7], 
   CircleBox[{1.8364951926362667`*^7, 0}, 1.3], 
   CircleBox[{1.8364951926362667`*^7, 0}, 1.3], 
   CircleBox[{-1.8364951926362667`*^7, 0}, 1.3], 
   CircleBox[{-1.8364951926362667`*^7, 0}, 1.3], 
   CircleBox[{1.8364951926362667`*^7, 1.15}, 1.15], 
   CircleBox[{1.8364951926362667`*^7, -1.15}, 1.15], 
   CircleBox[{-1.8364951926362667`*^7, 1.15}, 1.15], 
   CircleBox[{-1.8364951926362667`*^7, -1.15}, 1.15]},
  AspectRatio->Automatic,
  Axes->True,
  Frame->True,
  GridLines->Automatic,
  PlotRange->{{-4, 4}, {-3, 3}}]], "Output",
 CellChangeTimes->{3.5774504995558867`*^9, 3.577450631984519*^9, 
  3.5774507317466946`*^9, 3.5774516279238195`*^9, 3.5774516740687003`*^9, 
  3.5774520187825055`*^9, 3.577455119752071*^9, 3.5774553713337126`*^9, 
  3.577455500236739*^9, 3.577455581637682*^9, 3.577455683474661*^9, 
  3.577456051884908*^9, 3.577456355523842*^9, 3.5774564459884005`*^9, 
  3.57745657085102*^9, 3.577456632970329*^9, 3.5774566932956347`*^9, 
  3.5774571570844493`*^9}]
}, Open  ]]
},
WindowToolbars->"EditBar",
WindowSize->{1008, 964},
WindowMargins->{{Automatic, 108}, {0, Automatic}},
PrintingCopies->1,
PrintingPageRange->{Automatic, Automatic},
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, 5831, 115, 272, "Input"],
Cell[CellGroupData[{
Cell[6413, 139, 12468, 309, 672, "Input"],
Cell[18884, 450, 695, 13, 56, "Message"],
Cell[19582, 465, 783, 16, 33, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[20402, 486, 2297, 45, 152, "Input"],
Cell[22702, 533, 500, 7, 31, "Output"],
Cell[23205, 542, 516, 8, 31, "Output"],
Cell[23724, 552, 500, 7, 31, "Output"],
Cell[24227, 561, 504, 7, 31, "Output"],
Cell[24734, 570, 500, 7, 31, "Output"],
Cell[25237, 579, 503, 7, 31, "Output"]
}, Open  ]],
Cell[25755, 589, 2016, 33, 31, "Input"],
Cell[27774, 624, 859, 15, 31, "Input"],
Cell[CellGroupData[{
Cell[28658, 643, 5975, 143, 212, "Input"],
Cell[34636, 788, 2115, 44, 364, "Output"]
}, Open  ]],
Cell[36766, 835, 885, 15, 31, "Input"],
Cell[CellGroupData[{
Cell[37676, 854, 2307, 53, 92, "Input"],
Cell[39986, 909, 1323, 25, 293, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[41346, 939, 2538, 59, 92, "Input"],
Cell[43887, 1000, 1326, 25, 293, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)