(* 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[     45493,       1074]
NotebookOptionsPosition[     43654,       1011]
NotebookOutlinePosition[     44193,       1031]
CellTagsIndexPosition[     44150,       1028]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{"B2H6", " ", "with", " ", "spherical", " ", "bananas"}], ",", " ", 
    RowBox[{"but", " ", "not", " ", "TSM", "  ", "28.06", ".2012"}]}], " ", 
   "*)"}], "\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", ",", "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}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"z", "=", "5.0"}], ";", "  ", 
  RowBox[{"nc", "=", "8"}], ";", " ", 
  RowBox[{"pi", "=", 
   RowBox[{"-", "0.6"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"bohr", "=", "0.529177"}], ";", " ", 
  RowBox[{"rad", "=", "57.29578"}], ";"}], "\n", 
 RowBox[{
  RowBox[{"k", "=", 
   RowBox[{"{", 
    RowBox[{
    "1.0", ",", "1.0", ",", "1.3", ",", "1.3", ",", "1.3", ",", "1.3", ",", 
     "1.6", ",", "1.6"}], "}"}]}], ";"}], "\[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}}],

Cell[BoxData[
 RowBox[{"(*", " ", 
  RowBox[{
   RowBox[{"C", " ", "He", " ", "shell"}], ",", " ", 
   RowBox[{"i", "=", 
    RowBox[{"j", " ", "each"}]}]}], " ", "*)"}]], "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}}],

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"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"R4", "=", "R2"}], ";"}], "*)"}]}]}], "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}],

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", ",", "1", ",", 
     "1"}], "}"}]}], ";"}], "\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", ",", "pb", ",", 
     "pb"}], "}"}]}], ";"}], "\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"}]}], ";", " ", 
   RowBox[{"d6", "=", 
    RowBox[{"R4", "+", "pb"}]}], ";"}], "\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", ",", 
      RowBox[{"pi", "+", "d6"}], ",", 
      RowBox[{
       RowBox[{"-", "pi"}], "-", "d6"}]}], "}"}]}], ";"}], "        "}], "\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", ",", 
     RowBox[{"pi", "+", "R4"}], ",", 
     RowBox[{
      RowBox[{"-", "pi"}], "-", "R4"}]}], "}"}]}], ";"}], "\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}}],

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[{"i", " ", "\[NotEqual]", " ", "j"}], ",", " ", 
            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}}],

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[{
   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}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"t", " ", "=", " ", 
  RowBox[{"NMinimize", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"func", ",", 
      RowBox[{
       RowBox[{
        RowBox[{
         RowBox[{"2", "*", "d1"}], "-", "3.3"}], " ", "\[Equal]", "0"}], " ", 
       "&&", " ", 
       RowBox[{
        RowBox[{"R1", "+", "R2", "+", "p"}], " ", "\[Equal]", " ", "2.249"}], 
       " ", "&&", " ", 
       RowBox[{
        RowBox[{"d1", "/", 
         RowBox[{"(", 
          RowBox[{"pi", "+", "d6"}], ")"}]}], "\[Equal]", "0.89946"}], " ", "&&",
        " ", 
       RowBox[{"wi", "\[Equal]", "1.062"}]}]}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{"R1", ",", "0.31", ",", "0.332"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"R2", ",", "1.49", ",", "1.56"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"R3", ",", "0.7", ",", "0.9"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"R4", ",", "0.75", ",", "1.80"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"p", ",", "0.59", ",", "0.61"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"pb", ",", "0.78", ",", "0.80"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"wi", ",", "1.14", ",", "1.15"}], "}"}]}], "}"}]}], 
   "]"}]}]], "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}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", "54.07843290619366`"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"R1", "\[Rule]", "0.31420694646188174`"}], ",", 
     RowBox[{"R2", "\[Rule]", "1.4384650602864155`"}], ",", 
     RowBox[{"R3", "\[Rule]", "1.3357930535381182`"}], ",", 
     RowBox[{"R4", "\[Rule]", "1.8287287944529025`"}], ",", 
     RowBox[{"p", "\[Rule]", "0.49632799325170296`"}], ",", 
     RowBox[{"pb", "\[Rule]", "0.6057052042741818`"}], ",", 
     RowBox[{"wi", "\[Rule]", "1.062`"}]}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{
  3.57957650798874*^9, 3.579576869394575*^9, 3.579576933510688*^9, 
   3.5795769848035774`*^9, 3.5795771258122253`*^9, 3.5795773313425865`*^9, 
   3.5795773667390485`*^9, 3.579577400606708*^9, 3.5795774517903976`*^9, 
   3.5795775558269806`*^9, 3.5795776291003094`*^9, 3.5795776737007875`*^9, {
   3.57957792559443*^9, 3.5795779492908716`*^9}, 3.5795779886497407`*^9, 
   3.57957803343742*^9, 3.579579210085086*^9, 3.579579697383142*^9, 
   3.5795797319372025`*^9, 3.5795798842246704`*^9, 3.5795800145784993`*^9, 
   3.579580212121646*^9, 3.5795802717605505`*^9, 3.5795803440822783`*^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["48.20832518742991`"], "Output",
 CellChangeTimes->{
  3.57957650798874*^9, 3.5795768694101753`*^9, 3.5795769335262876`*^9, 
   3.5795769848191776`*^9, 3.5795771258278255`*^9, 3.5795773313581867`*^9, 
   3.5795773667546487`*^9, 3.579577400622308*^9, 3.579577451805998*^9, 
   3.5795775558269806`*^9, 3.5795776291159096`*^9, 3.5795776737007875`*^9, {
   3.57957792561003*^9, 3.5795779492908716`*^9}, 3.579577988665341*^9, 
   3.57957803343742*^9, 3.579579210100686*^9, 3.579579697383142*^9, 
   3.5795797319528027`*^9, 3.5795798842402706`*^9, 3.579580014594099*^9, 
   3.579580212137246*^9, 3.5795802717761507`*^9, 3.579580344097878*^9}],

Cell[BoxData[
 RowBox[{"-", "187.9279763777267`"}]], "Output",
 CellChangeTimes->{
  3.57957650798874*^9, 3.5795768694101753`*^9, 3.5795769335262876`*^9, 
   3.5795769848191776`*^9, 3.5795771258278255`*^9, 3.5795773313581867`*^9, 
   3.5795773667546487`*^9, 3.579577400622308*^9, 3.579577451805998*^9, 
   3.5795775558269806`*^9, 3.5795776291159096`*^9, 3.5795776737007875`*^9, {
   3.57957792561003*^9, 3.5795779492908716`*^9}, 3.579577988665341*^9, 
   3.57957803343742*^9, 3.579579210100686*^9, 3.579579697383142*^9, 
   3.5795797319528027`*^9, 3.5795798842402706`*^9, 3.579580014594099*^9, 
   3.579580212137246*^9, 3.5795802717761507`*^9, 3.579580344097878*^9}],

Cell[BoxData["32.25316376114499`"], "Output",
 CellChangeTimes->{
  3.57957650798874*^9, 3.5795768694101753`*^9, 3.5795769335262876`*^9, 
   3.5795769848191776`*^9, 3.5795771258278255`*^9, 3.5795773313581867`*^9, 
   3.5795773667546487`*^9, 3.579577400622308*^9, 3.579577451805998*^9, 
   3.5795775558269806`*^9, 3.5795776291159096`*^9, 3.5795776737007875`*^9, {
   3.57957792561003*^9, 3.5795779492908716`*^9}, 3.579577988665341*^9, 
   3.57957803343742*^9, 3.579579210100686*^9, 3.579579697383142*^9, 
   3.5795797319528027`*^9, 3.5795798842402706`*^9, 3.579580014594099*^9, 
   3.579580212137246*^9, 3.5795802717761507`*^9, 3.579580344097878*^9}],

Cell[BoxData["2.012931326855872`"], "Output",
 CellChangeTimes->{
  3.57957650798874*^9, 3.5795768694101753`*^9, 3.5795769335262876`*^9, 
   3.5795769848191776`*^9, 3.5795771258278255`*^9, 3.5795773313581867`*^9, 
   3.5795773667546487`*^9, 3.579577400622308*^9, 3.579577451805998*^9, 
   3.5795775558269806`*^9, 3.5795776291159096`*^9, 3.5795776737007875`*^9, {
   3.57957792561003*^9, 3.5795779492908716`*^9}, 3.579577988665341*^9, 
   3.57957803343742*^9, 3.579579210100686*^9, 3.579579697383142*^9, 
   3.5795797319528027`*^9, 3.5795798842402706`*^9, 3.579580014594099*^9, 
   3.579580212137246*^9, 3.5795802717761507`*^9, 3.579580344097878*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"2", "*", "d1", "*", "bohr"}], " ", "/.", " ", 
  RowBox[{"t", "[", 
   RowBox[{"[", "2", "]"}], "]"}], 
  "                                                     ", 
  RowBox[{"(*", " ", 
   RowBox[{"B", "-", 
    RowBox[{"B", " ", "1.753"}]}], " ", "*)"}]}]], "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}}],

Cell[BoxData["1.7462841`"], "Output",
 CellChangeTimes->{
  3.5795765080199404`*^9, 3.5795768694257755`*^9, 3.5795769335418873`*^9, 
   3.579576984834778*^9, 3.5795771258434258`*^9, 3.5795773313737864`*^9, 
   3.5795773667702484`*^9, 3.579577400637908*^9, 3.579577451821598*^9, 
   3.579577555842581*^9, 3.5795776291159096`*^9, 3.5795776737163877`*^9, {
   3.5795779256256304`*^9, 3.579577949306472*^9}, 3.579577988665341*^9, 
   3.5795780334530196`*^9, 3.5795792101162863`*^9, 3.5795796973987417`*^9, 
   3.579579731968403*^9, 3.5795798842558703`*^9, 3.5795800146096992`*^9, 
   3.5795802121528463`*^9, 3.579580271791751*^9, 3.5795803441134777`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"(", 
    RowBox[{"R1", "+", "R2", "+", "p"}], ")"}], "*", "bohr"}], " ", "/.", " ", 
  RowBox[{"t", "[", 
   RowBox[{"[", "2", "]"}], "]"}], "                                     ", 
  RowBox[{"(*", " ", 
   RowBox[{"B", "-", 
    RowBox[{"H2", " ", "1.19"}]}], " ", "*)"}]}]], "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}}],

