(* 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[     43189,       1210]
NotebookOptionsPosition[     41968,       1166]
NotebookOutlinePosition[     42467,       1186]
CellTagsIndexPosition[     42424,       1183]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"(*", " ", 
    RowBox[{"F2", " ", "identical", " ", "to", " ", 
     RowBox[{"F2_F2", ".", "PAS", ".", " ", "ES"}], " ", "23.05", ".2013", 
     " ", "after", " ", "pas"}], " ", "*)"}], "Text"], "\[IndentingNewLine]", 
  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", ",", "xw", ",", "yw", ",", "zw", ",", "xn", ",", "yn", ",", 
      "zn", ",", "ow", ",", "cn", ",", "rr", ",", "d1", ",", "d2", ",", "d3", 
      ",", "d4", ",", "d5", ",", "i", ",", "j", ",", "d", ",", "t", ",", "sq",
       ",", "Ekin", ",", "Epot", ",", "pi", ",", "p", ",", "cc", ",", "dv", 
      ",", "dw", ",", "dcw", ",", "dcc", ",", "tw", ",", "tv", ",", "w4", ",",
       "w5", ",", "q", ",", "x1", ",", "x2", ",", "x3", ",", "x4", ",", 
      "Ekin", ",", "Epot", ",", "func", ",", "nc", ",", "z"}], "]"}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"sq3", "=", 
     RowBox[{"Sqrt", "[", "3", "]"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"nc", "=", "9"}], ";", " ", 
    RowBox[{"z", "=", "9.0"}], ";"}], "\n", 
   RowBox[{
    RowBox[{"a", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
       RowBox[{"{", "nc", "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"b", "=", "a"}], ";", " ", 
    RowBox[{"c", "=", "a"}], ";"}], " ", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"g", "=", 
     RowBox[{"Table", "[", 
      RowBox[{"1.21283096", ",", 
       RowBox[{"{", "nc", "}"}]}], "]"}]}], ";"}], "  ", 
   RowBox[{"(*", " ", 
    RowBox[{"kin", " ", "param"}], " ", "*)"}], "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{"x1", "=", "1.1913455"}], ";", " ", 
     RowBox[{"x2", "=", "0.1593004"}], ";", " ", 
     RowBox[{"x3", "=", "0.8003616"}], ";", " ", 
     RowBox[{"x4", "=", "1.374017"}], ";"}], " ", "*)"}], " ", 
   RowBox[{"(*", " ", 
    RowBox[{"radii", " ", "and", " ", "angle"}], " ", "*)"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"g", "[", 
      RowBox[{"[", "1", "]"}], "]"}], "=", "2.35018989"}], ";", 
    RowBox[{
     RowBox[{"g", "[", 
      RowBox[{"[", "2", "]"}], "]"}], "=", "0.95875665"}], ";", 
    RowBox[{
     RowBox[{"g", "[", 
      RowBox[{"[", "6", "]"}], "]"}], "=", 
     RowBox[{"g", "[", 
      RowBox[{"[", "2", "]"}], "]"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"sig", "=", 
     RowBox[{"Table", "[", 
      RowBox[{"0.30820426", ",", 
       RowBox[{"{", "nc", "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"sig", "[", 
      RowBox[{"[", "1", "]"}], "]"}], "=", "0.31403448"}], ";", " ", 
    RowBox[{
     RowBox[{"sig", "[", 
      RowBox[{"[", "2", "]"}], "]"}], ":=", "0.30821367"}], ";", 
    RowBox[{
     RowBox[{"sig", "[", 
      RowBox[{"[", "6", "]"}], "]"}], "=", 
     RowBox[{"sig", "[", 
      RowBox[{"[", "2", "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", "\n", 
   RowBox[{
    RowBox[{"ss", "=", 
     RowBox[{"Sin", "[", "x4", "]"}]}], ";", 
    RowBox[{"cs", "=", 
     RowBox[{"Cos", "[", "x4", "]"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{
     RowBox[{"c", "[", 
      RowBox[{"[", 
       RowBox[{"1", ",", "1"}], "]"}], "]"}], "=", 
     RowBox[{"-", "2."}]}], ";", 
    RowBox[{
     RowBox[{"c", "[", 
      RowBox[{"[", 
       RowBox[{"1", ",", "2"}], "]"}], "]"}], "=", "0"}], ";", 
    RowBox[{
     RowBox[{"c", "[", 
      RowBox[{"[", 
       RowBox[{"1", ",", "3"}], "]"}], "]"}], "=", "x1"}], ";"}], "\n", 
   RowBox[{"(*", 
    RowBox[{"F", " ", "right"}], "*)"}], "\n", 
   RowBox[{
    RowBox[{
     RowBox[{"c", "[", 
      RowBox[{"[", 
       RowBox[{"2", ",", "1"}], "]"}], "]"}], "=", 
     RowBox[{"-", "2."}]}], ";", 
    RowBox[{
     RowBox[{"c", "[", 
      RowBox[{"[", 
       RowBox[{"2", ",", "2"}], "]"}], "]"}], "=", "z"}], ";", 
    RowBox[{
     RowBox[{"c", "[", 
      RowBox[{"[", 
       RowBox[{"2", ",", "3"}], "]"}], "]"}], "=", "x2"}], ";"}], 
   "\[IndentingNewLine]", "\n", 
   RowBox[{
    RowBox[{
     RowBox[{"c", "[", 
      RowBox[{"[", 
       RowBox[{"3", ",", "1"}], "]"}], "]"}], "=", 
     RowBox[{"-", "2."}]}], ";", 
    RowBox[{
     RowBox[{"c", "[", 
      RowBox[{"[", 
       RowBox[{"3", ",", "2"}], "]"}], "]"}], "=", "0"}], ";", 
    RowBox[{
     RowBox[{"c", "[", 
      RowBox[{"[", 
       RowBox[{"3", ",", "3"}], "]"}], "]"}], "=", "x3"}], ";"}], "\n", 
   RowBox[{
    RowBox[{
     RowBox[{"c", "[", 
      RowBox[{"[", 
       RowBox[{"4", ",", "1"}], "]"}], "]"}], "=", 
     RowBox[{"-", "2."}]}], ";", 
    RowBox[{
     RowBox[{"c", "[", 
      RowBox[{"[", 
       RowBox[{"4", ",", "2"}], "]"}], "]"}], "=", "0"}], ";", 
    RowBox[{
     RowBox[{"c", "[", 
      RowBox[{"[", 
       RowBox[{"4", ",", "3"}], "]"}], "]"}], "=", "x3"}], ";"}], "\n", 
   RowBox[{
    RowBox[{
     RowBox[{"c", "[", 
      RowBox[{"[", 
       RowBox[{"5", ",", "1"}], "]"}], "]"}], "=", 
     RowBox[{"-", "2."}]}], ";", 
    RowBox[{
     RowBox[{"c", "[", 
      RowBox[{"[", 
       RowBox[{"5", ",", "2"}], "]"}], "]"}], "=", "0"}], ";", 
    RowBox[{
     RowBox[{"c", "[", 
      RowBox[{"[", 
       RowBox[{"5", ",", "3"}], "]"}], "]"}], "=", "x3"}], ";"}], "\n", 
   RowBox[{"(*", 
    RowBox[{"F", " ", "left"}], "*)"}], "\n", 
   RowBox[{
    RowBox[{
     RowBox[{"c", "[", 
      RowBox[{"[", 
       RowBox[{"6", ",", "1"}], "]"}], "]"}], "=", 
     RowBox[{"-", "2."}]}], ";", 
    RowBox[{
     RowBox[{"c", "[", 
      RowBox[{"[", 
       RowBox[{"6", ",", "2"}], "]"}], "]"}], "=", "z"}], ";", 
    RowBox[{
     RowBox[{"c", "[", 
      RowBox[{"[", 
       RowBox[{"6", ",", "3"}], "]"}], "]"}], "=", "x2"}], ";"}], "\n", 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"c", "[", 
      RowBox[{"[", 
       RowBox[{"7", ",", "1"}], "]"}], "]"}], "=", 
     RowBox[{"-", "2."}]}], ";", 
    RowBox[{
     RowBox[{"c", "[", 
      RowBox[{"[", 
       RowBox[{"7", ",", "2"}], "]"}], "]"}], "=", "0"}], ";", 
    RowBox[{
     RowBox[{"c", "[", 
      RowBox[{"[", 
       RowBox[{"7", ",", "3"}], "]"}], "]"}], "=", "x3"}], ";"}], "\n", 
   RowBox[{
    RowBox[{
     RowBox[{"c", "[", 
      RowBox[{"[", 
       RowBox[{"8", ",", "1"}], "]"}], "]"}], "=", 
     RowBox[{"-", "2."}]}], ";", 
    RowBox[{
     RowBox[{"c", "[", 
      RowBox[{"[", 
       RowBox[{"8", ",", "2"}], "]"}], "]"}], "=", "0"}], ";", 
    RowBox[{
     RowBox[{"c", "[", 
      RowBox[{"[", 
       RowBox[{"8", ",", "3"}], "]"}], "]"}], "=", "x3"}], ";"}], "\n", 
   RowBox[{
    RowBox[{
     RowBox[{"c", "[", 
      RowBox[{"[", 
       RowBox[{"9", ",", "1"}], "]"}], "]"}], "=", 
     RowBox[{"-", "2."}]}], ";", 
    RowBox[{
     RowBox[{"c", "[", 
      RowBox[{"[", 
       RowBox[{"9", ",", "2"}], "]"}], "]"}], "=", "0"}], ";", 
    RowBox[{
     RowBox[{"c", "[", 
      RowBox[{"[", 
       RowBox[{"9", ",", "3"}], "]"}], "]"}], "=", "x3"}], ";"}], "\n", "\n", 
   RowBox[{"(*", 
    RowBox[{
    "coord", " ", "with", " ", "zero", " ", "value", " ", "not", " ", 
     "explicitely", " ", 
     RowBox[{"written", "!"}]}], "*)"}], "\n", 
   RowBox[{"(*", "nuclei", "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"b", "[", 
      RowBox[{"[", 
       RowBox[{"2", ",", "1"}], "]"}], "]"}], "=", 
     RowBox[{"x1", "+", "x2"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{
     RowBox[{"b", "[", 
      RowBox[{"[", 
       RowBox[{"6", ",", "1"}], "]"}], "]"}], "=", 
     RowBox[{"-", 
      RowBox[{"b", "[", 
       RowBox[{"[", 
        RowBox[{"2", ",", "1"}], "]"}], "]"}]}]}], ";"}], "\n", 
   RowBox[{"(*", "clouds", "*)"}], "\n", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"a", "[", 
      RowBox[{"[", 
       RowBox[{"2", ",", "1"}], "]"}], "]"}], "=", 
     RowBox[{"b", "[", 
      RowBox[{"[", 
       RowBox[{"2", ",", "1"}], "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
   "\n", 
   RowBox[{
    RowBox[{
     RowBox[{"a", "[", 
      RowBox[{"[", 
       RowBox[{"3", ",", "1"}], "]"}], "]"}], "=", 
     RowBox[{
      RowBox[{"b", "[", 
       RowBox[{"[", 
        RowBox[{"2", ",", "1"}], "]"}], "]"}], "+", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"x2", "+", "x3"}], ")"}], "*", "cs"}]}]}], ";", 
    RowBox[{
     RowBox[{"a", "[", 
      RowBox[{"[", 
       RowBox[{"3", ",", "2"}], "]"}], "]"}], "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"x2", "+", "x3"}], ")"}], "*", "ss"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{
     RowBox[{"a", "[", 
      RowBox[{"[", 
       RowBox[{"4", ",", "1"}], "]"}], "]"}], "=", 
     RowBox[{"a", "[", 
      RowBox[{"[", 
       RowBox[{"3", ",", "1"}], "]"}], "]"}]}], ";", 
    RowBox[{
     RowBox[{"a", "[", 
      RowBox[{"[", 
       RowBox[{"4", ",", "2"}], "]"}], "]"}], "=", 
     RowBox[{
      RowBox[{"-", "0.5"}], "*", 
      RowBox[{"a", "[", 
       RowBox[{"[", 
        RowBox[{"3", ",", "2"}], "]"}], "]"}]}]}], ";", 
    RowBox[{
     RowBox[{"a", "[", 
      RowBox[{"[", 
       RowBox[{"4", ",", "3"}], "]"}], "]"}], "=", 
     RowBox[{"0.5", "*", "sq3", "*", 
      RowBox[{"a", "[", 
       RowBox[{"[", 
        RowBox[{"3", ",", "2"}], "]"}], "]"}]}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{
     RowBox[{"a", "[", 
      RowBox[{"[", 
       RowBox[{"5", ",", "1"}], "]"}], "]"}], "=", 
     RowBox[{"a", "[", 
      RowBox[{"[", 
       RowBox[{"3", ",", "1"}], "]"}], "]"}]}], ";", 
    RowBox[{
     RowBox[{"a", "[", 
      RowBox[{"[", 
       RowBox[{"5", ",", "2"}], "]"}], "]"}], "=", 
     RowBox[{"a", "[", 
      RowBox[{"[", 
       RowBox[{"4", ",", "2"}], "]"}], "]"}]}], ";", 
    RowBox[{
     RowBox[{"a", "[", 
      RowBox[{"[", 
       RowBox[{"5", ",", "3"}], "]"}], "]"}], "=", 
     RowBox[{"-", 
      RowBox[{"a", "[", 
       RowBox[{"[", 
        RowBox[{"4", ",", "3"}], "]"}], "]"}]}]}], ";"}], "\n", 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"a", "[", 
      RowBox[{"[", 
       RowBox[{"6", ",", "1"}], "]"}], "]"}], "=", 
     RowBox[{"-", 
      RowBox[{"b", "[", 
       RowBox[{"[", 
        RowBox[{"2", ",", "1"}], "]"}], "]"}]}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{
     RowBox[{"a", "[", 
      RowBox[{"[", 
       RowBox[{"7", ",", "1"}], "]"}], "]"}], "=", 
     RowBox[{"-", 
      RowBox[{"a", "[", 
       RowBox[{"[", 
        RowBox[{"3", ",", "1"}], "]"}], "]"}]}]}], ";", 
    RowBox[{
     RowBox[{"a", "[", 
      RowBox[{"[", 
       RowBox[{"7", ",", "2"}], "]"}], "]"}], "=", 
     RowBox[{"-", 
      RowBox[{"a", "[", 
       RowBox[{"[", 
        RowBox[{"3", ",", "2"}], "]"}], "]"}]}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{
     RowBox[{"a", "[", 
      RowBox[{"[", 
       RowBox[{"8", ",", "1"}], "]"}], "]"}], "=", 
     RowBox[{"-", 
      RowBox[{"a", "[", 
       RowBox[{"[", 
        RowBox[{"3", ",", "1"}], "]"}], "]"}]}]}], ";", 
    RowBox[{
     RowBox[{"a", "[", 
      RowBox[{"[", 
       RowBox[{"8", ",", "2"}], "]"}], "]"}], "=", 
     RowBox[{"-", 
      RowBox[{"a", "[", 
       RowBox[{"[", 
        RowBox[{"4", ",", "2"}], "]"}], "]"}]}]}], ";", 
    RowBox[{
     RowBox[{"a", "[", 
      RowBox[{"[", 
       RowBox[{"8", ",", "3"}], "]"}], "]"}], "=", 
     RowBox[{"a", "[", 
      RowBox[{"[", 
       RowBox[{"4", ",", "3"}], "]"}], "]"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{
     RowBox[{"a", "[", 
      RowBox[{"[", 
       RowBox[{"9", ",", "1"}], "]"}], "]"}], "=", 
     RowBox[{"-", 
      RowBox[{"a", "[", 
       RowBox[{"[", 
        RowBox[{"3", ",", "1"}], "]"}], "]"}]}]}], ";", 
    RowBox[{
     RowBox[{"a", "[", 
      RowBox[{"[", 
       RowBox[{"9", ",", "2"}], "]"}], "]"}], "=", 
     RowBox[{"-", 
      RowBox[{"a", "[", 
       RowBox[{"[", 
        RowBox[{"4", ",", "2"}], "]"}], "]"}]}]}], ";", 
    RowBox[{
     RowBox[{"a", "[", 
      RowBox[{"[", 
       RowBox[{"9", ",", "3"}], "]"}], "]"}], "=", 
     RowBox[{"a", "[", 
      RowBox[{"[", 
       RowBox[{"5", ",", "3"}], "]"}], "]"}]}], ";"}], "\n", "\n", 
   RowBox[{
    RowBox[{"vnn", "=", 
     RowBox[{"z", "*", 
      RowBox[{"z", "/", 
       RowBox[{"(", 
        RowBox[{"2", "*", 
         RowBox[{"(", 
          RowBox[{"x1", "+", "x2"}], ")"}]}], ")"}]}]}]}], ";", "   ", 
    RowBox[{"vee", "=", "0.0"}], ";", 
    RowBox[{"vne", "=", "0.0"}], ";", 
    RowBox[{"ekin", "=", "0.0"}], ";"}], "\n", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{
      RowBox[{"WI", "/", "WJ"}], " ", 
      RowBox[{"EXCL", ".", "I"}]}], "=", "J"}], "*)"}], "\n", 
   RowBox[{"For", "[", " ", 
    RowBox[{
     RowBox[{"i", "=", "1"}], ",", 
     RowBox[{"i", "<", " ", "nc"}], " ", ",", 
     RowBox[{"i", "++"}], ",", "\n", 
     RowBox[{"For", "[", " ", 
      RowBox[{
       RowBox[{"j", "=", 
        RowBox[{"i", "+", "1"}]}], ",", 
       RowBox[{"j", "<", 
        RowBox[{"nc", "+", "1"}]}], ",", 
       RowBox[{"j", "++"}], ",", "\n", 
       RowBox[{
        RowBox[{"d", "=", 
         RowBox[{"Sqrt", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"(", 
             RowBox[{
              RowBox[{"a", "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", "1"}], "]"}], "]"}], "-", 
              RowBox[{"a", "[", 
               RowBox[{"[", 
                RowBox[{"j", ",", "1"}], "]"}], "]"}]}], ")"}], "^", "2"}], 
           "+", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{
              RowBox[{"a", "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", "2"}], "]"}], "]"}], "-", 
              RowBox[{"a", "[", 
               RowBox[{"[", 
                RowBox[{"j", ",", "2"}], "]"}], "]"}]}], ")"}], "^", "2"}], 
           "+", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{
              RowBox[{"a", "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", "3"}], "]"}], "]"}], "-", 
              RowBox[{"a", "[", 
               RowBox[{"[", 
                RowBox[{"j", ",", "3"}], "]"}], "]"}]}], ")"}], "^", "2"}]}], 
          "]"}]}], ";", "\n", 
        RowBox[{"If", "[", " ", 
         RowBox[{
          RowBox[{"d", ">", "0"}], " ", ",", 
          RowBox[{"vee", "=", 
           RowBox[{"vee", "+", 
            RowBox[{
             RowBox[{"c", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "1"}], "]"}], "]"}], "*", 
             RowBox[{
              RowBox[{"c", "[", 
               RowBox[{"[", 
                RowBox[{"j", ",", "1"}], "]"}], "]"}], "/", "d"}]}]}]}]}], 
         "]"}]}]}], "]"}]}], "]"}], "\n", "\n", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{
      RowBox[{"KJ", "/", "WI"}], " ", "I"}], "=", "J"}], "*)"}], 
   "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{"vne", "=", 
      RowBox[{"1.5", "*", 
       RowBox[{"c", "[", 
        RowBox[{"[", 
         RowBox[{"2", ",", "1"}], "]"}], "]"}], "*", 
       RowBox[{
        RowBox[{"c", "[", 
         RowBox[{"[", 
          RowBox[{"2", ",", "2"}], "]"}], "]"}], "/", 
        RowBox[{"c", "[", 
         RowBox[{"[", 
          RowBox[{"2", ",", "3"}], "]"}], "]"}]}]}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"vne", "=", 
      RowBox[{"vne", "+", 
       RowBox[{"1.5", "*", 
        RowBox[{"c", "[", 
         RowBox[{"[", 
          RowBox[{"6", ",", "1"}], "]"}], "]"}], "*", 
        RowBox[{
         RowBox[{"c", "[", 
          RowBox[{"[", 
           RowBox[{"6", ",", "2"}], "]"}], "]"}], "/", 
         RowBox[{"c", "[", 
          RowBox[{"[", 
           RowBox[{"6", ",", "3"}], "]"}], "]"}]}]}]}]}], ";"}], "*)"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"vne", "=", 
     RowBox[{
      RowBox[{"-", "2"}], "*", "3.0", "*", 
      RowBox[{"z", "/", "x2"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{
      RowBox[{"KJ", "/", "WI"}], " ", 
      RowBox[{"I", "<>", "J"}], " ", "inc"}], ";"}], "*)"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"For", "[", " ", 
     RowBox[{
      RowBox[{"i", "=", "1"}], ",", 
      RowBox[{"i", "<", 
       RowBox[{"nc", "+", "1"}]}], ",", 
      RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
      RowBox[{"For", "[", 
       RowBox[{
        RowBox[{"j", "=", "1"}], ",", 
        RowBox[{"j", "<", 
         RowBox[{"nc", "+", "1"}]}], ",", 
        RowBox[{"j", "++"}], ",", "\n", 
        RowBox[{
         RowBox[{"d", "=", 
          RowBox[{"Sqrt", "[", 
           RowBox[{
            RowBox[{
             RowBox[{"(", 
              RowBox[{
               RowBox[{"a", "[", 
                RowBox[{"[", 
                 RowBox[{"i", ",", "1"}], "]"}], "]"}], "-", 
               RowBox[{"b", "[", 
                RowBox[{"[", 
                 RowBox[{"j", ",", "1"}], "]"}], "]"}]}], ")"}], "^", "2"}], 
            "+", 
            RowBox[{
             RowBox[{"a", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "2"}], "]"}], "]"}], "^", "2"}], "+", 
            RowBox[{
             RowBox[{"a", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "3"}], "]"}], "]"}], "^", "2"}]}], "]"}]}], 
         ";", "\n", 
         RowBox[{"If", "[", " ", 
          RowBox[{
           RowBox[{"d", ">", "0"}], ",", 
           RowBox[{"vne", "=", 
            RowBox[{"vne", "+", 
             RowBox[{
              RowBox[{"c", "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", "1"}], "]"}], "]"}], "*", 
              RowBox[{
               RowBox[{"c", "[", 
                RowBox[{"[", 
                 RowBox[{"j", ",", "2"}], "]"}], "]"}], "/", "d"}]}]}]}]}], 
          "]"}]}]}], "]"}]}], "]"}], ";"}], "\n", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{
      RowBox[{
       RowBox[{"d", "=", 
        RowBox[{"Sqrt", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"(", 
            RowBox[{
             RowBox[{"a", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "1"}], "]"}], "]"}], "-", 
             RowBox[{"b", "[", 
              RowBox[{"[", 
               RowBox[{"6", ",", "1"}], "]"}], "]"}]}], ")"}], "^", "2"}], 
          "+", 
          RowBox[{
           RowBox[{"a", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "2"}], "]"}], "]"}], "^", "2"}], "+", 
          RowBox[{
           RowBox[{"a", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "3"}], "]"}], "]"}], "^", "2"}]}], "]"}]}], 
       ";", "\n", 
       RowBox[{"If", "[", " ", 
        RowBox[{
         RowBox[{"d", ">", "0"}], " ", ",", 
         RowBox[{"vne", "=", 
          RowBox[{"vne", "+", 
           RowBox[{
            RowBox[{"c", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "1"}], "]"}], "]"}], "*", 
            RowBox[{"z", "/", "d"}]}]}]}]}], "]"}]}], "]"}], ";"}], "*)"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"For", "[", " ", 
     RowBox[{
      RowBox[{"i", "=", "1"}], ",", 
      RowBox[{"i", "<", " ", 
       RowBox[{"nc", "+", "1"}]}], ",", 
      RowBox[{"i", "++"}], ",", "\n", "      ", 
      RowBox[{"ekin", "=", 
       RowBox[{"ekin", "-", 
        RowBox[{"1.125", "*", 
         RowBox[{"g", "[", 
          RowBox[{"[", "i", "]"}], "]"}], "*", 
         RowBox[{
          RowBox[{"c", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "1"}], "]"}], "]"}], "/", 
          RowBox[{
           RowBox[{"c", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "3"}], "]"}], "]"}], "^", "2"}]}]}]}]}]}], 
     "]"}], ";"}], "\n", 
   RowBox[{
    RowBox[{"For", "[", " ", 
     RowBox[{
      RowBox[{"i", "=", "1"}], ",", 
      RowBox[{"i", "<", " ", 
       RowBox[{"nc", "+", "1"}]}], ",", 
      RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
      RowBox[{"vee", "=", 
       RowBox[{"vee", "+", 
        RowBox[{"3.0", "*", 
         RowBox[{
          RowBox[{"sig", "[", 
           RowBox[{"[", "i", "]"}], "]"}], "/", 
          RowBox[{"c", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "3"}], "]"}], "]"}]}]}]}]}]}], "]"}], ";"}], 
   "\n", 
   RowBox[{
    RowBox[{"epot", "=", 
     RowBox[{"vne", "+", "vee", "+", "vnn"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"func", "=", 
     RowBox[{"ekin", "+", "epot"}]}], ";"}], "\n", "\n", 
   RowBox[{"t", "=", 
    RowBox[{"NMinimize", "[", 
     RowBox[{"func", ",", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"x1", ",", "1.190", ",", "1.192"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"x2", ",", "0.158", ",", "0.160"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"x3", ",", "0.79", ",", "0.81"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"x4", ",", "1.36", ",", "1.39"}], "}"}]}], "}"}]}], 
     "]"}]}]}]}]], "Input",
 CellChangeTimes->{{3.577942230718297*^9, 3.577942230718297*^9}, {
   3.577942264180356*^9, 3.577942265131957*^9}, {3.577942319014452*^9, 
   3.5779424054386044`*^9}, {3.577942440101865*^9, 3.5779424567158937`*^9}, {
   3.5779425225012093`*^9, 3.5779425465408516`*^9}, {3.5779425971161404`*^9, 
   3.577942597365741*^9}, {3.577942660701852*^9, 3.577942690669505*^9}, {
   3.5779427231487617`*^9, 3.577942868431817*^9}, {3.5779429100370903`*^9, 
   3.577943391641136*^9}, 3.577943422747591*^9, {3.577943618964735*^9, 
   3.5779436491819887`*^9}, {3.5779437045308857`*^9, 3.577943727384926*^9}, {
   3.577943802077857*^9, 3.5779438458203335`*^9}, {3.577943880187194*^9, 
   3.577943888221208*^9}, {3.577943985190978*^9, 3.577944388373687*^9}, {
   3.5779445035798893`*^9, 3.5779445527355757`*^9}, {3.5779445893488398`*^9, 
   3.577944718829067*^9}, {3.577944766611951*^9, 3.5779448631605206`*^9}, {
   3.5779449296478376`*^9, 3.5779449301314383`*^9}, {3.5779450076947746`*^9, 
   3.577945070765685*^9}, {3.57794511340056*^9, 3.577945117113367*^9}, {
   3.5779451860030875`*^9, 3.57794522158675*^9}, {3.577949900771303*^9, 
   3.577949901395304*^9}, {3.577950305217613*^9, 3.5779503055452137`*^9}, {
   3.577950404168587*^9, 3.5779504345418406`*^9}, {3.577959366380599*^9, 
   3.5779593666925993`*^9}, {3.577959431978714*^9, 3.5779594354419203`*^9}, 
   3.577959523098474*^9, {3.5779596062622204`*^9, 3.577959676743144*^9}, 
   3.577959795084952*^9, 3.5779598664082775`*^9, {3.577959911819957*^9, 
   3.5779599125219584`*^9}, {3.577959989804494*^9, 3.5779600126741343`*^9}, {
   3.5779605201742253`*^9, 3.5779606130099883`*^9}, {3.5779608636868286`*^9, 
   3.577960865153231*^9}, {3.5779609090829086`*^9, 3.5779609120001135`*^9}, 
   3.5779613402052655`*^9, {3.577966477736745*^9, 3.5779665028683896`*^9}, {
   3.577966969605609*^9, 3.5779670293069143`*^9}, {3.577967059570967*^9, 
   3.577967065608178*^9}, {3.5779672243384566`*^9, 3.577967240578085*^9}, {
   3.5779673600118947`*^9, 3.577967481785709*^9}, {3.577967692776079*^9, 
   3.5779677259105377`*^9}, {3.577967762071401*^9, 3.5779678096202846`*^9}, {
   3.577967920193279*^9, 3.5779679205832796`*^9}, {3.577968291948732*^9, 
   3.5779685312531524`*^9}, 3.577968612373295*^9, 3.5779687868908014`*^9, {
   3.5779688258596697`*^9, 3.577968838152491*^9}, {3.577968880288165*^9, 
   3.577968889726182*^9}, {3.5779689372438655`*^9, 3.5779689429378757`*^9}, {
   3.577974893307839*^9, 3.577974913806275*^9}, {3.5779749802311916`*^9, 
   3.5779749826491957`*^9}, 3.5779750560785246`*^9, 3.577975158898305*^9, {
   3.5779751952463694`*^9, 3.5779751967907715`*^9}, 3.5779752512816677`*^9, {
   3.5779753532278466`*^9, 3.5779753576114545`*^9}, {3.5779754239583707`*^9, 
   3.5779754248631725`*^9}, {3.5779754901648874`*^9, 
   3.5779755171841345`*^9}, {3.577975665259595*^9, 3.577975678067217*^9}, 
   3.5779759903485656`*^9, {3.577976309883927*^9, 3.5779763118651304`*^9}, {
   3.57797643688375*^9, 3.5779764986286583`*^9}, {3.577976617220067*^9, 
   3.5779766469069185`*^9}, {3.5779795513378515`*^9, 3.57797955583066*^9}, {
   3.577979610758356*^9, 3.5779796280743866`*^9}, {3.577979870670412*^9, 
   3.577979874539219*^9}, {3.5779802053222003`*^9, 3.5779802246974344`*^9}, {
   3.5779803659712825`*^9, 3.5779803728508945`*^9}, {3.577980417357773*^9, 
   3.5779804285585923`*^9}, {3.5779804859198933`*^9, 
   3.5779805058411283`*^9}, {3.577980577414054*^9, 3.5779806325289507`*^9}, {
   3.5779808443149223`*^9, 3.577980876482179*^9}, {3.5779809082594347`*^9, 
   3.5779809259030657`*^9}, {3.5779810200648313`*^9, 
   3.5779810469904785`*^9}, {3.577981174910703*^9, 3.577981193303136*^9}, {
   3.57798124718563*^9, 3.577981257216448*^9}, {3.5779815034004803`*^9, 
   3.5779815145389*^9}, 3.5779826478408904`*^9, {3.5779827432662582`*^9, 
   3.577982829284809*^9}, {3.577982965317048*^9, 3.577983009231125*^9}, {
   3.5780217722811317`*^9, 3.578021781672348*^9}, {3.5780246076329117`*^9, 
   3.578024624262541*^9}, {3.5780249676347437`*^9, 3.578024985808776*^9}, 
   3.5780250527796936`*^9, {3.5780252283580017`*^9, 3.5780253028169327`*^9}, {
   3.5780253387593956`*^9, 3.578025371488253*^9}, {3.5780255027000837`*^9, 
   3.578025515897707*^9}, {3.578025582665824*^9, 3.5780256468599367`*^9}, {
   3.578025716576459*^9, 3.5780257213500676`*^9}, {3.578025810270224*^9, 
   3.5780258121110272`*^9}, {3.578026140569604*^9, 3.578026141739606*^9}, {
   3.5780262592546124`*^9, 3.578026333105142*^9}, {3.5780263742736144`*^9, 
   3.578026385053233*^9}, 3.578064864045664*^9, {3.578064902281331*^9, 
   3.5780649961934958`*^9}, 3.5780650732576313`*^9, 3.578065104410886*^9, {
   3.578065135564141*^9, 3.5780652004446545`*^9}, {3.5780652326119113`*^9, 
   3.5780653812333727`*^9}, {3.5780656119265776`*^9, 3.578065683327903*^9}, {
   3.578065723388773*^9, 3.578065828174157*^9}, {3.5780660532825527`*^9, 
   3.57806630253979*^9}, {3.5782826346786675`*^9, 3.578282683865554*^9}, {
   3.578282719184016*^9, 3.5782829078507476`*^9}, {3.5782830043057165`*^9, 
   3.578283004898518*^9}, {3.5782831659844007`*^9, 3.5782831812568274`*^9}, {
   3.5782832506145496`*^9, 3.5782834066148233`*^9}, 3.5782834589685154`*^9, 
   3.5782835673419056`*^9, 3.578747458898316*^9, {3.5787474916427736`*^9, 
   3.578747519052022*^9}, {3.578747587816942*^9, 3.578747595960157*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", "157.59235032398846`"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"x1", "\[Rule]", "26220.7566254865`"}], ",", 
     RowBox[{"x2", "\[Rule]", "0.16545908060765693`"}], ",", 
     RowBox[{"x3", "\[Rule]", "75953.89838626009`"}], ",", 
     RowBox[{"x4", "\[Rule]", 
      RowBox[{"-", "75292.53646961565`"}]}]}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.578747614430589*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"ekin", " ", "/.", " ", 
  RowBox[{"t", "[", 
   RowBox[{"[", "2", "]"}], "]"}]}], "\[IndentingNewLine]", 
 RowBox[{"vne", " ", "/.", " ", 
  RowBox[{"t", "[", 
   RowBox[{"[", "2", "]"}], "]"}]}], "\[IndentingNewLine]", 
 RowBox[{"vee", " ", "/.", " ", 
  RowBox[{"t", "[", 
   RowBox[{"[", "2", "]"}], "]"}]}], "\[IndentingNewLine]", 
 RowBox[{"vnn", " ", "/.", " ", 
  RowBox[{"t", "[", 
   RowBox[{"[", "2", "]"}], "]"}]}], "\[IndentingNewLine]", 
 RowBox[{"vir", " ", "=", 
  RowBox[{
   RowBox[{
    RowBox[{"-", "epot"}], "/", "ekin"}], " ", "/.", " ", 
   RowBox[{"t", "[", 
    RowBox[{"[", "2", "]"}], "]"}]}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"2", "*", 
   RowBox[{"(", 
    RowBox[{"x1", "+", "x2"}], ")"}], "*", "0.529177"}], " ", "/.", " ", 
  RowBox[{"t", "[", 
   RowBox[{"[", "2", "]"}], "]"}]}]}], "Input",
 CellChangeTimes->{{3.577960629405617*^9, 3.577960709823758*^9}, {
  3.57796078464149*^9, 3.577960831893973*^9}, {3.5779801723437424`*^9, 
  3.5779801737945447`*^9}, {3.5779832465543423`*^9, 3.577983249752348*^9}, {
  3.577984706342506*^9, 3.5779847140489197`*^9}, {3.5779847533609886`*^9, 
  3.5779847556073923`*^9}, {3.578283501993391*^9, 3.5782835326318445`*^9}}],

Cell[BoxData["157.594003882888`"], "Output",
 CellChangeTimes->{3.578747614461789*^9}],

Cell[BoxData[
 RowBox[{"-", "326.36468063089825`"}]], "Output",
 CellChangeTimes->{3.578747614461789*^9}],

Cell[BoxData["11.1767818559144`"], "Output",
 CellChangeTimes->{3.578747614461789*^9}],

Cell[BoxData["0.0015445681074593921`"], "Output",
 CellChangeTimes->{3.578747614461789*^9}],

Cell[BoxData["1.9999895074757998`"], "Output",
 CellChangeTimes->{3.578747614461789*^9}],

Cell[BoxData["27751.017771889936`"], "Output",
 CellChangeTimes->{3.578747614461789*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"rr", "=", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
     "x1", ",", "x2", ",", "x3", ",", "x3", ",", "x3", ",", "x2", ",", "x3", 
      ",", "x3", ",", "x3"}], "}"}], "   ", "/.", " ", 
    RowBox[{"t", "[", 
     RowBox[{"[", "2", "]"}], "]"}]}]}], " ", " ", ";"}]], "Input",
 CellChangeTimes->{{3.577960992699055*^9, 3.5779610311999226`*^9}, {
   3.5779672624493237`*^9, 3.577967267160532*^9}, 3.577967381009532*^9, 
   3.5779676939772816`*^9, {3.5779678483239527`*^9, 3.57796785264516*^9}, 
   3.578283440201682*^9}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"plot0", "=", 
   RowBox[{
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"Graphics3D", "[", 
       RowBox[{"Sphere", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"a", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "1"}], "]"}], "]"}], ",", 
           RowBox[{"a", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "2"}], "]"}], "]"}], ",", 
           RowBox[{"a", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "3"}], "]"}], "]"}]}], "}"}], ",", 
         RowBox[{"rr", "[", 
          RowBox[{"[", "i", "]"}], "]"}]}], "]"}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "1"}], "}"}]}], "]"}], "     ", "/.", " ", 
    RowBox[{"t", "[", 
     RowBox[{"[", "2", "]"}], "]"}]}]}], " ", " ", " ", ";"}], "\n", 
 RowBox[{
  RowBox[{"plot1", "=", 
   RowBox[{
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"Graphics3D", "[", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"Darker", "[", 
          RowBox[{"Green", ",", "0.3"}], "]"}], ",", 
         RowBox[{"Sphere", "[", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{
             RowBox[{"a", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "1"}], "]"}], "]"}], ",", 
             RowBox[{"a", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "2"}], "]"}], "]"}], ",", 
             RowBox[{"a", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "3"}], "]"}], "]"}]}], "}"}], ",", 
           RowBox[{"rr", "[", 
            RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}], "}"}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "2", ",", "2"}], "}"}]}], "]"}], "      ", "/.", " ", 
    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[{"a", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "1"}], "]"}], "]"}], ",", 
             RowBox[{"a", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "2"}], "]"}], "]"}], ",", 
             RowBox[{"a", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "3"}], "]"}], "]"}]}], "}"}], ",", 
           RowBox[{"rr", "[", 
            RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}], "}"}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "3", ",", "5"}], "}"}]}], "]"}], "     ", "/.", " ", 
    RowBox[{"t", "[", 
     RowBox[{"[", "2", "]"}], "]"}]}]}], " ", "  ", 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"plot2a", "=", 
   RowBox[{
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"Graphics3D", "[", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"Darker", "[", 
          RowBox[{"Green", ",", "0.3"}], "]"}], ",", 
         RowBox[{"Sphere", "[", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{
             RowBox[{"a", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "1"}], "]"}], "]"}], ",", 
             RowBox[{"a", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "2"}], "]"}], "]"}], ",", 
             RowBox[{"a", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "3"}], "]"}], "]"}]}], "}"}], ",", 
           RowBox[{"rr", "[", 
            RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}], "}"}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "6", ",", "6"}], "}"}]}], "]"}], "      ", "/.", " ", 
    RowBox[{"t", "[", 
     RowBox[{"[", "2", "]"}], "]"}]}]}], " ", " ", " ", ";"}], "\n", 
 RowBox[{
  RowBox[{"plot2b", "=", 
   RowBox[{
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"Graphics3D", "[", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"Opacity", "[", "0.5", "]"}], ",", 
         RowBox[{"Lighter", "[", 
          RowBox[{"Green", ",", "0.4"}], "]"}], ",", 
         RowBox[{"Sphere", "[", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{
             RowBox[{"a", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "1"}], "]"}], "]"}], ",", 
             RowBox[{"a", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "2"}], "]"}], "]"}], ",", 
             RowBox[{"a", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "3"}], "]"}], "]"}]}], "}"}], ",", 
           RowBox[{"rr", "[", 
            RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}], "}"}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "7", ",", "nc"}], "}"}]}], "]"}], "     ", "/.", " ", 
    RowBox[{"t", "[", 
     RowBox[{"[", "2", "]"}], "]"}]}]}], " ", "  ", ";"}]}], "Input",
 CellChangeTimes->{{3.5779611162980723`*^9, 3.577961138418911*^9}, {
   3.5779611697593665`*^9, 3.5779612099294367`*^9}, {3.5779660815272493`*^9, 
   3.577966088110461*^9}, {3.577966305465643*^9, 3.5779664035586147`*^9}, {
   3.577967381009532*^9, 3.577967381025132*^9}, {3.5779676939772816`*^9, 
   3.5779676941800823`*^9}, {3.5779680122490406`*^9, 3.577968087534773*^9}, 
   3.5779682529174633`*^9, {3.578066366655903*^9, 3.5780663683719063`*^9}, {
   3.5782834402172823`*^9, 3.5782834402328825`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Show", "[", 
  RowBox[{
  "plot0", ",", "plot1", ",", "plot2", ",", "plot2a", ",", "plot2b", ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"AspectRatio", " ", "\[Rule]", " ", "Automatic"}], ",", 
     RowBox[{"Axes", " ", "->", " ", "True"}], ",", 
     RowBox[{"PlotRange", " ", "\[Rule]", " ", "All"}], ",", " ", 
     RowBox[{"AxesLabel", " ", "->", " ", 
      RowBox[{"{", 
       RowBox[{"xl", ",", " ", "yl", ",", " ", "zl"}], "}"}]}], ",", 
     RowBox[{"SphericalRegion", "\[Rule]", "True"}]}], "}"}]}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.5779681027604*^9, 3.577968117939226*^9}, 
   3.578066372459113*^9}],

