(* 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[     40785,       1106]
NotebookOptionsPosition[     39213,       1057]
NotebookOutlinePosition[     39722,       1076]
CellTagsIndexPosition[     39679,       1073]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{"(*", " ", 
  RowBox[{"CH4", "   ", "15.06", ".2012"}], " ", "*)"}]], "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5783693559648447`*^9, 3.5783693662920628`*^9}, 
   3.578369554833994*^9, {3.5783729062258806`*^9, 3.5783729062258806`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"Clear", "[", 
   RowBox[{
   "z", ",", "sig1", ",", "sig2", ",", "sig4", ",", "k1", ",", "k2", ",", 
    "k4", ",", "nc", ",", "R1", ",", "R2", ",", "R3", ",", "R4", ",", "w", 
    ",", "p", ",", "vee", ",", "vne", ",", "vnn", ",", "\n", "a1", ",", "a2", 
    ",", "a3", ",", "b1", ",", "b2", ",", "b3", ",", "c1", ",", "c2", ",", 
    "c3", ",", "cs", ",", "ss", ",", "d1", ",", "d2", ",", "d3", ",", "d4", 
    ",", "d5", ",", "pi", ",", "i", ",", "j", ",", "d", ",", "t"}], "]"}], 
  ";"}]], "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5783693559648447`*^9, 3.5783693662920628`*^9}, 
   3.578369554833994*^9, 3.578370271764453*^9, {3.5783729062258806`*^9, 
   3.5783729062258806`*^9}}],

Cell[BoxData[
 RowBox[{"(*", 
  RowBox[{
   RowBox[{"z", "=", "6.0"}], ";", " ", 
   RowBox[{"sig1", "=", "0.3"}], ";", " ", 
   RowBox[{"sig2", "=", "0.3"}], ";", " ", 
   RowBox[{"sig4", "=", "0.3"}], ";", " ", 
   RowBox[{"nc", "=", "5"}], ";", " ", 
   RowBox[{"pi", "=", "0.0"}], ";", "\n", 
   RowBox[{"k1", "=", "1.0"}], ";", " ", 
   RowBox[{"k2", "=", "1.2"}], ";", " ", 
   RowBox[{"k4", "=", "1.66"}], ";", " ", 
   RowBox[{"bohr", "=", "0.529177"}], ";"}], "*)"}]], "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5783693559648447`*^9, 3.5783693662920628`*^9}, 
   3.578369554833994*^9, 3.578370271764453*^9, {3.578370376315837*^9, 
   3.578370376315837*^9}, {3.5783729062258806`*^9, 3.5783729062258806`*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"z", "=", "6.0"}], ";"}], " "}], "\n", 
 RowBox[{
  RowBox[{"k1", "=", "1.04605075"}], ";", "  ", 
  RowBox[{"k2", "=", "1.20889711"}], ";", "  ", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"sig1", "=", "0.25273353"}], ";", "  ", 
    RowBox[{"sig2", "=", "0.32960771"}], ";"}], "*)"}], "\n", 
  RowBox[{"sig1", "=", "0.2650207"}], ";", " ", 
  RowBox[{"sig2", "=", "0.30451569"}], ";"}], "\n", 
 RowBox[{
  RowBox[{"nc", "=", "5"}], ";", " ", 
  RowBox[{"pi", "=", "0.0"}], ";", " ", 
  RowBox[{"bohr", "=", "0.529177"}], ";", 
  RowBox[{"sq", "=", 
   RowBox[{"Sqrt", "[", "3", "]"}]}], ";"}]}], "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5783693559648447`*^9, 3.5783693662920628`*^9}, 
   3.578369554833994*^9, 3.578370271764453*^9, {3.5783703674238214`*^9, 
   3.5783703954726706`*^9}, {3.578371324750703*^9, 3.5783713343291197`*^9}, 
   3.5783717590710654`*^9, {3.5783729062258806`*^9, 3.5783729062258806`*^9}, 
   3.5783735444386015`*^9, {3.5783737801394157`*^9, 3.578373791355835*^9}}],

