(* 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[     44062,       1277]
NotebookOptionsPosition[     42716,       1230]
NotebookOutlinePosition[     43216,       1250]
CellTagsIndexPosition[     43173,       1247]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell[BoxData[{
 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[{
   RowBox[{"b", "=", "a"}], ";", " ", 
   RowBox[{"c", "=", "a"}], ";"}], " "}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"g", "=", 
    RowBox[{"Table", "[", 
     RowBox[{"1.2", ",", 
      RowBox[{"{", "nc", "}"}]}], "]"}]}], ";"}], 
  " "}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"sig", "=", 
   RowBox[{"Table", "[", 
    RowBox[{"0.325", ",", 
     RowBox[{"{", "nc", "}"}]}], "]"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"x1", "=", "1.166"}], ";", " ", 
   RowBox[{"x2", "=", "0.1486"}], ";", " ", 
   RowBox[{"x3", "=", "0.96"}], ";", " ", 
   RowBox[{"x4", "=", "1.23"}], ";"}], "\n", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"sig1", "=", "0.31403448"}], ";", 
    RowBox[{"sig2", "=", "0.30821367"}], ";", 
    RowBox[{"sig3", "=", "0.30820426"}], ";", 
    RowBox[{"sig4", "=", "0.3"}], ";", 
    RowBox[{"sig5", "=", "0.3"}], ";", "\n", 
    RowBox[{"sig6", "=", "0.3"}], ";", 
    RowBox[{"sig7", "=", "0.3"}], ";"}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"For", "[", " ", 
     RowBox[{
      RowBox[{"i", "=", "1"}], ",", 
      RowBox[{"i", "<", "10"}], ",", 
      RowBox[{"i", "++"}], ",", " ", 
      RowBox[{
       RowBox[{"g", "[", 
        RowBox[{"[", "i", "]"}], "]"}], "=", "1.0"}]}], "]"}], ";"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"g", "[", 
     RowBox[{"[", "1", "]"}], "]"}], "=", "1.33"}], ";", 
   RowBox[{
    RowBox[{"g", "[", 
     RowBox[{"[", "2", "]"}], "]"}], "=", "0.8828"}], ";", 
   RowBox[{
    RowBox[{"g", "[", 
     RowBox[{"[", "6", "]"}], "]"}], "=", 
    RowBox[{"g", "[", 
     RowBox[{"[", "2", "]"}], "]"}]}], ";"}], "\n", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
     RowBox[{"sig", "[", 
      RowBox[{"[", "1", "]"}], "]"}], "=", "0.325"}], ";", 
    RowBox[{
     RowBox[{"sig", "[", 
      RowBox[{"[", "2", "]"}], "]"}], "=", "0.325"}], ";", 
    RowBox[{
     RowBox[{"sig", "[", 
      RowBox[{"[", "6", "]"}], "]"}], "=", 
     RowBox[{"sig", "[", 
      RowBox[{"[", "2", "]"}], "]"}]}], ";", 
    RowBox[{
     RowBox[{"sig", "[", 
      RowBox[{"[", "3", "]"}], "]"}], "=", "0.325"}], ";", "\n", 
    RowBox[{
     RowBox[{"sig", "[", 
      RowBox[{"[", "4", "]"}], "]"}], "=", 
     RowBox[{"sig", "[", 
      RowBox[{"[", "3", "]"}], "]"}]}], ";", 
    RowBox[{
     RowBox[{"sig", "[", 
      RowBox[{"[", "5", "]"}], "]"}], "=", 
     RowBox[{"sig", "[", 
      RowBox[{"[", "3", "]"}], "]"}]}], ";", 
    RowBox[{
     RowBox[{"sig", "[", 
      RowBox[{"[", "8", "]"}], "]"}], "=", 
     RowBox[{"sig", "[", 
      RowBox[{"[", "3", "]"}], "]"}]}], ";", 
    RowBox[{
     RowBox[{"sig", "[", 
      RowBox[{"[", "9", "]"}], "]"}], "=", 
     RowBox[{"sig", "[", 
      RowBox[{"[", "3", "]"}], "]"}]}], ";", 
    RowBox[{
     RowBox[{"sig", "[", 
      RowBox[{"[", "7", "]"}], "]"}], "=", 
     RowBox[{"sig", "[", 
      RowBox[{"[", "3", "]"}], "]"}]}], ";"}], "*)"}], 
  "\[IndentingNewLine]"}], "\n", 
 RowBox[{
  RowBox[{"ss", "=", 
   RowBox[{"Sin", "[", "x4", "]"}]}], ";", 
  RowBox[{"cs", "=", 
   RowBox[{"Cos", "[", "x4", "]"}]}], ";"}], "\n", 
 RowBox[{
  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[{
    RowBox[{"c", "[", 
     RowBox[{"[", 
      RowBox[{"2", ",", "1"}], "]"}], "]"}], "=", 
    RowBox[{"-", "2."}]}], ";", 
   RowBox[{
    RowBox[{"c", "[", 
     RowBox[{"[", 
      RowBox[{"2", ",", "2"}], "]"}], "]"}], "=", "z"}], ";", 
   RowBox[{
    RowBox[{"c", "[", 
     RowBox[{"[", 
      RowBox[{"2", ",", "3"}], "]"}], "]"}], "=", "x2"}], ";"}], 
  "\n"}], "\[IndentingNewLine]", 
 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[{
    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[{
    RowBox[{"c", "[", 
     RowBox[{"[", 
      RowBox[{"6", ",", "1"}], "]"}], "]"}], "=", 
    RowBox[{"-", "2."}]}], ";", 
   RowBox[{
    RowBox[{"c", "[", 
     RowBox[{"[", 
      RowBox[{"6", ",", "2"}], "]"}], "]"}], "=", "z"}], ";", 
   RowBox[{
    RowBox[{"c", "[", 
     RowBox[{"[", 
      RowBox[{"6", ",", "3"}], "]"}], "]"}], "=", "x2"}], ";"}], 
  "\[IndentingNewLine]"}], "\n", 
 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[{
    RowBox[{"c", "[", 
     RowBox[{"[", 
      RowBox[{"9", ",", "1"}], "]"}], "]"}], "=", 
    RowBox[{"-", "2."}]}], ";", 
   RowBox[{
    RowBox[{"c", "[", 
     RowBox[{"[", 
      RowBox[{"9", ",", "2"}], "]"}], "]"}], "=", "0"}], ";", 
   RowBox[{
    RowBox[{"c", "[", 
     RowBox[{"[", 
      RowBox[{"9", ",", "3"}], "]"}], "]"}], "=", "x3"}], ";"}], 
  "\[IndentingNewLine]", "\n", 
  RowBox[{"(*", 
   RowBox[{
   "coord", " ", "with", " ", "zero", " ", "value", " ", "not", " ", 
    "explicitely", " ", 
    RowBox[{"written", "!"}]}], "*)"}], "\n", 
  RowBox[{"(*", "nuclei", "*)"}]}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"b", "[", 
    RowBox[{"[", 
     RowBox[{"2", ",", "1"}], "]"}], "]"}], "=", 
   RowBox[{"x1", "+", "x2"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"b", "[", 
     RowBox[{"[", 
      RowBox[{"6", ",", "1"}], "]"}], "]"}], "=", 
    RowBox[{"-", 
     RowBox[{"b", "[", 
      RowBox[{"[", 
       RowBox[{"2", ",", "1"}], "]"}], "]"}]}]}], ";"}], "\n", 
  RowBox[{"(*", "clouds", "*)"}], "\[IndentingNewLine]"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"a", "[", 
     RowBox[{"[", 
      RowBox[{"2", ",", "1"}], "]"}], "]"}], "=", 
    RowBox[{"b", "[", 
     RowBox[{"[", 
      RowBox[{"2", ",", "1"}], "]"}], "]"}]}], ";"}], 
  "\n"}], "\[IndentingNewLine]", 
 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[{
    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"}], "]"}], "]"}]}]}], ";"}], 
  "\[IndentingNewLine]"}], "\n", 
 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[{
    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[{
   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[{
  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"}]}]}]}]}], 
        "]"}]}]}], "]"}]}], "]"}], "\[IndentingNewLine]", 
  "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
     RowBox[{"KJ", "/", "WI"}], " ", "I"}], "=", "J"}], "*)"}], "\n", 
  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"}], "]"}], "]"}]}]}]}]}], ";"}], 
   "*)"}]}], "\n", 
 RowBox[{
  RowBox[{"vne", "=", 
   RowBox[{
    RowBox[{"-", "2"}], "*", "3.0", "*", 
    RowBox[{"z", "/", "x2"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"vne", "\n", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
     RowBox[{"KJ", "/", "WI"}], " ", 
     RowBox[{"I", "<>", "J"}], " ", "inc"}], ";"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  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"}]}]}]}]}], 
         "]"}]}]}], "]"}]}], "]"}], ";"}], "\[IndentingNewLine]", 
  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"}]}]}]}]}], "]"}]}], "]"}], ";"}], 
   "*)"}]}], "\n", "vne", "\n", 
 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[{
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"t", "=", 
     RowBox[{"FindMinimum", "[", 
      RowBox[{"func", ",", 
       RowBox[{"{", 
        RowBox[{"x1", ",", "1.14", ",", "1.167"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"x2", ",", "0.14", ",", "0.156"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"x3", ",", "0.9", ",", "1.0"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"=", 
         RowBox[{",", "1.21", ",", "1.3"}]}], "}"}], ",", 
       RowBox[{"Method", "\[Rule]", "Automatic"}], ",", 
       RowBox[{"MaxIterations", "\[Rule]", "500"}]}], "]"}]}], ";"}], 
   "*)"}]}]}], "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}}],

