(* 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[     52200,       1318]
NotebookOptionsPosition[     50381,       1257]
NotebookOutlinePosition[     50921,       1277]
CellTagsIndexPosition[     50878,       1274]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"(*", " ", 
    RowBox[{
     RowBox[{
     "C2H4", " ", "derived", " ", "from", " ", "B2H6", " ", "not", " ", 
      "TSM"}], ",", " ", 
     RowBox[{
     "overlaps", " ", "exchange", " ", "corrected", " ", "ES", " ", "08.06", 
      ".2013"}]}], " ", "*)"}], "Text"], "\n", 
  RowBox[{
   RowBox[{
    RowBox[{"Clear", "[", 
     RowBox[{
     "z", ",", "nc", ",", "R1", ",", "R2", ",", "R3", ",", "R4", ",", "wi", 
      ",", "p", ",", "vee", ",", "vne", ",", "vnn", ",", "\[Sigma]", ",", 
      "a1", ",", "a2", ",", "a3", ",", "b1", ",", "b2", ",", "b3", ",", "cs", 
      ",", "ss", ",", "d1", ",", "d2", ",", "d3", ",", "d4", ",", "d5", ",", 
      "pi", ",", "pb", ",", "i", ",", "j", ",", "d", ",", "t", ",", "q", ",", 
      "r", ",", "r2", ",", "ow", ",", "cn", ",", "c", ",", "w", ",", "n", ",",
       "dwr", ",", "dwnr", ",", "dnr", ",", "c2"}], "]"}], ";"}], 
   "\n"}]}]], "Input",
 CellChangeTimes->{{3.5793386253714504`*^9, 3.5793386509182944`*^9}, {
   3.579338894319522*^9, 3.5793389424768066`*^9}, {3.5793389747532635`*^9, 
   3.5793389925684943`*^9}, {3.579339023487749*^9, 3.5793390381361747`*^9}, {
   3.5793390697886305`*^9, 3.5793390710990324`*^9}, {3.579339124373126*^9, 
   3.5793392373953247`*^9}, {3.579339355503132*^9, 3.579339448806896*^9}, {
   3.5793399816258316`*^9, 3.579339981657032*^9}, {3.579340107627253*^9, 
   3.5793401136020637`*^9}, 3.5793409787327833`*^9, {3.57934431334344*^9, 
   3.5793443261510625`*^9}, {3.5793443573823175`*^9, 3.579344360112322*^9}, 
   3.5793443911407766`*^9, 3.5793444273796406`*^9, 3.5793444951773596`*^9, {
   3.5796661192565126`*^9, 3.579666120005314*^9}, {3.5796719837948136`*^9, 
   3.579672026024088*^9}, {3.5797526579487133`*^9, 3.5797526633151226`*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"z", "=", "6.0"}], ";", "  ", 
  RowBox[{"nc", "=", "8"}], ";", " ", 
  RowBox[{"pi", "=", 
   RowBox[{"-", "0.2"}]}], ";", " ", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"R3", "=", "0.6"}], ";"}], "*)"}], "\n", 
  RowBox[{"bohr", "=", "0.529177"}], ";", " ", 
  RowBox[{"rad", "=", "57.29578"}], ";"}], "\n", 
 RowBox[{
  RowBox[{"k", "=", 
   RowBox[{"{", 
    RowBox[{
    "1.0135", ",", "1.0135", ",", "1.3", ",", "1.3", ",", "1.3", ",", "1.3", 
     ",", "1.0", ",", "1.0"}], "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[Sigma]", "=", 
   RowBox[{"{", 
    RowBox[{
    "0.3", ",", "0.3", ",", "0.3", ",", "0.3", ",", "0.3", ",", "0.3", ",", 
     "0.3", ",", "0.3"}], "}"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.5793386253714504`*^9, 3.5793386509182944`*^9}, {
   3.579338894319522*^9, 3.5793389424768066`*^9}, {3.5793389747532635`*^9, 
   3.5793389925684943`*^9}, {3.579339023487749*^9, 3.5793390381361747`*^9}, {
   3.5793390697886305`*^9, 3.5793390710990324`*^9}, {3.579339124373126*^9, 
   3.5793392373953247`*^9}, {3.579339355503132*^9, 3.579339448806896*^9}, {
   3.5793399816258316`*^9, 3.579339981657032*^9}, {3.579340107627253*^9, 
   3.5793401136020637`*^9}, {3.5793402017266183`*^9, 3.579340301675994*^9}, {
   3.5793406867626705`*^9, 3.5793406954674854`*^9}, 3.579349891502223*^9, 
   3.579350146796671*^9, 3.579350269615687*^9, 3.579350853930313*^9, 
   3.5793508955043864`*^9, {3.57935093144685*^9, 3.579350935253256*^9}, {
   3.5793509795105343`*^9, 3.579351010772989*^9}, {3.579351043299046*^9, 
   3.5793510457950506`*^9}, {3.579364450837943*^9, 3.5793644593711576`*^9}, {
   3.5793651957548513`*^9, 3.57936521787569*^9}, 3.579365250464147*^9, 
   3.579365312505456*^9, {3.579365342831909*^9, 3.5793653574179354`*^9}, 
   3.5795764865231023`*^9, {3.5795768362601166`*^9, 3.579576863872165*^9}, {
   3.5795769265998755`*^9, 3.5795769273174767`*^9}, 3.5795769796867685`*^9, {
   3.5795773238545732`*^9, 3.579577325055775*^9}, 3.5795773626830416`*^9, 
   3.579577395583499*^9, 3.57957744712599*^9, {3.57957755000817*^9, 
   3.579577550725772*^9}, 3.579577624435901*^9, {3.5795779205712214`*^9, 
   3.579577944501663*^9}, {3.5795779836109324`*^9, 3.579577984032133*^9}, {
   3.579578023250602*^9, 3.579578023562602*^9}, {3.5796716344166*^9, 
   3.5796716368034043`*^9}, {3.5796716693138614`*^9, 
   3.5796716747894707`*^9}, {3.5796718550945873`*^9, 3.57967189630986*^9}, {
   3.579752617794243*^9, 3.5797526185118437`*^9}, {3.5797559420064864`*^9, 
   3.5797559428488874`*^9}, {3.579755975312545*^9, 3.5797559835649595`*^9}, 
   3.579756073935918*^9, 3.5798451728604136`*^9, {3.579845254401757*^9, 
   3.5798452577245626`*^9}, {3.5798452892834177`*^9, 3.57984536419475*^9}, 
   3.579845415019639*^9}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"dwr", "=", 
   RowBox[{"Table", "[", 
    RowBox[{"0", ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "nc"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"j", ",", "1", ",", "nc"}], "}"}]}], "]"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"dnr", "=", "dwr"}], ";", " ", 
  RowBox[{"dwnr", "=", "dwr"}], ";"}], "\n", 
 RowBox[{
  RowBox[{"c2", "=", 
   RowBox[{"Table", "[", 
    RowBox[{"0", ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "nc"}], "}"}]}], "]"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.5793386253714504`*^9, 3.5793386509182944`*^9}, {
   3.579338894319522*^9, 3.5793389424768066`*^9}, {3.5793389747532635`*^9, 
   3.5793389925684943`*^9}, {3.579339023487749*^9, 3.5793390381361747`*^9}, {
   3.5793390697886305`*^9, 3.5793390710990324`*^9}, {3.579339124373126*^9, 
   3.5793392373953247`*^9}, {3.579339355503132*^9, 3.579339448806896*^9}, {
   3.5793399816258316`*^9, 3.579339981657032*^9}, {3.579340107627253*^9, 
   3.579340120341275*^9}, 3.5793407052955027`*^9, 3.5793437072199755`*^9, {
   3.5793454193073826`*^9, 3.579345425344593*^9}, 3.5793503939479055`*^9, {
   3.5797544213782153`*^9, 3.5797544223766174`*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"ow", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "2"}], ",", 
     RowBox[{"-", "2"}], ",", 
     RowBox[{"-", "2"}], ",", 
     RowBox[{"-", "2"}], ",", 
     RowBox[{"-", "2"}], ",", 
     RowBox[{"-", "2"}], ",", 
     RowBox[{"-", "2"}], ",", 
     RowBox[{"-", "2"}]}], "}"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"cn", "=", 
   RowBox[{"{", 
    RowBox[{
    "z", ",", "z", ",", "1", ",", "1", ",", "1", ",", "1", ",", "0", ",", 
     "0"}], "}"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"c", "=", 
   RowBox[{"{", 
    RowBox[{
    "R1", ",", "R1", ",", "R2", ",", "R2", ",", "R2", ",", "R2", ",", "R4", 
     ",", "R4"}], "}"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"ec", "=", 
   RowBox[{"{", 
    RowBox[{
    "0", ",", "0", ",", "p", ",", "p", ",", "p", ",", "p", ",", "0", ",", 
     "0"}], "}"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"cs", "=", 
   RowBox[{"Cos", "[", "wi", "]"}]}], ";", "  ", 
  RowBox[{"ss", "=", 
   RowBox[{"Sin", "[", "wi", "]"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"d1", "=", 
    RowBox[{"R1", "+", "R3"}]}], ";", " ", 
   RowBox[{"d2", "=", 
    RowBox[{"d1", "+", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"R1", "+", "R2", "+", "p"}], ")"}], "*", "cs"}]}]}], ";", " ", 
   RowBox[{"d3", "=", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"R1", "+", "R2", "+", "p"}], ")"}], "*", "ss"}]}], ";"}], 
  " "}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"d4", "=", 
    RowBox[{"d1", "+", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"R1", "+", "R2"}], ")"}], "*", "cs"}]}]}], ";", "  ", 
   RowBox[{"d5", "=", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"R1", "+", "R2"}], ")"}], "*", "ss"}]}], ";"}], 
  " "}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"r", "=", 
    RowBox[{"pi", "+", "R4"}]}], ";", " ", 
   RowBox[{"r2", "=", 
    RowBox[{"2", "*", "r"}]}], ";"}], "\n", 
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{
     RowBox[{"Nuclear", " ", "positions"}], ";", " ", "1"}], ",", "2", ",", 
    RowBox[{"3", " ", "=", " ", "x"}], ",", "y", ",", "z"}], " ", 
   "*)"}]}], "\n", 
 RowBox[{
  RowBox[{"b1", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "d1"}], ",", "d1", ",", 
     RowBox[{"-", "d2"}], ",", 
     RowBox[{"-", "d2"}], ",", "d2", ",", "d2", ",", "0", ",", "0"}], "}"}]}],
   ";"}], "\n", 
 RowBox[{
  RowBox[{"b2", "=", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "d3", ",", 
     RowBox[{"-", "d3"}], ",", 
     RowBox[{"-", "d3"}], ",", "d3", ",", "0", ",", "0"}], "}"}]}], 
  ";"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"b3", "=", 
    RowBox[{"{", 
     RowBox[{
     "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", 
      "0"}], "}"}]}], ";"}], "        "}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"n", "=", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        RowBox[{"b1", "[", 
         RowBox[{"[", "i", "]"}], "]"}], ",", 
        RowBox[{"b2", "[", 
         RowBox[{"[", "i", "]"}], "]"}], ",", 
        RowBox[{"b3", "[", 
         RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "nc"}], "}"}]}], "]"}]}], ";"}], "\n", 
  RowBox[{"(*", " ", 
   RowBox[{"Cloud", " ", "positions"}], " ", "*)"}]}], "\n", 
 RowBox[{
  RowBox[{"a1", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "d1"}], ",", "d1", ",", 
     RowBox[{"-", "d4"}], ",", 
     RowBox[{"-", "d4"}], ",", "d4", ",", "d4", ",", "0", ",", "0"}], "}"}]}],
   ";"}], "\n", 
 RowBox[{
  RowBox[{"a2", "=", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "d5", ",", 
     RowBox[{"-", "d5"}], ",", 
     RowBox[{"-", "d5"}], ",", "d5", ",", "0", ",", "0"}], "}"}]}], 
  ";"}], "\n", 
 RowBox[{
  RowBox[{"a3", "=", 
   RowBox[{"{", 
    RowBox[{
    "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "r", ",", 
     RowBox[{"-", "r"}]}], "}"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"w", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"a1", "[", 
        RowBox[{"[", "i", "]"}], "]"}], ",", 
       RowBox[{"a2", "[", 
        RowBox[{"[", "i", "]"}], "]"}], ",", 
       RowBox[{"a3", "[", 
        RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "nc"}], "}"}]}], "]"}]}], 
  ";"}], "\n"}], "Input",
 CellChangeTimes->{{3.5793386253714504`*^9, 3.5793386509182944`*^9}, {
   3.579338894319522*^9, 3.5793389424768066`*^9}, {3.5793389747532635`*^9, 
   3.5793389925684943`*^9}, {3.579339023487749*^9, 3.5793390381361747`*^9}, {
   3.5793390697886305`*^9, 3.5793390710990324`*^9}, {3.579339124373126*^9, 
   3.5793392373953247`*^9}, {3.579339355503132*^9, 3.579339448806896*^9}, {
   3.5793399816258316`*^9, 3.579339981657032*^9}, {3.579340107627253*^9, 
   3.579340120341275*^9}, {3.5793404189258*^9, 3.57934045884627*^9}, {
   3.5793407200843287`*^9, 3.579340724171536*^9}, 3.5793435809845533`*^9, 
   3.579343698998761*^9, {3.5793440842570376`*^9, 3.5793441144118905`*^9}, {
   3.5793445016201706`*^9, 3.5793445051301765`*^9}, {3.579345216756627*^9, 
   3.5793452263194437`*^9}, 3.579365706499748*^9, {3.5795771074665933`*^9, 
   3.579577115750208*^9}, {3.5795792011306705`*^9, 3.5795792047654767`*^9}, {
   3.579579999789673*^9, 3.579580007496087*^9}, {3.5796657404098473`*^9, 
   3.579665809705169*^9}, {3.5796660974944744`*^9, 3.5796661080400934`*^9}, {
   3.579752631054266*^9, 3.579752634891873*^9}, {3.579753342415516*^9, 
   3.579753363085552*^9}, {3.5797537262697897`*^9, 3.579753730247797*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"For", "[", 
    RowBox[{
     RowBox[{"i", " ", "=", " ", "1"}], ",", " ", 
     RowBox[{"i", " ", "\[LessEqual]", " ", "nc"}], ",", " ", 
     RowBox[{"i", "++"}], ",", "\n", "     ", 
     RowBox[{
      RowBox[{
       RowBox[{"c2", "[", 
        RowBox[{"[", "i", "]"}], "]"}], " ", "=", " ", 
       RowBox[{
        RowBox[{"k", "[", 
         RowBox[{"[", "i", "]"}], "]"}], "/", 
        RowBox[{
         RowBox[{"c", "[", 
          RowBox[{"[", "i", "]"}], "]"}], "^", "2"}]}]}], ";", "\n", "     ", 
      RowBox[{"For", "[", 
       RowBox[{
        RowBox[{"j", " ", "=", " ", "1"}], ",", " ", 
        RowBox[{"j", " ", "\[LessEqual]", " ", "nc"}], ",", " ", 
        RowBox[{"j", "++"}], ",", "\n", "         ", 
        RowBox[{
         RowBox[{"d", " ", "=", " ", 
          RowBox[{"Norm", "[", 
           RowBox[{
            RowBox[{"w", "[", 
             RowBox[{"[", "i", "]"}], "]"}], " ", "-", " ", 
            RowBox[{"w", "[", 
             RowBox[{"[", "j", "]"}], "]"}]}], "]"}]}], ";", "\n", 
         "         ", 
         RowBox[{
          RowBox[{"dwr", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "j"}], "]"}], "]"}], " ", "=", " ", 
          RowBox[{"If", "[", 
           RowBox[{
            RowBox[{"i", " ", "\[NotEqual]", " ", "j"}], ",", " ", 
            RowBox[{"1", "/", "d"}], ",", 
            RowBox[{"1.5", "*", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"\[Sigma]", "[", 
                RowBox[{"[", "i", "]"}], "]"}], "/", 
               RowBox[{"c", "[", 
                RowBox[{"[", "i", "]"}], "]"}]}], ")"}]}]}], "]"}]}], ";", 
         "\n", "        ", 
         RowBox[{"d", "=", " ", 
          RowBox[{"Norm", "[", 
           RowBox[{
            RowBox[{"n", "[", 
             RowBox[{"[", "i", "]"}], "]"}], " ", "-", " ", 
            RowBox[{"n", "[", 
             RowBox[{"[", "j", "]"}], "]"}]}], "]"}]}], ";", "\n", "        ", 
         RowBox[{
          RowBox[{"dnr", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "j"}], "]"}], "]"}], " ", "=", " ", 
          RowBox[{"If", "[", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{
              RowBox[{"i", " ", "\[NotEqual]", " ", "j"}], " ", "&&", " ", 
              RowBox[{"d", "\[NotEqual]", "0"}]}], ")"}], ",", " ", 
            RowBox[{"1", "/", "d"}], ",", " ", "0"}], "]"}]}], ";", "\n", 
         "        ", 
         RowBox[{"d", "=", " ", 
          RowBox[{"Norm", "[", 
           RowBox[{
            RowBox[{"w", "[", 
             RowBox[{"[", "i", "]"}], "]"}], " ", "-", " ", 
            RowBox[{"n", "[", 
             RowBox[{"[", "j", "]"}], "]"}]}], "]"}]}], ";", "\n", "        ", 
         RowBox[{
          RowBox[{"dwnr", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "j"}], "]"}], "]"}], " ", "=", " ", 
          RowBox[{"If", "[", 
           RowBox[{
            RowBox[{"i", " ", "\[NotEqual]", " ", "j"}], ",", " ", 
            RowBox[{"1", "/", "d"}], ",", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{"1.5", " ", "-", " ", 
               RowBox[{"0.5", "*", 
                RowBox[{
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{"ec", "[", 
                    RowBox[{"[", "i", "]"}], "]"}], "/", 
                   RowBox[{"c", "[", 
                    RowBox[{"[", "i", "]"}], "]"}]}], ")"}], "^", "2"}]}]}], 
              ")"}], "/", 
             RowBox[{"c", "[", 
              RowBox[{"[", "i", "]"}], "]"}]}]}], "]"}]}]}]}], "]"}]}]}], 
    "]"}], ";"}], " "}]], "Input",
 CellChangeTimes->{{3.5793386253714504`*^9, 3.5793386509182944`*^9}, {
  3.579338894319522*^9, 3.5793389424768066`*^9}, {3.5793389747532635`*^9, 
  3.5793389925684943`*^9}, {3.579339023487749*^9, 3.5793390381361747`*^9}, {
  3.5793390697886305`*^9, 3.5793390710990324`*^9}, {3.579339124373126*^9, 
  3.5793392373953247`*^9}, {3.579339355503132*^9, 3.579339448806896*^9}, {
  3.579339981657032*^9, 3.5793399816726317`*^9}, {3.579341222280411*^9, 
  3.5793412545412674`*^9}, {3.5797541569577513`*^9, 3.5797541794997907`*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"Ekin", " ", "=", " ", 
    RowBox[{
     RowBox[{"-", 
      RowBox[{"(", 
       RowBox[{"9", "/", "8"}], ")"}]}], "*", 
     RowBox[{"ow", " ", ".", " ", "c2"}]}]}], ";"}], " "}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"vee", " ", "=", " ", 
    RowBox[{"0.5", "*", 
     RowBox[{"ow", " ", ".", " ", 
      RowBox[{"(", 
       RowBox[{"ow", " ", ".", " ", "dwr"}], ")"}]}]}]}], ";"}], " "}], "\n", 
 RowBox[{
  RowBox[{"delvee", "=", 
   RowBox[{"If", "[", 
    RowBox[{
     RowBox[{"pi", "<", "0"}], ",", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"4", "/", "R4"}], ")"}], "*", 
      RowBox[{"(", 
       RowBox[{"1.2", "-", 
        RowBox[{"0.5", "*", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"r2", "/", "R4"}], ")"}], "^", "2"}]}], "+", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"3", "/", "16"}], ")"}], "*", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"r2", "/", "R4"}], ")"}], "^", "3"}]}], "-", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"1", "/", "160"}], ")"}], "*", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"r2", "/", "R4"}], ")"}], "^", "5"}]}]}], ")"}]}], ",", 
     "0"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"vee", "=", 
   RowBox[{"vee", "+", 
    RowBox[{"If", "[", 
     RowBox[{
      RowBox[{"pi", "<", "0"}], ",", 
      RowBox[{
       RowBox[{
        RowBox[{"-", "4"}], "/", "r2"}], "+", "delvee", "-", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"2", "/", "R4"}], ")"}], "*", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{
           RowBox[{"-", "0.5"}], "*", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{"r2", "/", "R4"}], ")"}], "^", "2"}]}], "+", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{"3", "/", "16"}], ")"}], "*", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{"r2", "/", "R4"}], ")"}], "^", "3"}]}], "-", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{"1", "/", "160"}], ")"}], "*", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{"r2", "/", "R4"}], ")"}], "^", "5"}]}]}], ")"}]}]}], ",",
       "0"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"vnn", " ", "=", " ", 
    RowBox[{"0.5", "*", 
     RowBox[{"cn", " ", ".", " ", 
      RowBox[{"(", 
       RowBox[{"cn", " ", ".", " ", "dnr"}], ")"}]}]}]}], ";"}], " "}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"vne", " ", "=", " ", 
    RowBox[{"cn", " ", ".", " ", 
     RowBox[{"(", 
      RowBox[{"ow", " ", ".", " ", "dwnr"}], ")"}]}]}], ";"}], " "}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"Epot", " ", "=", " ", 
    RowBox[{"vne", " ", "+", " ", "vee", " ", "+", " ", "vnn"}]}], ";"}], 
  " "}], "\n", 
 RowBox[{
  RowBox[{"vir", "=", 
   RowBox[{
    RowBox[{"-", "Epot"}], "/", "Ekin"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"func", " ", "=", " ", 
    RowBox[{"Ekin", " ", "+", " ", "Epot"}]}], ";"}], " "}]}], "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5793386253714504`*^9, 3.5793386509182944`*^9}, {
   3.579338894319522*^9, 3.5793389424768066`*^9}, {3.5793389747532635`*^9, 
   3.5793389925684943`*^9}, {3.579339023487749*^9, 3.5793390381361747`*^9}, {
   3.5793390697886305`*^9, 3.5793390710990324`*^9}, {3.579339124373126*^9, 
   3.5793392373953247`*^9}, {3.579339355503132*^9, 3.579339506199397*^9}, {
   3.579339540051456*^9, 3.5793395413774586`*^9}, {3.5793395728427134`*^9, 
   3.579339683743308*^9}, {3.5793399816726317`*^9, 3.5793400106886826`*^9}, 
   3.579340927377493*^9, 3.5793453744885035`*^9, {3.579666680935499*^9, 
   3.5796668175293393`*^9}, {3.5796668656242237`*^9, 
   3.5796669598951893`*^9}, {3.5796671065042467`*^9, 
   3.5796671286562853`*^9}, {3.5798427365417347`*^9, 3.579842744840949*^9}, {
   3.5798427861186213`*^9, 3.5798429198576565`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"t", " ", "=", " ", 
  RowBox[{"NMinimize", "[", 
   RowBox[{"func", ",", 
    RowBox[{"(*", 
     RowBox[{
      RowBox[{
       RowBox[{
        RowBox[{
         RowBox[{
          RowBox[{"2", "*", "d1"}], "-", "2.53"}], " ", "\[Equal]", "0"}], 
        " ", "&&", " ", 
        RowBox[{
         RowBox[{"R1", "+", "R2", "+", "p"}], " ", "\[Equal]", " ", "2.054"}],
         " ", "&&", " ", 
        RowBox[{"wi", "\[Equal]", "1.0245"}]}], "}"}], ","}], "*)"}], 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{"R1", ",", "0.25", ",", "0.28"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"R2", ",", "1.15", ",", "1.25"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"R3", ",", "0.5", ",", "1.2"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"R4", ",", "1.2", ",", "1.90"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"p", ",", "0.49", ",", "0.61"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"wi", ",", "1.00", ",", "1.1"}], "}"}]}], "}"}]}], 
   "]"}]}]], "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5793386253714504`*^9, 3.5793386509182944`*^9}, {
   3.579338894319522*^9, 3.5793389424768066`*^9}, {3.5793389747532635`*^9, 
   3.5793389925684943`*^9}, {3.579339023487749*^9, 3.5793390381361747`*^9}, {
   3.5793390697886305`*^9, 3.5793390710990324`*^9}, {3.579339124373126*^9, 
   3.5793392373953247`*^9}, {3.579339355503132*^9, 3.579339506199397*^9}, {
   3.579339540051456*^9, 3.5793395413774586`*^9}, {3.5793395728427134`*^9, 
   3.579339683743308*^9}, {3.5793399816726317`*^9, 3.5793400106886826`*^9}, 
   3.579340927377493*^9, {3.5793433833366065`*^9, 3.5793433837422075`*^9}, {
   3.579344512243789*^9, 3.57934455821707*^9}, 3.579344675232876*^9, {
   3.579344794885086*^9, 3.5793447997678947`*^9}, {3.5793452439006743`*^9, 
   3.579345253806692*^9}, {3.5793453515408635`*^9, 3.579345356361272*^9}, 
   3.579345443300225*^9, {3.579350403417122*^9, 3.579350409579133*^9}, {
   3.579365415216037*^9, 3.5793654535453043`*^9}, {3.5794458283765016`*^9, 
   3.579445830700906*^9}, {3.5794458648805656`*^9, 3.5794458940682173`*^9}, {
   3.579446370758254*^9, 3.5794464147659316`*^9}, {3.5794466902156153`*^9, 
   3.579446734426093*^9}, {3.5794468710979333`*^9, 3.579446890207967*^9}, 
   3.579576260743906*^9, {3.5795763810981174`*^9, 3.5795764136085744`*^9}, 
   3.579576825995299*^9, {3.579580198721223*^9, 3.5795802058504353`*^9}, {
   3.579662269871752*^9, 3.5796623591663084`*^9}, {3.5796659620706367`*^9, 
   3.579665967827047*^9}, 3.5797533970156116`*^9, {3.5797534604765224`*^9, 
   3.5797534623953266`*^9}, {3.579753493345781*^9, 3.579753526199438*^9}, {
   3.5797540484595604`*^9, 3.579754087022828*^9}, 3.579754226268673*^9, {
   3.5797543136444263`*^9, 3.5797543547660985`*^9}, 3.5797546678742485`*^9, 
   3.5797547213511424`*^9, {3.5797554161295624`*^9, 3.579755419233968*^9}, {
   3.579755449997222*^9, 3.579755464645648*^9}, 3.579755999414587*^9, {
   3.579845424676056*^9, 3.5798454285292625`*^9}}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"NMinimize", "::", "cvdiv"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Failed to converge to a solution. The function may be \
