(* 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[     61929,       1488]
NotebookOptionsPosition[     60256,       1435]
NotebookOutlinePosition[     60624,       1451]
CellTagsIndexPosition[     60581,       1448]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{"(*", " ", 
  RowBox[{"Ethene", ",", " ", 
   RowBox[{"20.01", ".2013"}]}], " ", "*)"}]], "Input",
 CellChangeTimes->{{3.5675987354944196`*^9, 3.567599332289068*^9}, {
   3.567599364347124*^9, 3.567599388854767*^9}, {3.567599446824469*^9, 
   3.567599671277663*^9}, {3.5675997743158445`*^9, 3.567599879038828*^9}, {
   3.567599920722101*^9, 3.567600053883935*^9}, {3.5676001118380375`*^9, 
   3.5676001705877404`*^9}, {3.5676002126766143`*^9, 
   3.5676003454796476`*^9}, {3.5676004083477583`*^9, 3.567600442527418*^9}, {
   3.5676004745698743`*^9, 3.567600484132691*^9}, {3.5676005145995445`*^9, 
   3.5676005320091753`*^9}, {3.5676005791056576`*^9, 
   3.5676012559440465`*^9}, {3.5676013203097596`*^9, 3.5676013219321623`*^9}, 
   3.5676013524770164`*^9, {3.5676850872857094`*^9, 3.5676850878785105`*^9}, 
   3.567685750380474*^9}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"Clear", "[", 
   RowBox[{
   "Cr", ",", "CC", ",", "CH", ",", "Cpi", ",", "nc", ",", "nn", ",", "xn", 
    ",", "yn", ",", "zn", ",", "xc", ",", "yc", ",", "zc", ",", "cx", ",", 
    "chx", ",", "chy", ",", "\[IndentingNewLine]", "cchx", ",", "cchy", ",", 
    "Ekin", ",", "vne", ",", "vee", ",", "vnn", ",", "Epot", ",", "w", ",", 
    "csw", ",", "ssw", ",", "p", ",", "sq", ",", "sn", ",", "sr", ",", 
    "\[IndentingNewLine]", "k1", ",", "k2", ",", "k3", ",", "k4", ",", "s1", 
    ",", "s2", ",", "s3", ",", "s4", ",", "Z"}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"nc", "=", "9"}], ";", "  ", 
  RowBox[{"nn", "=", "6"}], ";", " ", 
  RowBox[{"Z", "=", "6.0"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"bohr", "=", "0.529177"}], ";", 
  RowBox[{"rad", "=", "57.29578"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"k1", "=", "1.08"}], ";", " ", 
  RowBox[{"k2", "=", "1.05"}], ";", " ", 
  RowBox[{"k3", "=", "1.25"}], ";", " ", 
  RowBox[{"k4", "=", "1.1"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"s1", "=", "0.3"}], ";", 
   RowBox[{"s2", "=", "0.3"}], ";", " ", 
   RowBox[{"s3", "=", "0.3"}], ";"}], "  "}]}], "Input",
 CellChangeTimes->{{3.5675987354944196`*^9, 3.567599332289068*^9}, {
   3.567599364347124*^9, 3.567599388854767*^9}, {3.567599446824469*^9, 
   3.567599671277663*^9}, {3.5675997743158445`*^9, 3.567599879038828*^9}, {
   3.567599920722101*^9, 3.567600053883935*^9}, {3.5676001118380375`*^9, 
   3.5676001705877404`*^9}, {3.5676002126766143`*^9, 
   3.5676003454796476`*^9}, {3.5676004083477583`*^9, 3.567600442527418*^9}, {
   3.5676004745698743`*^9, 3.567600484132691*^9}, {3.5676005145995445`*^9, 
   3.5676005320091753`*^9}, {3.5676005791056576`*^9, 
   3.5676012559440465`*^9}, {3.5676013203097596`*^9, 3.5676013219321623`*^9}, 
   3.5676013524770164`*^9, {3.5676015878502293`*^9, 3.56760173054368*^9}, {
   3.5676125832717414`*^9, 3.5676126027561755`*^9}, {3.5676126533938646`*^9, 
   3.5676126797423105`*^9}, {3.5676127131263695`*^9, 
   3.5676127364484105`*^9}, {3.567613357173501*^9, 3.567613363585112*^9}, {
   3.5676168815728908`*^9, 3.567616936672188*^9}, {3.5676611367108335`*^9, 
   3.567661137740435*^9}, {3.5676854112046785`*^9, 3.567685413638283*^9}, 
   3.5676857759021187`*^9, {3.5676860335109715`*^9, 3.5676860647422266`*^9}, {
   3.567687447325855*^9, 3.5676874507266607`*^9}, {3.567698526156251*^9, 
   3.567698526967453*^9}}],

