(* 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[     43645,       1151]
NotebookOptionsPosition[     41818,       1088]
NotebookOutlinePosition[     42356,       1108]
CellTagsIndexPosition[     42313,       1105]
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", "=", "0.0"}], ";"}], "\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}}],

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

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[{"2", "*", "d1"}], "-", "3.3"}], " ", "\[Equal]", "0"}]}], 
     "}"}], ",", 
    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}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", "53.485273664874335`"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"R1", "\[Rule]", "0.31424347442785283`"}], ",", 
     RowBox[{"R2", "\[Rule]", "1.4630429475211477`"}], ",", 
     RowBox[{"R3", "\[Rule]", "1.335756525572147`"}], ",", 
     RowBox[{"R4", "\[Rule]", "1.9929963867556548`"}], ",", 
     RowBox[{"p", "\[Rule]", "0.5894782671639446`"}], ",", 
     RowBox[{"pb", "\[Rule]", "0.8428666495739683`"}], ",", 
     RowBox[{"wi", "\[Rule]", "1.1618438662572672`"}]}], "}"}]}], 
  "}"}]], "Output",
 CellChangeTimes->{3.579365490798169*^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["44.905379362998026`"], "Output",
 CellChangeTimes->{3.5793654908137693`*^9}],

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

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

Cell[BoxData["2.012043315900975`"], "Output",
 CellChangeTimes->{3.5793654908293695`*^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.5793654908293695`*^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.252437437890721`"], "Output",
 CellChangeTimes->{3.5793654908449697`*^9}]
}, Open  ]],