unbounded. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", ButtonStyle->\\\"Link\
\\\", ButtonFrame->None, ButtonData:>\\\"paclet:ref/message/NMinimize/cvdiv\\\
\", ButtonNote -> \\\"NMinimize::cvdiv\\\"]\\)\"\>"}]], "Message", "MSG",
 CellChangeTimes->{
  3.579843018340629*^9, 3.579845178445223*^9, 3.5798452682545815`*^9, {
   3.5798452984874344`*^9, 3.5798453388603053`*^9}, 3.5798453737575665`*^9, 
   3.5798454359860754`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", "6.456618614949853`*^10"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"R1", "\[Rule]", "0.29054649214032524`"}], ",", 
     RowBox[{"R2", "\[Rule]", "1.195815860774469`"}], ",", 
     RowBox[{"R3", "\[Rule]", 
      RowBox[{"-", "0.29054649229520507`"}]}], ",", 
     RowBox[{"R4", "\[Rule]", "2.1856604075996877`"}], ",", 
     RowBox[{"p", "\[Rule]", "0.6017027314352494`"}], ",", 
     RowBox[{"wi", "\[Rule]", "0.9836430254160703`"}]}], "}"}]}], 
  "}"}]], "Output",
 CellChangeTimes->{
  3.5798430183562293`*^9, 3.579845178445223*^9, 3.5798452682545815`*^9, {
   3.5798452984874344`*^9, 3.5798453388603053`*^9}, 3.5798453737575665`*^9, 
   3.5798454359860754`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"vee", " ", "/.", " ", 
  RowBox[{"t", "[", 
   RowBox[{"[", "2", "]"}], "]"}]}], "\n", 
 RowBox[{"vne", " ", "/.", " ", 
  RowBox[{"t", "[", 
   RowBox[{"[", "2", "]"}], "]"}]}], "\n", 
 RowBox[{"vnn", " ", "/.", " ", 
  RowBox[{"t", "[", 
   RowBox[{"[", "2", "]"}], "]"}]}], "\n", 
 RowBox[{"vir", " ", "/.", " ", 
  RowBox[{"t", "[", 
   RowBox[{"[", "2", "]"}], "]"}]}]}], "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5793386253714504`*^9, 3.5793386509182944`*^9}, {
   3.579338894319522*^9, 3.5793389424768066`*^9}, {3.5793389747532635`*^9, 
   3.5793389925684943`*^9}, {3.579339023487749*^9, 3.5793390381361747`*^9}, {
   3.5793390697886305`*^9, 3.5793390710990324`*^9}, {3.579339124373126*^9, 
   3.5793392373953247`*^9}, {3.579339355503132*^9, 3.579339506199397*^9}, {
   3.579339540051456*^9, 3.5793395413774586`*^9}, {3.5793395728427134`*^9, 
   3.579339683743308*^9}, {3.5793399816726317`*^9, 3.5793400106886826`*^9}, 
   3.579340927377493*^9, {3.5793433833366065`*^9, 3.5793433837422075`*^9}, {
   3.579344512243789*^9, 3.57934455821707*^9}, 3.579344675232876*^9, {
   3.579344794885086*^9, 3.5793447997678947`*^9}, {3.5793452439006743`*^9, 
   3.579345253806692*^9}, 3.5793453515408635`*^9, {3.5793504919472775`*^9, 
   3.5793505047393*^9}, {3.5793505769830265`*^9, 3.5793506236271086`*^9}, {
   3.579350654125162*^9, 3.579350657697569*^9}, {3.579350819657053*^9, 
   3.5793508278626676`*^9}, {3.5793645321452856`*^9, 3.579364610285823*^9}, {
   3.579364652530697*^9, 3.5793647063663917`*^9}}],