Cell[BoxData[
 RowBox[{"-", "363.3916554508748`"}]], "Output",
 CellChangeTimes->{3.5780250940105658`*^9, 3.5780253898182855`*^9, 
  3.578025658341557*^9, 3.578025732098487*^9, 3.578025830097859*^9, 
  3.5780261506160216`*^9, 3.5780264335537186`*^9}],

Cell[BoxData[
 RowBox[{"-", "535.8887104005684`"}]], "Output",
 CellChangeTimes->{3.5780250940105658`*^9, 3.5780253898182855`*^9, 
  3.578025658341557*^9, 3.578025732098487*^9, 3.578025830097859*^9, 
  3.5780261506160216`*^9, 3.5780264335537186`*^9}],

Cell[BoxData[
 RowBox[{"-", "213.3416411181068`"}]], "Output",
 CellChangeTimes->{3.5780250940105658`*^9, 3.5780253898182855`*^9, 
  3.578025658341557*^9, 3.578025732098487*^9, 3.578025830097859*^9, 
  3.5780261506160216`*^9, 3.5780264335537186`*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"(*", 
  RowBox[{"MatrixForm", "[", "a", "]"}], "*)"}]], "Input",
 CellChangeTimes->{{3.5779595463893147`*^9, 3.5779595552345304`*^9}, {
   3.5779597259924307`*^9, 3.577959732669242*^9}, 3.5779607530670347`*^9, 
   3.577961309317211*^9, 3.577983470867136*^9}],