Cell[BoxData[
 RowBox[{"(*", " ", 
  RowBox[{
   RowBox[{"C", " ", "He", " ", "shell"}], ",", " ", 
   RowBox[{"i", "=", 
    RowBox[{"j", " ", "each"}]}]}], " ", "*)"}]], "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5783693559648447`*^9, 3.5783693662920628`*^9}, {
   3.578369554833994*^9, 3.5783695604344034`*^9}, {3.5783729062258806`*^9, 
   3.5783729062258806`*^9}, 3.578373808328665*^9}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"Ekin", " ", "=", " ", 
   RowBox[{"2.25", "*", 
    RowBox[{"k1", "/", 
     RowBox[{"R1", "^", "2"}]}]}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"vee", "=", 
   RowBox[{"3.0", "*", 
    RowBox[{"sig1", "/", "R1"}]}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"vne", "=", 
    RowBox[{
     RowBox[{"-", "3.0"}], "*", 
     RowBox[{"z", "/", "R1"}]}]}], ";"}], "\n", 
  RowBox[{"(*", " ", 
   RowBox[{"bonding", " ", "pairs"}], " ", "*)"}]}], "\n", 
 RowBox[{
  RowBox[{"Ekin", " ", "=", " ", 
   RowBox[{"Ekin", " ", "+", " ", 
    RowBox[{"2.25", "*", 
     RowBox[{"(", 
      RowBox[{"4", "*", 
       RowBox[{"k2", "/", 
        RowBox[{"R2", "^", "2"}]}]}], ")"}]}]}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"vee", "=", 
   RowBox[{"vee", "+", 
    RowBox[{"3.0", "*", 
     RowBox[{"(", 
      RowBox[{"4", "*", 
       RowBox[{"sig2", "/", "R2"}]}], ")"}]}]}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"c1", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "2"}], ",", 
     RowBox[{"-", "2"}], ",", 
     RowBox[{"-", "2"}], ",", 
     RowBox[{"-", "2"}], ",", 
     RowBox[{"-", "2"}]}], "}"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"c2", "=", 
   RowBox[{"{", 
    RowBox[{"z", ",", "1", ",", "1", ",", "1", ",", "1"}], "}"}]}], 
  ";"}], "\n", 
 RowBox[{
  RowBox[{"c3", "=", 
   RowBox[{"{", 
    RowBox[{"R1", ",", "R2", ",", "R2", ",", "R2", ",", "R2"}], "}"}]}], 
  ";"}], "\n", 
 RowBox[{
  RowBox[{"d2", "=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"R1", "+", "R2", "+", "p"}], ")"}], "/", "sq"}]}], ";", "  ", 
  RowBox[{"d4", "=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"R1", "+", "R2"}], ")"}], "/", "sq"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"b1", "=", 
   RowBox[{"{", 
    RowBox[{"0", ",", "d2", ",", 
     RowBox[{"-", "d2"}], ",", "d2", ",", 
     RowBox[{"-", "d2"}]}], "}"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"b2", "=", 
   RowBox[{"{", 
    RowBox[{"0", ",", "d2", ",", 
     RowBox[{"-", "d2"}], ",", 
     RowBox[{"-", "d2"}], ",", "d2"}], "}"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"b3", "=", 
    RowBox[{"{", 
     RowBox[{"0", ",", "d2", ",", "d2", ",", 
      RowBox[{"-", "d2"}], ",", 
      RowBox[{"-", "d2"}]}], "}"}]}], ";"}], "        "}], "\n", 
 RowBox[{
  RowBox[{"a1", "=", 
   RowBox[{"{", 
    RowBox[{"0", ",", "d4", ",", 
     RowBox[{"-", "d4"}], ",", "d4", ",", 
     RowBox[{"-", "d4"}]}], "}"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"a2", "=", 
   RowBox[{"{", 
    RowBox[{"0", ",", "d4", ",", 
     RowBox[{"-", "d4"}], ",", 
     RowBox[{"-", "d4"}], ",", "d4"}], "}"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"a3", "=", 
   RowBox[{"{", 
    RowBox[{"0", ",", "d4", ",", "d4", ",", 
     RowBox[{"-", "d4"}], ",", 
     RowBox[{"-", "d4"}]}], "}"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"vne", "=", 
    RowBox[{"vne", "-", 
     RowBox[{"4", "*", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"3", "-", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"p", "/", "R2"}], ")"}], "^", "2"}]}], ")"}], "/", 
       "R2"}]}]}]}], ";"}], "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"vne", "=", 
     RowBox[{"vne", "-", 
      RowBox[{"4", "*", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"3", "-", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{"p", "-", "1"}], ")"}], "*", 
             RowBox[{"(", 
              RowBox[{"1", "+", 
               RowBox[{"R1", "/", "R2"}]}], ")"}]}], ")"}], "^", "2"}]}], 
         ")"}], "/", "R2"}]}]}]}], ";"}], "*)"}]}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"vnn", " ", "=", " ", "0.0"}], ";"}], "\n", 
  RowBox[{"(*", "ww", "*)"}]}], "\n", 
 RowBox[{
  RowBox[{"For", "[", 
   RowBox[{
    RowBox[{"i", " ", "=", " ", "1"}], ",", " ", 
    RowBox[{"i", " ", "<", " ", "nc"}], ",", " ", 
    RowBox[{"i", "++"}], ",", "\n", "  ", 
    RowBox[{"For", "[", 
     RowBox[{
      RowBox[{"j", " ", "=", " ", 
       RowBox[{"i", "+", "1"}]}], ",", " ", 
      RowBox[{"j", " ", "<", " ", 
       RowBox[{"nc", "+", "1"}]}], ",", " ", 
      RowBox[{"j", "++"}], ",", "\n", 
      RowBox[{"vee", " ", "=", 
       RowBox[{"vee", " ", "+", " ", 
        RowBox[{
         RowBox[{"c1", "[", 
          RowBox[{"[", "i", "]"}], "]"}], "*", 
         RowBox[{
          RowBox[{"c1", "[", 
           RowBox[{"[", "j", "]"}], "]"}], "/", 
          RowBox[{"Sqrt", "[", 
           RowBox[{
            RowBox[{
             RowBox[{"(", 
              RowBox[{
               RowBox[{"a1", "[", 
                RowBox[{"[", "i", "]"}], "]"}], "-", 
               RowBox[{"a1", "[", 
                RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}], "+", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{
               RowBox[{"a2", "[", 
                RowBox[{"[", "i", "]"}], "]"}], "-", 
               RowBox[{"a2", "[", 
                RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}], "+", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{
               RowBox[{"a3", "[", 
                RowBox[{"[", "i", "]"}], "]"}], "-", 
               RowBox[{"a3", "[", 
                RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}]}], 
           "]"}]}]}]}]}]}], "]"}]}], "]"}], "\n", "\n", 
  RowBox[{"(*", "nn", "*)"}]}], "\n", 
 RowBox[{
  RowBox[{"For", "[", 
   RowBox[{
    RowBox[{"i", " ", "=", " ", "1"}], ",", " ", 
    RowBox[{"i", " ", "<", " ", "nc"}], ",", " ", 
    RowBox[{"i", "++"}], ",", "\n", "  ", 
    RowBox[{"For", "[", 
     RowBox[{
      RowBox[{"j", " ", "=", " ", 
       RowBox[{"i", "+", "1"}]}], ",", " ", 
      RowBox[{"j", " ", "<", " ", 
       RowBox[{"nc", "+", "1"}]}], ",", " ", 
      RowBox[{"j", "++"}], ",", 
      RowBox[{"vnn", " ", "=", " ", 
       RowBox[{"vnn", " ", "+", " ", 
        RowBox[{
         RowBox[{"c2", "[", 
          RowBox[{"[", "i", "]"}], "]"}], "*", 
         RowBox[{
          RowBox[{"c2", "[", 
           RowBox[{"[", "j", "]"}], "]"}], "/", 
          RowBox[{"Sqrt", "[", 
           RowBox[{
            RowBox[{
             RowBox[{"(", 
              RowBox[{
               RowBox[{"b1", "[", 
                RowBox[{"[", "i", "]"}], "]"}], "-", 
               RowBox[{"b1", "[", 
                RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}], "+", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{
               RowBox[{"b2", "[", 
                RowBox[{"[", "i", "]"}], "]"}], "-", 
               RowBox[{"b2", "[", 
                RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}], "+", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{
               RowBox[{"b3", "[", 
                RowBox[{"[", "i", "]"}], "]"}], "-", 
               RowBox[{"b3", "[", 
                RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}]}], 
           "]"}]}]}]}]}]}], "]"}]}], "]"}], "\n", 
  RowBox[{"(*", "nw", "*)"}]}], "\n", 
 RowBox[{
  RowBox[{"For", "[", 
   RowBox[{
    RowBox[{"i", " ", "=", " ", "1"}], ",", " ", 
    RowBox[{"i", " ", "<", " ", 
     RowBox[{"nc", "+", "1"}]}], ",", " ", 
    RowBox[{"i", "++"}], ",", "\n", "  ", 
    RowBox[{"For", "[", 
     RowBox[{
      RowBox[{"j", " ", "=", " ", "1"}], ",", " ", 
      RowBox[{"j", " ", "<", " ", 
       RowBox[{"nc", "+", "1"}]}], ",", " ", 
      RowBox[{"j", "++"}], ",", 
      RowBox[{
       RowBox[{"d", "=", 
        RowBox[{"Sqrt", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"(", 
            RowBox[{
             RowBox[{"a1", "[", 
              RowBox[{"[", "i", "]"}], "]"}], "-", 
             RowBox[{"b1", "[", 
              RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}], "+", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{
             RowBox[{"a2", "[", 
              RowBox[{"[", "i", "]"}], "]"}], "-", 
             RowBox[{"b2", "[", 
              RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}], "+", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{
             RowBox[{"a3", "[", 
              RowBox[{"[", "i", "]"}], "]"}], "-", 
             RowBox[{"b3", "[", 
              RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}]}], "]"}]}], 
       ";", "\n", "    ", 
       RowBox[{"If", "[", 
        RowBox[{
         RowBox[{"d", "\[NotEqual]", "0"}], ",", 
         RowBox[{"vne", " ", "=", " ", 
          RowBox[{"vne", " ", "+", " ", 
           RowBox[{
            RowBox[{"c1", "[", 
             RowBox[{"[", "i", "]"}], "]"}], "*", 
            RowBox[{
             RowBox[{"c2", "[", 
              RowBox[{"[", "j", "]"}], "]"}], "/", "d"}]}]}]}]}], "]"}]}]}], 
     "]"}]}], "]"}], "      ", "\n"}], "\n", 
 RowBox[{
  RowBox[{"Epot", "=", 
   RowBox[{"vne", "+", "vee", "+", "vnn"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"func", "=", 
   RowBox[{"Ekin", "+", "Epot"}]}], ";"}], "\n", 
 RowBox[{"t", " ", "=", " ", 
  RowBox[{"FindMinimum", "[", 
   RowBox[{"func", ",", 
    RowBox[{"{", 
     RowBox[{"R1", ",", "0.267"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"R2", ",", "1.244"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"p", ",", "0.65"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"Method", " ", "->", " ", "Automatic"}], "}"}], ",", " ", 
    RowBox[{"{", 
     RowBox[{"MaxIterations", " ", "->", " ", "500"}], "}"}]}], 
   "]"}]}]}], "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5783693559648447`*^9, 3.5783693662920628`*^9}, {
   3.578369554833994*^9, 3.5783695604344034`*^9}, 3.5783701337978106`*^9, {
   3.578371289853441*^9, 3.578371302442663*^9}, {3.57837134613834*^9, 
   3.57837135733916*^9}, {3.578371392314421*^9, 3.5783716741757164`*^9}, {
   3.578371712130583*^9, 3.578371724501405*^9}, {3.5783726153945694`*^9, 
   3.5783726153945694`*^9}, 3.5783726553930397`*^9, 3.578372749586005*^9, {
   3.5783728157925215`*^9, 3.5783728573041945`*^9}, {3.578373248116081*^9, 
   3.5783732528896894`*^9}, {3.5783732892065535`*^9, 
   3.5783733600150776`*^9}, {3.578373605637509*^9, 3.5783736465095806`*^9}, 
   3.578373878918789*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"FindMinimum", "::", "lstol"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"The line search decreased the step size to within the \
tolerance specified by AccuracyGoal and PrecisionGoal but was unable to find \
a sufficient decrease in the function. You may need more than \
\[NoBreak]\\!\\(MachinePrecision\\)\[NoBreak] digits of working precision to \
meet these tolerances. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \
ButtonStyle->\\\"Link\\\", ButtonFrame->None, \
ButtonData:>\\\"paclet:ref/message/FindMinimum/lstol\\\", ButtonNote -> \
\\\"FindMinimum::lstol\\\"]\\)\"\>"}]], "Message", "MSG",
 CellChangeTimes->{3.5783738912116103`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", "31.442057252568695`"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"R1", "\[Rule]", "0.2735975221316634`"}], ",", 
     RowBox[{"R2", "\[Rule]", "5.653787030318295`*^10"}], ",", 
     RowBox[{"p", "\[Rule]", "1.2570011015068634`*^10"}]}], "}"}]}], 
  "}"}]], "Output",
 CellChangeTimes->{3.5783738912116103`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"vne", " ", "/.", " ", 
  RowBox[{"t", "[", 
   RowBox[{"[", "2", "]"}], "]"}]}], "\n", 
 RowBox[{"vee", " ", "/.", " ", 
  RowBox[{"t", "[", 
   RowBox[{"[", "2", "]"}], "]"}]}], "\n", 
 RowBox[{"vnn", " ", "/.", " ", 
  RowBox[{"t", "[", 
   RowBox[{"[", "2", "]"}], "]"}]}], "\n", 
 RowBox[{"vir", "=", 
  RowBox[{
   RowBox[{
    RowBox[{"-", "Epot"}], "/", "Ekin"}], " ", "/.", " ", 
   RowBox[{"t", "[", 
    RowBox[{"[", "2", "]"}], "]"}]}]}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"(", 
    RowBox[{"R1", "+", "R2", "+", "p"}], ")"}], "*", "bohr"}], " ", "/.", " ", 
  RowBox[{"t", "[", 
   RowBox[{"[", "2", "]"}], "]"}]}]}], "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5783693559648447`*^9, 3.5783693662920628`*^9}, {
   3.578369554833994*^9, 3.5783695604344034`*^9}, 3.5783696029912786`*^9, {
   3.5783698599861298`*^9, 3.5783699151166267`*^9}, {3.5783699654735155`*^9, 
   3.5783699759255333`*^9}, 3.578370096825746*^9, {3.5783723049696245`*^9, 
   3.578372305827626*^9}, {3.5783729062258806`*^9, 3.5783729062258806`*^9}, {
   3.578373477623684*^9, 3.578373497186118*^9}, {3.5783736524999914`*^9, 
   3.5783736582720013`*^9}, 3.578373849497137*^9}],

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

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