Cell[BoxData["1.190119073`"], "Output",
 CellChangeTimes->{
  3.5795765080199404`*^9, 3.5795768694257755`*^9, 3.5795769335418873`*^9, 
   3.579576984850378*^9, 3.5795771258434258`*^9, 3.5795773313737864`*^9, 
   3.5795773667858486`*^9, 3.579577400653508*^9, 3.579577451821598*^9, 
   3.579577555858181*^9, 3.5795776291315093`*^9, 3.579577673731988*^9, {
   3.57957792564123*^9, 3.579577949322072*^9}, 3.5795779886809406`*^9, 
   3.5795780334686193`*^9, 3.5795792101162863`*^9, 3.579579697414342*^9, 
   3.579579731984003*^9, 3.57957988427147*^9, 3.5795800146252995`*^9, 
   3.5795802121684465`*^9, 3.579580271807351*^9, 3.579580344129078*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"Sqrt", "[", 
    RowBox[{
     RowBox[{"d1", "^", "2"}], "+", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"pi", "+", "d6"}], ")"}], "^", "2"}]}], "]"}], "*", "bohr"}], 
  " ", "/.", " ", 
  RowBox[{"t", "[", 
   RowBox[{"[", "2", "]"}], "]"}], "    ", 
  RowBox[{"(*", " ", 
   RowBox[{"B", "-", 
    RowBox[{"H1", "  ", "1.31"}]}], " ", "*)"}], " "}]], "Input",
 CellChangeTimes->{{3.5793649976657033`*^9, 3.5793650337017665`*^9}, {
  3.579446074950535*^9, 3.5794460764637375`*^9}, {3.579580323677442*^9, 
  3.5795803312434554`*^9}}],

