(* 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[     41869,       1094]
NotebookOptionsPosition[     40398,       1046]
NotebookOutlinePosition[     40937,       1066]
CellTagsIndexPosition[     40894,       1063]
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.2"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"bohr", "=", "0.529177"}], ";", " ", 
  RowBox[{"rad", "=", "57.29578"}], ";"}], "\n", 
 RowBox[{
  RowBox[{"k", "=", 
   RowBox[{"{", 
    RowBox[{
    "1.0", ",", "1.0", ",", "1.31", ",", "1.31", ",", "1.31", ",", "1.31", 
     ",", "3.0", ",", "3.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}}],

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[{"(*", 
    RowBox[{"Sqrt", "[", 
     RowBox[{
      RowBox[{"R1", "^", "2"}], "+", 
      RowBox[{"2", "*", "R1", "*", "R4"}]}], "]"}], "*)"}], ";", " ", 
   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}}],

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

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", "52.60805971652283`"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"R1", "\[Rule]", "0.3137393360956287`"}], ",", 
     RowBox[{"R2", "\[Rule]", "1.4357366244058762`"}], ",", 
     RowBox[{"R3", "\[Rule]", "1.558674140616056`"}], ",", 
     RowBox[{"R4", "\[Rule]", "3.223334654981386`"}], ",", 
     RowBox[{"p", "\[Rule]", "0.608978765258157`"}], ",", 
     RowBox[{"pb", "\[Rule]", "1.24819953584896`"}], ",", 
     RowBox[{"wi", "\[Rule]", "1.3501946616782374`"}]}], "}"}]}], 
  "}"}]], "Output",
 CellChangeTimes->{
  3.579345382600518*^9, 3.5793454545010443`*^9, 3.5793500527441063`*^9, 
   3.5793501522722807`*^9, 3.5793502752784967`*^9, 3.5793504157255435`*^9, 
   3.5793505176249228`*^9, 3.5793506655911827`*^9, {3.5793507235452843`*^9, 
   3.579350735385705*^9}, {3.5793508330262766`*^9, 3.5793508582203207`*^9}, 
   3.579350901057996*^9, 3.579350940463665*^9, 3.579350987076547*^9, 
   3.5793510186198025`*^9, 3.5793510499290576`*^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", "]"}], "]"}]}], "\n", 
 RowBox[{
  RowBox[{"2", "*", "d1", "*", "bohr"}], " ", "/.", " ", 
  RowBox[{"t", "[", 
   RowBox[{"[", "2", "]"}], "]"}]}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"(", 
    RowBox[{"R1", "+", "R2", "+", "p"}], ")"}], "*", "bohr"}], " ", "/.", " ", 
  RowBox[{"t", "[", 
   RowBox[{"[", "2", "]"}], "]"}]}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"Sqrt", "[", 
    RowBox[{
     RowBox[{"d1", "^", "2"}], "+", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"b3", "[", 
        RowBox[{"[", "7", "]"}], "]"}], ")"}], "^", "2"}]}], "]"}], "*", 
   "bohr"}], " ", "/.", " ", 
  RowBox[{"t", "[", 
   RowBox[{"[", "2", "]"}], "]"}]}], "\n", 
 RowBox[{
  RowBox[{"2", "*", "wi", "*", "rad"}], " ", "/.", " ", 
  RowBox[{"t", "[", 
   RowBox[{"[", "2", "]"}], "]"}]}], "\n", 
 RowBox[{
  RowBox[{"2", "*", 
   RowBox[{"ArcTan", "[", 
    RowBox[{"d1", "/", 
     RowBox[{"b3", "[", 
      RowBox[{"[", "7", "]"}], "]"}]}], "]"}], "*", "rad"}], " ", "/.", " ", 
  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}}],