Cell[BoxData["4.0044976182414616`*^-10"], "Output",
 CellChangeTimes->{3.57837389122721*^9}],

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

Cell[BoxData["3.6570301312493195`*^10"], "Output",
 CellChangeTimes->{3.57837389122721*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"plot1", "=", 
   RowBox[{
    RowBox[{"Graphics", "[", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"Circle", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"a1", "[", 
            RowBox[{"[", "1", "]"}], "]"}], ",", 
           RowBox[{"a2", "[", 
            RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", "R1"}], "]"}], ",", 
       "      ", 
       RowBox[{"Circle", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"a1", "[", 
            RowBox[{"[", "2", "]"}], "]"}], ",", 
           RowBox[{"a2", "[", 
            RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", "R1"}], "]"}], ",", 
       RowBox[{"Circle", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"a1", "[", 
            RowBox[{"[", "3", "]"}], "]"}], ",", 
           RowBox[{"a2", "[", 
            RowBox[{"[", "3", "]"}], "]"}]}], "}"}], ",", "R2"}], "]"}], ",", 
       RowBox[{"Circle", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"a1", "[", 
            RowBox[{"[", "4", "]"}], "]"}], ",", 
           RowBox[{"a2", "[", 
            RowBox[{"[", "4", "]"}], "]"}]}], "}"}], ",", "R2"}], "]"}], ",", 
       RowBox[{"Circle", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"a1", "[", 
            RowBox[{"[", "5", "]"}], "]"}], ",", 
           RowBox[{"a2", "[", 
            RowBox[{"[", "5", "]"}], "]"}]}], "}"}], ",", "R2"}], "]"}], ",", 
       RowBox[{"Circle", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"a1", "[", 
            RowBox[{"[", "6", "]"}], "]"}], ",", 
           RowBox[{"a2", "[", 
            RowBox[{"[", "6", "]"}], "]"}]}], "}"}], ",", "R2"}], "]"}], ",", 
       RowBox[{"Circle", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"a1", "[", 
            RowBox[{"[", "7", "]"}], "]"}], ",", 
           RowBox[{"a2", "[", 
            RowBox[{"[", "7", "]"}], "]"}]}], "}"}], ",", "R4"}], "]"}], ",", 
       RowBox[{"Circle", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"a1", "[", 
            RowBox[{"[", "8", "]"}], "]"}], ",", 
           RowBox[{"a2", "[", 
            RowBox[{"[", "8", "]"}], "]"}]}], "}"}], ",", "R4"}], "]"}], ",", 
       RowBox[{"Disk", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"b1", "[", 
            RowBox[{"[", "1", "]"}], "]"}], ",", 
           RowBox[{"b2", "[", 
            RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], 
       ",", " ", 
       RowBox[{"Disk", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"b1", "[", 
            RowBox[{"[", "2", "]"}], "]"}], ",", 
           RowBox[{"b2", "[", 
            RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], 
       ",", " ", 
       RowBox[{"Disk", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"b1", "[", 
            RowBox[{"[", "3", "]"}], "]"}], ",", 
           RowBox[{"b2", "[", 
            RowBox[{"[", "3", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], 
       ",", " ", 
       RowBox[{"Disk", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"b1", "[", 
            RowBox[{"[", "4", "]"}], "]"}], ",", 
           RowBox[{"b2", "[", 
            RowBox[{"[", "4", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], 
       ",", 
       RowBox[{"Disk", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"b1", "[", 
            RowBox[{"[", "5", "]"}], "]"}], ",", 
           RowBox[{"b2", "[", 
            RowBox[{"[", "5", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], 
       ",", 
       RowBox[{"Disk", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"b1", "[", 
            RowBox[{"[", "6", "]"}], "]"}], ",", 
           RowBox[{"b2", "[", 
            RowBox[{"[", "6", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}]}], 
      "}"}], " ", "]"}], " ", "/.", " ", 
    RowBox[{"t", "[", 
     RowBox[{"[", "2", "]"}], "]"}]}]}], "\n"}], "\n", 
 RowBox[{
  RowBox[{"Show", "[", 
   RowBox[{"plot1", ",", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"AspectRatio", " ", "\[Rule]", " ", "Automatic"}], ",", 
      RowBox[{"Axes", " ", "->", " ", "True"}], ",", 
      RowBox[{"GridLines", " ", "->", " ", "Automatic"}], ",", " ", 
      RowBox[{"PlotRange", " ", "\[Rule]", " ", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"-", "4"}], ",", "4"}], "}"}], ",", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"-", "3"}], ",", "3"}], "}"}]}], "}"}]}], ",", " ", 
      RowBox[{"Frame", " ", "->", " ", "True"}]}], "}"}]}], "]"}], 
  "\n"}], "\n", 
 RowBox[{
  RowBox[{"plot2", "=", 
   RowBox[{
    RowBox[{"Graphics", "[", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"Circle", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"a1", "[", 
            RowBox[{"[", "1", "]"}], "]"}], ",", 
           RowBox[{"a2", "[", 
            RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", "R1"}], "]"}], ",", 
       "      ", 
       RowBox[{"Circle", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"a1", "[", 
            RowBox[{"[", "2", "]"}], "]"}], ",", 
           RowBox[{"a3", "[", 
            RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", "R1"}], "]"}], ",", 
       RowBox[{"Circle", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"a1", "[", 
            RowBox[{"[", "3", "]"}], "]"}], ",", 
           RowBox[{"a3", "[", 
            RowBox[{"[", "3", "]"}], "]"}]}], "}"}], ",", "R2"}], "]"}], ",", 
       RowBox[{"Circle", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"a1", "[", 
            RowBox[{"[", "4", "]"}], "]"}], ",", 
           RowBox[{"a3", "[", 
            RowBox[{"[", "4", "]"}], "]"}]}], "}"}], ",", "R2"}], "]"}], ",", 
       RowBox[{"Circle", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"a1", "[", 
            RowBox[{"[", "5", "]"}], "]"}], ",", 
           RowBox[{"a3", "[", 
            RowBox[{"[", "5", "]"}], "]"}]}], "}"}], ",", "R2"}], "]"}], ",", 
       RowBox[{"Circle", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"a1", "[", 
            RowBox[{"[", "6", "]"}], "]"}], ",", 
           RowBox[{"a3", "[", 
            RowBox[{"[", "6", "]"}], "]"}]}], "}"}], ",", "R2"}], "]"}], ",", 
       RowBox[{"Circle", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"a1", "[", 
            RowBox[{"[", "7", "]"}], "]"}], ",", 
           RowBox[{"a3", "[", 
            RowBox[{"[", "7", "]"}], "]"}]}], "}"}], ",", "R4"}], "]"}], ",", 
       RowBox[{"Circle", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"a1", "[", 
            RowBox[{"[", "8", "]"}], "]"}], ",", 
           RowBox[{"a3", "[", 
            RowBox[{"[", "8", "]"}], "]"}]}], "}"}], ",", "R4"}], "]"}], ",", 
       RowBox[{"Disk", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"b1", "[", 
            RowBox[{"[", "1", "]"}], "]"}], ",", 
           RowBox[{"b3", "[", 
            RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], 
       ",", " ", 
       RowBox[{"Disk", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"b1", "[", 
            RowBox[{"[", "2", "]"}], "]"}], ",", 
           RowBox[{"b3", "[", 
            RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], 
       ",", " ", 
       RowBox[{"Disk", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"b1", "[", 
            RowBox[{"[", "3", "]"}], "]"}], ",", 
           RowBox[{"b3", "[", 
            RowBox[{"[", "3", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], 
       ",", " ", 
       RowBox[{"Disk", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"b1", "[", 
            RowBox[{"[", "4", "]"}], "]"}], ",", 
           RowBox[{"b3", "[", 
            RowBox[{"[", "4", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], 
       ",", 
       RowBox[{"Disk", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"b1", "[", 
            RowBox[{"[", "5", "]"}], "]"}], ",", 
           RowBox[{"b3", "[", 
            RowBox[{"[", "5", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], 
       ",", 
       RowBox[{"Disk", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"b1", "[", 
            RowBox[{"[", "6", "]"}], "]"}], ",", 
           RowBox[{"b3", "[", 
            RowBox[{"[", "6", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}]}], 
      "}"}], " ", "]"}], " ", "/.", " ", 
    RowBox[{"t", "[", 
     RowBox[{"[", "2", "]"}], "]"}]}]}], "\n"}], "\n", 
 RowBox[{
  RowBox[{"Show", "[", 
   RowBox[{"plot2", ",", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"AspectRatio", " ", "\[Rule]", " ", "Automatic"}], ",", 
      RowBox[{"Axes", " ", "->", " ", "True"}], ",", 
      RowBox[{"GridLines", " ", "->", " ", "Automatic"}], ",", " ", 
      RowBox[{"PlotRange", " ", "\[Rule]", " ", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"-", "4"}], ",", "4"}], "}"}], ",", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"-", "3"}], ",", "3"}], "}"}]}], "}"}]}], ",", " ", 
      RowBox[{"Frame", " ", "->", " ", "True"}]}], "}"}]}], "]"}], 
  "\n"}], "\n"}], "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5783729062414804`*^9, 3.5783729062414804`*^9}}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Part", "::", "partw"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Part \[NoBreak]\\!\\(6\\)\[NoBreak] of \[NoBreak]\\!\\({0, \
\\(R1 + R2\\)\\/\\@3, \\(\\(-\\(\\(\\(R1 + R2\\)\\/\\@3\\)\\)\\)\\), \\(R1 + \
R2\\)\\/\\@3, \\(\\(-\\(\\(\\(R1 + R2\\)\\/\\@3\\)\\)\\)\\)}\\)\[NoBreak] \
does not exist. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \
ButtonStyle->\\\"Link\\\", ButtonFrame->None, \
ButtonData:>\\\"paclet:ref/message/General/partw\\\", ButtonNote -> \
\\\"Part::partw\\\"]\\)\"\>"}]], "Message", "MSG",
 CellChangeTimes->{3.578373891305211*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Part", "::", "partw"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Part \[NoBreak]\\!\\(6\\)\[NoBreak] of \[NoBreak]\\!\\({0, \
\\(R1 + R2\\)\\/\\@3, \\(\\(-\\(\\(\\(R1 + R2\\)\\/\\@3\\)\\)\\)\\), \
\\(\\(-\\(\\(\\(R1 + R2\\)\\/\\@3\\)\\)\\)\\), \\(R1 + R2\\)\\/\\@3}\\)\
\[NoBreak] does not exist. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \
ButtonStyle->\\\"Link\\\", ButtonFrame->None, \
ButtonData:>\\\"paclet:ref/message/General/partw\\\", ButtonNote -> \
\\\"Part::partw\\\"]\\)\"\>"}]], "Message", "MSG",
 CellChangeTimes->{3.5783738913520107`*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Part", "::", "partw"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Part \[NoBreak]\\!\\(7\\)\[NoBreak] of \[NoBreak]\\!\\({0, \
\\(R1 + R2\\)\\/\\@3, \\(\\(-\\(\\(\\(R1 + R2\\)\\/\\@3\\)\\)\\)\\), \\(R1 + \
R2\\)\\/\\@3, \\(\\(-\\(\\(\\(R1 + R2\\)\\/\\@3\\)\\)\\)\\)}\\)\[NoBreak] \
does not exist. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \
ButtonStyle->\\\"Link\\\", ButtonFrame->None, \
ButtonData:>\\\"paclet:ref/message/General/partw\\\", ButtonNote -> \
\\\"Part::partw\\\"]\\)\"\>"}]], "Message", "MSG",
 CellChangeTimes->{3.578373891398811*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"General", "::", "stop"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Further output of \[NoBreak]\\!\\(\\*StyleBox[\\(Part :: \
partw\\), \\\"MessageName\\\"]\\)\[NoBreak] will be suppressed during this \
calculation. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \
ButtonStyle->\\\"Link\\\", ButtonFrame->None, \
ButtonData:>\\\"paclet:ref/message/General/stop\\\", ButtonNote -> \
\\\"General::stop\\\"]\\)\"\>"}]], "Message", "MSG",
 CellChangeTimes->{3.578373891430011*^9}],

Cell[BoxData[
 GraphicsBox[{CircleBox[{0, 0}, 0.2735975221316634], 
   CircleBox[{3.264215463910879*^10, 3.264215463910879*^10}, 
    0.2735975221316634], 
   CircleBox[{-3.264215463910879*^10, -3.264215463910879*^10}, 
    5.653787030318295*^10], 
   CircleBox[{3.264215463910879*^10, -3.264215463910879*^10}, 
    5.653787030318295*^10], 
   CircleBox[{-3.264215463910879*^10, 3.264215463910879*^10}, 
    5.653787030318295*^10], 
   CircleBox[{
    Part[{0, 3.264215463910879*^10, -3.264215463910879*^10, 
       3.264215463910879*^10, -3.264215463910879*^10}, 6], 
     Part[{0, 3.264215463910879*^10, -3.264215463910879*^10, \
-3.264215463910879*^10, 3.264215463910879*^10}, 6]}, 5.653787030318295*^10], 
   CircleBox[{
    Part[{0, 3.264215463910879*^10, -3.264215463910879*^10, 
       3.264215463910879*^10, -3.264215463910879*^10}, 7], 
     Part[{0, 3.264215463910879*^10, -3.264215463910879*^10, \
-3.264215463910879*^10, 3.264215463910879*^10}, 7]}, $CellContext`R4], 
   CircleBox[{
    Part[{0, 3.264215463910879*^10, -3.264215463910879*^10, 
       3.264215463910879*^10, -3.264215463910879*^10}, 8], 
     Part[{0, 3.264215463910879*^10, -3.264215463910879*^10, \
-3.264215463910879*^10, 3.264215463910879*^10}, 8]}, $CellContext`R4], 
   DiskBox[{0, 0}, 0.08], 
   DiskBox[{3.989945388237523*^10, 3.989945388237523*^10}, 0.08], 
   DiskBox[{-3.989945388237523*^10, -3.989945388237523*^10}, 0.08], 
   DiskBox[{3.989945388237523*^10, -3.989945388237523*^10}, 0.08], 
   DiskBox[{-3.989945388237523*^10, 3.989945388237523*^10}, 0.08], 
   DiskBox[{Part[{
      0, 3.989945388237523*^10, -3.989945388237523*^10, 
       3.989945388237523*^10, -3.989945388237523*^10}, 6], 
     Part[{0, 3.989945388237523*^10, -3.989945388237523*^10, \
-3.989945388237523*^10, 3.989945388237523*^10}, 6]}, 0.08]}]], "Output",
 CellChangeTimes->{3.578373891430011*^9}],

Cell[BoxData[
 GraphicsBox[{CircleBox[{0, 0}, 0.2735975221316634], 
   CircleBox[{3.264215463910879*^10, 3.264215463910879*^10}, 
    0.2735975221316634], 
   CircleBox[{-3.264215463910879*^10, -3.264215463910879*^10}, 
    5.653787030318295*^10], 
   CircleBox[{3.264215463910879*^10, -3.264215463910879*^10}, 
    5.653787030318295*^10], 
   CircleBox[{-3.264215463910879*^10, 3.264215463910879*^10}, 
    5.653787030318295*^10], 
   CircleBox[{
    Part[{0, 3.264215463910879*^10, -3.264215463910879*^10, 
       3.264215463910879*^10, -3.264215463910879*^10}, 6], 
     Part[{0, 3.264215463910879*^10, -3.264215463910879*^10, \
-3.264215463910879*^10, 3.264215463910879*^10}, 6]}, 5.653787030318295*^10], 
   CircleBox[{
    Part[{0, 3.264215463910879*^10, -3.264215463910879*^10, 
       3.264215463910879*^10, -3.264215463910879*^10}, 7], 
     Part[{0, 3.264215463910879*^10, -3.264215463910879*^10, \
-3.264215463910879*^10, 3.264215463910879*^10}, 7]}, $CellContext`R4], 
   CircleBox[{
    Part[{0, 3.264215463910879*^10, -3.264215463910879*^10, 
       3.264215463910879*^10, -3.264215463910879*^10}, 8], 
     Part[{0, 3.264215463910879*^10, -3.264215463910879*^10, \
-3.264215463910879*^10, 3.264215463910879*^10}, 8]}, $CellContext`R4], 
   DiskBox[{0, 0}, 0.08], 
   DiskBox[{3.989945388237523*^10, 3.989945388237523*^10}, 0.08], 
   DiskBox[{-3.989945388237523*^10, -3.989945388237523*^10}, 0.08], 
   DiskBox[{3.989945388237523*^10, -3.989945388237523*^10}, 0.08], 
   DiskBox[{-3.989945388237523*^10, 3.989945388237523*^10}, 0.08], 
   DiskBox[{Part[{
      0, 3.989945388237523*^10, -3.989945388237523*^10, 
       3.989945388237523*^10, -3.989945388237523*^10}, 6], 
     Part[{0, 3.989945388237523*^10, -3.989945388237523*^10, \
-3.989945388237523*^10, 3.989945388237523*^10}, 6]}, 0.08]},
  AspectRatio->Automatic,
  Axes->True,
  Frame->True,
  GridLines->Automatic,
  PlotRange->{{-4, 4}, {-3, 3}}]], "Output",
 CellChangeTimes->{3.578373891430011*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Part", "::", "partw"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Part \[NoBreak]\\!\\(6\\)\[NoBreak] of \[NoBreak]\\!\\({0, \
\\(R1 + R2\\)\\/\\@3, \\(\\(-\\(\\(\\(R1 + R2\\)\\/\\@3\\)\\)\\)\\), \\(R1 + \
R2\\)\\/\\@3, \\(\\(-\\(\\(\\(R1 + R2\\)\\/\\@3\\)\\)\\)\\)}\\)\[NoBreak] \
does not exist. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \
ButtonStyle->\\\"Link\\\", ButtonFrame->None, \
ButtonData:>\\\"paclet:ref/message/General/partw\\\", ButtonNote -> \
\\\"Part::partw\\\"]\\)\"\>"}]], "Message", "MSG",
 CellChangeTimes->{3.5783738914924107`*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Part", "::", "partw"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Part \[NoBreak]\\!\\(6\\)\[NoBreak] of \[NoBreak]\\!\\({0, \
\\(R1 + R2\\)\\/\\@3, \\(R1 + R2\\)\\/\\@3, \\(\\(-\\(\\(\\(R1 + R2\\)\\/\\@3\
\\)\\)\\)\\), \\(\\(-\\(\\(\\(R1 + R2\\)\\/\\@3\\)\\)\\)\\)}\\)\[NoBreak] \
does not exist. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \
ButtonStyle->\\\"Link\\\", ButtonFrame->None, \
ButtonData:>\\\"paclet:ref/message/General/partw\\\", ButtonNote -> \
\\\"Part::partw\\\"]\\)\"\>"}]], "Message", "MSG",
 CellChangeTimes->{3.5783738915392113`*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"Part", "::", "partw"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Part \[NoBreak]\\!\\(7\\)\[NoBreak] of \[NoBreak]\\!\\({0, \
\\(R1 + R2\\)\\/\\@3, \\(\\(-\\(\\(\\(R1 + R2\\)\\/\\@3\\)\\)\\)\\), \\(R1 + \
R2\\)\\/\\@3, \\(\\(-\\(\\(\\(R1 + R2\\)\\/\\@3\\)\\)\\)\\)}\\)\[NoBreak] \
does not exist. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \
ButtonStyle->\\\"Link\\\", ButtonFrame->None, \
ButtonData:>\\\"paclet:ref/message/General/partw\\\", ButtonNote -> \
\\\"Part::partw\\\"]\\)\"\>"}]], "Message", "MSG",
 CellChangeTimes->{3.578373891586011*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"General", "::", "stop"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"Further output of \[NoBreak]\\!\\(\\*StyleBox[\\(Part :: \
partw\\), \\\"MessageName\\\"]\\)\[NoBreak] will be suppressed during this \
calculation. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \
ButtonStyle->\\\"Link\\\", ButtonFrame->None, \
ButtonData:>\\\"paclet:ref/message/General/stop\\\", ButtonNote -> \
\\\"General::stop\\\"]\\)\"\>"}]], "Message", "MSG",
 CellChangeTimes->{3.5783738916172113`*^9}],