Cell[BoxData["1.3056468630426206`"], "Output",
 CellChangeTimes->{
  3.5795765080355406`*^9, 3.5795768694413753`*^9, 3.5795769335574875`*^9, 
   3.579576984850378*^9, 3.579577125859026*^9, 3.579577331389386*^9, 
   3.5795773667858486`*^9, 3.579577400653508*^9, 3.5795774518371983`*^9, 
   3.579577555858181*^9, 3.5795776291471095`*^9, 3.579577673731988*^9, {
   3.57957792564123*^9, 3.5795779493376718`*^9}, 3.579577988696541*^9, 
   3.5795780334842196`*^9, 3.5795792101318865`*^9, 3.579579697429942*^9, 
   3.579579731984003*^9, 3.57957988427147*^9, 3.579580014640899*^9, 
   3.5795802121840467`*^9, 3.579580271807351*^9, 3.579580344144678*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"2", "*", "wi", "*", "rad"}], " ", "/.", " ", 
  RowBox[{"t", "[", 
   RowBox[{"[", "2", "]"}], "]"}], 
  "                                                         ", 
  RowBox[{"(*", " ", 
   RowBox[{"H2BH2", " ", "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}}],

Cell[BoxData["121.69623672`"], "Output",
 CellChangeTimes->{
  3.5795765080511403`*^9, 3.579576869456975*^9, 3.5795769335730877`*^9, 
   3.579576984865978*^9, 3.5795771258746257`*^9, 3.5795773314049864`*^9, 
   3.579577366801449*^9, 3.579577400669108*^9, 3.579577451852798*^9, 
   3.5795775558737807`*^9, 3.5795776291627097`*^9, 3.5795776737475877`*^9, {
   3.5795779256568303`*^9, 3.5795779493376718`*^9}, 3.579577988712141*^9, 
   3.5795780334842196`*^9, 3.5795792101474867`*^9, 3.5795796974455423`*^9, 
   3.579579731999603*^9, 3.5795798842870703`*^9, 3.579580014640899*^9, 
   3.5795802121996465`*^9, 3.5795802718229513`*^9, 3.5795803441602783`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"2", "*", 
   RowBox[{"ArcTan", "[", 
    RowBox[{"d1", "/", 
     RowBox[{"(", 
      RowBox[{"pi", "+", "d6"}], ")"}]}], "]"}], "*", "rad"}], " ", "/.", " ", 
  RowBox[{"t", "[", 
   RowBox[{"[", "2", "]"}], "]"}], "             ", 
  RowBox[{"(*", " ", 
   RowBox[{"BH1B", " ", "83.94"}], " ", "*)"}], " "}]], "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.5793647690629015`*^9}, {3.5793650698938303`*^9, 3.579365073965437*^9}, {
   3.579446086369755*^9, 3.579446107585792*^9}, {3.579580259296129*^9, 
   3.5795802657857404`*^9}}],

Cell[BoxData["83.940228831732`"], "Output",
 CellChangeTimes->{
  3.5795765080511403`*^9, 3.579576869456975*^9, 3.5795769335730877`*^9, 
   3.579576984881578*^9, 3.5795771258746257`*^9, 3.5795773314049864`*^9, 
   3.579577366817049*^9, 3.579577400684708*^9, 3.579577451852798*^9, 
   3.5795775558893805`*^9, 3.5795776291627097`*^9, 3.579577673763188*^9, {
   3.5795779256724305`*^9, 3.579577949353272*^9}, 3.579577988712141*^9, 
   3.5795780334998198`*^9, 3.5795792101630864`*^9, 3.5795796974611425`*^9, 
   3.5795797320152025`*^9, 3.5795798843026705`*^9, 3.5795800146564994`*^9, 
   3.579580212215246*^9, 3.5795802718385515`*^9, 3.5795803441602783`*^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[{"Show", "[", 
  RowBox[{"plot1", ",", "plot1a", ",", 
   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"}]}], "}"}]}], "]"}]}], "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}}],

Cell[BoxData[
 GraphicsBox[{CircleBox[{-1.65, 0}, 0.31420694646188174`], 
   CircleBox[{1.65, 0}, 0.31420694646188174`], 
   CircleBox[{-2.503772795064379, 1.53066363961749}, 1.4384650602864155`], 
   CircleBox[{-2.503772795064379, -1.53066363961749}, 1.4384650602864155`], 
   CircleBox[{2.503772795064379, -1.53066363961749}, 1.4384650602864155`], 
   CircleBox[{2.503772795064379, 1.53066363961749}, 1.4384650602864155`], 
   CircleBox[{0, 0}, 1.8287287944529025`], 
   CircleBox[{0, 0}, 1.8287287944529025`], DiskBox[{-1.65, 0}, 0.08], 
   DiskBox[{1.65, 0}, 0.08], 
   DiskBox[{-2.7455472608147504`, 1.9641225010984673`}, 0.08], 
   DiskBox[{-2.7455472608147504`, -1.9641225010984673`}, 0.08], 
   DiskBox[{2.7455472608147504`, -1.9641225010984673`}, 0.08], 
   DiskBox[{2.7455472608147504`, 1.9641225010984673`}, 0.08]},
  AspectRatio->Automatic,
  Axes->True,
  Frame->True,
  GridLines->Automatic,
  PlotRange->{{-4, 4}, {-4, 4}}]], "Output",
 CellChangeTimes->{
  3.5795765080667405`*^9, 3.579576869472575*^9, 3.579576933588688*^9, 
   3.5795769848971777`*^9, 3.5795771259058256`*^9, 3.5795773314205866`*^9, 
   3.5795773668326488`*^9, 3.579577400700308*^9, 3.5795774518683977`*^9, 
   3.5795775559049807`*^9, 3.5795776291783094`*^9, 3.5795776737787876`*^9, {
   3.5795779256880302`*^9, 3.5795779493688717`*^9}, 3.5795779887277412`*^9, 
   3.5795780335154195`*^9, 3.579579210178686*^9, 3.5795796974767427`*^9, 
   3.5795797320308027`*^9, 3.5795798843026705`*^9, 3.579580014672099*^9, 
   3.5795802122308464`*^9, 3.5795802718541512`*^9, 3.579580344175878*^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", ",", "nc"}], "}"}]}], "]"}], " ", "/.", 
    RowBox[{"t", "[", 
     RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"Show", "[", 
  RowBox[{"plot2", ",", "plot2a", ",", 
   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"}]}], "}"}]}], "]"}]}], "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}}],

Cell[BoxData[
 GraphicsBox[{CircleBox[{-1.65, 0}, 0.31420694646188174`], 
   CircleBox[{1.65, 0}, 0.31420694646188174`], 
   CircleBox[{-2.503772795064379, 0}, 1.4384650602864155`], 
   CircleBox[{-2.503772795064379, 0}, 1.4384650602864155`], 
   CircleBox[{2.503772795064379, 0}, 1.4384650602864155`], 
   CircleBox[{2.503772795064379, 0}, 1.4384650602864155`], 
   CircleBox[{0, 1.2287287944529024`}, 1.8287287944529025`], 
   CircleBox[{0, -1.2287287944529024`}, 1.8287287944529025`], 
   DiskBox[{-1.65, 0}, 0.08], DiskBox[{1.65, 0}, 0.08], 
   DiskBox[{-2.7455472608147504`, 0}, 0.08], 
   DiskBox[{-2.7455472608147504`, 0}, 0.08], 
   DiskBox[{2.7455472608147504`, 0}, 0.08], 
   DiskBox[{2.7455472608147504`, 0}, 0.08], 
   DiskBox[{0, 1.8344339987270843`}, 0.08], 
   DiskBox[{0, -1.8344339987270843`}, 0.08]},
  AspectRatio->Automatic,
  Axes->True,
  Frame->True,
  GridLines->Automatic,
  PlotRange->{{-4, 4}, {-4, 4}}]], "Output",
 CellChangeTimes->{
  3.5795765080823402`*^9, 3.5795768694881754`*^9, 3.579576933604288*^9, 
   3.579576984912778*^9, 3.579577125921426*^9, 3.5795773314361863`*^9, 
   3.5795773668482485`*^9, 3.579577400715908*^9, 3.579577451883998*^9, 
   3.579577555920581*^9, 3.5795776291939096`*^9, 3.579577673794388*^9, {
   3.5795779257036304`*^9, 3.579577949384472*^9}, 3.579577988743341*^9, 
   3.5795780335310197`*^9, 3.5795792101942863`*^9, 3.5795796974923425`*^9, 
   3.579579732046403*^9, 3.5795798843338704`*^9, 3.5795800146876993`*^9, 
   3.5795802122464466`*^9, 3.579580271869751*^9, 3.579580344207078*^9}]
}, Open  ]]
},
WindowToolbars->"EditBar",
WindowSize->{774, 826},
WindowMargins->{{559, Automatic}, {69, Automatic}},
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, 1517, 26, 112, "Input"],
Cell[2077, 48, 2221, 41, 92, "Input"],
Cell[4301, 91, 652, 12, 31, "Input"],
Cell[4956, 105, 1255, 30, 92, "Input"],
Cell[6214, 137, 5309, 158, 372, "Input"],
Cell[11526, 297, 4017, 100, 192, "Input"],
Cell[15546, 399, 1832, 50, 152, "Input"],
Cell[CellGroupData[{
Cell[17403, 453, 2632, 55, 112, "Input"],
Cell[20038, 510, 1161, 21, 52, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[21236, 536, 1542, 28, 92, "Input"],
Cell[22781, 566, 649, 9, 31, "Output"],
Cell[23433, 577, 666, 10, 31, "Output"],
Cell[24102, 589, 649, 9, 31, "Output"],
Cell[24754, 600, 649, 9, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[25440, 614, 1506, 25, 31, "Input"],
Cell[26949, 641, 650, 9, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[27636, 655, 1588, 27, 31, "Input"],
Cell[29227, 684, 642, 9, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[29906, 698, 587, 17, 31, "Input"],
Cell[30496, 717, 645, 9, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[31178, 731, 1563, 25, 31, "Input"],
Cell[32744, 758, 653, 9, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[33434, 772, 1687, 29, 31, "Input"],
Cell[35124, 803, 654, 9, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[35815, 817, 2297, 63, 112, "Input"],
Cell[38115, 882, 1567, 27, 379, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[39719, 914, 2369, 63, 112, "Input"],
Cell[42091, 979, 1547, 29, 411, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)