Cell[BoxData["1.2913237280696339`*^10"], "Output",
 CellChangeTimes->{
  3.5798430183718295`*^9, 3.579845178460823*^9, 3.579845268270181*^9, {
   3.5798452984874344`*^9, 3.579845338875905*^9}, 3.5798453737731667`*^9, 
   3.579845436017276*^9}],

Cell[BoxData[
 RowBox[{"-", "7.747942353867004`*^10"}]], "Output",
 CellChangeTimes->{
  3.5798430183718295`*^9, 3.579845178460823*^9, 3.579845268270181*^9, {
   3.5798452984874344`*^9, 3.579845338875905*^9}, 3.5798453737731667`*^9, 
   3.579845436017276*^9}],

Cell[BoxData[
 RowBox[{"0.5`", " ", 
  RowBox[{"(", 
   RowBox[{
    RowBox[{"12.`", " ", 
     RowBox[{"(", 
      RowBox[{"0.`", "\[VeryThinSpace]", "+", 
       FractionBox["10.`", 
        RowBox[{"Abs", "[", 
         RowBox[{"0.2`", "\[VeryThinSpace]", "-", "R4"}], "]"}]]}], ")"}]}], 
    "+", 
    FractionBox["60.`", 
     RowBox[{"Abs", "[", 
      RowBox[{"0.2`", "\[VeryThinSpace]", "-", "R4"}], "]"}]]}], 
   ")"}]}]], "Output",
 CellChangeTimes->{
  3.5798430183718295`*^9, 3.579845178460823*^9, 3.579845268270181*^9, {
   3.5798452984874344`*^9, 3.579845338875905*^9}, 3.5798453737731667`*^9, 
   3.579845436017276*^9}],

Cell[BoxData[
 RowBox[{"0.01583525734439699`", " ", 
  RowBox[{"(", 
   RowBox[{"6.4566186257973694`*^10", "-", 
    RowBox[{"0.5`", " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"12.`", " ", 
        RowBox[{"(", 
         RowBox[{"0.`", "\[VeryThinSpace]", "+", 
          FractionBox["10.`", 
           RowBox[{"Abs", "[", 
            RowBox[{"0.2`", "\[VeryThinSpace]", "-", "R4"}], "]"}]]}], 
         ")"}]}], "+", 
       FractionBox["60.`", 
        RowBox[{"Abs", "[", 
         RowBox[{"0.2`", "\[VeryThinSpace]", "-", "R4"}], "]"}]]}], ")"}]}]}],
    ")"}]}]], "Output",
 CellChangeTimes->{
  3.5798430183718295`*^9, 3.579845178460823*^9, 3.579845268270181*^9, {
   3.5798452984874344`*^9, 3.579845338875905*^9}, 3.5798453737731667`*^9, 
   3.579845436017276*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"2", "*", "d1", "*", "bohr"}], " ", "/.", " ", 
  RowBox[{"t", "[", 
   RowBox[{"[", "2", "]"}], "]"}], 
  "                                                     ", 
  RowBox[{"(*", " ", 
   RowBox[{"C", "-", 
    RowBox[{"C", " ", "1.33"}]}], " ", "*)"}]}]], "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5793386253714504`*^9, 3.5793386509182944`*^9}, {
   3.579338894319522*^9, 3.5793389424768066`*^9}, {3.5793389747532635`*^9, 
   3.5793389925684943`*^9}, {3.579339023487749*^9, 3.5793390381361747`*^9}, {
   3.5793390697886305`*^9, 3.5793390710990324`*^9}, {3.579339124373126*^9, 
   3.5793392373953247`*^9}, {3.579339355503132*^9, 3.579339506199397*^9}, {
   3.579339540051456*^9, 3.5793395413774586`*^9}, {3.5793395728427134`*^9, 
   3.579339683743308*^9}, {3.5793399816726317`*^9, 3.5793400106886826`*^9}, 
   3.579340927377493*^9, {3.5793433833366065`*^9, 3.5793433837422075`*^9}, {
   3.579344512243789*^9, 3.57934455821707*^9}, 3.579344675232876*^9, {
   3.579344794885086*^9, 3.5793447997678947`*^9}, {3.5793452439006743`*^9, 
   3.579345253806692*^9}, 3.5793453515408635`*^9, {3.5793504919472775`*^9, 
   3.5793505047393*^9}, {3.5793505769830265`*^9, 3.5793506236271086`*^9}, {
   3.579350654125162*^9, 3.579350657697569*^9}, {3.579350819657053*^9, 
   3.5793508278626676`*^9}, {3.5793645321452856`*^9, 3.579364610285823*^9}, {
   3.579364652530697*^9, 3.5793647063663917`*^9}, 3.579364778485318*^9, {
   3.579364848888242*^9, 3.579364851321846*^9}, {3.5797535852299423`*^9, 
   3.579753589348349*^9}, {3.579753639970438*^9, 3.5797536430280437`*^9}}],

Cell[BoxData[
 RowBox[{"-", "1.639176888511562`*^-10"}]], "Output",
 CellChangeTimes->{
  3.5798430184030294`*^9, 3.5798451784764233`*^9, 3.579845268285781*^9, {
   3.579845298503034*^9, 3.5798453388915052`*^9}, 3.5798453737887664`*^9, 
   3.5798454360484757`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"(", 
    RowBox[{"R1", "+", "R2", "+", "p"}], ")"}], "*", "bohr"}], " ", "/.", " ", 
  RowBox[{"t", "[", 
   RowBox[{"[", "2", "]"}], "]"}], "                                     ", 
  RowBox[{"(*", " ", 
   RowBox[{"C", "-", 
    RowBox[{"H", " ", "1.1"}]}], " ", "*)"}]}]], "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5793386253714504`*^9, 3.5793386509182944`*^9}, {
   3.579338894319522*^9, 3.5793389424768066`*^9}, {3.5793389747532635`*^9, 
   3.5793389925684943`*^9}, {3.579339023487749*^9, 3.5793390381361747`*^9}, {
   3.5793390697886305`*^9, 3.5793390710990324`*^9}, {3.579339124373126*^9, 
   3.5793392373953247`*^9}, {3.579339355503132*^9, 3.579339506199397*^9}, {
   3.579339540051456*^9, 3.5793395413774586`*^9}, {3.5793395728427134`*^9, 
   3.579339683743308*^9}, {3.5793399816726317`*^9, 3.5793400106886826`*^9}, 
   3.579340927377493*^9, {3.5793433833366065`*^9, 3.5793433837422075`*^9}, {
   3.579344512243789*^9, 3.57934455821707*^9}, 3.579344675232876*^9, {
   3.579344794885086*^9, 3.5793447997678947`*^9}, {3.5793452439006743`*^9, 
   3.579345253806692*^9}, 3.5793453515408635`*^9, {3.5793504919472775`*^9, 
   3.5793505047393*^9}, {3.5793505769830265`*^9, 3.5793506236271086`*^9}, {
   3.579350654125162*^9, 3.579350657697569*^9}, {3.579350819657053*^9, 
   3.5793508278626676`*^9}, {3.5793645321452856`*^9, 3.579364610285823*^9}, {
   3.579364652530697*^9, 3.579364711623601*^9}, 3.5793647825569253`*^9, {
   3.579364856142255*^9, 3.5793648764846907`*^9}, {3.579365039161776*^9, 
   3.57936504123658*^9}, {3.579753593700757*^9, 3.579753597694364*^9}, 
   3.579753647942052*^9}],