Cell[CellGroupData[{

Cell[BoxData[{"ekin", "\[IndentingNewLine]", "vne", "\[IndentingNewLine]", \
"vee", "\[IndentingNewLine]", "vnn", "\[IndentingNewLine]", 
 RowBox[{"vir", " ", "=", 
  RowBox[{
   RowBox[{"-", "epot"}], "/", "ekin"}]}], "\[IndentingNewLine]", 
 RowBox[{"2", "*", 
  RowBox[{"(", 
   RowBox[{"x1", "+", "x2"}], ")"}], "*", "0.529177"}]}], "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}}],

Cell[BoxData["199.68171989110044`"], "Output",
 CellChangeTimes->{3.578025094072966*^9, 3.5780253898650856`*^9, 
  3.578025658388357*^9, 3.5780257321296864`*^9, 3.578025830144659*^9, 
  3.578026150662822*^9, 3.5780264336005187`*^9}],

Cell[BoxData[
 RowBox[{"-", "535.8887104005684`"}]], "Output",
 CellChangeTimes->{3.578025094072966*^9, 3.5780253898650856`*^9, 
  3.578025658388357*^9, 3.5780257321296864`*^9, 3.578025830144659*^9, 
  3.578026150662822*^9, 3.5780264336005187`*^9}],

Cell[BoxData["92.05749909469296`"], "Output",
 CellChangeTimes->{3.578025094072966*^9, 3.5780253898650856`*^9, 
  3.578025658388357*^9, 3.5780257321296864`*^9, 3.578025830144659*^9, 
  3.578026150662822*^9, 3.5780264336005187`*^9}],

Cell[BoxData["30.80785029666819`"], "Output",
 CellChangeTimes->{3.578025094072966*^9, 3.5780253898650856`*^9, 
  3.578025658388357*^9, 3.5780257321296864`*^9, 3.578025830144659*^9, 
  3.578026150662822*^9, 3.5780264336005187`*^9}],

Cell[BoxData["2.068408471413688`"], "Output",
 CellChangeTimes->{3.578025094072966*^9, 3.5780253898650856`*^9, 
  3.578025658388357*^9, 3.5780257321296864`*^9, 3.578025830144659*^9, 
  3.578026150662822*^9, 3.578026433616119*^9}],

Cell[BoxData["1.3913121684`"], "Output",
 CellChangeTimes->{3.578025094072966*^9, 3.5780253898650856`*^9, 
  3.578025658388357*^9, 3.5780257321296864`*^9, 3.578025830144659*^9, 
  3.578026150662822*^9, 3.578026433616119*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"rr", "=", 
   RowBox[{"{", 
    RowBox[{
    "x1", ",", "x2", ",", "x3", ",", "x3", ",", "x3", ",", "x2", ",", "x3", 
     ",", "x3", ",", "x3"}], "}"}]}], "  ", 
  RowBox[{"(*", " ", 
   RowBox[{"/.", " ", 
    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}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"plot0", "=", 
   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[{"(*", " ", 
   RowBox[{"/.", " ", 
    RowBox[{"t", "[", 
     RowBox[{"[", "2", "]"}], "]"}]}], " ", "*)"}], " ", " ", ";"}], "\n", 
 RowBox[{
  RowBox[{"plot1", "=", 
   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[{"(*", " ", 
   RowBox[{"/.", " ", 
    RowBox[{"t", "[", 
     RowBox[{"[", "2", "]"}], "]"}]}], " ", "*)"}], " ", " ", ";"}], "\n", 
 RowBox[{
  RowBox[{"plot2", "=", 
   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[{"(*", " ", 
   RowBox[{"/.", " ", 
    RowBox[{"t", "[", 
     RowBox[{"[", "2", "]"}], "]"}]}], " ", "*)"}], "  ", 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"plot2a", "=", 
   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[{"(*", " ", 
   RowBox[{"/.", " ", 
    RowBox[{"t", "[", 
     RowBox[{"[", "2", "]"}], "]"}]}], " ", "*)"}], " ", " ", ";"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"plot2b", "=", 
    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[{"(*", " ", 
    RowBox[{"/.", " ", 
     RowBox[{"t", "[", 
      RowBox[{"[", "2", "]"}], "]"}]}], " ", "*)"}], "  ", ";"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]"}], "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}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  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"}]}], "}"}]}], "]"}], 
  "\n"}]], "Input",
 CellChangeTimes->{{3.5779681027604*^9, 3.577968117939226*^9}}],

Cell[BoxData[
 Graphics3DBox[{SphereBox[{0, 0, 0}, 1.166], 
   {RGBColor[0., 0.7, 0.], SphereBox[{1.3146, 0, 0}, 0.1486]}, 
   {RGBColor[0.4, 1., 0.4], Opacity[0.5], 
    SphereBox[{1.6851359442902236`, 1.0448430858214799`, 0}, 0.96]}, 
   {RGBColor[0.4, 1., 0.4], Opacity[0.5], 
    SphereBox[{1.6851359442902236`, -0.5224215429107399, 0.904860655289926}, 
     0.96]}, 
   {RGBColor[0.4, 1., 0.4], Opacity[0.5], 
    SphereBox[{1.6851359442902236`, -0.5224215429107399, -0.904860655289926}, 
     0.96]}, 
   {RGBColor[0., 0.7, 0.], SphereBox[{-1.3146, 0, 0}, 0.1486]}, 
   {RGBColor[0.4, 1., 0.4], Opacity[0.5], 
    SphereBox[{-1.6851359442902236`, -1.0448430858214799`, 0}, 0.96]}, 
   {RGBColor[0.4, 1., 0.4], Opacity[0.5], 
    SphereBox[{-1.6851359442902236`, 0.5224215429107399, 0.904860655289926}, 
     0.96]}, 
   {RGBColor[0.4, 1., 0.4], Opacity[0.5], 
    SphereBox[{-1.6851359442902236`, 0.5224215429107399, -0.904860655289926}, 
     0.96]}},
  AspectRatio->Automatic,
  Axes->True,
  AxesLabel->{
    FormBox["xl", TraditionalForm], 
    FormBox["yl", TraditionalForm], 
    FormBox["zl", TraditionalForm]},
  PlotRange->All,
  SphericalRegion->True]], "Output",
 CellChangeTimes->{3.578025094088566*^9, 3.5780253898962855`*^9, 
  3.578025658419557*^9, 3.578025732160887*^9, 3.5780258301758585`*^9, 
  3.5780261506940217`*^9, 3.5780264336317186`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"plot3", "=", 
    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[{"(*", " ", 
    RowBox[{"/.", " ", 
     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",
 CellChangeTimes->{{3.5779612308334737`*^9, 3.5779612392574883`*^9}, 
   3.577967381025132*^9, 3.577967694382882*^9}],

Cell[BoxData[
 GraphicsBox[{CircleBox[{0, 0}, 1.166], CircleBox[{1.3146, 0}, 0.1486], 
   CircleBox[{1.6851359442902236`, 1.0448430858214799`}, 0.96], 
   CircleBox[{1.6851359442902236`, -0.5224215429107399}, 0.96], 
   CircleBox[{1.6851359442902236`, -0.5224215429107399}, 0.96], 
   CircleBox[{-1.3146, 0}, 0.1486], 
   CircleBox[{-1.6851359442902236`, -1.0448430858214799`}, 0.96], 
   CircleBox[{-1.6851359442902236`, 0.5224215429107399}, 0.96], 
   CircleBox[{-1.6851359442902236`, 0.5224215429107399}, 0.96]},
  AspectRatio->Automatic,
  Axes->True,
  Frame->True,
  GridLines->Automatic,
  PlotRange->{{-4, 4}, {-3, 3}}]], "Output",
 CellChangeTimes->{3.578025094104166*^9, 3.5780253898962855`*^9, 
  3.5780256584351573`*^9, 3.578025732160887*^9, 3.5780258301758585`*^9, 
  3.578026150709622*^9, 3.578026433647319*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"plot4", "=", 
   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[{"(*", " ", 
   RowBox[{"/.", " ", 
    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}],

Cell[BoxData[
 GraphicsBox[{CircleBox[{0, 0}, 1.166], CircleBox[{0, 0}, 0.1486], 
   CircleBox[{1.0448430858214799`, 0}, 0.96], 
   CircleBox[{-0.5224215429107399, 0.904860655289926}, 0.96], 
   CircleBox[{-0.5224215429107399, -0.904860655289926}, 0.96]},
  AspectRatio->Automatic,
  Axes->True,
  Frame->True,
  GridLines->Automatic,
  PlotRange->{{-4, 4}, {-3, 3}}]], "Output",
 CellChangeTimes->{3.578025094119766*^9, 3.5780253899118853`*^9, 
  3.5780256584351573`*^9, 3.578025732176487*^9, 3.5780258301914587`*^9, 
  3.578026150725222*^9, 3.578026433647319*^9}]
}, Open  ]]
},
WindowSize->{703, 817},
WindowMargins->{{Automatic, 191}, {Automatic, 44}},
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, 26094, 776, 1652, "Input"],
Cell[26676, 800, 250, 4, 31, "Output"],
Cell[26929, 806, 250, 4, 31, "Output"],
Cell[27182, 812, 250, 4, 31, "Output"]
}, Open  ]],
Cell[27447, 819, 282, 5, 31, "Input"],
Cell[CellGroupData[{
Cell[27754, 828, 664, 12, 132, "Input"],
Cell[28421, 842, 232, 3, 31, "Output"],
Cell[28656, 847, 248, 4, 31, "Output"],
Cell[28907, 853, 231, 3, 31, "Output"],
Cell[29141, 858, 231, 3, 31, "Output"],
Cell[29375, 863, 229, 3, 31, "Output"],
Cell[29607, 868, 224, 3, 31, "Output"]
}, Open  ]],
Cell[29846, 874, 556, 13, 31, "Input"],
Cell[30405, 889, 5560, 157, 452, "Input"],
Cell[CellGroupData[{
Cell[35990, 1050, 653, 15, 92, "Input"],
Cell[36646, 1067, 1367, 30, 368, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[38050, 1102, 1640, 45, 112, "Input"],
Cell[39693, 1149, 825, 16, 293, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[40555, 1170, 1577, 43, 92, "Input"],
Cell[42135, 1215, 565, 12, 293, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)