Cell[BoxData[
 Graphics3DBox[{SphereBox[{0, 0, 0}, 26220.7566254865], 
   {RGBColor[0., 0.7, 0.], 
    SphereBox[{26220.922084567108`, 0, 0}, 0.16545908060765693`]}, 
   {RGBColor[0.4, 1., 0.4], Opacity[0.5], 
    SphereBox[{58837.95263054443, -68594.08963595095, 0}, 75953.89838626009]}, 
   {RGBColor[0.4, 1., 0.4], Opacity[0.5], 
    SphereBox[{58837.95263054443, 34297.044817975475`, -59404.224174200404`}, 
     75953.89838626009]}, 
   {RGBColor[0.4, 1., 0.4], Opacity[0.5], 
    SphereBox[{58837.95263054443, 34297.044817975475`, 59404.224174200404`}, 
     75953.89838626009]}, 
   {RGBColor[0., 0.7, 0.], 
    SphereBox[{-26220.922084567108`, 0, 0}, 0.16545908060765693`]}, 
   {RGBColor[0.4, 1., 0.4], Opacity[0.5], 
    SphereBox[{-58837.95263054443, 68594.08963595095, 0}, 75953.89838626009]}, 
   {RGBColor[0.4, 1., 0.4], Opacity[0.5], 
    SphereBox[{-58837.95263054443, -34297.044817975475`, \
-59404.224174200404`}, 75953.89838626009]}, 
   {RGBColor[0.4, 1., 0.4], Opacity[0.5], 
    SphereBox[{-58837.95263054443, -34297.044817975475`, 59404.224174200404`},
      75953.89838626009]}},
  AspectRatio->Automatic,
  Axes->True,
  AxesLabel->{
    FormBox["xl", TraditionalForm], 
    FormBox["yl", TraditionalForm], 
    FormBox["zl", TraditionalForm]},
  PlotRange->All,
  SphericalRegion->True]], "Output",
 CellChangeTimes->{3.5787476144929895`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"plot3", "=", 
   RowBox[{
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"Graphics", "[", 
       RowBox[{"{", 
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"a", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "1"}], "]"}], "]"}], ",", 
            RowBox[{"a", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "2"}], "]"}], "]"}]}], "}"}], ",", 
          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"}]}], "}"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.5779612308334737`*^9, 3.5779612392574883`*^9}, 
   3.577967381025132*^9, 3.577967694382882*^9, 3.5780663764995203`*^9, 
   3.578283440435683*^9}],

Cell[BoxData[
 GraphicsBox[{CircleBox[{0, 0}, 26220.7566254865], 
   CircleBox[{26220.922084567108`, 0}, 0.16545908060765693`], 
   CircleBox[{58837.95263054443, -68594.08963595095}, 75953.89838626009], 
   CircleBox[{58837.95263054443, 34297.044817975475`}, 75953.89838626009], 
   CircleBox[{58837.95263054443, 34297.044817975475`}, 75953.89838626009], 
   CircleBox[{-26220.922084567108`, 0}, 0.16545908060765693`], 
   CircleBox[{-58837.95263054443, 68594.08963595095}, 75953.89838626009], 
   CircleBox[{-58837.95263054443, -34297.044817975475`}, 75953.89838626009], 
   CircleBox[{-58837.95263054443, -34297.044817975475`}, 75953.89838626009]},
  AspectRatio->Automatic,
  Axes->True,
  Frame->True,
  GridLines->Automatic,
  PlotRange->{{-4, 4}, {-3, 3}}]], "Output",
 CellChangeTimes->{3.5787476145085897`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"plot4", "=", 
   RowBox[{
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"Graphics", "[", 
       RowBox[{"{", 
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"a", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "2"}], "]"}], "]"}], ",", 
            RowBox[{"a", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "3"}], "]"}], "]"}]}], "}"}], ",", 
          RowBox[{"rr", "[", 
           RowBox[{"[", "i", "]"}], "]"}]}], "]"}], "}"}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "5"}], "}"}]}], "]"}], "     ", "/.", " ", 
    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",
 CellChangeTimes->{{3.5780221720630336`*^9, 3.5780221920154686`*^9}, {
   3.578022273104411*^9, 3.578022316909288*^9}, 3.5780223753469906`*^9, 
   3.578283440638483*^9}],

Cell[BoxData[
 GraphicsBox[{CircleBox[{0, 0}, 26220.7566254865], 
   CircleBox[{0, 0}, 0.16545908060765693`], 
   CircleBox[{-68594.08963595095, 0}, 75953.89838626009], 
   CircleBox[{34297.044817975475`, -59404.224174200404`}, 75953.89838626009], 
   CircleBox[{34297.044817975475`, 59404.224174200404`}, 75953.89838626009]},
  AspectRatio->Automatic,
  Axes->True,
  Frame->True,
  GridLines->Automatic,
  PlotRange->{{-4, 4}, {-3, 3}}]], "Output",
 CellChangeTimes->{3.5787476145085897`*^9}]
}, Open  ]]
},
WindowSize->{824, 817},
WindowMargins->{{501, Automatic}, {Automatic, 0}},
PrintingCopies->1,
PrintingPageRange->{32000, 32000},
PrintingOptions->{"Magnification"->1.,
"PaperOrientation"->"Portrait",
"PaperSize"->{595.2, 841.8}},
FrontEndVersion->"9.0 for Microsoft Windows (64-bit) (January 25, 2013)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[579, 22, 26352, 729, 1492, "Input"],
Cell[26934, 753, 444, 11, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[27415, 769, 1228, 29, 132, "Input"],
Cell[28646, 800, 86, 1, 31, "Output"],
Cell[28735, 803, 105, 2, 31, "Output"],
Cell[28843, 807, 86, 1, 31, "Output"],
Cell[28932, 810, 91, 1, 31, "Output"],
Cell[29026, 813, 88, 1, 31, "Output"],
Cell[29117, 816, 88, 1, 31, "Output"]
}, Open  ]],
Cell[29220, 820, 557, 13, 31, "Input"],
Cell[29780, 835, 5543, 151, 332, "Input"],
Cell[CellGroupData[{
Cell[35348, 990, 649, 15, 72, "Input"],
Cell[36000, 1007, 1368, 30, 365, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[37405, 1042, 1600, 43, 72, "Input"],
Cell[39008, 1087, 818, 15, 293, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[39863, 1107, 1592, 43, 72, "Input"],
Cell[41458, 1152, 494, 11, 325, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)