Cell[BoxData["1.104956017141103`"], "Output",
 CellChangeTimes->{
  3.5798430184186296`*^9, 3.5798451784920235`*^9, 3.579845268301381*^9, {
   3.579845298518634*^9, 3.5798453389071054`*^9}, 3.579845373804366*^9, 
   3.5798454360640755`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"2", "*", "wi", "*", "rad"}], " ", "/.", " ", 
  RowBox[{"t", "[", 
   RowBox[{"[", "2", "]"}], "]"}], 
  "                                                         ", 
  RowBox[{"(*", " ", 
   RowBox[{"HCH", " ", "121.7"}], " ", "*)"}]}]], "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5793386253714504`*^9, 3.5793386509182944`*^9}, {
   3.579338894319522*^9, 3.5793389424768066`*^9}, {3.5793389747532635`*^9, 
   3.5793389925684943`*^9}, {3.579339023487749*^9, 3.5793390381361747`*^9}, {
   3.5793390697886305`*^9, 3.5793390710990324`*^9}, {3.579339124373126*^9, 
   3.5793392373953247`*^9}, {3.579339355503132*^9, 3.579339506199397*^9}, {
   3.579339540051456*^9, 3.5793395413774586`*^9}, {3.5793395728427134`*^9, 
   3.579339683743308*^9}, {3.5793399816726317`*^9, 3.5793400106886826`*^9}, 
   3.579340927377493*^9, {3.5793433833366065`*^9, 3.5793433837422075`*^9}, {
   3.579344512243789*^9, 3.57934455821707*^9}, 3.579344675232876*^9, {
   3.579344794885086*^9, 3.5793447997678947`*^9}, {3.5793452439006743`*^9, 
   3.579345253806692*^9}, 3.5793453515408635`*^9, {3.5793504919472775`*^9, 
   3.5793505047393*^9}, {3.5793505769830265`*^9, 3.5793506236271086`*^9}, {
   3.579350654125162*^9, 3.579350657697569*^9}, {3.579350819657053*^9, 
   3.5793508278626676`*^9}, {3.5793645321452856`*^9, 3.579364610285823*^9}, {
   3.579364652530697*^9, 3.579364717567211*^9}, {3.579364755553278*^9, 
   3.57936478491253*^9}, {3.579364865486671*^9, 3.5793648811646986`*^9}, {
   3.579365054418603*^9, 3.5793650555730047`*^9}, {3.579753603076373*^9, 
   3.5797536115627885`*^9}}],