Cell[BoxData["41.019857202335345`"], "Output",
 CellChangeTimes->{
  3.579345382616118*^9, 3.5793454545166445`*^9, 3.5793500527597065`*^9, 
   3.5793501522722807`*^9, 3.579350275294097*^9, 3.5793504157411437`*^9, 
   3.579350517640523*^9, 3.579350665606783*^9, {3.5793507235452843`*^9, 
   3.579350735401305*^9}, {3.579350833041877*^9, 3.579350858235921*^9}, 
   3.579350901073596*^9, 3.579350940479265*^9, 3.579350987076547*^9, 
   3.5793510186354027`*^9, 3.5793510499446573`*^9}],

Cell[BoxData[
 RowBox[{"-", "172.4786008470017`"}]], "Output",
 CellChangeTimes->{
  3.579345382616118*^9, 3.5793454545166445`*^9, 3.5793500527597065`*^9, 
   3.5793501522722807`*^9, 3.579350275294097*^9, 3.5793504157411437`*^9, 
   3.579350517640523*^9, 3.579350665606783*^9, {3.5793507235452843`*^9, 
   3.579350735401305*^9}, {3.579350833041877*^9, 3.579350858235921*^9}, 
   3.579350901073596*^9, 3.579350940479265*^9, 3.579350987076547*^9, 
   3.5793510186354027`*^9, 3.5793510499446573`*^9}],

Cell[BoxData["26.115094273409756`"], "Output",
 CellChangeTimes->{
  3.579345382616118*^9, 3.5793454545166445`*^9, 3.5793500527597065`*^9, 
   3.5793501522722807`*^9, 3.579350275294097*^9, 3.5793504157411437`*^9, 
   3.579350517640523*^9, 3.579350665606783*^9, {3.5793507235452843`*^9, 
   3.579350735401305*^9}, {3.579350833041877*^9, 3.579350858235921*^9}, 
   3.579350901073596*^9, 3.579350940479265*^9, 3.579350987076547*^9, 
   3.5793510186354027`*^9, 3.5793510499446573`*^9}],

Cell[BoxData["1.9975817102066002`"], "Output",
 CellChangeTimes->{
  3.579345382616118*^9, 3.5793454545166445`*^9, 3.5793500527597065`*^9, 
   3.5793501522722807`*^9, 3.579350275294097*^9, 3.5793504157411437`*^9, 
   3.579350517640523*^9, 3.579350665606783*^9, {3.5793507235452843`*^9, 
   3.579350735401305*^9}, {3.579350833041877*^9, 3.579350858235921*^9}, 
   3.579350901073596*^9, 3.579350940479265*^9, 3.579350987076547*^9, 
   3.5793510186354027`*^9, 3.5793510499446573`*^9}],

Cell[BoxData["1.9816762927317184`"], "Output",
 CellChangeTimes->{
  3.579345382616118*^9, 3.5793454545166445`*^9, 3.5793500527597065`*^9, 
   3.5793501522722807`*^9, 3.579350275294097*^9, 3.5793504157411437`*^9, 
   3.579350517640523*^9, 3.579350665606783*^9, {3.5793507235452843`*^9, 
   3.579350735401305*^9}, {3.579350833041877*^9, 3.579350858235921*^9}, 
   3.579350901073596*^9, 3.579350940479265*^9, 3.579350987076547*^9, 
   3.5793510186354027`*^9, 3.5793510499446573`*^9}],

Cell[BoxData["1.2480399964133209`"], "Output",
 CellChangeTimes->{
  3.579345382616118*^9, 3.5793454545166445`*^9, 3.5793500527597065`*^9, 
   3.5793501522722807`*^9, 3.579350275294097*^9, 3.5793504157411437`*^9, 
   3.579350517640523*^9, 3.579350665606783*^9, {3.5793507235452843`*^9, 
   3.579350735401305*^9}, {3.579350833041877*^9, 3.579350858235921*^9}, 
   3.579350901073596*^9, 3.579350940479265*^9, 3.579350987076547*^9, 
   3.5793510186354027`*^9, 3.5793510499446573`*^9}],

Cell[BoxData["2.565310716484888`"], "Output",
 CellChangeTimes->{
  3.579345382616118*^9, 3.5793454545166445`*^9, 3.5793500527597065`*^9, 
   3.5793501522722807`*^9, 3.579350275294097*^9, 3.5793504157411437`*^9, 
   3.579350517640523*^9, 3.579350665606783*^9, {3.5793507235452843`*^9, 
   3.579350735401305*^9}, {3.579350833041877*^9, 3.579350858235921*^9}, 
   3.579350901073596*^9, 3.579350940479265*^9, 3.579350987076547*^9, 
   3.5793510186354027`*^9, 3.5793510499446573`*^9}],

Cell[BoxData["154.72091258538146`"], "Output",
 CellChangeTimes->{
  3.579345382616118*^9, 3.5793454545166445`*^9, 3.5793500527597065`*^9, 
   3.5793501522722807`*^9, 3.579350275294097*^9, 3.5793504157411437`*^9, 
   3.579350517640523*^9, 3.579350665606783*^9, {3.5793507235452843`*^9, 
   3.579350735401305*^9}, {3.579350833041877*^9, 3.579350858235921*^9}, 
   3.579350901073596*^9, 3.579350940479265*^9, 3.579350987076547*^9, 
   3.5793510186354027`*^9, 3.5793510499446573`*^9}],

Cell[BoxData["45.442092587967934`"], "Output",
 CellChangeTimes->{
  3.579345382616118*^9, 3.5793454545166445`*^9, 3.5793500527597065`*^9, 
   3.5793501522722807`*^9, 3.579350275294097*^9, 3.5793504157411437`*^9, 
   3.579350517640523*^9, 3.579350665606783*^9, {3.5793507235452843`*^9, 
   3.579350735401305*^9}, {3.579350833041877*^9, 3.579350858235921*^9}, 
   3.579350901073596*^9, 3.579350940479265*^9, 3.579350987076547*^9, 
   3.5793510186354027`*^9, 3.5793510499446573`*^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.8724134767116847`, 0}, 0.3137393360956287], 
   CircleBox[{1.8724134767116847`, 0}, 0.3137393360956287], 
   CircleBox[{-2.2552281153406697`, 1.707079110300405}, 1.4357366244058762`], 
   CircleBox[{-2.2552281153406697`, -1.707079110300405}, 1.4357366244058762`],
    CircleBox[{2.2552281153406697`, -1.707079110300405}, 1.4357366244058762`],
    CircleBox[{2.2552281153406697`, 1.707079110300405}, 1.4357366244058762`], 
   CircleBox[{0, 0}, 3.223334654981386], CircleBox[{0, 0}, 3.223334654981386],
    DiskBox[{-1.8724134767116847`, 0}, 0.08], 
   DiskBox[{1.8724134767116847`, 0}, 0.08], 
   DiskBox[{-2.388482867745612, 2.301299866835253}, 0.08], 
   DiskBox[{-2.388482867745612, -2.301299866835253}, 0.08], 
   DiskBox[{2.388482867745612, -2.301299866835253}, 0.08], 
   DiskBox[{2.388482867745612, 2.301299866835253}, 0.08]},
  AspectRatio->Automatic,
  Axes->True,
  Frame->True,
  GridLines->Automatic,
  PlotRange->{{-4, 4}, {-4, 4}}]], "Output",
 CellChangeTimes->{
  3.5793453826317177`*^9, 3.5793454545478444`*^9, 3.579350052790906*^9, 
   3.579350152303481*^9, 3.579350275309697*^9, 3.5793504157723436`*^9, 
   3.579350517671723*^9, 3.5793506656223826`*^9, {3.5793507235764847`*^9, 
   3.5793507354169054`*^9}, {3.5793508330730767`*^9, 3.579350858251521*^9}, 
   3.5793509011047964`*^9, 3.5793509405104656`*^9, 3.579350987107747*^9, 
   3.579351018666603*^9, 3.5793510499602575`*^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.8724134767116847`, 0}, 0.3137393360956287], 
   CircleBox[{1.8724134767116847`, 0}, 0.3137393360956287], 
   CircleBox[{-2.2552281153406697`, 0}, 1.4357366244058762`], 
   CircleBox[{-2.2552281153406697`, 0}, 1.4357366244058762`], 
   CircleBox[{2.2552281153406697`, 0}, 1.4357366244058762`], 
   CircleBox[{2.2552281153406697`, 0}, 1.4357366244058762`], 
   CircleBox[{0, 3.0233346549813858`}, 3.223334654981386], 
   CircleBox[{0, -3.0233346549813858`}, 3.223334654981386], 
   DiskBox[{-1.8724134767116847`, 0}, 0.08], 
   DiskBox[{1.8724134767116847`, 0}, 0.08], 
   DiskBox[{-2.388482867745612, 0}, 0.08], 
   DiskBox[{-2.388482867745612, 0}, 0.08], 
   DiskBox[{2.388482867745612, 0}, 0.08], 
   DiskBox[{2.388482867745612, 0}, 0.08], 
   DiskBox[{0, 4.471534190830345}, 0.08], 
   DiskBox[{0, -4.471534190830345}, 0.08]},
  AspectRatio->Automatic,
  Axes->True,
  Frame->True,
  GridLines->Automatic,
  PlotRange->{{-4, 4}, {-6, 6}}]], "Output",
 CellChangeTimes->{
  3.579345382662918*^9, 3.579345454563444*^9, 3.579350052806506*^9, 
   3.5793501523190813`*^9, 3.5793502753252974`*^9, 3.579350415787944*^9, 
   3.5793505176873226`*^9, 3.5793506656535826`*^9, {3.5793507235920844`*^9, 
   3.579350735432505*^9}, {3.5793508330886765`*^9, 3.5793508582671213`*^9}, 
   3.5793509011047964`*^9, 3.5793509405260653`*^9, 3.5793509871233473`*^9, 
   3.579351018682203*^9, 3.5793510499758577`*^9}]
}, Open  ]]
},
WindowToolbars->"EditBar",
WindowSize->{774, 826},
WindowMargins->{{Automatic, 503}, {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, 1532, 32, 92, "Input"],
Cell[3612, 82, 652, 12, 31, "Input"],
Cell[4267, 96, 1255, 30, 92, "Input"],
Cell[5525, 128, 5241, 160, 372, "Input"],
Cell[10769, 290, 4017, 100, 192, "Input"],
Cell[14789, 392, 1832, 50, 152, "Input"],
Cell[CellGroupData[{
Cell[16646, 446, 1655, 33, 52, "Input"],
Cell[18304, 481, 1007, 20, 52, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[19348, 506, 2466, 61, 192, "Input"],
Cell[21817, 569, 481, 7, 31, "Output"],
Cell[22301, 578, 497, 8, 31, "Output"],
Cell[22801, 588, 481, 7, 31, "Output"],
Cell[23285, 597, 481, 7, 31, "Output"],
Cell[23769, 606, 481, 7, 31, "Output"],
Cell[24253, 615, 481, 7, 31, "Output"],
Cell[24737, 624, 480, 7, 31, "Output"],
Cell[25220, 633, 481, 7, 31, "Output"],
Cell[25704, 642, 481, 7, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[26222, 654, 5286, 154, 352, "Input"],
Cell[31511, 810, 1439, 25, 379, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[32987, 840, 5955, 173, 392, "Input"],
Cell[38945, 1015, 1437, 28, 479, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)