Cell[CellGroupData[{

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

Cell[BoxData["1.736202169967467`"], "Output",
 CellChangeTimes->{3.5793654908605695`*^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["133.13750111085162`"], "Output",
 CellChangeTimes->{3.5793654908761697`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"2", "*", 
   RowBox[{"ArcTan", "[", 
    RowBox[{"d1", "/", 
     RowBox[{"b3", "[", 
      RowBox[{"[", "7", "]"}], "]"}]}], "]"}], "*", "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}}],

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

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"plot1", "=", 
    RowBox[{
     RowBox[{"Graphics", "[", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"a1", "[", 
             RowBox[{"[", "1", "]"}], "]"}], ",", 
            RowBox[{"a2", "[", 
             RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", "R1"}], "]"}], ",",
         "\n", 
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"a1", "[", 
             RowBox[{"[", "2", "]"}], "]"}], ",", 
            RowBox[{"a2", "[", 
             RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", "R1"}], "]"}], ",",
         "\n", 
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"a1", "[", 
             RowBox[{"[", "3", "]"}], "]"}], ",", 
            RowBox[{"a2", "[", 
             RowBox[{"[", "3", "]"}], "]"}]}], "}"}], ",", "R2"}], "]"}], ",",
         "\n", 
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"a1", "[", 
             RowBox[{"[", "4", "]"}], "]"}], ",", 
            RowBox[{"a2", "[", 
             RowBox[{"[", "4", "]"}], "]"}]}], "}"}], ",", "R2"}], "]"}], ",",
         "\n", 
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"a1", "[", 
             RowBox[{"[", "5", "]"}], "]"}], ",", 
            RowBox[{"a2", "[", 
             RowBox[{"[", "5", "]"}], "]"}]}], "}"}], ",", "R2"}], "]"}], ",",
         "\n", 
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"a1", "[", 
             RowBox[{"[", "6", "]"}], "]"}], ",", 
            RowBox[{"a2", "[", 
             RowBox[{"[", "6", "]"}], "]"}]}], "}"}], ",", "R2"}], "]"}], ",",
         "\n", 
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"a1", "[", 
             RowBox[{"[", "7", "]"}], "]"}], ",", 
            RowBox[{"a2", "[", 
             RowBox[{"[", "7", "]"}], "]"}]}], "}"}], ",", "R4"}], "]"}], ",",
         "\n", 
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"a1", "[", 
             RowBox[{"[", "8", "]"}], "]"}], ",", 
            RowBox[{"a2", "[", 
             RowBox[{"[", "8", "]"}], "]"}]}], "}"}], ",", "R4"}], "]"}], ",",
         "\n", 
        RowBox[{"Disk", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"b1", "[", 
             RowBox[{"[", "1", "]"}], "]"}], ",", 
            RowBox[{"b2", "[", 
             RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], 
        ",", " ", "\n", 
        RowBox[{"Disk", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"b1", "[", 
             RowBox[{"[", "2", "]"}], "]"}], ",", 
            RowBox[{"b2", "[", 
             RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], 
        ",", " ", "\n", 
        RowBox[{"Disk", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"b1", "[", 
             RowBox[{"[", "3", "]"}], "]"}], ",", 
            RowBox[{"b2", "[", 
             RowBox[{"[", "3", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], 
        ",", " ", "\n", 
        RowBox[{"Disk", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"b1", "[", 
             RowBox[{"[", "4", "]"}], "]"}], ",", 
            RowBox[{"b2", "[", 
             RowBox[{"[", "4", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], 
        ",", "\n", 
        RowBox[{"Disk", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"b1", "[", 
             RowBox[{"[", "5", "]"}], "]"}], ",", 
            RowBox[{"b2", "[", 
             RowBox[{"[", "5", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], 
        ",", "\n", 
        RowBox[{"Disk", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"b1", "[", 
             RowBox[{"[", "6", "]"}], "]"}], ",", 
            RowBox[{"b2", "[", 
             RowBox[{"[", "6", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}]}], 
       "}"}], " ", "]"}], " ", "/.", " ", 
     RowBox[{"t", "[", 
      RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], "\n"}], "\n", 
 RowBox[{"Show", "[", 
  RowBox[{"plot1", ",", 
   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}}],

Cell[BoxData[
 GraphicsBox[{CircleBox[{-1.65, 0}, 0.31424347442785283`], 
   CircleBox[{1.65, 0}, 0.31424347442785283`], 
   CircleBox[{-2.35673508283811, 1.630727613162324}, 1.4630429475211477`], 
   CircleBox[{-2.35673508283811, -1.630727613162324}, 1.4630429475211477`], 
   CircleBox[{2.35673508283811, -1.630727613162324}, 1.4630429475211477`], 
   CircleBox[{2.35673508283811, 1.630727613162324}, 1.4630429475211477`], 
   CircleBox[{0, 0}, 1.9929963867556548`], 
   CircleBox[{0, 0}, 1.9929963867556548`], DiskBox[{-1.65, 0}, 0.08], 
   DiskBox[{1.65, 0}, 0.08], 
   DiskBox[{-2.591140166245274, 2.1715962518644467`}, 0.08], 
   DiskBox[{-2.591140166245274, -2.1715962518644467`}, 0.08], 
   DiskBox[{2.591140166245274, -2.1715962518644467`}, 0.08], 
   DiskBox[{2.591140166245274, 2.1715962518644467`}, 0.08]},
  AspectRatio->Automatic,
  Axes->True,
  Frame->True,
  GridLines->Automatic,
  PlotRange->{{-4, 4}, {-4, 4}}]], "Output",
 CellChangeTimes->{3.5793654909073696`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"plot2", "=", 
    RowBox[{
     RowBox[{"Graphics", "[", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"a1", "[", 
             RowBox[{"[", "1", "]"}], "]"}], ",", 
            RowBox[{"a2", "[", 
             RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", "R1"}], "]"}], ",",
         "      ", "\n", 
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"a1", "[", 
             RowBox[{"[", "2", "]"}], "]"}], ",", 
            RowBox[{"a3", "[", 
             RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", "R1"}], "]"}], ",",
         "\n", 
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"a1", "[", 
             RowBox[{"[", "3", "]"}], "]"}], ",", 
            RowBox[{"a3", "[", 
             RowBox[{"[", "3", "]"}], "]"}]}], "}"}], ",", "R2"}], "]"}], ",",
         "\n", 
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"a1", "[", 
             RowBox[{"[", "4", "]"}], "]"}], ",", 
            RowBox[{"a3", "[", 
             RowBox[{"[", "4", "]"}], "]"}]}], "}"}], ",", "R2"}], "]"}], ",",
         "\n", 
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"a1", "[", 
             RowBox[{"[", "5", "]"}], "]"}], ",", 
            RowBox[{"a3", "[", 
             RowBox[{"[", "5", "]"}], "]"}]}], "}"}], ",", "R2"}], "]"}], ",",
         "\n", 
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"a1", "[", 
             RowBox[{"[", "6", "]"}], "]"}], ",", 
            RowBox[{"a3", "[", 
             RowBox[{"[", "6", "]"}], "]"}]}], "}"}], ",", "R2"}], "]"}], ",",
         "\n", 
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"a1", "[", 
             RowBox[{"[", "7", "]"}], "]"}], ",", 
            RowBox[{"a3", "[", 
             RowBox[{"[", "7", "]"}], "]"}]}], "}"}], ",", "R4"}], "]"}], ",",
         "\n", 
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"a1", "[", 
             RowBox[{"[", "8", "]"}], "]"}], ",", 
            RowBox[{"a3", "[", 
             RowBox[{"[", "8", "]"}], "]"}]}], "}"}], ",", "R4"}], "]"}], ",",
         "\n", 
        RowBox[{"Disk", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"b1", "[", 
             RowBox[{"[", "1", "]"}], "]"}], ",", 
            RowBox[{"b3", "[", 
             RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], 
        ",", " ", "\n", 
        RowBox[{"Disk", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"b1", "[", 
             RowBox[{"[", "2", "]"}], "]"}], ",", 
            RowBox[{"b3", "[", 
             RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], 
        ",", " ", "\n", 
        RowBox[{"Disk", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"b1", "[", 
             RowBox[{"[", "3", "]"}], "]"}], ",", 
            RowBox[{"b3", "[", 
             RowBox[{"[", "3", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], 
        ",", " ", "\n", 
        RowBox[{"Disk", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"b1", "[", 
             RowBox[{"[", "4", "]"}], "]"}], ",", 
            RowBox[{"b3", "[", 
             RowBox[{"[", "4", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], 
        ",", "\n", 
        RowBox[{"Disk", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"b1", "[", 
             RowBox[{"[", "5", "]"}], "]"}], ",", 
            RowBox[{"b3", "[", 
             RowBox[{"[", "5", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], 
        ",", "\n", 
        RowBox[{"Disk", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"b1", "[", 
             RowBox[{"[", "6", "]"}], "]"}], ",", 
            RowBox[{"b3", "[", 
             RowBox[{"[", "6", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], 
        ",", "\n", 
        RowBox[{"Disk", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"b1", "[", 
             RowBox[{"[", "7", "]"}], "]"}], ",", 
            RowBox[{"b3", "[", 
             RowBox[{"[", "7", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], 
        ",", "\n", 
        RowBox[{"Disk", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"b1", "[", 
             RowBox[{"[", "8", "]"}], "]"}], ",", 
            RowBox[{"b3", "[", 
             RowBox[{"[", "8", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}]}], 
       "}"}], " ", "]"}], " ", "/.", " ", 
     RowBox[{"t", "[", 
      RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], "\n"}], "\n", 
 RowBox[{"Show", "[", 
  RowBox[{"plot2", ",", 
   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[{"-", "6"}], ",", "6"}], "}"}]}], "}"}]}], ",", " ", 
     RowBox[{"Frame", " ", "->", " ", "True"}]}], "}"}]}], "]"}]}], "Input",
 CellChangeTimes->{{3.579339709701754*^9, 3.5793398772928486`*^9}, {
  3.579339981735032*^9, 3.5793399817506323`*^9}, {3.5793507174300737`*^9, 
  3.579350731251698*^9}}],

Cell[BoxData[
 GraphicsBox[{CircleBox[{-1.65, 0}, 0.31424347442785283`], 
   CircleBox[{1.65, 0}, 0.31424347442785283`], 
   CircleBox[{-2.35673508283811, 0}, 1.4630429475211477`], 
   CircleBox[{-2.35673508283811, 0}, 1.4630429475211477`], 
   CircleBox[{2.35673508283811, 0}, 1.4630429475211477`], 
   CircleBox[{2.35673508283811, 0}, 1.4630429475211477`], 
   CircleBox[{0, 1.9929963867556548`}, 1.9929963867556548`], 
   CircleBox[{0, -1.9929963867556548`}, 1.9929963867556548`], 
   DiskBox[{-1.65, 0}, 0.08], DiskBox[{1.65, 0}, 0.08], 
   DiskBox[{-2.591140166245274, 0}, 0.08], 
   DiskBox[{-2.591140166245274, 0}, 0.08], 
   DiskBox[{2.591140166245274, 0}, 0.08], 
   DiskBox[{2.591140166245274, 0}, 0.08], 
   DiskBox[{0, 2.835863036329623}, 0.08], 
   DiskBox[{0, -2.835863036329623}, 0.08]},
  AspectRatio->Automatic,
  Axes->True,
  Frame->True,
  GridLines->Automatic,
  PlotRange->{{-4, 4}, {-6, 6}}]], "Output",
 CellChangeTimes->{3.57936549092297*^9}]
}, Open  ]]
},
WindowToolbars->"EditBar",
WindowSize->{774, 826},
WindowMargins->{{Automatic, 310}, {Automatic, 9}},
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, 1704, 33, 92, "Input"],
Cell[3784, 83, 652, 12, 31, "Input"],
Cell[4439, 97, 1255, 30, 92, "Input"],
Cell[5697, 129, 5106, 155, 372, "Input"],
Cell[10806, 286, 4017, 100, 192, "Input"],
Cell[14826, 388, 1832, 50, 152, "Input"],
Cell[CellGroupData[{
Cell[16683, 442, 1863, 40, 72, "Input"],
Cell[18549, 484, 610, 14, 52, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[19196, 503, 1542, 28, 92, "Input"],
Cell[20741, 533, 90, 1, 31, "Output"],
Cell[20834, 536, 106, 2, 31, "Output"],
Cell[20943, 540, 89, 1, 31, "Output"],
Cell[21035, 543, 89, 1, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[21161, 549, 1506, 25, 31, "Input"],
Cell[22670, 576, 81, 1, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[22788, 582, 1588, 27, 31, "Input"],
Cell[24379, 611, 89, 1, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[24505, 617, 519, 16, 31, "Input"],
Cell[25027, 635, 89, 1, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[25153, 641, 1563, 25, 31, "Input"],
Cell[26719, 668, 90, 1, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[26846, 674, 1589, 27, 31, "Input"],
Cell[28438, 703, 90, 1, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[28565, 709, 5286, 154, 352, "Input"],
Cell[33854, 865, 986, 19, 379, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[34877, 889, 5955, 173, 392, "Input"],
Cell[40835, 1064, 967, 21, 447, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)