Cell[BoxData["112.71718876554715`"], "Output",
 CellChangeTimes->{
  3.5798430184498296`*^9, 3.5798451784920235`*^9, 3.5798452683169813`*^9, {
   3.579845298534234*^9, 3.579845338922705*^9}, 3.5798453738199663`*^9, 
   3.579845436095276*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"plot1", "=", 
   RowBox[{
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"Graphics", "[", 
       RowBox[{"{", 
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"a1", "[", 
             RowBox[{"[", "i", "]"}], "]"}], ",", 
            RowBox[{"a2", "[", 
             RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", 
          RowBox[{"c", "[", 
           RowBox[{"[", "i", "]"}], "]"}]}], "]"}], "}"}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "nc"}], "}"}]}], "]"}], " ", "/.", 
    RowBox[{"t", "[", 
     RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"plot1a", "=", 
   RowBox[{
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"Graphics", "[", 
       RowBox[{"{", 
        RowBox[{"Disk", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"b1", "[", 
             RowBox[{"[", "i", "]"}], "]"}], ",", 
            RowBox[{"b2", "[", 
             RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], 
        "}"}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "6"}], "}"}]}], "]"}], " ", "/.", 
    RowBox[{"t", "[", 
     RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"plot1b", "=", 
   RowBox[{
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"Graphics", "[", 
       RowBox[{"{", 
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"a1", "[", 
             RowBox[{"[", "i", "]"}], "]"}], ",", 
            RowBox[{"a2", "[", 
             RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], 
        "}"}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "nc"}], "}"}]}], "]"}], " ", "/.", 
    RowBox[{"t", "[", 
     RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"Show", "[", 
  RowBox[{"plot1", ",", "plot1a", ",", "plot1b", ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"AspectRatio", " ", "\[Rule]", " ", "Automatic"}], ",", 
     RowBox[{"Axes", " ", "->", " ", "True"}], ",", 
     RowBox[{"GridLines", " ", "->", " ", "Automatic"}], ",", " ", "\n", 
     RowBox[{"PlotRange", " ", "\[Rule]", " ", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{
          RowBox[{"-", "4"}], ",", "4"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"-", "4"}], ",", "4"}], "}"}]}], "}"}]}], ",", " ", 
     RowBox[{"Frame", " ", "->", " ", "True"}]}], "}"}]}], 
  "]"}], "\[IndentingNewLine]", 
 StyleBox[
  RowBox[{"xy", " ", "projection"}], "Text"]}], "Input",
 CellChangeTimes->{{3.579339709701754*^9, 3.5793398596336174`*^9}, {
  3.5793399817038317`*^9, 3.579339981735032*^9}, {3.579579406271031*^9, 
  3.5795795916617565`*^9}, {3.5795796337662306`*^9, 3.579579644499049*^9}, {
  3.5795797243243895`*^9, 3.5795797255879917`*^9}, {3.5795798997466974`*^9, 
  3.5795799348935595`*^9}, {3.5796626781556687`*^9, 3.579662707670921*^9}, {
  3.5796717667828326`*^9, 3.5796717759712486`*^9}}],

Cell[BoxData[
 GraphicsBox[{CircleBox[{1.5487983118234183`*^-10, 0}, 0.29054649214032524`], 
   CircleBox[{-1.5487983118234183`*^-10, 0}, 0.29054649214032524`], 
   CircleBox[{-0.8234348354678755, 1.2374280244302411`}, 1.195815860774469], 
   CircleBox[{-0.8234348354678755, -1.2374280244302411`}, 1.195815860774469], 
   CircleBox[{0.8234348354678755, -1.2374280244302411`}, 1.195815860774469], 
   CircleBox[{0.8234348354678755, 1.2374280244302411`}, 1.195815860774469], 
   CircleBox[{0, 0}, 2.1856604075996877`], 
   CircleBox[{0, 0}, 2.1856604075996877`], 
   DiskBox[{1.5487983118234183`*^-10, 0}, 0.08], 
   DiskBox[{-1.5487983118234183`*^-10, 0}, 0.08], 
   DiskBox[{-1.1567741378133178`, 1.7383582456472222`}, 0.08], 
   DiskBox[{-1.1567741378133178`, -1.7383582456472222`}, 0.08], 
   DiskBox[{1.1567741378133178`, -1.7383582456472222`}, 0.08], 
   DiskBox[{1.1567741378133178`, 1.7383582456472222`}, 0.08], 
   CircleBox[{1.5487983118234183`*^-10, 0}, 0.08], 
   CircleBox[{-1.5487983118234183`*^-10, 0}, 0.08], 
   CircleBox[{-0.8234348354678755, 1.2374280244302411`}, 0.08], 
   CircleBox[{-0.8234348354678755, -1.2374280244302411`}, 0.08], 
   CircleBox[{0.8234348354678755, -1.2374280244302411`}, 0.08], 
   CircleBox[{0.8234348354678755, 1.2374280244302411`}, 0.08], 
   CircleBox[{0, 0}, 0.08], CircleBox[{0, 0}, 0.08]},
  AspectRatio->Automatic,
  Axes->True,
  Frame->True,
  GridLines->Automatic,
  PlotRange->{{-4, 4}, {-4, 4}}]], "Output",
 CellChangeTimes->{
  3.5798430184810295`*^9, 3.579845178523224*^9, 3.5798452683325815`*^9, {
   3.5798452985498343`*^9, 3.579845338938305*^9}, 3.5798453738355665`*^9, 
   3.579845436110876*^9}],

Cell[BoxData[
 RowBox[{"projection", " ", "xy"}]], "Output",
 CellChangeTimes->{
  3.5798430184810295`*^9, 3.579845178523224*^9, 3.5798452683325815`*^9, {
   3.5798452985498343`*^9, 3.579845338938305*^9}, 3.5798453738355665`*^9, 
   3.579845436110876*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"plot2", "=", 
   RowBox[{
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"Graphics", "[", 
       RowBox[{"{", 
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"a1", "[", 
             RowBox[{"[", "i", "]"}], "]"}], ",", 
            RowBox[{"a3", "[", 
             RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", 
          RowBox[{"c", "[", 
           RowBox[{"[", "i", "]"}], "]"}]}], "]"}], "}"}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "nc"}], "}"}]}], "]"}], " ", "/.", 
    RowBox[{"t", "[", 
     RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"plot2a", "=", 
   RowBox[{
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"Graphics", "[", 
       RowBox[{"{", 
        RowBox[{"Disk", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"b1", "[", 
             RowBox[{"[", "i", "]"}], "]"}], ",", 
            RowBox[{"b3", "[", 
             RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], 
        "}"}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", 
        RowBox[{"nc", "-", "2"}]}], "}"}]}], "]"}], " ", "/.", 
    RowBox[{"t", "[", 
     RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"plot2b", "=", 
   RowBox[{
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"Graphics", "[", 
       RowBox[{"{", 
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"a1", "[", 
             RowBox[{"[", "i", "]"}], "]"}], ",", 
            RowBox[{"a3", "[", 
             RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], 
        "}"}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "nc"}], "}"}]}], "]"}], " ", "/.", 
    RowBox[{"t", "[", 
     RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"Show", "[", 
  RowBox[{"plot2", ",", "plot2a", ",", "plot2b", ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"AspectRatio", " ", "\[Rule]", " ", "Automatic"}], ",", 
     RowBox[{"Axes", " ", "->", " ", "True"}], ",", 
     RowBox[{"GridLines", " ", "->", " ", "Automatic"}], ",", " ", "\n", 
     RowBox[{"PlotRange", " ", "\[Rule]", " ", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{
          RowBox[{"-", "4"}], ",", "4"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"-", "4"}], ",", "4"}], "}"}]}], "}"}]}], ",", " ", 
     RowBox[{"Frame", " ", "->", " ", "True"}]}], "}"}]}], 
  "]"}], "\[IndentingNewLine]", 
 StyleBox[
  RowBox[{"xz", " ", "projection"}], "Text"]}], "Input",
 CellChangeTimes->{{3.579339709701754*^9, 3.5793398772928486`*^9}, {
  3.579339981735032*^9, 3.5793399817506323`*^9}, {3.5793507174300737`*^9, 
  3.579350731251698*^9}, {3.5795776651831727`*^9, 3.579577668599579*^9}, {
  3.579579759954852*^9, 3.5795797795172863`*^9}, {3.579579823166163*^9, 
  3.57957986711144*^9}, {3.57957991244512*^9, 3.57957991831073*^9}, {
  3.5796624612996883`*^9, 3.579662624101574*^9}, {3.5796717987004886`*^9, 
  3.5796718037080975`*^9}, {3.579753282761011*^9, 3.5797532835098124`*^9}}],