Cell[BoxData[
 RowBox[{"(*", " ", 
  RowBox[{
   RowBox[{"6", " ", 
    RowBox[{"Variable", ":", " ", "Cr"}]}], ",", "CC", ",", "CH", ",", "Cpi", 
   ",", "p", ",", "w"}], " ", "*)"}]], "Input",
 CellChangeTimes->{{3.5675987354944196`*^9, 3.567599332289068*^9}, {
   3.567599364347124*^9, 3.567599388854767*^9}, {3.567599446824469*^9, 
   3.567599671277663*^9}, {3.5675997743158445`*^9, 3.567599879038828*^9}, {
   3.567599920722101*^9, 3.567600053883935*^9}, {3.5676001118380375`*^9, 
   3.5676001705877404`*^9}, {3.5676002126766143`*^9, 
   3.5676003454796476`*^9}, {3.5676004083477583`*^9, 3.567600442527418*^9}, {
   3.5676004745698743`*^9, 3.567600484132691*^9}, {3.5676005145995445`*^9, 
   3.5676005320091753`*^9}, {3.5676005791056576`*^9, 
   3.5676012559440465`*^9}, {3.5676013203097596`*^9, 
   3.5676013219321623`*^9}, {3.5676013524770164`*^9, 3.5676013644734373`*^9}, 
   3.56761275877205*^9, 3.5676857814557285`*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"Cr", "=", "0.26"}], ";", " ", 
   RowBox[{"CC", "=", "1.0"}], ";", " ", 
   RowBox[{"CH", "=", "1.22"}], ";", " ", 
   RowBox[{"Cpi", "=", "1.4"}], ";", " ", 
   RowBox[{"w", "=", "1.02"}], ";", " ", 
   RowBox[{"p", "=", "0.56"}], ";"}], "\[IndentingNewLine]"}]], "Input",
 CellChangeTimes->{{3.5676156777655764`*^9, 3.567615748683301*^9}, {
   3.567615781240558*^9, 3.5676157870125685`*^9}, {3.5676162471041765`*^9, 
   3.567616264030206*^9}, 3.567616341141142*^9, {3.567616615342823*^9, 
   3.567616620444032*^9}, {3.567616774650303*^9, 3.567616775617505*^9}, {
   3.5676171308145285`*^9, 3.5676171359469376`*^9}, 3.5676857870561385`*^9, {
   3.567700521524556*^9, 3.567700527187366*^9}, {3.5677006006790953`*^9, 
   3.567700601334296*^9}, {3.5677011684888926`*^9, 3.5677011693000937`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"C", " ", "He", " ", "shell"}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"Ekin", "=", 
     RowBox[{"2.", "*", "2.25", "*", 
      RowBox[{"k1", "/", 
       RowBox[{"Cr", "^", "2"}]}]}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"vee", "=", 
     RowBox[{"2.", "*", "3.0", "*", 
      RowBox[{"s1", "/", "Cr"}]}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"vne", "=", 
     RowBox[{
      RowBox[{"-", "2."}], "*", "3.0", "*", 
      RowBox[{"Z", "/", "Cr"}]}]}], ";"}], "\[IndentingNewLine]"}]}]], "Input",
 CellChangeTimes->{{3.5675987354944196`*^9, 3.567599332289068*^9}, {
  3.567599364347124*^9, 3.567599388854767*^9}, {3.567599446824469*^9, 
  3.567599671277663*^9}, {3.5675997743158445`*^9, 3.567599879038828*^9}, {
  3.567599920722101*^9, 3.567600053883935*^9}, {3.5676001118380375`*^9, 
  3.5676001705877404`*^9}, {3.5676002126766143`*^9, 3.5676003454796476`*^9}, {
  3.5676004083477583`*^9, 3.567600442527418*^9}, {3.5676004745698743`*^9, 
  3.567600484132691*^9}, {3.5676005145995445`*^9, 3.5676005320091753`*^9}, {
  3.5676005791056576`*^9, 3.5676012559440465`*^9}, {3.5676013203097596`*^9, 
  3.5676013219321623`*^9}, {3.5676013524770164`*^9, 3.5676013644734373`*^9}, {
  3.5676917901504855`*^9, 3.567691799978503*^9}}],

Cell[BoxData[
 RowBox[{"(*", 
  RowBox[{"bonding", " ", "pairs"}], " ", "*)"}]], "Input",
 CellChangeTimes->{{3.5675987354944196`*^9, 3.567599332289068*^9}, {
  3.567599364347124*^9, 3.567599388854767*^9}, {3.567599446824469*^9, 
  3.567599671277663*^9}, {3.5675997743158445`*^9, 3.567599879038828*^9}, {
  3.567599920722101*^9, 3.567600053883935*^9}, {3.5676001118380375`*^9, 
  3.5676001705877404`*^9}, {3.5676002126766143`*^9, 3.5676003454796476`*^9}, {
  3.5676004083477583`*^9, 3.567600442527418*^9}, {3.5676004745698743`*^9, 
  3.567600484132691*^9}, {3.5676005145995445`*^9, 3.5676005320091753`*^9}, {
  3.5676005791056576`*^9, 3.5676012559440465`*^9}, {3.5676013203097596`*^9, 
  3.5676013219321623`*^9}, {3.5676013524770164`*^9, 3.567601370448248*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"Ekin", "=", 
   RowBox[{"Ekin", "+", " ", 
    RowBox[{"2.25", "*", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"k2", "/", 
        RowBox[{"CC", "^", "2"}]}], "+", 
       RowBox[{"4.0", "*", 
        RowBox[{"k3", "/", 
         RowBox[{"CH", "^", "2"}]}]}], "+", 
       RowBox[{"k4", "/", 
        RowBox[{"Cpi", "^", "2"}]}]}], ")"}]}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"vee", "=", 
   RowBox[{"vee", "+", 
    RowBox[{"3.0", "*", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"s2", "/", "CC"}], "+", 
       RowBox[{"4.0", "*", 
        RowBox[{"s3", "/", "CH"}]}]}], ")"}]}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"vne", "=", 
   RowBox[{"vne", "-", 
    RowBox[{"4.0", "*", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"3.0", "-", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"p", "/", "CH"}], ")"}], "^", "2"}]}], ")"}], "/", 
      "CH"}]}]}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.5675987354944196`*^9, 3.567599332289068*^9}, {
   3.567599364347124*^9, 3.567599388854767*^9}, {3.567599446824469*^9, 
   3.567599671277663*^9}, {3.5675997743158445`*^9, 3.567599879038828*^9}, {
   3.567599920722101*^9, 3.567600053883935*^9}, {3.5676001118380375`*^9, 
   3.5676001705877404`*^9}, {3.5676002126766143`*^9, 
   3.5676003454796476`*^9}, {3.5676004083477583`*^9, 3.567600442527418*^9}, {
   3.5676004745698743`*^9, 3.567600484132691*^9}, {3.5676005145995445`*^9, 
   3.5676005320091753`*^9}, {3.5676005791056576`*^9, 
   3.5676012559440465`*^9}, {3.5676013203097596`*^9, 
   3.5676013219321623`*^9}, {3.5676013524770164`*^9, 3.567601370448248*^9}, {
   3.5676046193892126`*^9, 3.5676046713841043`*^9}, {3.567611708656205*^9, 
   3.567611709685807*^9}, {3.567611936042205*^9, 3.5676119371030064`*^9}, {
   3.567612182787838*^9, 3.5676121931306562`*^9}, 3.5676857932649493`*^9, {
   3.5676918030049086`*^9, 3.5676918121465244`*^9}, {3.5676919908760386`*^9, 
   3.567691999190853*^9}, {3.567699705705523*^9, 3.5676997106819315`*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"cn", "=", 
   RowBox[{"{", 
    RowBox[{"Z", ",", "Z", ",", "1.", ",", "1.", ",", "1.", ",", "1."}], 
    "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"cc", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "2."}], ",", 
     RowBox[{"-", "2."}], ",", 
     RowBox[{"-", "2."}], ",", 
     RowBox[{"-", "2."}], ",", 
     RowBox[{"-", "2."}], ",", 
     RowBox[{"-", "2."}], ",", 
     RowBox[{"-", "2."}], ",", 
     RowBox[{"-", "1."}], ",", 
     RowBox[{"-", "1."}]}], "}"}]}], ";"}], "\[IndentingNewLine]"}], "Input",
 CellChangeTimes->{{3.5675987354944196`*^9, 3.567599332289068*^9}, {
  3.567599364347124*^9, 3.567599388854767*^9}, {3.567599446824469*^9, 
  3.567599671277663*^9}, {3.5675997743158445`*^9, 3.567599879038828*^9}, {
  3.567599920722101*^9, 3.567600053883935*^9}, {3.5676001118380375`*^9, 
  3.5676001705877404`*^9}, {3.5676002126766143`*^9, 3.5676003454796476`*^9}, {
  3.5676004083477583`*^9, 3.567600442527418*^9}, {3.5676004745698743`*^9, 
  3.567600484132691*^9}, {3.5676005145995445`*^9, 3.5676005320091753`*^9}, {
  3.5676005791056576`*^9, 3.5676012559440465`*^9}, {3.5676013203097596`*^9, 
  3.5676013219321623`*^9}, {3.5676013524770164`*^9, 3.567601375861457*^9}, {
  3.5676918198373375`*^9, 3.5676918343609633`*^9}, {3.5676984128844523`*^9, 
  3.5676984356604924`*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"csw", "=", 
   RowBox[{"Cos", "[", "w", "]"}]}], ";", "    ", 
  RowBox[{"ssw", "=", 
   RowBox[{"Sin", "[", "w", "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"cx", "=", 
   RowBox[{"Cr", "+", "CC"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"chx", "=", "   ", 
   RowBox[{"cx", "+", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"Cr", "+", "CH", "+", "p"}], ")"}], "*", "csw"}]}]}], ";", 
  "          ", 
  RowBox[{"chy", "=", "  ", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"Cr", "+", "CH", "+", "p"}], ")"}], "*", "ssw"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"cchx", "=", 
   RowBox[{"cx", "+", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"Cr", "+", "CH"}], ")"}], "*", "csw"}]}]}], ";", 
  "                 ", 
  RowBox[{"cchy", "=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"Cr", "+", "CH"}], ")"}], "*", "ssw"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"ccpi", "=", 
   RowBox[{
    RowBox[{"CC", "/", "3."}], "+", "Cpi"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.5675987354944196`*^9, 3.567599332289068*^9}, {
   3.567599364347124*^9, 3.567599388854767*^9}, {3.567599446824469*^9, 
   3.567599671277663*^9}, {3.5675997743158445`*^9, 3.567599879038828*^9}, {
   3.567599920722101*^9, 3.567600053883935*^9}, {3.5676001118380375`*^9, 
   3.5676001705877404`*^9}, {3.5676002126766143`*^9, 
   3.5676003454796476`*^9}, {3.5676004083477583`*^9, 3.567600442527418*^9}, {
   3.5676004745698743`*^9, 3.567600484132691*^9}, {3.5676005145995445`*^9, 
   3.5676005320091753`*^9}, {3.5676005791056576`*^9, 
   3.5676012559440465`*^9}, {3.5676013203097596`*^9, 
   3.5676013219321623`*^9}, {3.5676013524770164`*^9, 3.567601384847073*^9}, {
   3.567601484843249*^9, 3.5676015176657066`*^9}, 3.5676019519548693`*^9, {
   3.5676050017146845`*^9, 3.567605007845495*^9}, {3.5676052437491093`*^9, 
   3.5676052486163177`*^9}, {3.567660817518673*^9, 3.5676608349127035`*^9}, 
   3.5676855136344585`*^9, {3.567685802016565*^9, 3.567685806400172*^9}, {
   3.567698556061504*^9, 3.5676985705851293`*^9}, {3.567698658257283*^9, 
   3.5676986617984896`*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"xn", "=", 
   RowBox[{"{", 
    RowBox[{"cx", ",", 
     RowBox[{"-", "cx"}], ",", "chx", ",", "   ", "chx", ",", 
     RowBox[{"-", "chx"}], ",", 
     RowBox[{"-", "chx"}]}], "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"yn", "=", 
   RowBox[{"{", 
    RowBox[{"0", ",", "       ", "0", ",", "chy", ",", 
     RowBox[{"-", "chy"}], ",", "   ", "chy", ",", 
     RowBox[{"-", "chy"}]}], "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"zn", "=", 
   RowBox[{"{", 
    RowBox[{"0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0"}], "}"}]}], 
  ";"}], "\[IndentingNewLine]"}], "Input",
 CellChangeTimes->{{3.5675987354944196`*^9, 3.567599332289068*^9}, {
  3.567599364347124*^9, 3.567599388854767*^9}, {3.567599446824469*^9, 
  3.567599671277663*^9}, {3.5675997743158445`*^9, 3.567599879038828*^9}, {
  3.567599920722101*^9, 3.567600053883935*^9}, {3.5676001118380375`*^9, 
  3.5676001705877404`*^9}, {3.5676002126766143`*^9, 3.5676003454796476`*^9}, {
  3.5676004083477583`*^9, 3.567600442527418*^9}, {3.5676004745698743`*^9, 
  3.567600484132691*^9}, {3.5676005145995445`*^9, 3.5676005320091753`*^9}, {
  3.5676005791056576`*^9, 3.5676012559440465`*^9}, {3.5676013203097596`*^9, 
  3.5676013219321623`*^9}, {3.5676013524770164`*^9, 3.5676013901354823`*^9}, {
  3.5676048225327697`*^9, 3.567604838756798*^9}, {3.5676114185892954`*^9, 
  3.5676114417709365`*^9}, {3.5676854961000276`*^9, 3.5676855041184416`*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"xc", "=", 
   RowBox[{"{", 
    RowBox[{"cx", ",", 
     RowBox[{"-", "cx"}], ",", "0", ",", "cchx", ",", "   ", "cchx", ",", 
     RowBox[{"-", "cchx"}], ",", 
     RowBox[{"-", "cchx"}], ",", "     ", "0", ",", "          ", "0"}], 
    "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"yc", "=", 
   RowBox[{"{", 
    RowBox[{"0", ",", "       ", "0", ",", "0", ",", "cchy", ",", 
     RowBox[{"-", "cchy"}], ",", "   ", "cchy", ",", 
     RowBox[{"-", "cchy"}], ",", "      ", "0", ",", "         ", "0"}], 
    "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"zc", "=", 
   RowBox[{"{", 
    RowBox[{
    "0", ",", "       ", "0", ",", "0", ",", "       ", "0", ",", 
     "          ", "0", ",", "         ", "0", ",", "          ", "0", ",", 
     "ccpi", ",", 
     RowBox[{"-", "ccpi"}]}], "}"}]}], ";"}], "\[IndentingNewLine]"}], "Input",
 CellChangeTimes->{{3.5675987354944196`*^9, 3.567599332289068*^9}, {
  3.567599364347124*^9, 3.567599388854767*^9}, {3.567599446824469*^9, 
  3.567599671277663*^9}, {3.5675997743158445`*^9, 3.567599879038828*^9}, {
  3.567599920722101*^9, 3.567600053883935*^9}, {3.5676001118380375`*^9, 
  3.5676001705877404`*^9}, {3.5676002126766143`*^9, 3.5676003454796476`*^9}, {
  3.5676004083477583`*^9, 3.567600442527418*^9}, {3.5676004745698743`*^9, 
  3.567600484132691*^9}, {3.5676005145995445`*^9, 3.5676005320091753`*^9}, {
  3.5676005791056576`*^9, 3.5676012559440465`*^9}, {3.5676013203097596`*^9, 
  3.5676013219321623`*^9}, {3.5676013524770164`*^9, 3.5676013952366915`*^9}, {
  3.567604906585717*^9, 3.5676049391429744`*^9}, {3.5676607907646255`*^9, 
  3.567660802589446*^9}, {3.567685532853692*^9, 3.5676856314926653`*^9}, {
  3.5676984574849305`*^9, 3.567698503707812*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", "ww", "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"For", "[", 
    RowBox[{
     RowBox[{"i", "=", "1"}], ",", 
     RowBox[{"i", "<", "nc"}], ",", 
     RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
     RowBox[{"For", "[", 
      RowBox[{
       RowBox[{"j", "=", 
        RowBox[{"i", "+", "1"}]}], ",", 
       RowBox[{"j", "<", 
        RowBox[{"nc", "+", "1"}]}], ",", 
       RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"sq", "=", 
         RowBox[{"Sqrt", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"(", 
             RowBox[{
              RowBox[{"xc", "[", 
               RowBox[{"[", "i", "]"}], "]"}], "-", 
              RowBox[{"xc", "[", 
               RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}], "+", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{
              RowBox[{"yc", "[", 
               RowBox[{"[", "i", "]"}], "]"}], "-", 
              RowBox[{"yc", "[", 
               RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}], "+", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{
              RowBox[{"zc", "[", 
               RowBox[{"[", "i", "]"}], "]"}], "-", 
              RowBox[{"zc", "[", 
               RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}]}], "]"}]}],
         ";", "\[IndentingNewLine]", 
        RowBox[{"If", "[", 
         RowBox[{
          RowBox[{"sq", "\[NotEqual]", "0"}], ",", 
          RowBox[{"vee", "=", 
           RowBox[{"vee", "+", 
            RowBox[{
             RowBox[{"cc", "[", 
              RowBox[{"[", "i", "]"}], "]"}], "*", 
             RowBox[{
              RowBox[{"cc", "[", 
               RowBox[{"[", "j", "]"}], "]"}], "/", "sq"}]}]}]}], ",", 
          RowBox[{"Play", "[", 
           RowBox[{
            RowBox[{"Sin", "[", 
             RowBox[{"2", "Pi", " ", "440", " ", "u"}], "]"}], ",", 
            RowBox[{"{", 
             RowBox[{"u", ",", "0", ",", "0.004"}], "}"}]}], "]"}]}], 
         "]"}]}]}], "]"}]}], "]"}], ";"}]}]], "Input",
 CellChangeTimes->{{3.5675987354944196`*^9, 3.567599332289068*^9}, {
   3.567599364347124*^9, 3.567599388854767*^9}, {3.567599446824469*^9, 
   3.567599671277663*^9}, {3.5675997743158445`*^9, 3.567599879038828*^9}, {
   3.567599920722101*^9, 3.567600053883935*^9}, {3.5676001118380375`*^9, 
   3.5676001705877404`*^9}, {3.5676002126766143`*^9, 
   3.5676003454796476`*^9}, {3.5676004083477583`*^9, 3.567600442527418*^9}, {
   3.5676004745698743`*^9, 3.567600484132691*^9}, {3.5676005145995445`*^9, 
   3.5676005320091753`*^9}, {3.5676005791056576`*^9, 
   3.5676012559440465`*^9}, {3.5676013203097596`*^9, 
   3.5676013219321623`*^9}, {3.5676013524770164`*^9, 
   3.5676014051115084`*^9}, {3.5676177595268326`*^9, 3.567617780103269*^9}, {
   3.5676178274649525`*^9, 3.5676178737346334`*^9}, {3.5676608539135365`*^9, 
   3.5676608549587383`*^9}, {3.567686936378157*^9, 3.5676870257195144`*^9}, {
   3.5676986954009485`*^9, 3.5676986962901497`*^9}, 3.567701041489069*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", "nn", "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"vnn", "=", "0.0"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"For", "[", 
     RowBox[{
      RowBox[{"i", "=", "1"}], ",", 
      RowBox[{"i", "<", "nn"}], ",", 
      RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
      RowBox[{"For", "[", 
       RowBox[{
        RowBox[{"j", "=", 
         RowBox[{"i", "+", "1"}]}], ",", 
        RowBox[{"j", "<", 
         RowBox[{"nn", "+", "1"}]}], ",", 
        RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{"sn", "=", 
          RowBox[{"Sqrt", "[", 
           RowBox[{
            RowBox[{
             RowBox[{"(", 
              RowBox[{
               RowBox[{"xn", "[", 
                RowBox[{"[", "i", "]"}], "]"}], "-", 
               RowBox[{"xn", "[", 
                RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}], "+", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{
               RowBox[{"yn", "[", 
                RowBox[{"[", "i", "]"}], "]"}], "-", 
               RowBox[{"yn", "[", 
                RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}], "+", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{
               RowBox[{"zn", "[", 
                RowBox[{"[", "i", "]"}], "]"}], "-", 
               RowBox[{"zn", "[", 
                RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}]}], 
           "]"}]}], ";", "\[IndentingNewLine]", 
         RowBox[{"If", "[", 
          RowBox[{
           RowBox[{"sn", "\[NotEqual]", "0"}], ",", 
           RowBox[{"vnn", "=", 
            RowBox[{"vnn", "+", 
             RowBox[{
              RowBox[{"cn", "[", 
               RowBox[{"[", "i", "]"}], "]"}], "*", 
              RowBox[{
               RowBox[{"cn", "[", 
                RowBox[{"[", "j", "]"}], "]"}], "/", "sn"}]}]}]}], ",", 
           RowBox[{"Play", "[", 
            RowBox[{
             RowBox[{"Sin", "[", 
              RowBox[{"2", "Pi", " ", "340", " ", "u"}], "]"}], ",", 
             RowBox[{"{", 
              RowBox[{"u", ",", "0", ",", "0.004"}], "}"}]}], "]"}]}], 
          "]"}]}]}], "]"}]}], "]"}], ";"}], 
   "\[IndentingNewLine]"}]}]], "Input",
 CellChangeTimes->{{3.5675987354944196`*^9, 3.567599332289068*^9}, {
   3.567599364347124*^9, 3.567599388854767*^9}, {3.567599446824469*^9, 
   3.567599671277663*^9}, {3.5675997743158445`*^9, 3.567599879038828*^9}, {
   3.567599920722101*^9, 3.567600053883935*^9}, {3.5676001118380375`*^9, 
   3.5676001705877404`*^9}, {3.5676002126766143`*^9, 
   3.5676003454796476`*^9}, {3.5676004083477583`*^9, 3.567600442527418*^9}, {
   3.5676004745698743`*^9, 3.567600484132691*^9}, {3.5676005145995445`*^9, 
   3.5676005320091753`*^9}, {3.5676005791056576`*^9, 
   3.5676012559440465`*^9}, {3.5676013203097596`*^9, 
   3.5676013219321623`*^9}, {3.5676013524770164`*^9, 3.56760141183512*^9}, {
   3.5676636515600505`*^9, 3.5676636621992693`*^9}, {3.5676872211098576`*^9, 
   3.5676872966451902`*^9}, 3.5676928599223647`*^9, 3.567699729963566*^9, 
   3.56769980080329*^9, {3.567700354276662*^9, 3.5677003607506733`*^9}, {
   3.567700402621147*^9, 3.567700403572749*^9}, 3.5677004786400805`*^9, {
   3.5677008744439754`*^9, 3.5677009568901205`*^9}, {3.5677010667611136`*^9, 
   3.5677010809259386`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", "nw", "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"For", "[", 
     RowBox[{
      RowBox[{"i", "=", "1"}], ",", 
      RowBox[{"i", "<", 
       RowBox[{"nc", "+", "1"}]}], ",", 
      RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
      RowBox[{"For", "[", 
       RowBox[{
        RowBox[{"j", "=", "1"}], ",", 
        RowBox[{"j", "<", 
         RowBox[{"nn", "+", "1"}]}], ",", 
        RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{"sr", "=", 
          RowBox[{"Sqrt", "[", 
           RowBox[{
            RowBox[{
             RowBox[{"(", 
              RowBox[{
               RowBox[{"xc", "[", 
                RowBox[{"[", "i", "]"}], "]"}], "-", 
               RowBox[{"xn", "[", 
                RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}], "+", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{
               RowBox[{"yc", "[", 
                RowBox[{"[", "i", "]"}], "]"}], "-", 
               RowBox[{"yn", "[", 
                RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}], "+", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{
               RowBox[{"zc", "[", 
                RowBox[{"[", "i", "]"}], "]"}], "-", 
               RowBox[{"zn", "[", 
                RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}]}], 
           "]"}]}], ";", "\[IndentingNewLine]", 
         RowBox[{"If", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"(", 
             RowBox[{"i", "\[NotEqual]", "j"}], ")"}], " ", "&&", " ", 
            RowBox[{"(", 
             RowBox[{"sr", "\[NotEqual]", "0"}], ")"}]}], ",", 
           RowBox[{"vne", "=", 
            RowBox[{"vne", "+", 
             RowBox[{
              RowBox[{"cc", "[", 
               RowBox[{"[", "i", "]"}], "]"}], "*", 
              RowBox[{
               RowBox[{"cn", "[", 
                RowBox[{"[", "j", "]"}], "]"}], "/", "sr"}]}]}]}], ",", 
           RowBox[{"Play", "[", 
            RowBox[{
             RowBox[{"Sin", "[", 
              RowBox[{"2", "Pi", " ", "540", " ", "u"}], "]"}], ",", 
             RowBox[{"{", 
              RowBox[{"u", ",", "0", ",", "0.004"}], "}"}]}], "]"}]}], 
          "]"}]}]}], "]"}]}], "]"}], ";"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{
      RowBox[{"{", 
       RowBox[{"Method", "\[Rule]", "\"\<Newton\>\""}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"MaxIterations", "\[Rule]", "800"}], "}"}]}], "]"}], 
    "*)"}]}]}]], "Input",
 CellChangeTimes->{{3.5675987354944196`*^9, 3.567599332289068*^9}, {
   3.567599364347124*^9, 3.567599388854767*^9}, {3.567599446824469*^9, 
   3.567599671277663*^9}, {3.5675997743158445`*^9, 3.567599879038828*^9}, {
   3.567599920722101*^9, 3.567600053883935*^9}, {3.5676001118380375`*^9, 
   3.5676001705877404`*^9}, {3.5676002126766143`*^9, 
   3.5676003454796476`*^9}, {3.5676004083477583`*^9, 3.567600442527418*^9}, {
   3.5676004745698743`*^9, 3.567600484132691*^9}, {3.5676005145995445`*^9, 
   3.5676005320091753`*^9}, {3.5676005791056576`*^9, 
   3.5676012559440465`*^9}, {3.5676013203097596`*^9, 
   3.5676013219321623`*^9}, {3.5676013524770164`*^9, 3.567601421335537*^9}, {
   3.5676156065414515`*^9, 3.5676156341635*^9}, 3.567663665537675*^9, {
   3.5676849746847115`*^9, 3.5676849977571526`*^9}, 3.5676858120005827`*^9, {
   3.5676873370648613`*^9, 3.5676874232706127`*^9}, 3.567692864898773*^9, {
   3.567698958729411*^9, 3.5676989761858416`*^9}, {3.567699079504823*^9, 
   3.567699113762483*^9}, {3.5677010555758944`*^9, 3.567701074202327*^9}}],

Cell[BoxData[
 RowBox[{"(*", 
  RowBox[{
   RowBox[{"Cr", "=", "0.276"}], ";", " ", 
   RowBox[{"CC", "=", "1.1"}], ";", " ", 
   RowBox[{"CH", "=", "1.25"}], ";", " ", 
   RowBox[{"Cpi", "=", "1.5"}], ";", " ", 
   RowBox[{"w", "=", "1.02"}], ";", " ", 
   RowBox[{"p", "=", "0.56"}], ";"}], "*)"}]], "Input",
 CellChangeTimes->{{3.5675987354944196`*^9, 3.567599332289068*^9}, {
   3.567599364347124*^9, 3.567599388854767*^9}, {3.567599446824469*^9, 
   3.567599671277663*^9}, {3.5675997743158445`*^9, 3.567599879038828*^9}, {
   3.567599920722101*^9, 3.567600053883935*^9}, {3.5676001118380375`*^9, 
   3.5676001705877404`*^9}, {3.5676002126766143`*^9, 
   3.5676003454796476`*^9}, {3.5676004083477583`*^9, 3.567600442527418*^9}, {
   3.5676004745698743`*^9, 3.567600484132691*^9}, {3.5676005145995445`*^9, 
   3.5676005320091753`*^9}, {3.5676005791056576`*^9, 
   3.5676012559440465`*^9}, {3.5676013203097596`*^9, 
   3.5676013219321623`*^9}, {3.5676013524770164`*^9, 3.567601421335537*^9}, {
   3.5676156065414515`*^9, 3.5676156341635*^9}, 3.567663665537675*^9, {
   3.5676849746847115`*^9, 3.5676849977571526`*^9}, 3.5676858120005827`*^9, {
   3.5676873370648613`*^9, 3.5676874232706127`*^9}, 3.567692864898773*^9, {
   3.567698958729411*^9, 3.5676989761858416`*^9}, {3.567699079504823*^9, 
   3.567699136288923*^9}, {3.5677006772596292`*^9, 3.56770068318764*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"Epot", "=", 
   RowBox[{"vne", "+", "vee", "+", "vnn"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"func", "=", 
    RowBox[{"Ekin", "+", "Epot"}]}], ";"}], "\n", 
  RowBox[{"(*", 
   RowBox[{"t", "=", 
    RowBox[{"FindMinimum", "[", 
     RowBox[{"func", ",", 
      RowBox[{"{", 
       RowBox[{"Cr", ",", "0.276"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"CC", ",", "1.1"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"CH", ",", "1.25"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"Cpi", ",", "1.5"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"w", ",", "1.1"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"p", ",", "0.5"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"AccuracyGoal", "\[Rule]", "6"}], ",", 
        RowBox[{"PrecisionGoal", "\[Rule]", "6"}], ",", 
        RowBox[{"Method", "\[Rule]", "\"\<Newton\>\""}], ",", 
        RowBox[{"MaxIterations", "\[Rule]", "800"}]}], "}"}]}], "]"}]}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"t", "=", 
   RowBox[{"{", 
    RowBox[{"func", ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"Cr", "\[Rule]", "0.26"}], ",", 
       RowBox[{"CC", "\[Rule]", "1.0"}], ",", 
       RowBox[{"CH", "\[Rule]", "1.22"}], ",", 
       RowBox[{"Cpi", "\[Rule]", "1.4"}], ",", 
       RowBox[{"w", "\[Rule]", "1.02"}], ",", 
       RowBox[{"p", "\[Rule]", "0.56"}]}], "}"}]}], "}"}]}], 
  ";"}], "\[IndentingNewLine]", "func"}], "Input",
 CellChangeTimes->{{3.5675987354944196`*^9, 3.567599332289068*^9}, {
   3.567599364347124*^9, 3.567599388854767*^9}, {3.567599446824469*^9, 
   3.567599671277663*^9}, {3.5675997743158445`*^9, 3.567599879038828*^9}, {
   3.567599920722101*^9, 3.567600053883935*^9}, {3.5676001118380375`*^9, 
   3.5676001705877404`*^9}, {3.5676002126766143`*^9, 
   3.5676003454796476`*^9}, {3.5676004083477583`*^9, 3.567600442527418*^9}, {
   3.5676004745698743`*^9, 3.567600484132691*^9}, {3.5676005145995445`*^9, 
   3.5676005320091753`*^9}, {3.5676005791056576`*^9, 
   3.5676012559440465`*^9}, {3.5676013203097596`*^9, 
   3.5676013219321623`*^9}, {3.5676013524770164`*^9, 
   3.5676014283711495`*^9}, {3.567604737731021*^9, 3.5676047599610596`*^9}, {
   3.567611568287159*^9, 3.5676115840743866`*^9}, {3.5676124995619946`*^9, 
   3.5676125004199963`*^9}, {3.5676147541247544`*^9, 3.567614755169956*^9}, {
   3.567615587415818*^9, 3.5676155981798368`*^9}, {3.5676156384067073`*^9, 
   3.5676156583747425`*^9}, {3.5676158998007665`*^9, 
   3.5676159408288383`*^9}, {3.567615989126523*^9, 3.567616164252431*^9}, {
   3.5676162788658323`*^9, 3.567616283686241*^9}, {3.567616477313781*^9, 
   3.5676165242074633`*^9}, {3.5676165629735312`*^9, 
   3.5676165647363343`*^9}, {3.5676166087284117`*^9, 
   3.5676166312860513`*^9}, {3.5676169704150467`*^9, 3.567617073671628*^9}, {
   3.567617106774886*^9, 3.5676171121256957`*^9}, 3.567617182341419*^9, {
   3.567617294022015*^9, 3.567617336079689*^9}, {3.567660899793217*^9, 
   3.5676609125384398`*^9}, {3.5676609439100943`*^9, 3.567660949073704*^9}, {
   3.5676610469482756`*^9, 3.567661068710314*^9}, 3.567661166756486*^9, {
   3.5676850168515854`*^9, 3.5676850349164176`*^9}, {3.5676852124135294`*^9, 
   3.567685248761593*^9}, {3.5676856839243574`*^9, 3.5676856880427647`*^9}, {
   3.5676858185057936`*^9, 3.567685824714605*^9}, 3.567686114506314*^9, {
   3.567686148249173*^9, 3.5676861890588446`*^9}, {3.567686447551299*^9, 
   3.5676864707953396`*^9}, {3.567686506925003*^9, 3.567686565409506*^9}, {
   3.567686809066334*^9, 3.567686826366764*^9}, {3.5676874811935143`*^9, 
   3.5676874818643155`*^9}, {3.5676921395598993`*^9, 
   3.5676921509791193`*^9}, {3.567692753654978*^9, 3.5676927950730505`*^9}, {
   3.567698206184089*^9, 3.5676982494429655`*^9}, 3.567698778985895*^9, {
   3.5676988631324434`*^9, 3.567698879060071*^9}, {3.5676990016918864`*^9, 
   3.567699072718811*^9}, {3.5677005516170087`*^9, 3.5677006099767113`*^9}, {
   3.567701125339216*^9, 3.5677011286464224`*^9}, {3.567701158676475*^9, 
   3.567701160096078*^9}}],