Cell[BoxData[
 GraphicsBox[{CircleBox[{0, 0}, 0.2735975221316634], 
   CircleBox[{3.264215463910879*^10, 3.264215463910879*^10}, 
    0.2735975221316634], 
   CircleBox[{-3.264215463910879*^10, 3.264215463910879*^10}, 
    5.653787030318295*^10], 
   CircleBox[{3.264215463910879*^10, -3.264215463910879*^10}, 
    5.653787030318295*^10], 
   CircleBox[{-3.264215463910879*^10, -3.264215463910879*^10}, 
    5.653787030318295*^10], 
   CircleBox[{
    Part[{0, 3.264215463910879*^10, -3.264215463910879*^10, 
       3.264215463910879*^10, -3.264215463910879*^10}, 6], 
     Part[{0, 3.264215463910879*^10, 
       3.264215463910879*^10, -3.264215463910879*^10, -3.264215463910879*^10},
       6]}, 5.653787030318295*^10], 
   CircleBox[{
    Part[{0, 3.264215463910879*^10, -3.264215463910879*^10, 
       3.264215463910879*^10, -3.264215463910879*^10}, 7], 
     Part[{0, 3.264215463910879*^10, 
       3.264215463910879*^10, -3.264215463910879*^10, -3.264215463910879*^10},
       7]}, $CellContext`R4], 
   CircleBox[{
    Part[{0, 3.264215463910879*^10, -3.264215463910879*^10, 
       3.264215463910879*^10, -3.264215463910879*^10}, 8], 
     Part[{0, 3.264215463910879*^10, 
       3.264215463910879*^10, -3.264215463910879*^10, -3.264215463910879*^10},
       8]}, $CellContext`R4], DiskBox[{0, 0}, 0.08], 
   DiskBox[{3.989945388237523*^10, 3.989945388237523*^10}, 0.08], 
   DiskBox[{-3.989945388237523*^10, 3.989945388237523*^10}, 0.08], 
   DiskBox[{3.989945388237523*^10, -3.989945388237523*^10}, 0.08], 
   DiskBox[{-3.989945388237523*^10, -3.989945388237523*^10}, 0.08], 
   DiskBox[{Part[{
      0, 3.989945388237523*^10, -3.989945388237523*^10, 
       3.989945388237523*^10, -3.989945388237523*^10}, 6], 
     Part[{0, 3.989945388237523*^10, 
       3.989945388237523*^10, -3.989945388237523*^10, -3.989945388237523*^10},
       6]}, 0.08]}]], "Output",
 CellChangeTimes->{3.5783738916172113`*^9}],

Cell[BoxData[
 GraphicsBox[{CircleBox[{0, 0}, 0.2735975221316634], 
   CircleBox[{3.264215463910879*^10, 3.264215463910879*^10}, 
    0.2735975221316634], 
   CircleBox[{-3.264215463910879*^10, 3.264215463910879*^10}, 
    5.653787030318295*^10], 
   CircleBox[{3.264215463910879*^10, -3.264215463910879*^10}, 
    5.653787030318295*^10], 
   CircleBox[{-3.264215463910879*^10, -3.264215463910879*^10}, 
    5.653787030318295*^10], 
   CircleBox[{
    Part[{0, 3.264215463910879*^10, -3.264215463910879*^10, 
       3.264215463910879*^10, -3.264215463910879*^10}, 6], 
     Part[{0, 3.264215463910879*^10, 
       3.264215463910879*^10, -3.264215463910879*^10, -3.264215463910879*^10},
       6]}, 5.653787030318295*^10], 
   CircleBox[{
    Part[{0, 3.264215463910879*^10, -3.264215463910879*^10, 
       3.264215463910879*^10, -3.264215463910879*^10}, 7], 
     Part[{0, 3.264215463910879*^10, 
       3.264215463910879*^10, -3.264215463910879*^10, -3.264215463910879*^10},
       7]}, $CellContext`R4], 
   CircleBox[{
    Part[{0, 3.264215463910879*^10, -3.264215463910879*^10, 
       3.264215463910879*^10, -3.264215463910879*^10}, 8], 
     Part[{0, 3.264215463910879*^10, 
       3.264215463910879*^10, -3.264215463910879*^10, -3.264215463910879*^10},
       8]}, $CellContext`R4], DiskBox[{0, 0}, 0.08], 
   DiskBox[{3.989945388237523*^10, 3.989945388237523*^10}, 0.08], 
   DiskBox[{-3.989945388237523*^10, 3.989945388237523*^10}, 0.08], 
   DiskBox[{3.989945388237523*^10, -3.989945388237523*^10}, 0.08], 
   DiskBox[{-3.989945388237523*^10, -3.989945388237523*^10}, 0.08], 
   DiskBox[{Part[{
      0, 3.989945388237523*^10, -3.989945388237523*^10, 
       3.989945388237523*^10, -3.989945388237523*^10}, 6], 
     Part[{0, 3.989945388237523*^10, 
       3.989945388237523*^10, -3.989945388237523*^10, -3.989945388237523*^10},
       6]}, 0.08]},
  AspectRatio->Automatic,
  Axes->True,
  Frame->True,
  GridLines->Automatic,
  PlotRange->{{-4, 4}, {-3, 3}}]], "Output",
 CellChangeTimes->{3.5783738916172113`*^9}]
}, Open  ]]
},
WindowToolbars->"EditBar",
WindowSize->{890, 964},
WindowMargins->{{3, Automatic}, {Automatic, 6}},
PrintingCopies->1,
PrintingPageRange->{Automatic, Automatic},
PrivateNotebookOptions->{"VersionedStylesheet"->{"Default.nb"[8.] -> False}},
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, 272, 5, 31, "Input"],
Cell[832, 27, 735, 14, 52, "Input"],
Cell[1570, 43, 729, 16, 72, "Input"],
Cell[2302, 61, 1051, 24, 112, "Input"],
Cell[3356, 87, 403, 9, 31, "Input"],
Cell[CellGroupData[{
Cell[3784, 100, 10130, 304, 852, "Input"],
Cell[13917, 406, 697, 13, 56, "Message"],
Cell[14617, 421, 375, 10, 33, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[15029, 436, 1197, 29, 112, "Input"],
Cell[16229, 467, 103, 2, 31, "Output"],
Cell[16335, 471, 87, 1, 31, "Output"],
Cell[16425, 474, 92, 1, 31, "Output"],
Cell[16520, 477, 87, 1, 31, "Output"],
Cell[16610, 480, 91, 1, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[16738, 486, 9931, 294, 472, "Input"],
Cell[26672, 782, 608, 12, 42, "Message"],
Cell[27283, 796, 610, 12, 42, "Message"],
Cell[27896, 810, 608, 12, 42, "Message"],
Cell[28507, 824, 528, 11, 21, "Message"],
Cell[29038, 837, 1864, 35, 374, "Output"],
Cell[30905, 874, 1976, 40, 293, "Output"],
Cell[32884, 916, 610, 12, 42, "Message"],
Cell[33497, 930, 610, 12, 42, "Message"],
Cell[34110, 944, 608, 12, 42, "Message"],
Cell[34721, 958, 530, 11, 21, "Message"],
Cell[35254, 971, 1914, 38, 374, "Output"],
Cell[37171, 1011, 2026, 43, 293, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)