Cell[BoxData[
 GraphicsBox[{CircleBox[{1.5487983118234183`*^-10, 0}, 0.29054649214032524`], 
   CircleBox[{-1.5487983118234183`*^-10, 0}, 0.29054649214032524`], 
   CircleBox[{-0.8234348354678755, 0}, 1.195815860774469], 
   CircleBox[{-0.8234348354678755, 0}, 1.195815860774469], 
   CircleBox[{0.8234348354678755, 0}, 1.195815860774469], 
   CircleBox[{0.8234348354678755, 0}, 1.195815860774469], 
   CircleBox[{0, 1.9856604075996878`}, 2.1856604075996877`], 
   CircleBox[{0, -1.9856604075996878`}, 2.1856604075996877`], 
   DiskBox[{1.5487983118234183`*^-10, 0}, 0.08], 
   DiskBox[{-1.5487983118234183`*^-10, 0}, 0.08], 
   DiskBox[{-1.1567741378133178`, 0}, 0.08], 
   DiskBox[{-1.1567741378133178`, 0}, 0.08], 
   DiskBox[{1.1567741378133178`, 0}, 0.08], 
   DiskBox[{1.1567741378133178`, 0}, 0.08], 
   CircleBox[{1.5487983118234183`*^-10, 0}, 0.08], 
   CircleBox[{-1.5487983118234183`*^-10, 0}, 0.08], 
   CircleBox[{-0.8234348354678755, 0}, 0.08], 
   CircleBox[{-0.8234348354678755, 0}, 0.08], 
   CircleBox[{0.8234348354678755, 0}, 0.08], 
   CircleBox[{0.8234348354678755, 0}, 0.08], 
   CircleBox[{0, 1.9856604075996878`}, 0.08], 
   CircleBox[{0, -1.9856604075996878`}, 0.08]},
  AspectRatio->Automatic,
  Axes->True,
  Frame->True,
  GridLines->Automatic,
  PlotRange->{{-4, 4}, {-4, 4}}]], "Output",
 CellChangeTimes->{
  3.57984301851223*^9, 3.5798451785388236`*^9, 3.579845268363781*^9, {
   3.5798452985810347`*^9, 3.5798453389695053`*^9}, 3.579845373866767*^9, 
   3.579845436142076*^9}],

Cell[BoxData[
 RowBox[{"projection", " ", "xz"}]], "Output",
 CellChangeTimes->{
  3.57984301851223*^9, 3.5798451785388236`*^9, 3.579845268363781*^9, {
   3.5798452985810347`*^9, 3.5798453389695053`*^9}, 3.579845373866767*^9, 
   3.579845436142076*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"plot3", "=", 
   RowBox[{
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"Graphics3D", "[", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"Opacity", "[", "0.4", "]"}], ",", 
         RowBox[{"Sphere", "[", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{
             RowBox[{"a1", "[", 
              RowBox[{"[", "i", "]"}], "]"}], ",", 
             RowBox[{"a2", "[", 
              RowBox[{"[", "i", "]"}], "]"}], ",", 
             RowBox[{"a3", "[", 
              RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", 
           RowBox[{"c", "[", 
            RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}], "}"}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", 
        RowBox[{"nc", "-", "2"}]}], "}"}]}], "]"}], " ", "/.", 
    RowBox[{"t", "[", 
     RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"plot3a", "=", 
   RowBox[{
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"Graphics3D", "[", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"Darker", "[", 
          RowBox[{"Red", ",", "1.0"}], "]"}], ",", 
         RowBox[{"Sphere", "[", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{
             RowBox[{"b1", "[", 
              RowBox[{"[", "i", "]"}], "]"}], ",", 
             RowBox[{"b2", "[", 
              RowBox[{"[", "i", "]"}], "]"}], ",", 
             RowBox[{"b3", "[", 
              RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}]}],
         "}"}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", 
        RowBox[{"nc", "-", "2"}]}], "}"}]}], "]"}], " ", "/.", 
    RowBox[{"t", "[", 
     RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"plot3b", "=", 
   RowBox[{
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"Graphics3D", "[", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"Opacity", "[", "0.6", "]"}], ",", 
         RowBox[{"Lighter", "[", 
          RowBox[{"Blue", ",", "0.7"}], "]"}], ",", 
         RowBox[{"Sphere", "[", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{
             RowBox[{"a1", "[", 
              RowBox[{"[", "i", "]"}], "]"}], ",", 
             RowBox[{"a2", "[", 
              RowBox[{"[", "i", "]"}], "]"}], ",", 
             RowBox[{"a3", "[", 
              RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", 
           RowBox[{"c", "[", 
            RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}], "}"}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", 
        RowBox[{"nc", "-", "1"}], ",", "nc"}], "}"}]}], "]"}], " ", "/.", 
    RowBox[{"t", "[", 
     RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"plot3c", "=", 
    RowBox[{
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"Graphics3D", "[", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"Opacity", "[", "0.6", "]"}], ",", 
          RowBox[{"Lighter", "[", 
           RowBox[{"Magenta", ",", "0.7"}], "]"}], ",", 
          RowBox[{"Sphere", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{"a1", "[", 
               RowBox[{"[", "i", "]"}], "]"}], ",", 
              RowBox[{"a2", "[", 
               RowBox[{"[", "i", "]"}], "]"}], ",", 
              RowBox[{"a3", "[", 
               RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", 
            RowBox[{"c", "[", 
             RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}], "}"}], "]"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "2"}], "}"}]}], "]"}], " ", "/.", 
     RowBox[{"t", "[", 
      RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{"Show", "[", 
  RowBox[{"plot3", ",", "plot3a", ",", "plot3b", ",", "plot3c", ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"AspectRatio", " ", "\[Rule]", " ", "Automatic"}], ",", 
     RowBox[{"Axes", " ", "->", " ", "True"}], ",", 
     RowBox[{"SphericalRegion", "\[Rule]", "True"}]}], "}"}]}], 
  "]"}], "\[IndentingNewLine]"}], "Input",
 CellChangeTimes->{{3.579683275360855*^9, 3.5796833411461706`*^9}, {
   3.5796833720342245`*^9, 3.5796834557283716`*^9}, {3.579683488488429*^9, 
   3.5796835178788805`*^9}, {3.5796835721981764`*^9, 
   3.5796835803413906`*^9}, {3.579683651695916*^9, 3.579683739929671*^9}, {
   3.5796837704589243`*^9, 3.5796838321258326`*^9}, {3.5796838663210926`*^9, 
   3.5796839172551823`*^9}, 3.5796839627448626`*^9, {3.579753264040978*^9, 
   3.579753264820979*^9}}],

Cell[BoxData[
 Graphics3DBox[{
   {Opacity[0.4], 
    SphereBox[{1.5487983118234183`*^-10, 0, 0}, 0.29054649214032524`]}, 
   {Opacity[0.4], 
    SphereBox[{-1.5487983118234183`*^-10, 0, 0}, 0.29054649214032524`]}, 
   {Opacity[0.4], 
    SphereBox[{-0.8234348354678755, 1.2374280244302411`, 0}, 
     1.195815860774469]}, 
   {Opacity[0.4], 
    SphereBox[{-0.8234348354678755, -1.2374280244302411`, 0}, 
     1.195815860774469]}, 
   {Opacity[0.4], 
    SphereBox[{0.8234348354678755, -1.2374280244302411`, 0}, 
     1.195815860774469]}, 
   {Opacity[0.4], 
    SphereBox[{0.8234348354678755, 1.2374280244302411`, 0}, 
     1.195815860774469]}, 
   {RGBColor[0., 0., 0.], SphereBox[{1.5487983118234183`*^-10, 0, 0}, 0.08]}, 
   {RGBColor[0., 0., 0.], SphereBox[{-1.5487983118234183`*^-10, 0, 0}, 0.08]}, 
   {RGBColor[0., 0., 0.], 
    SphereBox[{-1.1567741378133178`, 1.7383582456472222`, 0}, 0.08]}, 
   {RGBColor[0., 0., 0.], 
    SphereBox[{-1.1567741378133178`, -1.7383582456472222`, 0}, 0.08]}, 
   {RGBColor[0., 0., 0.], 
    SphereBox[{1.1567741378133178`, -1.7383582456472222`, 0}, 0.08]}, 
   {RGBColor[0., 0., 0.], 
    SphereBox[{1.1567741378133178`, 1.7383582456472222`, 0}, 0.08]}, 
   {RGBColor[0.7, 0.7, 1.], Opacity[0.6], 
    SphereBox[{0, 0, 1.9856604075996878`}, 2.1856604075996877`]}, 
   {RGBColor[0.7, 0.7, 1.], Opacity[0.6], 
    SphereBox[{0, 0, -1.9856604075996878`}, 2.1856604075996877`]}, 
   {RGBColor[1., 0.7, 1.], Opacity[0.6], 
    SphereBox[{1.5487983118234183`*^-10, 0, 0}, 0.29054649214032524`]}, 
   {RGBColor[1., 0.7, 1.], Opacity[0.6], 
    SphereBox[{-1.5487983118234183`*^-10, 0, 0}, 0.29054649214032524`]}},
  AspectRatio->Automatic,
  Axes->True,
  SphericalRegion->True]], "Output",
 CellChangeTimes->{
  3.57984301888663*^9, 3.5798451785544233`*^9, 3.579845268379381*^9, {
   3.5798452985966344`*^9, 3.5798453389851055`*^9}, 3.5798453738823667`*^9, 
   3.5798454361576757`*^9}]
}, Open  ]]
},
WindowToolbars->"EditBar",
WindowSize->{969, 1002},
WindowMargins->{{Automatic, 219}, {Automatic, 11}},
PrintingCopies->1,
PrintingPageRange->{Automatic, Automatic},
PrivateNotebookOptions->{"VersionedStylesheet"->{"Default.nb"[8.] -> False}},
CellLabelAutoDelete->True,
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, 1778, 33, 92, "Input"],
Cell[2338, 55, 2777, 50, 92, "Input"],
Cell[5118, 107, 1195, 26, 72, "Input"],
Cell[6316, 135, 5552, 164, 392, "Input"],
Cell[11871, 301, 4192, 103, 192, "Input"],
Cell[16066, 406, 3931, 113, 192, "Input"],
Cell[CellGroupData[{
Cell[20022, 523, 2999, 57, 52, "Input"],
Cell[23024, 582, 608, 12, 21, "Message"],
Cell[23635, 596, 732, 17, 55, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[24404, 618, 1542, 28, 92, "Input"],
Cell[25949, 648, 243, 4, 31, "Output"],
Cell[26195, 654, 259, 5, 31, "Output"],
Cell[26457, 661, 634, 18, 48, "Output"],
Cell[27094, 681, 785, 21, 48, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[27916, 707, 1605, 26, 31, "Input"],
Cell[29524, 735, 264, 5, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[29825, 745, 1658, 28, 31, "Input"],
Cell[31486, 775, 240, 4, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[31763, 784, 1613, 26, 31, "Input"],
Cell[33379, 812, 241, 4, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[33657, 821, 3139, 87, 132, "Input"],
Cell[36799, 910, 1656, 30, 379, "Output"],
Cell[38458, 942, 255, 5, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[38750, 952, 3291, 89, 132, "Input"],
Cell[42044, 1043, 1509, 31, 379, "Output"],
Cell[43556, 1076, 252, 5, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[43845, 1086, 4594, 124, 252, "Input"],
Cell[48442, 1212, 1923, 42, 419, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)