Cell[BoxData[
 RowBox[{"-", "86.55501699978875`"}]], "Output",
 CellChangeTimes->{3.5677011359940352`*^9, 3.5677011746977034`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"Ekin", "/.", 
  RowBox[{"t", "[", 
   RowBox[{"[", "2", "]"}], "]"}]}], "\n", 
 RowBox[{"vne", "/.", 
  RowBox[{"t", "[", 
   RowBox[{"[", "2", "]"}], "]"}]}], "\n", 
 RowBox[{"vee", "/.", 
  RowBox[{"t", "[", 
   RowBox[{"[", "2", "]"}], "]"}]}], "\n", 
 RowBox[{"vnn", "/.", 
  RowBox[{"t", "[", 
   RowBox[{"[", "2", "]"}], "]"}]}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"-", "Epot"}], "/", "Ekin"}], "/.", 
  RowBox[{"t", "[", 
   RowBox[{"[", "2", "]"}], "]"}]}], "\n", 
 RowBox[{
  RowBox[{"2", "*", "w", "*", "rad"}], "/.", 
  RowBox[{"t", "[", 
   RowBox[{"[", "2", "]"}], "]"}]}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"(", 
    RowBox[{
     RowBox[{"2", "*", "Cr"}], "+", 
     RowBox[{"2", "*", "CC"}]}], ")"}], "*", "bohr"}], "/.", 
  RowBox[{"t", "[", 
   RowBox[{"[", "2", "]"}], "]"}]}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"(", 
    RowBox[{"Cr", "+", "CH", "+", "p"}], ")"}], "*", "bohr"}], "/.", 
  RowBox[{"t", "[", 
   RowBox[{"[", "2", "]"}], "]"}]}]}], "Input",
 CellChangeTimes->{{3.5675987354944196`*^9, 3.567599332289068*^9}, {
   3.567599364347124*^9, 3.567599388854767*^9}, {3.567599446824469*^9, 
   3.567599671277663*^9}, {3.5675997743158445`*^9, 3.567599879038828*^9}, {
   3.567599920722101*^9, 3.567600053883935*^9}, {3.5676001118380375`*^9, 
   3.5676001705877404`*^9}, {3.5676002126766143`*^9, 
   3.5676003454796476`*^9}, {3.5676004083477583`*^9, 3.567600442527418*^9}, {
   3.5676004745698743`*^9, 3.567600484132691*^9}, {3.5676005145995445`*^9, 
   3.5676005320091753`*^9}, {3.5676005791056576`*^9, 
   3.5676012559440465`*^9}, {3.5676013203097596`*^9, 
   3.5676013219321623`*^9}, {3.5676013524770164`*^9, 
   3.5676014283711495`*^9}, {3.5676166536252904`*^9, 3.567616654654892*^9}, 
   3.567685840720233*^9}],

Cell[BoxData["83.07719825532563`"], "Output",
 CellChangeTimes->{
  3.567692193972795*^9, {3.5676927717354097`*^9, 3.5676928003302603`*^9}, 
   3.5676982222677174`*^9, 3.567698257196179*^9, 3.567698749704644*^9, 
   3.5676987856939073`*^9, 3.5676988896056895`*^9, 3.5676991890638156`*^9, 
   3.567699602355342*^9, 3.5676997489175987`*^9, 3.5677004218247805`*^9, 
   3.567700488608498*^9, 3.5677006238451357`*^9, 3.5677006895212507`*^9, 
   3.567700977716157*^9, {3.567701092235958*^9, 3.5677011049031806`*^9}, 
   3.567701136009635*^9, 3.567701174713303*^9}],

Cell[BoxData[
 RowBox[{"-", "262.3448159960881`"}]], "Output",
 CellChangeTimes->{
  3.567692193972795*^9, {3.5676927717354097`*^9, 3.5676928003302603`*^9}, 
   3.5676982222677174`*^9, 3.567698257196179*^9, 3.567698749704644*^9, 
   3.5676987856939073`*^9, 3.5676988896056895`*^9, 3.5676991890638156`*^9, 
   3.567699602355342*^9, 3.5676997489175987`*^9, 3.5677004218247805`*^9, 
   3.567700488608498*^9, 3.5677006238451357`*^9, 3.5677006895212507`*^9, 
   3.567700977716157*^9, {3.567701092235958*^9, 3.5677011049031806`*^9}, 
   3.567701136009635*^9, 3.567701174713303*^9}],

Cell[BoxData["59.29291505160059`"], "Output",
 CellChangeTimes->{
  3.567692193972795*^9, {3.5676927717354097`*^9, 3.5676928003302603`*^9}, 
   3.5676982222677174`*^9, 3.567698257196179*^9, 3.567698749704644*^9, 
   3.5676987856939073`*^9, 3.5676988896056895`*^9, 3.5676991890638156`*^9, 
   3.567699602355342*^9, 3.5676997489175987`*^9, 3.5677004218247805`*^9, 
   3.567700488608498*^9, 3.5677006238451357`*^9, 3.5677006895212507`*^9, 
   3.567700977716157*^9, {3.567701092235958*^9, 3.5677011049031806`*^9}, 
   3.567701136009635*^9, 3.567701174713303*^9}],

Cell[BoxData["33.41968568937313`"], "Output",
 CellChangeTimes->{
  3.567692193972795*^9, {3.5676927717354097`*^9, 3.5676928003302603`*^9}, 
   3.5676982222677174`*^9, 3.567698257196179*^9, 3.567698749704644*^9, 
   3.5676987856939073`*^9, 3.5676988896056895`*^9, 3.5676991890638156`*^9, 
   3.567699602355342*^9, 3.5676997489175987`*^9, 3.5677004218247805`*^9, 
   3.567700488608498*^9, 3.5677006238451357`*^9, 3.5677006895212507`*^9, 
   3.567700977716157*^9, {3.567701092235958*^9, 3.5677011049031806`*^9}, 
   3.567701136009635*^9, 3.567701174713303*^9}],

Cell[BoxData["2.0418624943726984`"], "Output",
 CellChangeTimes->{
  3.567692193972795*^9, {3.5676927717354097`*^9, 3.5676928003302603`*^9}, 
   3.5676982222677174`*^9, 3.567698257196179*^9, 3.567698749704644*^9, 
   3.5676987856939073`*^9, 3.5676988896056895`*^9, 3.5676991890638156`*^9, 
   3.567699602355342*^9, 3.5676997489175987`*^9, 3.5677004218247805`*^9, 
   3.567700488608498*^9, 3.5677006238451357`*^9, 3.5677006895212507`*^9, 
   3.567700977716157*^9, {3.567701092235958*^9, 3.5677011049031806`*^9}, 
   3.567701136009635*^9, 3.5677011747289033`*^9}],

Cell[BoxData["116.8833912`"], "Output",
 CellChangeTimes->{
  3.567692193972795*^9, {3.5676927717354097`*^9, 3.5676928003302603`*^9}, 
   3.5676982222677174`*^9, 3.567698257196179*^9, 3.567698749704644*^9, 
   3.5676987856939073`*^9, 3.5676988896056895`*^9, 3.5676991890638156`*^9, 
   3.567699602355342*^9, 3.5676997489175987`*^9, 3.5677004218247805`*^9, 
   3.567700488608498*^9, 3.5677006238451357`*^9, 3.5677006895212507`*^9, 
   3.567700977716157*^9, {3.567701092235958*^9, 3.5677011049031806`*^9}, 
   3.567701136009635*^9, 3.5677011747289033`*^9}],

Cell[BoxData["1.33352604`"], "Output",
 CellChangeTimes->{
  3.567692193972795*^9, {3.5676927717354097`*^9, 3.5676928003302603`*^9}, 
   3.5676982222677174`*^9, 3.567698257196179*^9, 3.567698749704644*^9, 
   3.5676987856939073`*^9, 3.5676988896056895`*^9, 3.5676991890638156`*^9, 
   3.567699602355342*^9, 3.5676997489175987`*^9, 3.5677004218247805`*^9, 
   3.567700488608498*^9, 3.5677006238451357`*^9, 3.5677006895212507`*^9, 
   3.567700977716157*^9, {3.567701092235958*^9, 3.5677011049031806`*^9}, 
   3.567701136009635*^9, 3.5677011747289033`*^9}],

Cell[BoxData["1.0795210800000001`"], "Output",
 CellChangeTimes->{
  3.567692193972795*^9, {3.5676927717354097`*^9, 3.5676928003302603`*^9}, 
   3.5676982222677174`*^9, 3.567698257196179*^9, 3.567698749704644*^9, 
   3.5676987856939073`*^9, 3.5676988896056895`*^9, 3.5676991890638156`*^9, 
   3.567699602355342*^9, 3.5676997489175987`*^9, 3.5677004218247805`*^9, 
   3.567700488608498*^9, 3.5677006238451357`*^9, 3.5677006895212507`*^9, 
   3.567700977716157*^9, {3.567701092235958*^9, 3.5677011049031806`*^9}, 
   3.567701136009635*^9, 3.5677011747289033`*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"plot1", "=", 
    RowBox[{
     RowBox[{"Graphics", "[", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"xc", "[", 
             RowBox[{"[", "1", "]"}], "]"}], ",", 
            RowBox[{"yc", "[", 
             RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", "Cr"}], "]"}], ",",
         "\[IndentingNewLine]", 
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"xc", "[", 
             RowBox[{"[", "2", "]"}], "]"}], ",", 
            RowBox[{"yc", "[", 
             RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", "Cr"}], "]"}], ",",
         "\[IndentingNewLine]", 
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"xc", "[", 
             RowBox[{"[", "3", "]"}], "]"}], ",", 
            RowBox[{"yc", "[", 
             RowBox[{"[", "3", "]"}], "]"}]}], "}"}], ",", "CC"}], "]"}], ",",
         "\[IndentingNewLine]", 
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"xc", "[", 
             RowBox[{"[", "4", "]"}], "]"}], ",", 
            RowBox[{"yc", "[", 
             RowBox[{"[", "4", "]"}], "]"}]}], "}"}], ",", "CH"}], "]"}], ",",
         "\n", "            ", 
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"xc", "[", 
             RowBox[{"[", "5", "]"}], "]"}], ",", 
            RowBox[{"yc", "[", 
             RowBox[{"[", "5", "]"}], "]"}]}], "}"}], ",", "CH"}], "]"}], ",",
         "\[IndentingNewLine]", 
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"xc", "[", 
             RowBox[{"[", "6", "]"}], "]"}], ",", 
            RowBox[{"yc", "[", 
             RowBox[{"[", "6", "]"}], "]"}]}], "}"}], ",", "CH"}], "]"}], ",",
         "\[IndentingNewLine]", 
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"xc", "[", 
             RowBox[{"[", "7", "]"}], "]"}], ",", 
            RowBox[{"yc", "[", 
             RowBox[{"[", "7", "]"}], "]"}]}], "}"}], ",", "CH"}], "]"}], ",",
         "\[IndentingNewLine]", 
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"xc", "[", 
             RowBox[{"[", "8", "]"}], "]"}], ",", 
            RowBox[{"yc", "[", 
             RowBox[{"[", "8", "]"}], "]"}]}], "}"}], ",", "Cpi"}], "]"}], 
        ",", "\[IndentingNewLine]", 
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"xc", "[", 
             RowBox[{"[", "9", "]"}], "]"}], ",", 
            RowBox[{"yc", "[", 
             RowBox[{"[", "9", "]"}], "]"}]}], "}"}], ",", "Cpi"}], "]"}], 
        ",", "\[IndentingNewLine]", 
        RowBox[{"Disk", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"xn", "[", 
             RowBox[{"[", "1", "]"}], "]"}], ",", 
            RowBox[{"yn", "[", 
             RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], 
        ",", " ", "\[IndentingNewLine]", 
        RowBox[{"Disk", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"xn", "[", 
             RowBox[{"[", "2", "]"}], "]"}], ",", 
            RowBox[{"yn", "[", 
             RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], 
        ",", " ", "\[IndentingNewLine]", 
        RowBox[{"Disk", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"xn", "[", 
             RowBox[{"[", "3", "]"}], "]"}], ",", 
            RowBox[{"yn", "[", 
             RowBox[{"[", "3", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], 
        ",", "\[IndentingNewLine]", 
        RowBox[{"Disk", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"xn", "[", 
             RowBox[{"[", "4", "]"}], "]"}], ",", 
            RowBox[{"yn", "[", 
             RowBox[{"[", "4", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], 
        ",", " ", "\[IndentingNewLine]", 
        RowBox[{"Disk", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"xn", "[", 
             RowBox[{"[", "5", "]"}], "]"}], ",", 
            RowBox[{"yn", "[", 
             RowBox[{"[", "5", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], 
        ",", "\[IndentingNewLine]", 
        RowBox[{"Disk", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"xn", "[", 
             RowBox[{"[", "6", "]"}], "]"}], ",", 
            RowBox[{"yn", "[", 
             RowBox[{"[", "6", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], 
        ",", "\n", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"Thickness", "[", "0.008", "]"}], ",", 
          RowBox[{"Line", "[", 
           RowBox[{"{", 
            RowBox[{
             RowBox[{"{", 
              RowBox[{
               RowBox[{"xn", "[", 
                RowBox[{"[", "2", "]"}], "]"}], ",", 
               RowBox[{"yn", "[", 
                RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", 
             RowBox[{"{", 
              RowBox[{
               RowBox[{"xn", "[", 
                RowBox[{"[", "1", "]"}], "]"}], ",", 
               RowBox[{"yn", "[", 
                RowBox[{"[", "1", "]"}], "]"}]}], "}"}]}], "}"}], "]"}], ",", 
          "\[IndentingNewLine]", 
          RowBox[{"Line", "[", 
           RowBox[{"{", 
            RowBox[{
             RowBox[{"{", 
              RowBox[{
               RowBox[{"xn", "[", 
                RowBox[{"[", "1", "]"}], "]"}], ",", 
               RowBox[{"yn", "[", 
                RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", 
             RowBox[{"{", 
              RowBox[{
               RowBox[{"xn", "[", 
                RowBox[{"[", "3", "]"}], "]"}], ",", 
               RowBox[{"yn", "[", 
                RowBox[{"[", "3", "]"}], "]"}]}], "}"}]}], "}"}], "]"}], ",", 
          "\[IndentingNewLine]", 
          RowBox[{"Line", "[", 
           RowBox[{"{", 
            RowBox[{
             RowBox[{"{", 
              RowBox[{
               RowBox[{"xn", "[", 
                RowBox[{"[", "1", "]"}], "]"}], ",", 
               RowBox[{"yn", "[", 
                RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", 
             RowBox[{"{", 
              RowBox[{
               RowBox[{"xn", "[", 
                RowBox[{"[", "4", "]"}], "]"}], ",", 
               RowBox[{"yn", "[", 
                RowBox[{"[", "4", "]"}], "]"}]}], "}"}]}], "}"}], "]"}], ",", 
          "\[IndentingNewLine]", 
          RowBox[{"Line", "[", 
           RowBox[{"{", 
            RowBox[{
             RowBox[{"{", 
              RowBox[{
               RowBox[{"xn", "[", 
                RowBox[{"[", "2", "]"}], "]"}], ",", 
               RowBox[{"yn", "[", 
                RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", 
             RowBox[{"{", 
              RowBox[{
               RowBox[{"xn", "[", 
                RowBox[{"[", "5", "]"}], "]"}], ",", 
               RowBox[{"yn", "[", 
                RowBox[{"[", "5", "]"}], "]"}]}], "}"}]}], "}"}], "]"}], ",", 
          "\[IndentingNewLine]", 
          RowBox[{"Line", "[", 
           RowBox[{"{", 
            RowBox[{
             RowBox[{"{", 
              RowBox[{
               RowBox[{"xn", "[", 
                RowBox[{"[", "2", "]"}], "]"}], ",", 
               RowBox[{"yn", "[", 
                RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", 
             RowBox[{"{", 
              RowBox[{
               RowBox[{"xn", "[", 
                RowBox[{"[", "6", "]"}], "]"}], ",", 
               RowBox[{"yn", "[", 
                RowBox[{"[", "6", "]"}], "]"}]}], "}"}]}], "}"}], "]"}]}], 
         "}"}]}], "}"}], "]"}], " ", "/.", " ", 
     RowBox[{"t", "[", 
      RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], "\n"}]], "Input",
 CellChangeTimes->{{3.567661245973425*^9, 3.567661380742062*^9}, {
   3.567661415998124*^9, 3.5676614516285863`*^9}, {3.5676614877114496`*^9, 
   3.567661525089115*^9}, {3.5676615670531893`*^9, 3.5676615719827976`*^9}, {
   3.5676616025120516`*^9, 3.567661696112216*^9}, {3.567663489803366*^9, 
   3.5676635070881968`*^9}, {3.567663544450262*^9, 3.567663576336718*^9}, {
   3.567663626366006*^9, 3.567663629704412*^9}, {3.567692320255017*^9, 
   3.567692403044362*^9}, 3.5676924445872355`*^9, {3.567692490482516*^9, 
   3.56769258937109*^9}, {3.5676926267487555`*^9, 3.567692699054882*^9}, {
   3.567698728192206*^9, 3.567698729923809*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{"xy", "-", "plane"}], " ", "*)"}], "\n", 
  RowBox[{"Show", "[", 
   RowBox[{"plot1", ",", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"AspectRatio", " ", "\[Rule]", " ", "Automatic"}], ",", 
      RowBox[{"Axes", " ", "->", " ", "True"}], ",", 
      RowBox[{"GridLines", " ", "->", " ", "Automatic"}], ",", " ", 
      RowBox[{"PlotRange", " ", "\[Rule]", " ", "Automatic"}], ",", " ", 
      RowBox[{"Frame", " ", "->", " ", "True"}]}], "}"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.567661245973425*^9, 3.567661380742062*^9}, {
   3.567661415998124*^9, 3.5676614516285863`*^9}, {3.5676614877114496`*^9, 
   3.567661525089115*^9}, {3.5676615670531893`*^9, 3.5676615719827976`*^9}, {
   3.5676616025120516`*^9, 3.567661696112216*^9}, {3.567663489803366*^9, 
   3.5676635070881968`*^9}, {3.567663544450262*^9, 3.567663576336718*^9}, {
   3.567663626366006*^9, 3.567663629704412*^9}, {3.567692320255017*^9, 
   3.567692403044362*^9}, 3.5676924445872355`*^9}],

Cell[BoxData[
 GraphicsBox[{CircleBox[{1.26, 0}, 0.26], CircleBox[{-1.26, 0}, 0.26], 
   CircleBox[{0, 0}], 
   CircleBox[{2.0345816078524415`, 1.2611198724850572`}, 1.22], 
   CircleBox[{2.0345816078524415`, -1.2611198724850572`}, 1.22], 
   CircleBox[{-2.0345816078524415`, 1.2611198724850572`}, 1.22], 
   CircleBox[{-2.0345816078524415`, -1.2611198724850572`}, 1.22], 
   CircleBox[{0, 0}, 1.4], CircleBox[{0, 0}, 1.4], DiskBox[{1.26, 0}, 0.08], 
   DiskBox[{-1.26, 0}, 0.08], 
   DiskBox[{2.327666540553365, 1.7383003647767006`}, 0.08], 
   DiskBox[{2.327666540553365, -1.7383003647767006`}, 0.08], 
   DiskBox[{-2.327666540553365, 1.7383003647767006`}, 0.08], 
   DiskBox[{-2.327666540553365, -1.7383003647767006`}, 0.08], 
   {Thickness[0.008], LineBox[{{-1.26, 0}, {1.26, 0}}], 
    LineBox[{{1.26, 0}, {2.327666540553365, 1.7383003647767006`}}], 
    LineBox[{{1.26, 0}, {2.327666540553365, -1.7383003647767006`}}], 
    LineBox[{{-1.26, 0}, {-2.327666540553365, 1.7383003647767006`}}], 
    LineBox[{{-1.26, 0}, {-2.327666540553365, -1.7383003647767006`}}]}},
  AspectRatio->Automatic,
  Axes->True,
  Frame->True,
  GridLines->Automatic,
  PlotRange->Automatic]], "Output",
 CellChangeTimes->{
  3.567692194003995*^9, {3.5676927717822104`*^9, 3.5676928003614597`*^9}, 
   3.567698222298918*^9, 3.567698257227379*^9, 3.5676987497358437`*^9, 
   3.5676987857407074`*^9, 3.5676988896368895`*^9, 3.5676991890950155`*^9, 
   3.5676996023865414`*^9, 3.567699748948799*^9, 3.567700421855981*^9, 
   3.567700488655298*^9, 3.5677006238763356`*^9, 3.567700689552451*^9, 
   3.567700977762957*^9, {3.5677010922671585`*^9, 3.5677011049499807`*^9}, 
   3.567701136056435*^9, 3.567701174744503*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"plot2", "=", 
   RowBox[{
    RowBox[{"Graphics", "[", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"Circle", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"xc", "[", 
            RowBox[{"[", "1", "]"}], "]"}], ",", 
           RowBox[{"zc", "[", 
            RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", "Cr"}], "]"}], ",", 
       "\[IndentingNewLine]", 
       RowBox[{"Circle", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"xc", "[", 
            RowBox[{"[", "2", "]"}], "]"}], ",", 
           RowBox[{"zc", "[", 
            RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", "Cr"}], "]"}], ",", 
       "\[IndentingNewLine]", 
       RowBox[{"Circle", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"xc", "[", 
            RowBox[{"[", "3", "]"}], "]"}], ",", 
           RowBox[{"zc", "[", 
            RowBox[{"[", "3", "]"}], "]"}]}], "}"}], ",", "CC"}], "]"}], ",", 
       "\[IndentingNewLine]", 
       RowBox[{"Circle", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"xc", "[", 
            RowBox[{"[", "4", "]"}], "]"}], ",", 
           RowBox[{"zc", "[", 
            RowBox[{"[", "4", "]"}], "]"}]}], "}"}], ",", "CH"}], "]"}], ",", 
       "\n", "            ", 
       RowBox[{"Circle", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"xc", "[", 
            RowBox[{"[", "5", "]"}], "]"}], ",", 
           RowBox[{"zc", "[", 
            RowBox[{"[", "5", "]"}], "]"}]}], "}"}], ",", "CH"}], "]"}], ",", 
       "\[IndentingNewLine]", 
       RowBox[{"Circle", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"xc", "[", 
            RowBox[{"[", "6", "]"}], "]"}], ",", 
           RowBox[{"zc", "[", 
            RowBox[{"[", "6", "]"}], "]"}]}], "}"}], ",", "CH"}], "]"}], ",", 
       "\[IndentingNewLine]", 
       RowBox[{"Circle", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"xc", "[", 
            RowBox[{"[", "7", "]"}], "]"}], ",", 
           RowBox[{"zc", "[", 
            RowBox[{"[", "7", "]"}], "]"}]}], "}"}], ",", "CH"}], "]"}], ",", 
       "\[IndentingNewLine]", 
       RowBox[{"Circle", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"xc", "[", 
            RowBox[{"[", "8", "]"}], "]"}], ",", 
           RowBox[{"zc", "[", 
            RowBox[{"[", "8", "]"}], "]"}]}], "}"}], ",", "Cpi"}], "]"}], ",",
        "\[IndentingNewLine]", 
       RowBox[{"Circle", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"xc", "[", 
            RowBox[{"[", "9", "]"}], "]"}], ",", 
           RowBox[{"zc", "[", 
            RowBox[{"[", "9", "]"}], "]"}]}], "}"}], ",", "Cpi"}], "]"}], ",",
        "\[IndentingNewLine]", 
       RowBox[{"Disk", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"xn", "[", 
            RowBox[{"[", "1", "]"}], "]"}], ",", 
           RowBox[{"zn", "[", 
            RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], 
       ",", " ", "\[IndentingNewLine]", 
       RowBox[{"Disk", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"xn", "[", 
            RowBox[{"[", "2", "]"}], "]"}], ",", 
           RowBox[{"zn", "[", 
            RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], 
       ",", " ", "\[IndentingNewLine]", 
       RowBox[{"Disk", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"xn", "[", 
            RowBox[{"[", "3", "]"}], "]"}], ",", 
           RowBox[{"zn", "[", 
            RowBox[{"[", "3", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], 
       ",", "\[IndentingNewLine]", 
       RowBox[{"Disk", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"xn", "[", 
            RowBox[{"[", "4", "]"}], "]"}], ",", 
           RowBox[{"zn", "[", 
            RowBox[{"[", "4", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], 
       ",", " ", "\[IndentingNewLine]", 
       RowBox[{"Disk", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"xn", "[", 
            RowBox[{"[", "5", "]"}], "]"}], ",", 
           RowBox[{"zn", "[", 
            RowBox[{"[", "5", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], 
       ",", "\[IndentingNewLine]", 
       RowBox[{"Disk", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"xn", "[", 
            RowBox[{"[", "6", "]"}], "]"}], ",", 
           RowBox[{"zn", "[", 
            RowBox[{"[", "6", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], 
       ",", "\n", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"Thickness", "[", "0.008", "]"}], ",", "\n", 
         RowBox[{"Line", "[", 
          RowBox[{"{", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{"xn", "[", 
               RowBox[{"[", "2", "]"}], "]"}], ",", 
              RowBox[{"zn", "[", 
               RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", 
            RowBox[{"{", 
             RowBox[{
              RowBox[{"xn", "[", 
               RowBox[{"[", "1", "]"}], "]"}], ",", 
              RowBox[{"zn", "[", 
               RowBox[{"[", "1", "]"}], "]"}]}], "}"}]}], "}"}], "]"}], ",", 
         "\n", 
         RowBox[{"Line", "[", 
          RowBox[{"{", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{"xn", "[", 
               RowBox[{"[", "1", "]"}], "]"}], ",", 
              RowBox[{"zn", "[", 
               RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", 
            RowBox[{"{", 
             RowBox[{
              RowBox[{"xn", "[", 
               RowBox[{"[", "3", "]"}], "]"}], ",", 
              RowBox[{"zn", "[", 
               RowBox[{"[", "3", "]"}], "]"}]}], "}"}]}], "}"}], "]"}], ",", 
         "\n", 
         RowBox[{"Line", "[", 
          RowBox[{"{", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{"xn", "[", 
               RowBox[{"[", "1", "]"}], "]"}], ",", 
              RowBox[{"zn", "[", 
               RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", 
            RowBox[{"{", 
             RowBox[{
              RowBox[{"xn", "[", 
               RowBox[{"[", "4", "]"}], "]"}], ",", 
              RowBox[{"zn", "[", 
               RowBox[{"[", "4", "]"}], "]"}]}], "}"}]}], "}"}], "]"}], ",", 
         "\n", 
         RowBox[{"Line", "[", 
          RowBox[{"{", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{"xn", "[", 
               RowBox[{"[", "2", "]"}], "]"}], ",", 
              RowBox[{"zn", "[", 
               RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", 
            RowBox[{"{", 
             RowBox[{
              RowBox[{"xn", "[", 
               RowBox[{"[", "5", "]"}], "]"}], ",", 
              RowBox[{"zn", "[", 
               RowBox[{"[", "5", "]"}], "]"}]}], "}"}]}], "}"}], "]"}], ",", 
         RowBox[{"Line", "[", 
          RowBox[{"{", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{"xn", "[", 
               RowBox[{"[", "2", "]"}], "]"}], ",", 
              RowBox[{"zn", "[", 
               RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", 
            RowBox[{"{", 
             RowBox[{
              RowBox[{"xn", "[", 
               RowBox[{"[", "6", "]"}], "]"}], ",", 
              RowBox[{"zn", "[", 
               RowBox[{"[", "6", "]"}], "]"}]}], "}"}]}], "}"}], "]"}]}], 
        "}"}]}], "}"}], " ", "]"}], " ", "/.", " ", 
    RowBox[{"t", "[", 
     RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}]], "Input",
 CellChangeTimes->{{3.5676995778008986`*^9, 3.5676995871765146`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{"xz", "-", "plane"}], " ", "*)"}], "\n", 
  RowBox[{"Show", "[", 
   RowBox[{"plot2", ",", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"AspectRatio", " ", "\[Rule]", " ", "Automatic"}], ",", 
      RowBox[{"Axes", " ", "->", " ", "True"}], ",", 
      RowBox[{"GridLines", " ", "->", " ", "Automatic"}], ",", " ", 
      RowBox[{"PlotRange", " ", "\[Rule]", " ", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"-", "3.5"}], ",", "3.5"}], "}"}], ",", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"-", "3.5"}], ",", "3.5"}], "}"}]}], "}"}]}], ",", " ", 
      RowBox[{"Frame", " ", "->", " ", "True"}]}], "}"}]}], "]"}]}]], "Input"],

Cell[BoxData[
 GraphicsBox[{CircleBox[{1.26, 0}, 0.26], CircleBox[{-1.26, 0}, 0.26], 
   CircleBox[{0, 0}], CircleBox[{2.0345816078524415`, 0}, 1.22], 
   CircleBox[{2.0345816078524415`, 0}, 1.22], 
   CircleBox[{-2.0345816078524415`, 0}, 1.22], 
   CircleBox[{-2.0345816078524415`, 0}, 1.22], 
   CircleBox[{0, 1.7333333333333332`}, 1.4], 
   CircleBox[{0, -1.7333333333333332`}, 1.4], DiskBox[{1.26, 0}, 0.08], 
   DiskBox[{-1.26, 0}, 0.08], DiskBox[{2.327666540553365, 0}, 0.08], 
   DiskBox[{2.327666540553365, 0}, 0.08], 
   DiskBox[{-2.327666540553365, 0}, 0.08], 
   DiskBox[{-2.327666540553365, 0}, 0.08], 
   {Thickness[0.008], LineBox[{{-1.26, 0}, {1.26, 0}}], 
    LineBox[{{1.26, 0}, {2.327666540553365, 0}}], 
    LineBox[{{1.26, 0}, {2.327666540553365, 0}}], 
    LineBox[{{-1.26, 0}, {-2.327666540553365, 0}}], 
    LineBox[{{-1.26, 0}, {-2.327666540553365, 0}}]}},
  AspectRatio->Automatic,
  Axes->True,
  Frame->True,
  GridLines->Automatic,
  PlotRange->{{-3.5, 3.5}, {-3.5, 3.5}}]], "Output",
 CellChangeTimes->{
  3.5676996024021416`*^9, 3.5676997489643993`*^9, 3.567700421871581*^9, 
   3.5677004886708984`*^9, 3.5677006238919353`*^9, 3.5677006895836515`*^9, 
   3.5677009777785573`*^9, {3.5677010922827587`*^9, 3.567701104965581*^9}, 
   3.5677011360720353`*^9, 3.5677011747757034`*^9}]
}, Open  ]]
},
WindowSize->{1296, 1002},
WindowMargins->{{235, Automatic}, {Automatic, -5}},
Magnification->1.5,
FrontEndVersion->"9.0 for Microsoft Windows (64-bit) (November 20, 2012)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[557, 20, 857, 14, 45, "Input"],
Cell[1417, 36, 2479, 46, 249, "Input"],
Cell[3899, 84, 929, 16, 45, "Input"],
Cell[4831, 102, 840, 15, 75, "Input"],
Cell[5674, 119, 1307, 28, 162, "Input"],
Cell[6984, 149, 761, 11, 45, "Input"],
Cell[7748, 162, 2049, 50, 104, "Input"],
Cell[9800, 214, 1355, 29, 104, "Input"],
Cell[11158, 245, 2151, 52, 162, "Input"],
Cell[13312, 299, 1477, 29, 133, "Input"],
Cell[14792, 330, 1782, 35, 133, "Input"],
Cell[16577, 367, 3083, 73, 162, "Input"],
Cell[19663, 442, 3380, 79, 220, "Input"],
Cell[23046, 523, 3684, 88, 220, "Input"],
Cell[26733, 613, 1369, 23, 45, "Input"],
Cell[CellGroupData[{
Cell[28127, 640, 4068, 80, 191, "Input"],
Cell[32198, 722, 130, 2, 63, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[32365, 729, 1778, 47, 249, "Input"],
Cell[34146, 778, 558, 8, 63, "Output"],
Cell[34707, 788, 575, 9, 63, "Output"],
Cell[35285, 799, 558, 8, 63, "Output"],
Cell[35846, 809, 558, 8, 63, "Output"],
Cell[36407, 819, 561, 8, 63, "Output"],
Cell[36971, 829, 554, 8, 63, "Output"],
Cell[37528, 839, 553, 8, 63, "Output"],
Cell[38084, 849, 561, 8, 63, "Output"]
}, Open  ]],
Cell[38660, 860, 8806, 236, 626, "Input"],
Cell[CellGroupData[{
Cell[47491, 1100, 1021, 19, 104, "Input"],
Cell[48515, 1121, 1695, 30, 461, "Output"]
}, Open  ]],
Cell[50225, 1154, 7918, 226, 626, "Input"],
Cell[CellGroupData[{
Cell[58168, 1384, 760, 20, 104, "Input"],
Cell[58931, 1406, 1309, 26, 580, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)