(* 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[    188353,       4043]
NotebookOptionsPosition[    182574,       3857]
NotebookOutlinePosition[    183108,       3878]
CellTagsIndexPosition[    183065,       3875]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{
    "Kimball", " ", "model", " ", "and", " ", "energy", " ", "derived", " ", 
     "from", " ", "experimental", " ", "molecule", " ", "data", " ", 
     "source"}], ",", " ", 
    RowBox[{"ES", " ", "Febuary", " ", "2013"}]}], " ", "*)"}], " ", 
  "\[IndentingNewLine]", 
  RowBox[{"(*", " ", 
   RowBox[{"sch\[ODoubleDot]n", ",", " ", 
    RowBox[{"aber", " ", 
     RowBox[{"falsch", "!"}], " ", 
     RowBox[{"check", "!"}], " ", "Vne"}]}], " ", "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"Clear", "[", 
    RowBox[{
    "tw", ",", "n2", ",", "n", ",", "nb", ",", "bc", ",", "bh", ",", "u", ",",
      "u1", ",", "u2", ",", "Vnn", ",", "Vne", ",", "Vee"}], "]"}], 
   ";"}]}]], "Input",
 CellChangeTimes->{{3.568278170270548*^9, 3.568278174310955*^9}, {
  3.5682782109554195`*^9, 3.5682782183342323`*^9}, {3.5688220426984916`*^9, 
  3.5688220506389055`*^9}, {3.568822092353379*^9, 3.568822134660653*^9}, {
  3.5688728241768503`*^9, 3.568872825034852*^9}, {3.5688958900531473`*^9, 
  3.5688958965115585`*^9}, {3.5689027504159966`*^9, 3.568902751695199*^9}, {
  3.5689109808408637`*^9, 3.568910992618884*^9}, {3.569079828786936*^9, 
  3.5690798397381554`*^9}, {3.569233977980655*^9, 3.5692340122851157`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"ChemicalData", "[", 
  RowBox[{"\"\<Ethane\>\"", ",", "\"\<MoleculePlot\>\""}], "]"}]], "Input",
 CellChangeTimes->{{3.568278170270548*^9, 3.568278174310955*^9}, {
   3.5682782109554195`*^9, 3.5682782183342323`*^9}, {3.5688220426984916`*^9, 
   3.5688220506389055`*^9}, {3.568822092353379*^9, 3.568822123693834*^9}, 
   3.568873215924739*^9, 3.5689108676002645`*^9, 3.5689997098237543`*^9, 
   3.5689998170115423`*^9}],

Cell[BoxData[
 Graphics3DBox[
  {AbsoluteThickness[3], EdgeForm[None], Specularity[
    GrayLevel[1], 100], GraphicsComplex3DBox[CompressedData["
1:eJwBYQGe/iFib1JlAgAAAA4AAAADAAAAvjCZKhhFUUASFD/G3BUqQKMjufyH
RD3ARUdy+Q9FUcDb+X5qvBQqwN21hHzQQz1AhetRuB4xYkCRD3o2q7ZEwLra
iv1lhz5Aak3zjlMAUkDBOSNKe6M2wLByaJHt1GDA7nw/NV6qVkA/NV66SYxe
QEw3iUFgdTrAObTIdr4zYsC2hHzQs7FEQGfV52ordj7ANqs+V1sBUsBvEoPA
yrE2QDEIrBxa1GBAMnctIR+kVsAv3SQGgY1ewLbz/dR4aTpA5IOezarTWkAa
ldQJaGIswGBxGw3gLeQ/FD/G3LWiUUBwXwfOGTETwCSX/5B+fVTA1lbsL7v3
U0DhC5OpgudQQHgtIR/03DvA3NeBc0bWWsB/jLlrCVksQKD4MeauJeO/PnlY
qDWjUcACK4cW2U4TQO2ePCzUfFRAPN9PjZf0U8A1XrpJDOhQwMpUwaik1jtA
G32fOw==
    "], {
     {RGBColor[0.65, 0.7, 0.7], SphereBox[3, 24.], SphereBox[4, 24.], 
      SphereBox[5, 24.], SphereBox[6, 24.], SphereBox[7, 24.], 
      SphereBox[8, 24.]}, 
     {RGBColor[0.4, 0.4, 0.4], SphereBox[1, 34.], SphereBox[2, 34.]}, 
     {RGBColor[0.65, 0.7, 0.7], CylinderBox[{9, 3}, 15.], 
      CylinderBox[{10, 4}, 15.], CylinderBox[{11, 5}, 15.], 
      CylinderBox[{12, 6}, 15.], CylinderBox[{13, 7}, 15.], 
      CylinderBox[{14, 8}, 15.]}, 
     {RGBColor[0.4, 0.4, 0.4], CylinderBox[{1, 2}, 15.], 
      CylinderBox[{1, 9}, 15.], CylinderBox[{1, 10}, 15.], 
      CylinderBox[{1, 11}, 15.], CylinderBox[{2, 12}, 15.], 
      CylinderBox[{2, 13}, 15.], CylinderBox[{2, 14}, 15.]}}]},
  Boxed->False,
  ImageSize->139.30525886879303`,
  Lighting->"Neutral",
  ViewPoint->{1.3, -2.4, -2.}]], "Output",
 CellChangeTimes->{{3.569080626931138*^9, 3.569080653607185*^9}, {
   3.5690807604673724`*^9, 3.5690808415563145`*^9}, 3.569080894955209*^9, {
   3.569080925749663*^9, 3.569080959991723*^9}, {3.569080993734582*^9, 
   3.569081008975809*^9}, 3.5690821724258523`*^9, 3.569088595284734*^9, 
   3.569172508393571*^9, 3.569176035265765*^9, 3.569176099537878*^9, 
   3.569176242371729*^9, 3.569176470163329*^9}]
}, Open  ]],

Cell[BoxData[{
 RowBox[{
  RowBox[{"t1", "=", 
   RowBox[{"0.01889726", "*", 
    RowBox[{"ChemicalData", "[", 
     RowBox[{"\"\<Ethane\>\"", ",", "\"\<AtomPositions\>\""}], "]"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"MatrixForm", "[", "t", "]"}], "*)"}]}]}], "Input",
 CellChangeTimes->{{3.567862912170579*^9, 3.567862949329844*^9}, {
   3.5678630095303497`*^9, 3.5678630107471523`*^9}, {3.5679573013548217`*^9, 
   3.5679573132576427`*^9}, 3.568114367533684*^9, {3.5682782383334675`*^9, 
   3.5682782435438766`*^9}, {3.568381606879155*^9, 3.5683816154123707`*^9}, 
   3.5688204675481253`*^9, {3.5688215945097046`*^9, 3.568821600578115*^9}, 
   3.568873215924739*^9, 3.5689108676002645`*^9, 3.5689997098393545`*^9, 
   3.5689998170271425`*^9}],

Cell[BoxData[
 RowBox[{"(*", 
  RowBox[{"ChemicalData", "\[NoBreak]", "[", 
   RowBox[{"\"\<Ethane\>\"", ",", "\"\<BondTally\>\""}], "]"}], 
  "*)"}]], "Input",
 CellChangeTimes->{{3.5678630319943895`*^9, 3.5678630956269007`*^9}, {
   3.5678632280087337`*^9, 3.567863228336334*^9}, {3.5678632621727934`*^9, 
   3.5678633509213495`*^9}, {3.5678644293824434`*^9, 3.567864451908883*^9}, {
   3.567864923404111*^9, 3.56786494544695*^9}, {3.5678650143210707`*^9, 
   3.5678650845367947`*^9}, 3.567869231084096*^9, 3.5678693210962543`*^9, {
   3.567869354027912*^9, 3.5678693684423375`*^9}, {3.5682782542610955`*^9, 
   3.5682782564138994`*^9}, {3.568821613806938*^9, 3.568821623151355*^9}, 
   3.5688732159403386`*^9, 3.5689108676158648`*^9, 3.5689997098393545`*^9, 
   3.5689998170271425`*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"u", "=", 
  RowBox[{"ChemicalData", "\[NoBreak]", "[", 
   RowBox[{"\"\<Ethane\>\"", ",", "\"\<VertexTypes\>\""}], "]"}]}]], "Input",
 CellChangeTimes->{{3.567869377256353*^9, 3.5678694101568108`*^9}, {
   3.5678695244894114`*^9, 3.56786953511303*^9}, 3.56786988243204*^9, 
   3.5678699252853155`*^9, {3.567870069179968*^9, 3.5678700766835814`*^9}, {
   3.5678703216664114`*^9, 3.567870325738019*^9}, 3.5679271251072817`*^9, {
   3.567942840466694*^9, 3.567942854381918*^9}, {3.567942932272855*^9, 
   3.567942950509287*^9}, {3.567943087867528*^9, 3.567943093733138*^9}, 
   3.567943185211699*^9, 3.5679432253505697`*^9, 3.5679433277647495`*^9, {
   3.5679434430645523`*^9, 3.5679434440941534`*^9}, {3.568278260766307*^9, 
   3.5682782624355097`*^9}, 3.5688732159403386`*^9, 3.5689108676158648`*^9, 
   3.5689997098393545`*^9, 3.5689998170271425`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"\<\"C\"\>", ",", "\<\"C\"\>", ",", "\<\"H\"\>", ",", "\<\"H\"\>", 
   ",", "\<\"H\"\>", ",", "\<\"H\"\>", ",", "\<\"H\"\>", ",", "\<\"H\"\>"}], 
  "}"}]], "Output",
 CellChangeTimes->{{3.5690806269467382`*^9, 3.5690806536227846`*^9}, {
   3.5690807604985723`*^9, 3.5690808415719147`*^9}, 3.569080894970809*^9, {
   3.5690809257652626`*^9, 3.5690809600073233`*^9}, {3.569080993750182*^9, 
   3.569081009007009*^9}, 3.5690821724414525`*^9, 3.5690885953159337`*^9, 
   3.569172508409171*^9, 3.569176035281365*^9, 3.5691760995534782`*^9, 
   3.569176242387329*^9, 3.569176470178929*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"n", "=", 
  RowBox[{"Length", "[", "u", "]"}]}]], "Input",
 CellChangeTimes->{{3.568277663590458*^9, 3.568277712652544*^9}, 
   3.5682785838808746`*^9, {3.56880682872497*^9, 3.5688068517038097`*^9}, 
   3.5688069726508226`*^9}],

Cell[BoxData["8"], "Output",
 CellChangeTimes->{{3.5690806269467382`*^9, 3.569080653638385*^9}, {
   3.5690807605141726`*^9, 3.569080841587515*^9}, 3.5690808949864087`*^9, {
   3.569080925780863*^9, 3.569080960022923*^9}, {3.5690809937657824`*^9, 
   3.5690810090226088`*^9}, 3.569082172457052*^9, 3.5690885953159337`*^9, 
   3.5691725084247713`*^9, 3.5691760352969646`*^9, 3.569176099569078*^9, 
   3.569176242402929*^9, 3.5691764701945286`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"un", "=", 
  RowBox[{"Table", "[", 
   RowBox[{
    RowBox[{"If", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"u", "[", 
        RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\<H\>\""}], ",", 
      "\"\<Z\>\"", ",", 
      RowBox[{"u", "[", 
       RowBox[{"[", "i", "]"}], "]"}]}], "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.568807971988578*^9, 3.5688080897687845`*^9}, {
   3.5688081769573374`*^9, 3.5688081865669546`*^9}, {3.568808274317109*^9, 
   3.5688082853151283`*^9}, {3.5688085365847692`*^9, 3.568808578330443*^9}, {
   3.568808822704872*^9, 3.5688088761037655`*^9}, {3.56880890727262*^9, 
   3.568808937365073*^9}, {3.568809097608555*^9, 3.5688091594002633`*^9}, 
   3.5688092418308077`*^9, {3.5688092753552666`*^9, 3.568809302577315*^9}, {
   3.5688096257786827`*^9, 3.568809703201618*^9}, {3.5688097456336927`*^9, 
   3.568809788268568*^9}, 3.568809850575077*^9, {3.568809946562046*^9, 
   3.5688099535664577`*^9}, {3.5688100003197403`*^9, 3.568810008556555*^9}, {
   3.568810075371472*^9, 3.56881007981748*^9}, {3.568813431920168*^9, 
   3.568813533086345*^9}, {3.5688135960324554`*^9, 3.5688136045344706`*^9}, {
   3.5688136407889347`*^9, 3.568813707494652*^9}, {3.568813812529636*^9, 
   3.5688138275524626`*^9}, {3.5688140611940727`*^9, 3.568814062738476*^9}, {
   3.568814121113778*^9, 3.5688141965555105`*^9}, {3.568814326675339*^9, 
   3.568814331074547*^9}, 3.568814589473401*^9, 3.5688146762563534`*^9, 
   3.56881472549004*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"\<\"C\"\>", ",", "\<\"C\"\>", ",", "\<\"Z\"\>", ",", "\<\"Z\"\>", 
   ",", "\<\"Z\"\>", ",", "\<\"Z\"\>", ",", "\<\"Z\"\>", ",", "\<\"Z\"\>"}], 
  "}"}]], "Output",
 CellChangeTimes->{{3.569080626962338*^9, 3.569080653653985*^9}, {
   3.5690807605297728`*^9, 3.569080841603115*^9}, 3.569080895002009*^9, {
   3.569080925796463*^9, 3.5690809600385227`*^9}, {3.5690809937813826`*^9, 
   3.5690810090226088`*^9}, 3.5690821724726524`*^9, 3.5690885953315334`*^9, 
   3.5691725084247713`*^9, 3.569176035312565*^9, 3.5691760995846777`*^9, 
   3.5691762424185286`*^9, 3.569176470210129*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"num", "=", 
   RowBox[{"Table", "[", 
    RowBox[{"i", ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.568807971988578*^9, 3.5688080897687845`*^9}, {
   3.5688081769573374`*^9, 3.5688081865669546`*^9}, {3.568808274317109*^9, 
   3.5688082853151283`*^9}, {3.5688085365847692`*^9, 3.568808578330443*^9}, {
   3.568808822704872*^9, 3.5688088761037655`*^9}, {3.56880890727262*^9, 
   3.568808937365073*^9}, {3.568809097608555*^9, 3.5688091594002633`*^9}, 
   3.5688092418308077`*^9, {3.5688092753552666`*^9, 3.568809302577315*^9}, {
   3.5688096257786827`*^9, 3.568809703201618*^9}, {3.5688097456336927`*^9, 
   3.568809788268568*^9}, 3.568809850575077*^9, {3.568809946562046*^9, 
   3.5688099535664577`*^9}, {3.5688100003197403`*^9, 3.568810008556555*^9}, {
   3.568810075371472*^9, 3.56881007981748*^9}, {3.568813431920168*^9, 
   3.568813533086345*^9}, {3.5688135960324554`*^9, 3.5688136045344706`*^9}, {
   3.5688136407889347`*^9, 3.568813707494652*^9}, {3.568813812529636*^9, 
   3.5688138275524626`*^9}, {3.5688140611940727`*^9, 3.568814062738476*^9}, {
   3.568814121113778*^9, 3.5688141965555105`*^9}, {3.568814326675339*^9, 
   3.568814331074547*^9}, 3.568814589473401*^9, {3.5688146762563534`*^9, 
   3.5688146829331646`*^9}, 3.5688147231032352`*^9, 3.5688216336189733`*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"um", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"un", "[", 
        RowBox[{"[", "i", "]"}], "]"}], ",", 
       RowBox[{"num", "[", 
        RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.568807971988578*^9, 3.5688080897687845`*^9}, {
   3.5688081769573374`*^9, 3.5688081865669546`*^9}, {3.568808274317109*^9, 
   3.5688082853151283`*^9}, {3.5688085365847692`*^9, 3.568808578330443*^9}, {
   3.568808822704872*^9, 3.5688088761037655`*^9}, {3.56880890727262*^9, 
   3.568808937365073*^9}, {3.568809097608555*^9, 3.5688091594002633`*^9}, 
   3.5688092418308077`*^9, {3.5688092753552666`*^9, 3.568809302577315*^9}, {
   3.5688096257786827`*^9, 3.568809703201618*^9}, {3.5688097456336927`*^9, 
   3.568809788268568*^9}, 3.568809850575077*^9, {3.568809946562046*^9, 
   3.5688099535664577`*^9}, {3.5688100003197403`*^9, 3.568810008556555*^9}, {
   3.568810075371472*^9, 3.56881007981748*^9}, {3.568813431920168*^9, 
   3.568813533086345*^9}, {3.5688135960324554`*^9, 3.5688136045344706`*^9}, {
   3.5688136407889347`*^9, 3.568813707494652*^9}, {3.568813812529636*^9, 
   3.5688138275524626`*^9}, {3.5688140611940727`*^9, 3.568814062738476*^9}, {
   3.568814121113778*^9, 3.5688141965555105`*^9}, {3.568814326675339*^9, 
   3.568814331074547*^9}, 3.568814589473401*^9, {3.5688146762563534`*^9, 
   3.5688146885023746`*^9}, 3.5688147203576307`*^9, 3.56882164316619*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"us", "=", 
  RowBox[{"Sort", "[", "um", "]"}]}]], "Input",
 CellChangeTimes->{{3.568807971988578*^9, 3.5688080897687845`*^9}, {
   3.5688081769573374`*^9, 3.5688081865669546`*^9}, {3.568808274317109*^9, 
   3.5688082853151283`*^9}, {3.5688085365847692`*^9, 3.568808578330443*^9}, {
   3.568808822704872*^9, 3.5688088761037655`*^9}, {3.56880890727262*^9, 
   3.568808937365073*^9}, {3.568809097608555*^9, 3.5688091594002633`*^9}, 
   3.5688092418308077`*^9, {3.5688092753552666`*^9, 3.568809302577315*^9}, {
   3.5688096257786827`*^9, 3.568809703201618*^9}, {3.5688097456336927`*^9, 
   3.568809788268568*^9}, 3.568809850575077*^9, {3.568809946562046*^9, 
   3.5688099535664577`*^9}, {3.5688100003197403`*^9, 3.568810008556555*^9}, {
   3.568810075371472*^9, 3.56881007981748*^9}, {3.568813431920168*^9, 
   3.568813533086345*^9}, {3.5688135960324554`*^9, 3.5688136045344706`*^9}, {
   3.5688136407889347`*^9, 3.568813707494652*^9}, {3.568813812529636*^9, 
   3.5688138275524626`*^9}, {3.5688140611940727`*^9, 3.568814062738476*^9}, {
   3.568814121113778*^9, 3.5688141965555105`*^9}, {3.568814326675339*^9, 
   3.568814331074547*^9}, 3.568814589473401*^9, {3.5688146762563534`*^9, 
   3.5688146945707855`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"\<\"C\"\>", ",", "1"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"\<\"C\"\>", ",", "2"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"\<\"Z\"\>", ",", "3"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"\<\"Z\"\>", ",", "4"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"\<\"Z\"\>", ",", "5"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"\<\"Z\"\>", ",", "6"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"\<\"Z\"\>", ",", "7"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"\<\"Z\"\>", ",", "8"}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{{3.5690806269779377`*^9, 3.5690806536695848`*^9}, {
   3.5690807605453725`*^9, 3.5690808416187153`*^9}, 3.5690808950176086`*^9, {
   3.5690809258120627`*^9, 3.569080960054123*^9}, {3.5690809937969823`*^9, 
   3.569081009053809*^9}, 3.5690821724882526`*^9, 3.569088595362734*^9, 
   3.5691725084403715`*^9, 3.569176035328165*^9, 3.569176099600278*^9, 
   3.5691762424341288`*^9, 3.569176470225729*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"num", "=", 
  RowBox[{"Table", "[", 
   RowBox[{
    RowBox[{"us", "[", 
     RowBox[{"[", 
      RowBox[{"i", ",", "2"}], "]"}], "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.56881016518083*^9, 3.5688102745058217`*^9}, {
   3.5688103430679426`*^9, 3.5688103543155622`*^9}, {3.5688124958093233`*^9, 
   3.568812496386524*^9}, {3.5688126345871673`*^9, 3.568812635491969*^9}, {
   3.568812673025635*^9, 3.5688126732128353`*^9}, {3.5688127492161684`*^9, 
   3.568812848541543*^9}, 3.5688129280080824`*^9, {3.568813062152718*^9, 
   3.5688131360032477`*^9}, 3.568813172678912*^9, 3.568821701869093*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "1", ",", "2", ",", "3", ",", "4", ",", "5", ",", "6", ",", "7", ",", "8"}],
   "}"}]], "Output",
 CellChangeTimes->{{3.5690806269779377`*^9, 3.5690806536851845`*^9}, {
   3.5690807605609727`*^9, 3.569080841634315*^9}, 3.569080895033209*^9, {
   3.5690809258276625`*^9, 3.569080960069723*^9}, {3.5690809938125825`*^9, 
   3.569081009053809*^9}, 3.569082172503853*^9, 3.569088595362734*^9, 
   3.5691725084559712`*^9, 3.5691760353437653`*^9, 3.569176099615878*^9, 
   3.569176242449729*^9, 3.569176470241329*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"For", "[", 
   RowBox[{
    RowBox[{"i", "=", "1"}], ",", 
    RowBox[{"i", "<", 
     RowBox[{"n", "+", "1"}]}], ",", 
    RowBox[{"i", "++"}], ",", 
    RowBox[{"If", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"us", "[", 
        RowBox[{"[", 
         RowBox[{"i", ",", "1"}], "]"}], "]"}], "\[Equal]", "\"\<Z\>\""}], 
      ",", 
      RowBox[{
       RowBox[{"us", "[", 
        RowBox[{"[", 
         RowBox[{"i", ",", "1"}], "]"}], "]"}], "=", "\"\<H\>\""}]}], "]"}]}],
    "]"}], ";"}]], "Input",
 CellChangeTimes->{{3.568277663590458*^9, 3.568277712652544*^9}, 
   3.5682785838808746`*^9, {3.56880682872497*^9, 3.568806862015428*^9}, 
   3.5688069650848093`*^9, 3.568807273668951*^9, 3.5688073087846127`*^9, {
   3.568807384553946*^9, 3.5688073893119545`*^9}, {3.56881489654434*^9, 
   3.5688149048591547`*^9}, {3.568821714349115*^9, 3.568821718951123*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"u1", "=", 
  RowBox[{"Table", "[", 
   RowBox[{
    RowBox[{"If", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"us", "[", 
        RowBox[{"[", 
         RowBox[{"i", ",", "1"}], "]"}], "]"}], "\[Equal]", "\"\<C\>\""}], 
      ",", "0.255", ",", 
      RowBox[{"If", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"us", "[", 
          RowBox[{"[", 
           RowBox[{"i", ",", "1"}], "]"}], "]"}], "\[Equal]", "\"\<H\>\""}], 
        ",", "0.1", ",", 
        RowBox[{"If", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"us", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "1"}], "]"}], "]"}], "\[Equal]", "\"\<N\>\""}],
           ",", "0.24", ",", 
          RowBox[{"If", "[", 
           RowBox[{
            RowBox[{
             RowBox[{"us", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "1"}], "]"}], "]"}], "\[Equal]", 
             "\"\<O\>\""}], ",", "0.178"}], "]"}]}], "]"}]}], "]"}]}], "]"}], 
    ",", 
    RowBox[{"{", 
     RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], 
   "]"}]}], "\[IndentingNewLine]", 
 RowBox[{"u2", "=", 
  RowBox[{"Table", "[", 
   RowBox[{
    RowBox[{"If", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"us", "[", 
        RowBox[{"[", 
         RowBox[{"i", ",", "1"}], "]"}], "]"}], "\[Equal]", "\"\<C\>\""}], 
      ",", "0.254", ",", 
      RowBox[{"If", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"us", "[", 
          RowBox[{"[", 
           RowBox[{"i", ",", "1"}], "]"}], "]"}], "\[Equal]", "\"\<H\>\""}], 
        ",", "0.6", ",", 
        RowBox[{"If", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"us", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "1"}], "]"}], "]"}], "\[Equal]", "\"\<N\>\""}],
           ",", "0.22522", ",", 
          RowBox[{"If", "[", 
           RowBox[{
            RowBox[{
             RowBox[{"us", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "1"}], "]"}], "]"}], "\[Equal]", 
             "\"\<O\>\""}], ",", "0.178"}], "]"}]}], "]"}]}], "]"}]}], "]"}], 
    ",", 
    RowBox[{"{", 
     RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}]}], "Input",
 CellChangeTimes->{{3.5679480369450207`*^9, 3.567948069112277*^9}, {
   3.5679486805553513`*^9, 3.5679487948723516`*^9}, 3.567950206581231*^9, {
   3.5681216077716007`*^9, 3.5681216839309344`*^9}, {3.568266380090639*^9, 
   3.5682664108850937`*^9}, {3.568266481350417*^9, 3.56826649971165*^9}, 
   3.5682795170433135`*^9, {3.5682795636873956`*^9, 3.5682797843809834`*^9}, {
   3.568279837187076*^9, 3.5682798459854913`*^9}, {3.568279877201146*^9, 
   3.568279880711152*^9}, {3.5682799570264864`*^9, 3.5682799599124913`*^9}, {
   3.5683770044498715`*^9, 3.5683771523693314`*^9}, 3.5688133425788107`*^9, {
   3.568814931192001*^9, 3.5688150075073347`*^9}, {3.5688765350957685`*^9, 
   3.5688765981198792`*^9}, {3.568876634452343*^9, 3.5688766352167444`*^9}, {
   3.568894358551657*^9, 3.5688944054453397`*^9}, {3.568894522929146*^9, 
   3.5688945259087515`*^9}, {3.5688946101644993`*^9, 
   3.5688946171221113`*^9}, {3.568894651067771*^9, 3.568894651956972*^9}, 
   3.5688946827982264`*^9, {3.568894734527918*^9, 3.568894735495119*^9}, {
   3.5688947726699843`*^9, 3.568894773496786*^9}, {3.5688948078948464`*^9, 
   3.568895097914956*^9}, {3.5689558039020233`*^9, 3.568955804619624*^9}, {
   3.5689559169086213`*^9, 3.5689559179070234`*^9}, {3.5689559496998787`*^9, 
   3.568955950589081*^9}, {3.568956090817727*^9, 3.568956118273775*^9}, 
   3.5689561556982408`*^9, {3.5690796452994137`*^9, 3.5690796760158677`*^9}, {
   3.569079718213942*^9, 3.5690797405219812`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "0.255`", ",", "0.255`", ",", "0.1`", ",", "0.1`", ",", "0.1`", ",", "0.1`",
    ",", "0.1`", ",", "0.1`"}], "}"}]], "Output",
 CellChangeTimes->{{3.5690806269779377`*^9, 3.5690806537007847`*^9}, {
   3.569080760576573*^9, 3.569080841649915*^9}, 3.569080895048809*^9, {
   3.5690809258432627`*^9, 3.5690809600853233`*^9}, {3.569080993828182*^9, 
   3.569081009085009*^9}, 3.5690821725194526`*^9, 3.569088595393934*^9, 
   3.569172508471571*^9, 3.5691760353593655`*^9, 3.5691760996314783`*^9, 
   3.5691762424653287`*^9, 3.569176470256929*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "0.254`", ",", "0.254`", ",", "0.6`", ",", "0.6`", ",", "0.6`", ",", "0.6`",
    ",", "0.6`", ",", "0.6`"}], "}"}]], "Output",
 CellChangeTimes->{{3.5690806269779377`*^9, 3.5690806537007847`*^9}, {
   3.569080760576573*^9, 3.569080841649915*^9}, 3.569080895048809*^9, {
   3.5690809258432627`*^9, 3.5690809600853233`*^9}, {3.569080993828182*^9, 
   3.569081009085009*^9}, 3.5690821725194526`*^9, 3.569088595393934*^9, 
   3.569172508471571*^9, 3.5691760353593655`*^9, 3.5691760996314783`*^9, 
   3.5691762424653287`*^9, 3.569176470256929*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"u", "=", 
  RowBox[{"Table", "[", 
   RowBox[{
    RowBox[{"us", "[", 
     RowBox[{"[", 
      RowBox[{"i", ",", "1"}], "]"}], "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.5679472243083935`*^9, 3.567947394894693*^9}, {
  3.567947453425996*^9, 3.567947471241227*^9}, {3.5679475415037503`*^9, 
  3.5679475448733563`*^9}, {3.5679488138107853`*^9, 3.567948820877598*^9}, {
  3.568820296244224*^9, 3.5688203477867146`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"\<\"C\"\>", ",", "\<\"C\"\>", ",", "\<\"H\"\>", ",", "\<\"H\"\>", 
   ",", "\<\"H\"\>", ",", "\<\"H\"\>", ",", "\<\"H\"\>", ",", "\<\"H\"\>"}], 
  "}"}]], "Output",
 CellChangeTimes->{{3.569080626993538*^9, 3.569080653716385*^9}, {
   3.5690807605921726`*^9, 3.569080841665515*^9}, 3.569080895064409*^9, {
   3.569080925858863*^9, 3.5690809601009235`*^9}, {3.5690809938437824`*^9, 
   3.5690810091006093`*^9}, 3.5690821725350523`*^9, 3.569088595393934*^9, 
   3.569172508487171*^9, 3.569176035374965*^9, 3.569176099647078*^9, 
   3.569176242480929*^9, 3.569176470256929*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"cn", "=", 
  RowBox[{"Table", "[", 
   RowBox[{
    RowBox[{"If", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"u", "[", 
        RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\<C\>\""}], ",", "6.0",
       ",", 
      RowBox[{"If", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"u", "[", 
          RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\<H\>\""}], ",", 
        "1.0", ",", 
        RowBox[{"If", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"u", "[", 
            RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\<N\>\""}], ",", 
          "7.0", ",", 
          RowBox[{"If", "[", 
           RowBox[{
            RowBox[{
             RowBox[{"u", "[", 
              RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\<O\>\""}], ",", 
            "8.0"}], "]"}]}], "]"}]}], "]"}]}], "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.5688224161787477`*^9, 3.5688224173799496`*^9}, {
  3.5688224905908785`*^9, 3.568822500434496*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "6.`", ",", "6.`", ",", "1.`", ",", "1.`", ",", "1.`", ",", "1.`", ",", 
   "1.`", ",", "1.`"}], "}"}]], "Output",
 CellChangeTimes->{{3.569080626993538*^9, 3.569080653731985*^9}, {
   3.569080760607773*^9, 3.569080841681115*^9}, 3.5690808950800085`*^9, {
   3.569080925874463*^9, 3.5690809601165233`*^9}, {3.5690809938593826`*^9, 
   3.569081009116209*^9}, 3.5690821725506525`*^9, 3.5690885954095335`*^9, 
   3.5691725085027714`*^9, 3.569176035390565*^9, 3.569176099662678*^9, 
   3.5691762424965286`*^9, 3.5691764702725286`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"bc", "=", "0"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"For", "[", 
   RowBox[{
    RowBox[{"i", "=", "1"}], ",", " ", 
    RowBox[{"i", "<", 
     RowBox[{"n", "+", "1"}]}], ",", 
    RowBox[{"i", "++"}], ",", 
    RowBox[{"If", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"u", "[", 
        RowBox[{"[", "i", "]"}], "]"}], "\[NotEqual]", "\"\<H\>\""}], ",", 
      " ", 
      RowBox[{"bc", "=", 
       RowBox[{"bc", "+", "1"}]}]}], "]"}]}], "]"}], 
  ";"}], "\[IndentingNewLine]", "bc", "\[IndentingNewLine]", 
 RowBox[{"bh", "=", 
  RowBox[{"n", "-", "bc"}]}]}], "Input",
 CellChangeTimes->{{3.5688281182607965`*^9, 3.5688281434392405`*^9}}],

Cell[BoxData["2"], "Output",
 CellChangeTimes->{{3.569080627009138*^9, 3.5690806537475853`*^9}, {
   3.569080760623373*^9, 3.5690808416967154`*^9}, 3.5690808950956087`*^9, {
   3.5690809258900633`*^9, 3.569080960132123*^9}, {3.5690809938749824`*^9, 
   3.569081009116209*^9}, 3.5690821725662527`*^9, 3.5690885954251337`*^9, 
   3.5691725085027714`*^9, 3.569176035406165*^9, 3.569176099678278*^9, 
   3.569176242512129*^9, 3.569176470288129*^9}],

Cell[BoxData["6"], "Output",
 CellChangeTimes->{{3.569080627009138*^9, 3.5690806537475853`*^9}, {
   3.569080760623373*^9, 3.5690808416967154`*^9}, 3.5690808950956087`*^9, {
   3.5690809258900633`*^9, 3.569080960132123*^9}, {3.5690809938749824`*^9, 
   3.569081009116209*^9}, 3.5690821725662527`*^9, 3.5690885954251337`*^9, 
   3.5691725085027714`*^9, 3.569176035406165*^9, 3.569176099678278*^9, 
   3.569176242512129*^9, 3.569176470288129*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"DeleteFile", "[", " ", "\"\<C:\\\\Kimb_Mathem\\\\tmp3\>\"", "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"For", "[", 
   RowBox[{
    RowBox[{"i", "=", "1"}], ",", 
    RowBox[{"i", "<", 
     RowBox[{"n", "+", "1"}]}], ",", 
    RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"u", "[", 
        RowBox[{"[", "i", "]"}], "]"}], ",", 
       RowBox[{"Extract", "[", 
        RowBox[{"t1", ",", 
         RowBox[{"{", "i", "}"}]}], "]"}]}], "}"}], ">>>", " ", 
     "\"\<C:\\\\Kimb_Mathem\\\\tmp3\>\""}]}], "]"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{"ChemicalData", "\[NoBreak]", "[", 
  RowBox[{"\"\<Ethane\>\"", ",", "\"\<MolecularFormulaString\>\""}], 
  "]"}]}], "Input",
 CellChangeTimes->{{3.567869377256353*^9, 3.5678694101568108`*^9}, {
   3.5678695244894114`*^9, 3.56786953511303*^9}, 3.56786988243204*^9, 
   3.5678699252853155`*^9, {3.567870069179968*^9, 3.5678700766835814`*^9}, {
   3.5678703216664114`*^9, 3.567870325738019*^9}, 3.5679271251072817`*^9, {
   3.567942840466694*^9, 3.567942854381918*^9}, {3.567942932272855*^9, 
   3.567942950509287*^9}, 3.5679430681802936`*^9, {3.568278272060727*^9, 
   3.56827827391713*^9}, 3.5688732159559383`*^9, 3.568910867631465*^9, 
   3.5689997098549547`*^9, 3.5689998170427427`*^9, {3.569059689202113*^9, 
   3.569059690730915*^9}}],

Cell[BoxData["\<\"CH3CH3\"\>"], "Output",
 CellChangeTimes->{{3.5690806270247383`*^9, 3.5690806537631855`*^9}, {
   3.5690807606389728`*^9, 3.5690808417279153`*^9}, 3.569080895111209*^9, {
   3.5690809259056635`*^9, 3.569080960147723*^9}, {3.569080993890582*^9, 
   3.5690810091474094`*^9}, 3.569082172581853*^9, 3.569088595440734*^9, 
   3.5691725085183716`*^9, 3.5691760354217653`*^9, 3.569176099693878*^9, 
   3.569176242527729*^9, 3.569176470303729*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"t", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"t1", "[", 
      RowBox[{"[", 
       RowBox[{"num", "[", 
        RowBox[{"[", "i", "]"}], "]"}], "]"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.567869377256353*^9, 3.567869394073182*^9}, {
   3.5678694248832364`*^9, 3.5678694526200852`*^9}, {3.5678694891241493`*^9, 
   3.5678694955045605`*^9}, {3.567869596046737*^9, 3.567869602645549*^9}, 
   3.56786971175214*^9, {3.5678697908910794`*^9, 3.567869798800293*^9}, {
   3.5679434548425727`*^9, 3.5679434621121855`*^9}, 3.568011875875622*^9, 
   3.56811435956207*^9, {3.5682782821071444`*^9, 3.5682782839479475`*^9}, {
   3.5688205455482616`*^9, 3.5688206262784038`*^9}, 3.5688208026055136`*^9, 
   3.5688732159559383`*^9, 3.568910867631465*^9, {3.5689830694736805`*^9, 
   3.568983072687286*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"Plot1", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"Graphics3D", "[", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"If", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"u", "[", 
            RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\<N\>\""}], ",", 
          RowBox[{"Lighter", "[", "Blue", "]"}], ",", 
          RowBox[{"If", "[", 
           RowBox[{
            RowBox[{
             RowBox[{"u", "[", 
              RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\<O\>\""}], ",", 
            RowBox[{"Lighter", "[", "Red", "]"}], ",", 
            RowBox[{"Lighter", "[", "Brown", "]"}]}], "]"}]}], "]"}], ",", 
        RowBox[{"Sphere", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"t", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "1"}], "]"}], "]"}], ",", 
            RowBox[{"t", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "2"}], "]"}], "]"}], ",", 
            RowBox[{"t", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "3"}], "]"}], "]"}]}], "}"}], ",", 
          RowBox[{"u1", "[", 
           RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}], "}"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.5679440895764875`*^9, 3.5679441000441055`*^9}, {
   3.5679441319773617`*^9, 3.5679441523041973`*^9}, {3.567944240210352*^9, 
   3.5679442813476243`*^9}, {3.567944323966899*^9, 3.567944345931738*^9}, {
   3.567944397567828*^9, 3.567944456941533*^9}, {3.567944491557994*^9, 
   3.5679445468444905`*^9}, {3.5679446505690727`*^9, 3.567944683875131*^9}, {
   3.567944736244423*^9, 3.5679448177233667`*^9}, {3.567945549707852*^9, 
   3.5679455499730525`*^9}, {3.5679456039959474`*^9, 3.567945606195551*^9}, {
   3.567945704460124*^9, 3.5679457164877453`*^9}, {3.5679457857362666`*^9, 
   3.5679458050023003`*^9}, {3.567946353670864*^9, 3.5679464127169676`*^9}, {
   3.5679466013368993`*^9, 3.5679466311173515`*^9}, {3.5679467590999765`*^9, 
   3.567946915739851*^9}, {3.5679469466123056`*^9, 3.567946946815106*^9}, {
   3.5679471385082426`*^9, 3.5679471497870626`*^9}, {3.5679483907224417`*^9, 
   3.5679484036080647`*^9}, {3.567948441937332*^9, 3.5679484421869326`*^9}, {
   3.567948496740228*^9, 3.56794849772303*^9}, {3.567948847912445*^9, 
   3.5679488949933276`*^9}, {3.568011728877214*^9, 3.5680117308213253`*^9}, {
   3.568011766325356*^9, 3.568011773436763*^9}, {3.5681173314452896`*^9, 
   3.5681173370612993`*^9}, {3.568117919800323*^9, 3.5681179678952074`*^9}, {
   3.56811809459863*^9, 3.568118105409449*^9}, {3.5681181818963833`*^9, 
   3.5681182389144835`*^9}, {3.568118274888147*^9, 3.5681182883977704`*^9}, {
   3.568118518513775*^9, 3.56811852739019*^9}, {3.5681188250543127`*^9, 
   3.568118862634779*^9}, {3.568118972396572*^9, 3.5681189862493963`*^9}, {
   3.568119020694257*^9, 3.568119047136303*^9}, {3.5681190903015785`*^9, 
   3.568119136930061*^9}, {3.5681191671005135`*^9, 3.568119202808976*^9}, {
   3.5681194788046613`*^9, 3.568119487603077*^9}, {3.5681195353235607`*^9, 
   3.5681195410331707`*^9}, 3.56811959201406*^9, {3.5681196426985493`*^9, 
   3.568119659702579*^9}, {3.568119726813897*^9, 3.56811974027672*^9}, {
   3.5682659184078283`*^9, 3.5682659267694435`*^9}, {3.568265972477524*^9, 
   3.568266002725977*^9}, 3.568277802649102*^9, {3.568820898608082*^9, 
   3.568820922429324*^9}, 3.5688217528499823`*^9, 3.568821789307246*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"Plot2", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"Graphics3D", "[", 
      RowBox[{"Text", "[", 
       RowBox[{"i", ",", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{
           RowBox[{"t", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "1"}], "]"}], "]"}], "-", "0.3"}], ",", 
          RowBox[{"t", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "2"}], "]"}], "]"}], ",", 
          RowBox[{"t", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "3"}], "]"}], "]"}]}], "}"}], " ", ",", 
        RowBox[{"{", 
         RowBox[{"1", ",", "0"}], "}"}]}], "]"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.5679440895764875`*^9, 3.5679441000441055`*^9}, {
   3.5679441319773617`*^9, 3.5679441523041973`*^9}, {3.567944240210352*^9, 
   3.5679442813476243`*^9}, {3.567944323966899*^9, 3.567944345931738*^9}, {
   3.567944397567828*^9, 3.567944456941533*^9}, {3.567944491557994*^9, 
   3.5679445468444905`*^9}, {3.5679446505690727`*^9, 3.567944683875131*^9}, {
   3.567944736244423*^9, 3.5679448177233667`*^9}, {3.567945549707852*^9, 
   3.5679455499730525`*^9}, {3.5679456039959474`*^9, 3.567945606195551*^9}, {
   3.567945704460124*^9, 3.5679457164877453`*^9}, {3.5679457857362666`*^9, 
   3.5679458050023003`*^9}, {3.567946353670864*^9, 3.5679464127169676`*^9}, {
   3.5679466013368993`*^9, 3.5679466311173515`*^9}, {3.5679467590999765`*^9, 
   3.567946915739851*^9}, {3.5679469466123056`*^9, 3.567946946815106*^9}, {
   3.5679471385082426`*^9, 3.5679471497870626`*^9}, {3.5679483907224417`*^9, 
   3.5679484036080647`*^9}, {3.567948441937332*^9, 3.5679484421869326`*^9}, {
   3.567948496740228*^9, 3.56794849772303*^9}, {3.567948847912445*^9, 
   3.5679488949933276`*^9}, {3.568011728877214*^9, 3.5680117308213253`*^9}, {
   3.568011766325356*^9, 3.568011773436763*^9}, {3.5681173314452896`*^9, 
   3.5681173370612993`*^9}, {3.568117919800323*^9, 3.5681179678952074`*^9}, {
   3.56811809459863*^9, 3.568118105409449*^9}, {3.5681181818963833`*^9, 
   3.5681182389144835`*^9}, {3.568118274888147*^9, 3.5681182883977704`*^9}, {
   3.568118518513775*^9, 3.56811852739019*^9}, {3.5681188250543127`*^9, 
   3.568118862634779*^9}, {3.568118972396572*^9, 3.5681189862493963`*^9}, {
   3.568119020694257*^9, 3.568119047136303*^9}, {3.5681190903015785`*^9, 
   3.568119136930061*^9}, {3.5681191671005135`*^9, 3.568119202808976*^9}, {
   3.5681194788046613`*^9, 3.568119487603077*^9}, {3.5681195353235607`*^9, 
   3.5681195410331707`*^9}, 3.56811959201406*^9, {3.5681196426985493`*^9, 
   3.568119659702579*^9}, {3.568119726813897*^9, 3.56811974027672*^9}, {
   3.5681199008322024`*^9, 3.568119920691037*^9}, {3.5681199536070952`*^9, 
   3.568119958521104*^9}, {3.568120000641178*^9, 3.568120067206495*^9}, {
   3.5681211528280015`*^9, 3.5681211643408213`*^9}, 3.5681213201694956`*^9, {
   3.568121503313817*^9, 3.568121508961027*^9}, {3.568266015439999*^9, 
   3.5682660461096525`*^9}, 3.568277806081108*^9, {3.568381589141924*^9, 
   3.568381589391525*^9}, 3.5688217602599955`*^9, {3.568821800305266*^9, 
   3.568821801397268*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"ChemicalData", "\[NoBreak]", "[", 
  RowBox[{"\"\<Ethane\>\"", ",", "\"\<CHStructureDiagram\>\""}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.5678630319943895`*^9, 3.5678630956269007`*^9}, {
   3.5678632280087337`*^9, 3.567863228336334*^9}, {3.5678632621727934`*^9, 
   3.5678633509213495`*^9}, {3.5678644293824434`*^9, 3.567864451908883*^9}, {
   3.567864923404111*^9, 3.56786494544695*^9}, {3.5678650143210707`*^9, 
   3.5678650845367947`*^9}, {3.567869231084096*^9, 3.5678692412397137`*^9}, 
   3.5678693168686466`*^9, {3.5678704450470285`*^9, 3.5678704621290584`*^9}, 
   3.5682783281584253`*^9, 3.5688732159715385`*^9, 3.568910867647065*^9, 
   3.5689997098549547`*^9, 3.5689998170427427`*^9}],

Cell[BoxData[
 GraphicsBox[{{
    {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[
      StyleBox["\<\"C\"\>",
       StripOnInput->False,
       ScriptSizeMultipliers->0.8,
       ScriptMinSize->8,
       FontSize->Scaled[0.1263477088948787],
       FontWeight->Plain], {200., 0.}], LineBox[{{220., 0.}, {250., 0.}}], 
     {GrayLevel[0], InsetBox[
       StyleBox["\<\"C\"\>",
        StripOnInput->False,
        ScriptSizeMultipliers->0.8,
        ScriptMinSize->8,
        FontSize->Scaled[0.1263477088948787],
        FontWeight->Plain], {300., 0.}], LineBox[{{250., 0.}, {280., 0.}}]}}, 
    {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[
      StyleBox["\<\"C\"\>",
       StripOnInput->False,
       ScriptSizeMultipliers->0.8,
       ScriptMinSize->8,
       FontSize->Scaled[0.1263477088948787],
       FontWeight->Plain], {200., 0.}], LineBox[{{200., 20.}, {200., 31.}}], 
     {GrayLevel[0], InsetBox[
       StyleBox["\<\"H\"\>",
        StripOnInput->False,
        ScriptSizeMultipliers->0.8,
        ScriptMinSize->8,
        FontSize->Scaled[0.1263477088948787],
        FontWeight->Plain], {200., 62.}], 
      LineBox[{{200., 31.}, {200., 42.}}]}}, 
    {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[
      StyleBox["\<\"C\"\>",
       StripOnInput->False,
       ScriptSizeMultipliers->0.8,
       ScriptMinSize->8,
       FontSize->Scaled[0.1263477088948787],
       FontWeight->Plain], {200., 0.}], LineBox[{{180., 0.}, {169., 0.}}], 
     {GrayLevel[0], InsetBox[
       StyleBox["\<\"H\"\>",
        StripOnInput->False,
        ScriptSizeMultipliers->0.8,
        ScriptMinSize->8,
        FontSize->Scaled[0.1263477088948787],
        FontWeight->Plain], {138., 0.}], LineBox[{{169., 0.}, {158., 0.}}]}}, 
    {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[
      StyleBox["\<\"C\"\>",
       StripOnInput->False,
       ScriptSizeMultipliers->0.8,
       ScriptMinSize->8,
       FontSize->Scaled[0.1263477088948787],
       FontWeight->Plain], {200., 0.}], LineBox[{{200., -20.}, {200., -31.}}], 
     {GrayLevel[0], InsetBox[
       StyleBox["\<\"H\"\>",
        StripOnInput->False,
        ScriptSizeMultipliers->0.8,
        ScriptMinSize->8,
        FontSize->Scaled[0.1263477088948787],
        FontWeight->Plain], {200., -62.}], 
      LineBox[{{200., -31.}, {200., -42.}}]}}, 
    {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[
      StyleBox["\<\"C\"\>",
       StripOnInput->False,
       ScriptSizeMultipliers->0.8,
       ScriptMinSize->8,
       FontSize->Scaled[0.1263477088948787],
       FontWeight->Plain], {300., 0.}], LineBox[{{300., -20.}, {300., -31.}}], 
     {GrayLevel[0], InsetBox[
       StyleBox["\<\"H\"\>",
        StripOnInput->False,
        ScriptSizeMultipliers->0.8,
        ScriptMinSize->8,
        FontSize->Scaled[0.1263477088948787],
        FontWeight->Plain], {300., -62.}], 
      LineBox[{{300., -31.}, {300., -42.}}]}}, 
    {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[
      StyleBox["\<\"C\"\>",
       StripOnInput->False,
       ScriptSizeMultipliers->0.8,
       ScriptMinSize->8,
       FontSize->Scaled[0.1263477088948787],
       FontWeight->Plain], {300., 0.}], LineBox[{{320., 0.}, {331., 0.}}], 
     {GrayLevel[0], InsetBox[
       StyleBox["\<\"H\"\>",
        StripOnInput->False,
        ScriptSizeMultipliers->0.8,
        ScriptMinSize->8,
        FontSize->Scaled[0.1263477088948787],
        FontWeight->Plain], {362., 0.}], LineBox[{{331., 0.}, {342., 0.}}]}}, 
    {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[
      StyleBox["\<\"C\"\>",
       StripOnInput->False,
       ScriptSizeMultipliers->0.8,
       ScriptMinSize->8,
       FontSize->Scaled[0.1263477088948787],
       FontWeight->Plain], {300., 0.}], LineBox[{{300., 20.}, {300., 31.}}], 
     {GrayLevel[0], InsetBox[
       StyleBox["\<\"H\"\>",
        StripOnInput->False,
        ScriptSizeMultipliers->0.8,
        ScriptMinSize->8,
        FontSize->Scaled[0.1263477088948787],
        FontWeight->Plain], {300., 62.}], 
      LineBox[{{300., 31.}, {300., 42.}}]}}}, {}, {}},
  ImageSize->{{0., 88.15713572731784}, {0., 58.82697336813459}},
  PlotRange->{{115.5, 380.}, {-92., 84.5}}]], "Output",
 CellChangeTimes->{{3.569080627040338*^9, 3.569080653809985*^9}, {
   3.569080760685773*^9, 3.569080841759115*^9}, 3.5690808951580095`*^9, {
   3.569080925952463*^9, 3.569080960178923*^9}, {3.5690809939373827`*^9, 
   3.569081009178609*^9}, 3.569082172613053*^9, 3.569088595487534*^9, 
   3.5691725085651717`*^9, 3.569176035468565*^9, 3.5691760997406783`*^9, 
   3.569176242574529*^9, 3.5691764703349295`*^9}]
}, Open  ]],

Cell[BoxData[{
 RowBox[{
  RowBox[{"x", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"t", "[", 
      RowBox[{"[", 
       RowBox[{"i", ",", "1"}], "]"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"y", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"t", "[", 
      RowBox[{"[", 
       RowBox[{"i", ",", "2"}], "]"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"z", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"t", "[", 
      RowBox[{"[", 
       RowBox[{"i", ",", "3"}], "]"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.5678630319943895`*^9, 3.5678630956269007`*^9}, {
  3.5678632280087337`*^9, 3.567863228336334*^9}, {3.5678632621727934`*^9, 
  3.5678633509213495`*^9}, {3.5678644293824434`*^9, 3.567864451908883*^9}, {
  3.567864923404111*^9, 3.56786494544695*^9}, {3.5678650143210707`*^9, 
  3.5678650845367947`*^9}, {3.567869231084096*^9, 3.567869278118179*^9}, {
  3.5682778102931156`*^9, 3.5682778156907253`*^9}}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"(*", " ", 
    RowBox[{"Distance", " ", "matrix", " ", "Bohr"}], " ", "*)"}], "Text"], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"Dij", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"If", "[", 
        RowBox[{
         RowBox[{"i", "\[GreaterEqual]", "j"}], ",", 
         RowBox[{"If", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"sq", "=", 
             RowBox[{"Sqrt", "[", 
              RowBox[{
               RowBox[{
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{"x", "[", 
                   RowBox[{"[", "i", "]"}], "]"}], "-", 
                  RowBox[{"x", "[", 
                   RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}], "+", 
               RowBox[{
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{"y", "[", 
                   RowBox[{"[", "i", "]"}], "]"}], "-", 
                  RowBox[{"y", "[", 
                   RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}], "+", 
               RowBox[{
                RowBox[{"(", 
                 RowBox[{
                  RowBox[{"z", "[", 
                   RowBox[{"[", "i", "]"}], "]"}], "-", 
                  RowBox[{"z", "[", 
                   RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}]}], 
              "]"}]}], ";", " ", 
            RowBox[{"sq", "<", "3.0"}]}], ",", "sq", ",", "0"}], "]"}], ",", 
         "0"}], "]"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "n"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"j", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{"MatrixForm", "[", "Dij", "]"}], "*)"}]}]}]], "Input",
 CellChangeTimes->{{3.567863535157673*^9, 3.5678635412104836`*^9}, {
   3.5678638182825704`*^9, 3.567863827283786*^9}, {3.5679503910983553`*^9, 
   3.5679504266352177`*^9}, {3.5679504653700857`*^9, 3.5679505048069553`*^9}, 
   3.567950534915008*^9, {3.567951107623214*^9, 3.567951136046464*^9}, 
   3.567951302545556*^9, {3.5679562986474605`*^9, 3.567956305355472*^9}, {
   3.5679563590975666`*^9, 3.567956370688387*^9}, {3.567956414930065*^9, 
   3.567956430467692*^9}, {3.567956484740187*^9, 3.567956500667815*^9}, {
   3.5679565456270943`*^9, 3.5679565868267665`*^9}, {3.5679569267045636`*^9, 
   3.5679569448941956`*^9}, {3.567956976562251*^9, 3.567957000711094*^9}, {
   3.5679570348595533`*^9, 3.5679570632672033`*^9}, {3.5679573263928657`*^9, 
   3.5679573286548696`*^9}, {3.5679573611029263`*^9, 
   3.5679573637393312`*^9}, {3.5679592820903006`*^9, 3.567959317127962*^9}, {
   3.567959397655304*^9, 3.5679594004945087`*^9}, {3.5679594704762316`*^9, 
   3.5679594724574347`*^9}, {3.567959655898157*^9, 3.567959673073787*^9}, 
   3.5681143121223865`*^9, 3.568121866388855*^9, 3.568220668221574*^9, {
   3.568220718968463*^9, 3.568220723742071*^9}, {3.568220862051914*^9, 
   3.5682208641111183`*^9}, {3.568220906761593*^9, 3.5682209083527956`*^9}, {
   3.5682209451844606`*^9, 3.5682209651680956`*^9}, {3.568221147782016*^9, 
   3.5682211695908546`*^9}, {3.568221237279373*^9, 3.568221253612602*^9}, {
   3.568221305420293*^9, 3.568221306247094*^9}, {3.5682215161142635`*^9, 
   3.568221546206716*^9}, 3.56822158247678*^9, {3.568277819185131*^9, 
   3.5682778284359474`*^9}, 3.56882183684053*^9, {3.5688219650415554`*^9, 
   3.568821965915157*^9}, {3.5688769979329815`*^9, 3.5688771619984694`*^9}, {
   3.568877204212144*^9, 3.5688772091729527`*^9}, {3.568877282493081*^9, 
   3.5688772843338842`*^9}, {3.568877374112042*^9, 3.568877378074449*^9}, {
   3.56887744125456*^9, 3.5688774664954042`*^9}, {3.56887758366721*^9, 
   3.568877604867647*^9}, {3.568877667220957*^9, 3.5688776710429635`*^9}, {
   3.5688778111156096`*^9, 3.56887782262843*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"For", "[", 
  RowBox[{
   RowBox[{"i", "=", "1"}], " ", ",", 
   RowBox[{"i", "<", 
    RowBox[{"n", "+", "1"}]}], ",", 
   RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
   RowBox[{"For", "[", 
    RowBox[{
     RowBox[{"j", "=", "1"}], ",", " ", 
     RowBox[{"j", "<", "i"}], ",", 
     RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", 
     RowBox[{"If", "[", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{
         RowBox[{
          RowBox[{"Dij", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "j"}], "]"}], "]"}], "\[NotEqual]", "0"}], " ", 
         "&&", " ", 
         RowBox[{
          RowBox[{"u", "[", 
           RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", " ", "\"\<H\>\""}]}], 
        ")"}], ",", 
       RowBox[{"If", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"u", "[", 
           RowBox[{"[", "j", "]"}], "]"}], "\[Equal]", "\"\<C\>\""}], ",", 
         RowBox[{
          RowBox[{"u2", "[", 
           RowBox[{"[", "i", "]"}], "]"}], "=", "0.5693"}], ",", 
         RowBox[{"If", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"u", "[", 
             RowBox[{"[", "j", "]"}], "]"}], "\[Equal]", "\"\<F\>\""}], ",", 
           RowBox[{
            RowBox[{"u2", "[", 
             RowBox[{"[", "i", "]"}], "]"}], "=", "0.5"}], ",", 
           RowBox[{"If", "[", 
            RowBox[{
             RowBox[{
              RowBox[{"u", "[", 
               RowBox[{"[", "j", "]"}], "]"}], "\[Equal]", "\"\<N\>\""}], ",", 
             RowBox[{
              RowBox[{"u2", "[", 
               RowBox[{"[", "i", "]"}], "]"}], "=", "0.5541"}], ",", 
             RowBox[{"If", "[", 
              RowBox[{
               RowBox[{
                RowBox[{"u", "[", 
                 RowBox[{"[", "j", "]"}], "]"}], "\[Equal]", "\"\<O\>\""}], 
               ",", 
               RowBox[{
                RowBox[{"u2", "[", 
                 RowBox[{"[", "i", "]"}], "]"}], "=", "0.5241"}]}], "]"}]}], 
            "]"}]}], "]"}]}], "]"}]}], "]"}]}], "]"}]}], 
  "]"}], "\[IndentingNewLine]", "u2"}], "Input",
 CellChangeTimes->{{3.5688800293771057`*^9, 3.5688801774057655`*^9}, {
   3.5688802290574565`*^9, 3.568880330208034*^9}, {3.5688803807053227`*^9, 
   3.5688804136525807`*^9}, 3.568880491434317*^9, 3.5688805648948464`*^9, {
   3.5688806280437574`*^9, 3.568880655234605*^9}, {3.568880712424306*^9, 
   3.568880716324312*^9}, {3.5688808087544746`*^9, 3.568880827786508*^9}, 
   3.5688808733385878`*^9, {3.568881358936241*^9, 3.5688813959239063`*^9}, {
   3.568881446561595*^9, 3.5688815849182377`*^9}, {3.5688816223739033`*^9, 
   3.568881642607139*^9}, 3.568881728953291*^9, {3.568881797624612*^9, 
   3.568881809870633*^9}, 3.568881863612727*^9, {3.568883098916497*^9, 
   3.5688831014749017`*^9}, {3.56888372124679*^9, 3.568883843239004*^9}, 
   3.568888211750622*^9, {3.5688882521546926`*^9, 3.5688885232519693`*^9}, {
   3.5688885890840845`*^9, 3.5688885895052853`*^9}, {3.568888675773437*^9, 
   3.5688886914514647`*^9}, {3.568888736285943*^9, 3.5688887717760057`*^9}, {
   3.5688888256429005`*^9, 3.5688888541441507`*^9}, {3.568892442649892*^9, 
   3.5688924499507046`*^9}, {3.568892758441247*^9, 3.5688927597360487`*^9}, {
   3.5688928234309607`*^9, 3.5688928379389863`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "0.254`", ",", "0.254`", ",", "0.5693`", ",", "0.5693`", ",", "0.5693`", 
   ",", "0.5693`", ",", "0.5693`", ",", "0.5693`"}], "}"}]], "Output",
 CellChangeTimes->{{3.5690806270559382`*^9, 3.5690806538255854`*^9}, {
   3.569080760701373*^9, 3.569080841790315*^9}, 3.5690808951736093`*^9, {
   3.5690809259680634`*^9, 3.569080960210123*^9}, {3.569080993952983*^9, 
   3.5690810092098093`*^9}, 3.569082172644253*^9, 3.569088595503134*^9, 
   3.5691725085807714`*^9, 3.569176035484165*^9, 3.569176099756278*^9, 
   3.569176242590129*^9, 3.569176470350529*^9}]
}, Open  ]],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"Fij", "=", 
    RowBox[{"Dij", "-", "u2"}]}], ";"}], "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"MatrixForm", "[", "Fij", "]"}], ";"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Gij", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"If", "[", 
      RowBox[{
       RowBox[{"i", "\[GreaterEqual]", "j"}], ",", 
       RowBox[{"If", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"Fij", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "j"}], "]"}], "]"}], ">", "0"}], ",", 
         RowBox[{"Fij", "[", 
          RowBox[{"[", 
           RowBox[{"i", ",", "j"}], "]"}], "]"}], ",", "0"}], "]"}], ",", 
       "\"\< \>\""}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "n"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"j", ",", "1", ",", "n"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"MatrixForm", "[", "Gij", "]"}], ";"}], "*)"}]}]}], "Input",
 CellChangeTimes->{{3.567958596890297*^9, 3.567958603270708*^9}, {
   3.5679587604565845`*^9, 3.5679587707682023`*^9}, {3.5679588272871017`*^9, 
   3.567958828254304*^9}, {3.567959580066424*^9, 3.5679595814704266`*^9}, {
   3.567959624994503*^9, 3.5679596300177116`*^9}, {3.5679597079086485`*^9, 
   3.567959902877791*^9}, {3.5679599961659546`*^9, 3.5679600446196404`*^9}, {
   3.5679609611992497`*^9, 3.567960974006872*^9}, 3.567961975154231*^9, {
   3.5679622637079377`*^9, 3.5679622945179915`*^9}, 3.5681142413606625`*^9, 
   3.568121878697276*^9, {3.5681219332193723`*^9, 3.56812197189184*^9}, {
   3.5682778370003624`*^9, 3.5682778397927675`*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"Hij", "=", 
   RowBox[{"Transpose", "[", "Gij", "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"MatrixForm", "[", "Hij", "]"}], ";"}], "*)"}]}]}], "Input",
 CellChangeTimes->{{3.5679601758938704`*^9, 3.5679601965327067`*^9}, {
   3.567960427147912*^9, 3.567960432217921*^9}, {3.5679605764089737`*^9, 
   3.5679605968294096`*^9}, {3.5679607389456596`*^9, 3.5679607643581038`*^9}, 
   3.5681142501278777`*^9, 3.5681218877140923`*^9}],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"Kij", "=", 
    RowBox[{"Hij", "-", "u2"}]}], ";"}], "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"MatrixForm", "[", "Kij", "]"}], ";"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Lij", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"If", "[", 
      RowBox[{
       RowBox[{"i", "\[LessEqual]", "j"}], ",", 
       RowBox[{"If", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"Kij", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "j"}], "]"}], "]"}], ">", "0"}], ",", 
         RowBox[{"Kij", "[", 
          RowBox[{"[", 
           RowBox[{"i", ",", "j"}], "]"}], "]"}], ",", "0"}], "]"}], ",", 
       "\"\< \>\""}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "n"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"j", ",", "1", ",", "n"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"MatrixForm", "[", "Lij", "]"}], ";"}], "*)"}]}]}], "Input",
 CellChangeTimes->{{3.5679602339571724`*^9, 3.567960257419614*^9}, {
   3.567960343406965*^9, 3.567960353500182*^9}, {3.5679604155258913`*^9, 
   3.567960421422702*^9}, {3.567960854323462*^9, 3.567960874197897*^9}, {
   3.5679611596939983`*^9, 3.567961170270817*^9}, {3.5679614000124207`*^9, 
   3.567961407016833*^9}, {3.5679615066074076`*^9, 3.567961509649413*^9}, 
   3.5679616124379935`*^9, 3.5679616766477065`*^9, {3.5679617424486217`*^9, 
   3.5679617456310277`*^9}, 3.567961801276325*^9, {3.5679618582788253`*^9, 
   3.5679618970760937`*^9}, {3.5679620279135237`*^9, 3.567962030113127*^9}, 
   3.567962107894864*^9, {3.567962149468937*^9, 3.5679621572533503`*^9}, 
   3.5679623389156694`*^9, {3.567962645815009*^9, 3.567962648888214*^9}, {
   3.567962691647889*^9, 3.567962693020692*^9}, {3.568048161486014*^9, 
   3.568048162827616*^9}, 3.5680482096900983`*^9, 3.5681142588638926`*^9, 
   3.568121899632513*^9, 3.5681219958534822`*^9, {3.568277844675576*^9, 
   3.5682778463603787`*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"Lij", "=", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"If", "[", 
       RowBox[{
        RowBox[{"i", "\[LessEqual]", "j"}], ",", 
        RowBox[{"If", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"Lij", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "j"}], "]"}], "]"}], ">", "0"}], ",", 
          RowBox[{"If", "[", 
           RowBox[{
            RowBox[{
             RowBox[{"u", "[", 
              RowBox[{"[", "j", "]"}], "]"}], "\[NotEqual]", "\"\<H\>\""}], 
            ",", 
            RowBox[{"0.5", "*", 
             RowBox[{"Lij", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "j"}], "]"}], "]"}]}], ",", 
            RowBox[{"Lij", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "j"}], "]"}], "]"}]}], "]"}], ",", "0"}], 
         "]"}], ",", "\"\< \>\""}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "n"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"j", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}], 
  "  "}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Qij", "=", "Lij"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"MatrixForm", "[", "Lij", "]"}], ";"}], "*)"}]}]}], "Input",
 CellChangeTimes->{{3.5679627617388124`*^9, 3.567962837461345*^9}, {
   3.567963118121438*^9, 3.567963118698639*^9}, {3.5679631520670977`*^9, 
   3.5679631541731014`*^9}, {3.5679631898191643`*^9, 
   3.5679632154656096`*^9}, {3.5679635823158536`*^9, 3.56796361468591*^9}, {
   3.567963691516045*^9, 3.567963693590849*^9}, {3.567963737957327*^9, 
   3.567963807299449*^9}, {3.568048114904332*^9, 3.5680481161835346`*^9}, 
   3.5680482221389203`*^9, {3.5680485003498087`*^9, 3.5680485123618298`*^9}, 
   3.5681142659931054`*^9, 3.5681219084621286`*^9, {3.568277848107582*^9, 
   3.5682778496519847`*^9}, {3.5683793122399254`*^9, 3.568379313019926*^9}, {
   3.568379456758579*^9, 3.56837947986222*^9}, 3.5683796000760307`*^9, {
   3.5683796473285136`*^9, 3.5683796831617765`*^9}, {3.568381744939398*^9, 
   3.5683817743298492`*^9}, {3.56846523720809*^9, 3.568465240218895*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"Nij", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"If", "[", 
      RowBox[{
       RowBox[{"i", "\[LessEqual]", "j"}], ",", 
       RowBox[{"If", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"Lij", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "j"}], "]"}], "]"}], ">", "0"}], ",", 
         RowBox[{"2.25", "/", 
          RowBox[{
           RowBox[{"Lij", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "j"}], "]"}], "]"}], "^", "2"}]}], ",", "0"}], 
        "]"}], ",", "\"\< \>\""}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "n"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"j", ",", "1", ",", "n"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"MatrixForm", "[", "Nij", "]"}], "*)"}]}]}], "Input",
 CellChangeTimes->{{3.5679628944794455`*^9, 3.5679629710755796`*^9}, {
   3.567963889963994*^9, 3.567963933800071*^9}, {3.5680478391738477`*^9, 
   3.5680478414358516`*^9}, 3.568048231171336*^9, {3.5680482823550262`*^9, 
   3.56804828449223*^9}, 3.5680483349271183`*^9, 3.568114273059918*^9, 
   3.5681219175101447`*^9, {3.568277851305588*^9, 3.568277857295998*^9}, {
   3.5688722635742655`*^9, 3.568872266616271*^9}, 3.5688723032763357`*^9, 
   3.5688725225503206`*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Tval", "=", 
  RowBox[{"Sum", "[", 
   RowBox[{
    RowBox[{"If", "[", 
     RowBox[{
      RowBox[{"i", "\[LessEqual]", "j"}], ",", 
      RowBox[{"Nij", "[", 
       RowBox[{"[", 
        RowBox[{"i", ",", "j"}], "]"}], "]"}], ",", "0"}], "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"i", ",", "1", ",", "n"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"j", ",", "1", ",", "n"}], "}"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.567963992456174*^9, 3.567964003890994*^9}, {
   3.5679640974287586`*^9, 3.5679642038209453`*^9}, {3.56796433180357*^9, 
   3.567964353487608*^9}, {3.567964403641696*^9, 3.567964453046983*^9}, {
   3.5679646160984693`*^9, 3.5679646168316708`*^9}, {3.567964688747797*^9, 
   3.567964762255126*^9}, {3.568277859480002*^9, 3.568277863598409*^9}, {
   3.568379511358675*^9, 3.5683795129030776`*^9}, 3.5688214977895346`*^9, {
   3.56882774467134*^9, 3.56882781892747*^9}}],

Cell[BoxData["9.842230699648079`"], "Output",
 CellChangeTimes->{{3.569080627087138*^9, 3.5690806538567853`*^9}, {
   3.5690807607325726`*^9, 3.569080841805915*^9}, 3.569080895204809*^9, {
   3.5690809259992633`*^9, 3.5690809602413235`*^9}, {3.5690809939841824`*^9, 
   3.5690810092410097`*^9}, 3.5690821726754527`*^9, 3.569088595534334*^9, 
   3.5691725085963717`*^9, 3.5691760355153656`*^9, 3.5691760997874784`*^9, 
   3.5691762426213293`*^9, 3.569176470381729*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Tcor", "=", 
  RowBox[{"Sum", "[", 
   RowBox[{
    RowBox[{"If", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"u", "[", 
        RowBox[{"[", "i", "]"}], "]"}], "\[NotEqual]", "\"\<H\>\""}], ",", 
      RowBox[{"2.25", "/", 
       RowBox[{
        RowBox[{"u2", "[", 
         RowBox[{"[", "i", "]"}], "]"}], "^", "2"}]}], ",", "0"}], "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}]], "Input",
 CellChangeTimes->{
  3.5688724058309155`*^9, {3.568872452646598*^9, 3.568872478480243*^9}, {
   3.5688764461756124`*^9, 3.5688764647552447`*^9}}],

Cell[BoxData["69.750139500279`"], "Output",
 CellChangeTimes->{{3.569080627087138*^9, 3.569080653872385*^9}, {
   3.5690807607481728`*^9, 3.569080841821515*^9}, 3.5690808952204094`*^9, {
   3.569080926014863*^9, 3.5690809602569237`*^9}, {3.5690809939997826`*^9, 
   3.56908100925661*^9}, 3.5690821726754527`*^9, 3.569088595549934*^9, 
   3.569172508611972*^9, 3.5691760355309653`*^9, 3.5691760998030787`*^9, 
   3.5691762426369295`*^9, 3.569176470381729*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Ttot", "=", 
  RowBox[{"Tval", "+", "Tcor", "  ", 
   RowBox[{"(*", " ", 
    RowBox[{"plus", " ", "\[Pi]", " ", "and", " ", "lone", " ", 
     RowBox[{"pairs", ":", " ", 
      RowBox[{"+", " ", "8.7616"}]}]}], " ", "*)"}]}]}]], "Input",
 CellChangeTimes->{
  3.5688724058309155`*^9, {3.568872452646598*^9, 3.568872478480243*^9}, {
   3.5688764461756124`*^9, 3.568876459248435*^9}, {3.568876733122516*^9, 
   3.5688767685033784`*^9}, {3.568895188410715*^9, 3.5688952349923964`*^9}}],

Cell[BoxData["79.59237019992709`"], "Output",
 CellChangeTimes->{{3.569080627087138*^9, 3.569080653887985*^9}, {
   3.569080760763773*^9, 3.5690808418371153`*^9}, 3.569080895236009*^9, {
   3.569080926030463*^9, 3.5690809602569237`*^9}, {3.569080994015383*^9, 
   3.5690810092722096`*^9}, 3.569082172691053*^9, 3.569088595565534*^9, 
   3.5691725086275716`*^9, 3.5691760355465655`*^9, 3.5691760998186784`*^9, 
   3.5691762426525292`*^9, 3.5691764703973293`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"TvalH", "=", 
  RowBox[{"Sum", "[", 
   RowBox[{
    RowBox[{"If", "[", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"i", "\[LessEqual]", "j"}], " ", "&&", " ", 
        RowBox[{
         RowBox[{"u", "[", 
          RowBox[{"[", "j", "]"}], "]"}], "\[Equal]", " ", "\"\<H\>\""}]}], 
       ")"}], ",", 
      RowBox[{"Nij", "[", 
       RowBox[{"[", 
        RowBox[{"i", ",", "j"}], "]"}], "]"}], ",", "0"}], "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"i", ",", "1", ",", "n"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"j", ",", "1", ",", "n"}], "}"}]}], "]"}]}]], "Input",
 CellChangeTimes->{
  3.5688724058309155`*^9, {3.568872452646598*^9, 3.568872478480243*^9}, {
   3.5688764461756124`*^9, 3.568876459248435*^9}, {3.568876733122516*^9, 
   3.5688767642445707`*^9}}],

Cell[BoxData["8.239905289478962`"], "Output",
 CellChangeTimes->{{3.5690806271027384`*^9, 3.5690806539035854`*^9}, {
   3.569080760779373*^9, 3.5690808418527155`*^9}, 3.569080895251609*^9, {
   3.5690809260460634`*^9, 3.5690809602725234`*^9}, {3.569080994030983*^9, 
   3.5690810092722096`*^9}, 3.5690821727066526`*^9, 3.569088595565534*^9, 
   3.5691725086431713`*^9, 3.5691760355621653`*^9, 3.5691760998186784`*^9, 
   3.569176242668129*^9, 3.5691764704129295`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"TvalC", "=", 
  RowBox[{"Tval", "-", "TvalH"}]}]], "Input",
 CellChangeTimes->{{3.5688726088808727`*^9, 3.568872610581275*^9}, {
  3.5688726952426243`*^9, 3.5688727146490583`*^9}}],

Cell[BoxData["1.602325410169117`"], "Output",
 CellChangeTimes->{{3.5690806271027384`*^9, 3.5690806539035854`*^9}, {
   3.569080760779373*^9, 3.5690808418683157`*^9}, 3.569080895267209*^9, {
   3.569080926061663*^9, 3.569080960288123*^9}, {3.5690809940465827`*^9, 
   3.5690810092878094`*^9}, 3.569082172722253*^9, 3.569088595581134*^9, 
   3.5691725086431713`*^9, 3.5691760355777655`*^9, 3.569176099834278*^9, 
   3.569176242683729*^9, 3.5691764704129295`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"Vnn", "=", "0.0"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"For", "[", " ", 
   RowBox[{
    RowBox[{"i", "=", "1"}], ",", 
    RowBox[{"i", "<", "n"}], ",", 
    RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
    RowBox[{"For", "[", 
     RowBox[{
      RowBox[{"j", "=", 
       RowBox[{"i", "+", "1"}]}], ",", 
      RowBox[{"j", "<", 
       RowBox[{"n", "+", "1"}]}], ",", 
      RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", 
      RowBox[{"Vnn", "=", 
       RowBox[{"Vnn", "+", 
        RowBox[{
         RowBox[{"cn", "[", 
          RowBox[{"[", "i", "]"}], "]"}], "*", 
         RowBox[{
          RowBox[{"cn", "[", 
           RowBox[{"[", "j", "]"}], "]"}], "/", 
          RowBox[{"Sqrt", "[", 
           RowBox[{
            RowBox[{
             RowBox[{"(", 
              RowBox[{
               RowBox[{"x", "[", 
                RowBox[{"[", "i", "]"}], "]"}], "-", 
               RowBox[{"x", "[", 
                RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}], "+", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{
               RowBox[{"y", "[", 
                RowBox[{"[", "i", "]"}], "]"}], "-", 
               RowBox[{"y", "[", 
                RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}], "+", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{
               RowBox[{"z", "[", 
                RowBox[{"[", "i", "]"}], "]"}], "-", 
               RowBox[{"z", "[", 
                RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}]}], 
           "]"}]}]}]}]}]}], "]"}]}], "]"}], 
  ";"}], "\[IndentingNewLine]", "Vnn"}], "Input"],

Cell[BoxData["41.956156943012274`"], "Output",
 CellChangeTimes->{{3.5690806271183386`*^9, 3.569080653919185*^9}, {
   3.5690807607949734`*^9, 3.569080841883916*^9}, 3.5690808952828093`*^9, {
   3.5690809260772634`*^9, 3.5690809603037233`*^9}, {3.569080994062183*^9, 
   3.5690810093034096`*^9}, 3.569082172737853*^9, 3.569088595596734*^9, 
   3.5691725086587715`*^9, 3.5691760355933657`*^9, 3.5691760998498783`*^9, 
   3.569176242683729*^9, 3.5691764704285293`*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"Mij", "=", 
     RowBox[{"Prepend", "[", 
      RowBox[{"Qij", ",", "u"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Oij", "=", 
     RowBox[{"Transpose", "[", "Mij", "]"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Pij", "=", 
     RowBox[{"Prepend", "[", 
      RowBox[{"Oij", ",", 
       RowBox[{"Prepend", "[", 
        RowBox[{"u", ",", "\"\< \>\""}], "]"}]}], "]"}]}], ";"}]}]}]], "Input",
 CellChangeTimes->{{3.568045828002715*^9, 3.56804588761042*^9}, 
   3.5680459828797874`*^9, {3.568046033938677*^9, 3.5680460518787084`*^9}, {
   3.568046989190755*^9, 3.5680470284404235`*^9}, {3.568047061028881*^9, 
   3.568047085567724*^9}, {3.5680471562982483`*^9, 3.568047169199471*^9}, {
   3.568047223393966*^9, 3.568047364059413*^9}, {3.568047417286707*^9, 
   3.5680474225127163`*^9}, {3.5680474583459787`*^9, 3.568047487861231*^9}, {
   3.5680475345521126`*^9, 3.5680475468605347`*^9}, {3.5680479168775845`*^9, 
   3.5680479311516094`*^9}, {3.568048066918648*^9, 3.568048067667449*^9}, {
   3.5680483816336*^9, 3.5680483930216208`*^9}, {3.568048424159275*^9, 
   3.568048450211321*^9}, {3.5680486326224413`*^9, 3.5680486370996494`*^9}, {
   3.5680487358166227`*^9, 3.5680487967191296`*^9}, 3.568048834377596*^9, 
   3.568048984418659*^9, {3.5680490255559316`*^9, 3.5680490790796256`*^9}, {
   3.5680491274553103`*^9, 3.5680491280949116`*^9}, {3.5680491782333994`*^9, 
   3.568049231944294*^9}, {3.568049361533722*^9, 3.5680493691933355`*^9}, {
   3.568049422171028*^9, 3.568049427490638*^9}, {3.568049476427923*^9, 
   3.5680495489368505`*^9}, {3.568049620634577*^9, 3.5680496274361887`*^9}, {
   3.5680496655314555`*^9, 3.568049693876705*^9}, 3.568049746573598*^9, {
   3.568049846460573*^9, 3.568049932759925*^9}, {3.568050009059659*^9, 
   3.5680500167348723`*^9}, {3.5680501572755194`*^9, 3.568050179911159*^9}, 
   3.5681087021497335`*^9}],

Cell[BoxData[
 RowBox[{"(*", 
  RowBox[{
   RowBox[{"MatrixForm", "[", "Pij", "]"}], ";"}], "*)"}]], "Input",
 CellChangeTimes->{{3.568045828002715*^9, 3.56804588761042*^9}, 
   3.5680459828797874`*^9, {3.568046033938677*^9, 3.5680460518787084`*^9}, {
   3.568046989190755*^9, 3.5680470284404235`*^9}, {3.568047061028881*^9, 
   3.568047085567724*^9}, {3.5680471562982483`*^9, 3.568047169199471*^9}, {
   3.568047223393966*^9, 3.568047364059413*^9}, {3.568047417286707*^9, 
   3.5680474225127163`*^9}, {3.5680474583459787`*^9, 3.568047487861231*^9}, {
   3.5680475345521126`*^9, 3.5680475468605347`*^9}, {3.5680479168775845`*^9, 
   3.5680479311516094`*^9}, {3.568048066918648*^9, 3.568048067667449*^9}, {
   3.5680483816336*^9, 3.5680483930216208`*^9}, {3.568048424159275*^9, 
   3.568048450211321*^9}, {3.5680486326224413`*^9, 3.5680486370996494`*^9}, {
   3.5680487358166227`*^9, 3.5680487967191296`*^9}, 3.568048834377596*^9, 
   3.568048984418659*^9, {3.5680490255559316`*^9, 3.5680490790796256`*^9}, {
   3.5680491274553103`*^9, 3.5680491280949116`*^9}, {3.5680491782333994`*^9, 
   3.568049231944294*^9}, {3.568049361533722*^9, 3.5680493691933355`*^9}, {
   3.568049422171028*^9, 3.568049427490638*^9}, {3.568049476427923*^9, 
   3.5680495489368505`*^9}, {3.568049620634577*^9, 3.5680496274361887`*^9}, {
   3.5680496655314555`*^9, 3.568049693876705*^9}, 3.568049746573598*^9, {
   3.568049846460573*^9, 3.568049932759925*^9}, {3.568050009059659*^9, 
   3.5680500167348723`*^9}, {3.5680501572755194`*^9, 3.568050179911159*^9}, {
   3.5681087021497335`*^9, 3.5681087223049684`*^9}, 3.568114284401138*^9, 
   3.5681220092695055`*^9}],

Cell[BoxData[
 RowBox[{"(*", " ", 
  RowBox[{"sigma", " ", "bond", " ", "cloud", " ", "radii"}], " ", 
  "*)"}]], "Input",
 CellChangeTimes->{{3.568045828002715*^9, 3.56804588761042*^9}, 
   3.5680459828797874`*^9, {3.568046033938677*^9, 3.5680460518787084`*^9}, {
   3.568046989190755*^9, 3.5680470284404235`*^9}, {3.568047061028881*^9, 
   3.568047085567724*^9}, {3.5680471562982483`*^9, 3.568047169199471*^9}, {
   3.568047223393966*^9, 3.568047364059413*^9}, {3.568047417286707*^9, 
   3.5680474225127163`*^9}, {3.5680474583459787`*^9, 3.568047487861231*^9}, {
   3.5680475345521126`*^9, 3.5680475468605347`*^9}, {3.5680479168775845`*^9, 
   3.5680479311516094`*^9}, {3.568048066918648*^9, 3.568048067667449*^9}, {
   3.5680483816336*^9, 3.5680483930216208`*^9}, {3.568048424159275*^9, 
   3.568048450211321*^9}, {3.5680486326224413`*^9, 3.5680486370996494`*^9}, {
   3.5680487358166227`*^9, 3.5680487967191296`*^9}, 3.568048834377596*^9, 
   3.568048984418659*^9, {3.5680490255559316`*^9, 3.5680490790796256`*^9}, {
   3.5680491274553103`*^9, 3.5680491280949116`*^9}, {3.5680491782333994`*^9, 
   3.568049231944294*^9}, {3.568049361533722*^9, 3.5680493691933355`*^9}, {
   3.568049422171028*^9, 3.568049427490638*^9}, {3.568049476427923*^9, 
   3.5680495489368505`*^9}, {3.568049620634577*^9, 3.5680496274361887`*^9}, {
   3.5680496655314555`*^9, 3.568049693876705*^9}, 3.568049746573598*^9, {
   3.568049846460573*^9, 3.568049932759925*^9}, {3.568050009059659*^9, 
   3.5680500167348723`*^9}, {3.5680501572755194`*^9, 3.568050179911159*^9}, {
   3.5681087021497335`*^9, 3.5681087223049684`*^9}, {3.5681133837179556`*^9, 
   3.5681133888191648`*^9}}],

Cell[BoxData[
 RowBox[{"(*", 
  RowBox[{"For", "[", 
   RowBox[{
    RowBox[{"i", "=", "2"}], ",", 
    RowBox[{"i", "<", 
     RowBox[{"n", "+", "2"}]}], ",", 
    RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
    RowBox[{"For", "[", 
     RowBox[{
      RowBox[{"j", "=", "2"}], ",", 
      RowBox[{"j", "\[LessEqual]", "i"}], ",", 
      RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", 
      RowBox[{"If", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"Pij", "[", 
          RowBox[{"[", 
           RowBox[{"i", ",", "j"}], "]"}], "]"}], ">", "0"}], ",", 
        RowBox[{"Print", "[", 
         RowBox[{
          RowBox[{"u", "[", 
           RowBox[{"[", 
            RowBox[{"i", "-", "1"}], "]"}], "]"}], ",", 
          RowBox[{"i", "-", "1"}], ",", "\"\< \>\"", ",", 
          RowBox[{"u", "[", 
           RowBox[{"[", 
            RowBox[{"j", "-", "1"}], "]"}], "]"}], ",", 
          RowBox[{"j", "-", "1"}], ",", "\"\<  \>\"", ",", 
          RowBox[{"Extract", "[", 
           RowBox[{"Pij", ",", 
            RowBox[{"{", 
             RowBox[{"i", ",", "j"}], "}"}]}], "]"}]}], "]"}]}], "]"}]}], 
     "]"}]}], "]"}], "*)"}]], "Input",
 CellChangeTimes->{{3.568108824079547*^9, 3.568109012824279*^9}, {
   3.568109558076436*^9, 3.568109588168889*^9}, {3.5681097523435774`*^9, 
   3.5681097637783976`*^9}, {3.5681100381828794`*^9, 
   3.5681101616726966`*^9}, {3.5681102015307665`*^9, 
   3.5681102288464146`*^9}, {3.568110263743676*^9, 3.568110267331682*^9}, {
   3.5681103163469687`*^9, 3.5681103268925867`*^9}, {3.5681106443375444`*^9, 
   3.5681106916368275`*^9}, 3.568110873002746*^9, {3.5681109388660617`*^9, 
   3.5681110268970165`*^9}, {3.568111092526332*^9, 3.568111123461186*^9}, {
   3.5681111884196997`*^9, 3.568111192756508*^9}, {3.5681112268113675`*^9, 
   3.568111366213212*^9}, {3.568111434946933*^9, 3.568111463354583*^9}, {
   3.568111544396725*^9, 3.568111547625931*^9}, {3.5681116856081734`*^9, 
   3.5681117131422215`*^9}, {3.568111864696488*^9, 3.5681119835998964`*^9}, 
   3.5681120218355637`*^9, {3.56811206543764*^9, 3.568112066389242*^9}, {
   3.568112332369709*^9, 3.5681123404505234`*^9}, {3.568112505935614*^9, 
   3.56811251481203*^9}, {3.5681125623609133`*^9, 3.568112569006525*^9}, {
   3.568112648067464*^9, 3.568112734210815*^9}, {3.568112764443668*^9, 
   3.5681127697320776`*^9}, {3.5681128018681335`*^9, 3.568112843426607*^9}, {
   3.568112884641879*^9, 3.5681128887602863`*^9}, {3.5681129663080225`*^9, 
   3.5681129699584293`*^9}, {3.568113133789917*^9, 3.568113227936082*^9}, {
   3.5681132946885996`*^9, 3.5681133044698167`*^9}, {3.568113502200164*^9, 
   3.568113507348173*^9}, 3.56811354556824*^9, {3.5681222886191964`*^9, 
   3.5681222896019983`*^9}, {3.568207424138049*^9, 3.5682074415164795`*^9}, 
   3.568277923134914*^9, 3.568883508322016*^9, 3.568883588350157*^9, 
   3.5688836441982546`*^9}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"DeleteFile", "[", " ", "\"\<C:\\\\Kimb_Mathem\\\\tmp4\>\"", "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"For", "[", 
  RowBox[{
   RowBox[{"i", "=", "2"}], ",", 
   RowBox[{"i", "<", 
    RowBox[{"n", "+", "2"}]}], ",", 
   RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
   RowBox[{"For", "[", 
    RowBox[{
     RowBox[{"j", "=", "2"}], ",", 
     RowBox[{"j", "\[LessEqual]", "i"}], ",", 
     RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", 
     RowBox[{"If", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"Pij", "[", 
         RowBox[{"[", 
          RowBox[{"i", ",", "j"}], "]"}], "]"}], ">", "0"}], ",", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{
          RowBox[{"i", "-", "1"}], ",", 
          RowBox[{"j", "-", "1"}], ",", 
          RowBox[{"Extract", "[", 
           RowBox[{"Pij", ",", 
            RowBox[{"{", 
             RowBox[{"i", ",", "j"}], "}"}]}], "]"}]}], "}"}], ">>>", " ", 
        "\"\<C:\\\\Kimb_Mathem\\\\tmp4\>\""}]}], "]"}]}], "]"}]}], 
  "]"}]}], "Input",
 CellChangeTimes->{{3.5681154444347754`*^9, 3.568115497147268*^9}, {
   3.568127214924899*^9, 3.568127265796589*^9}, {3.568128215276656*^9, 
   3.568128238536297*^9}, {3.568130718566253*^9, 3.5681307485339055`*^9}, 
   3.5681308175952272`*^9, {3.5681308916497574`*^9, 3.5681309458598523`*^9}, {
   3.5681309983383446`*^9, 3.5681310118947687`*^9}, {3.5681310723448744`*^9, 
   3.568131074357278*^9}, {3.5681315215476637`*^9, 3.568131584509374*^9}, {
   3.56813162789305*^9, 3.5681316427754765`*^9}, {3.5681317113063965`*^9, 
   3.568131715268804*^9}, {3.5681323203938665`*^9, 3.5681323260098763`*^9}, {
   3.568132387567585*^9, 3.568132440014877*^9}, {3.5681325121026034`*^9, 
   3.568132590945142*^9}, {3.56813274298981*^9, 3.5681327548302307`*^9}, {
   3.568132877867647*^9, 3.568132888335265*^9}, {3.568133052345154*^9, 
   3.568133072344389*^9}, {3.5681331740533676`*^9, 3.5681331951446047`*^9}, {
   3.568133267793932*^9, 3.5681332684803333`*^9}, {3.5681335079111547`*^9, 
   3.568133536833606*^9}, {3.5681336272201643`*^9, 3.5681336314945717`*^9}, {
   3.568133737994359*^9, 3.568133744000369*^9}, {3.568133978515581*^9, 
   3.5681339917912045`*^9}, {3.5681340324916763`*^9, 
   3.5681340418204927`*^9}, {3.568134244262048*^9, 3.5681342454944506`*^9}, {
   3.5681343354754086`*^9, 3.568134336099409*^9}, {3.5681345828762426`*^9, 
   3.568134583952645*^9}, {3.5681348406666956`*^9, 3.5681348474683075`*^9}, 
   3.568134889151581*^9, {3.5681352914606876`*^9, 3.568135296265496*^9}, {
   3.568135518877887*^9, 3.568135519595488*^9}, 3.5681357593679094`*^9, {
   3.5681787035360527`*^9, 3.568178728449296*^9}, 3.5682779128076954`*^9}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"a", "=", 
  RowBox[{
  "ReadList", "[", "\"\<C:\\\\Kimb_Mathem\\\\tmp4\>\"", 
   "]"}]}], "\[IndentingNewLine]", 
 RowBox[{"nb", "=", 
  RowBox[{"Length", "[", "a", "]"}]}], "\[IndentingNewLine]", 
 RowBox[{"cw", "=", 
  RowBox[{"nb", "-", "bh"}]}]}], "Input",
 CellChangeTimes->{{3.5681343270045934`*^9, 3.568134331481801*^9}, 
   3.5681343823222904`*^9, {3.5681345934218616`*^9, 3.5681345986322703`*^9}, {
   3.568134635479535*^9, 3.5681346521559644`*^9}, {3.5681351451792307`*^9, 
   3.568135145428831*^9}, {3.5681351781888885`*^9, 3.5681352349417877`*^9}, {
   3.568135340163973*^9, 3.568135340335573*^9}, {3.568135523495495*^9, 
   3.5681355244314966`*^9}, {3.5681357674799232`*^9, 
   3.5681357771675406`*^9}, {3.568136284433632*^9, 3.5681362987700567`*^9}, {
   3.568136395053426*^9, 3.5681363965510283`*^9}, {3.568381170343589*^9, 
   3.568381179844005*^9}, {3.5689533222924633`*^9, 3.568953376658559*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"2", ",", "1", ",", "1.1849933116283475`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"3", ",", "1", ",", "1.2799931097314476`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"4", ",", "1", ",", "1.2799901618360665`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"5", ",", "1", ",", "1.2799286426226604`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"6", ",", "2", ",", "1.2799743090211648`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"7", ",", "2", ",", "1.2800520427696944`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"8", ",", "2", ",", "1.2799875536973593`"}], "}"}]}], 
  "}"}]], "Output",
 CellChangeTimes->{{3.5690806271339383`*^9, 3.569080653950385*^9}, {
   3.5690807608261733`*^9, 3.5690808419151154`*^9}, 3.569080895314009*^9, {
   3.5690809261084633`*^9, 3.5690809603349237`*^9}, {3.569080994093383*^9, 
   3.5690810093346095`*^9}, 3.5690821727690535`*^9, 3.5690885956279345`*^9, 
   3.5691725087055717`*^9, 3.5691760356245656`*^9, 3.5691760998810787`*^9, 
   3.569176242714929*^9, 3.569176470444129*^9}],

Cell[BoxData["7"], "Output",
 CellChangeTimes->{{3.5690806271339383`*^9, 3.569080653950385*^9}, {
   3.5690807608261733`*^9, 3.5690808419151154`*^9}, 3.569080895314009*^9, {
   3.5690809261084633`*^9, 3.5690809603349237`*^9}, {3.569080994093383*^9, 
   3.5690810093346095`*^9}, 3.5690821727690535`*^9, 3.5690885956279345`*^9, 
   3.5691725087055717`*^9, 3.5691760356245656`*^9, 3.5691760998810787`*^9, 
   3.569176242714929*^9, 3.569176470459729*^9}],

Cell[BoxData["1"], "Output",
 CellChangeTimes->{{3.5690806271339383`*^9, 3.569080653950385*^9}, {
   3.5690807608261733`*^9, 3.5690808419151154`*^9}, 3.569080895314009*^9, {
   3.5690809261084633`*^9, 3.5690809603349237`*^9}, {3.569080994093383*^9, 
   3.5690810093346095`*^9}, 3.5690821727690535`*^9, 3.5690885956279345`*^9, 
   3.5691725087055717`*^9, 3.5691760356245656`*^9, 3.5691760998810787`*^9, 
   3.569176242714929*^9, 3.569176470459729*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"wxc", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"If", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"u", "[", 
         RowBox[{"[", 
          RowBox[{"a", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], 
        "\[NotEqual]", "\"\<H\>\""}], ",", 
       RowBox[{
        RowBox[{
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"a", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "3"}], "]"}], "]"}], "+", 
            RowBox[{"u2", "[", 
             RowBox[{"[", 
              RowBox[{"a", "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}]}], ")"}], 
          "/", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"u2", "[", 
             RowBox[{"[", 
              RowBox[{"a", "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "+", 
            RowBox[{"2", "*", 
             RowBox[{"a", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "3"}], "]"}], "]"}]}], "+", 
            RowBox[{"u2", "[", 
             RowBox[{"[", 
              RowBox[{"a", "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], 
           ")"}]}], "*", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"x", "[", 
            RowBox[{"[", 
             RowBox[{"a", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}], "-", 
           RowBox[{"x", "[", 
            RowBox[{"[", 
             RowBox[{"a", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}]}], ")"}]}],
         "+", 
        RowBox[{"x", "[", 
         RowBox[{"[", 
          RowBox[{"a", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}]}], ",", "0"}],
       "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "nb"}], "}"}]}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.5681365021632137`*^9, 3.5681365171392403`*^9}, {
   3.568136560351316*^9, 3.5681365863097615`*^9}, {3.5681369135895367`*^9, 
   3.5681369265999594`*^9}, {3.568137821916732*^9, 3.5681378378443604`*^9}, {
   3.5682072492149415`*^9, 3.5682072529901485`*^9}, 3.568209802830227*^9, {
   3.5682216503836985`*^9, 3.5682217583046885`*^9}, {3.5682217912831464`*^9, 
   3.568221843465238*^9}, {3.5682219760654707`*^9, 3.5682221133457117`*^9}, {
   3.568222159615393*^9, 3.5682221693498106`*^9}, {3.5682222565227633`*^9, 
   3.5682222578487663`*^9}, 3.5682224736283445`*^9, 3.5682227227607822`*^9, 
   3.5682227558640404`*^9, {3.568222818810151*^9, 3.5682228486062036`*^9}, {
   3.5682230057764797`*^9, 3.5682230079916835`*^9}, {3.568223041952943*^9, 
   3.5682230548229656`*^9}, 3.5682232619601297`*^9, 3.5682233000709963`*^9, {
   3.568223536645412*^9, 3.5682235381118145`*^9}, {3.568223591151908*^9, 
   3.568223591526308*^9}, {3.5682236491528096`*^9, 3.568223657483224*^9}, {
   3.568223713222122*^9, 3.5682237532361927`*^9}, {3.568223800379475*^9, 
   3.568223827648323*^9}, {3.5682238887380304`*^9, 3.5682238893152313`*^9}, {
   3.5682240287170763`*^9, 3.5682240612431335`*^9}, 3.5682241088076167`*^9, {
   3.568224198772975*^9, 3.5682242396450467`*^9}, {3.568224282170721*^9, 
   3.5682243185031853`*^9}, {3.568224437890195*^9, 3.568224439590598*^9}, {
   3.568224572253231*^9, 3.5682245867144566`*^9}, {3.5682247081762695`*^9, 
   3.568224733557514*^9}, {3.56822480520844*^9, 3.5682248433505073`*^9}, {
   3.568224884206979*^9, 3.5682248959069996`*^9}, {3.568224982830352*^9, 
   3.5682250047639904`*^9}, {3.56822503870965*^9, 3.5682250803305235`*^9}, {
   3.5682251223413973`*^9, 3.568225141451431*^9}, {3.568225234177994*^9, 
   3.5682252929588966`*^9}, {3.5682253775890455`*^9, 3.568225482140429*^9}, {
   3.56822551694409*^9, 3.5682255297205124`*^9}, {3.568225573197789*^9, 
   3.5682255809666023`*^9}, 3.568225660698343*^9, 3.5682257301964645`*^9, {
   3.568225765078126*^9, 3.5682257828933573`*^9}, {3.568225927026843*^9, 
   3.568225927416844*^9}, {3.568267062420238*^9, 3.568267081358671*^9}, {
   3.568267188328059*^9, 3.5682672830358257`*^9}, {3.5682676367040462`*^9, 
   3.568267747979042*^9}, 3.5682678404716043`*^9, 3.568267882030078*^9, {
   3.5682679252577534`*^9, 3.568267996612279*^9}, {3.5682680488879704`*^9, 
   3.568268061133992*^9}, {3.568268100914062*^9, 3.5682681150944867`*^9}, {
   3.5682681672609787`*^9, 3.5682681709113846`*^9}, {3.568268272904364*^9, 
   3.5682683182536435`*^9}, {3.5682685202739983`*^9, 3.568268599194537*^9}, {
   3.5682686348873997`*^9, 3.568268637773405*^9}, {3.5682686944951043`*^9, 
   3.568268872273017*^9}, {3.5682689420207396`*^9, 3.568269016947671*^9}, {
   3.5682695040741262`*^9, 3.568269601153097*^9}, {3.5682696652380095`*^9, 
   3.5682696728352227`*^9}, {3.5682697489633565`*^9, 
   3.5682698308635006`*^9}, {3.568271156491429*^9, 3.568271160063835*^9}, {
   3.568271244927984*^9, 3.568271369010602*^9}, {3.5682716250070515`*^9, 
   3.5682717312276382`*^9}, {3.5682717843301315`*^9, 
   3.5682718576814604`*^9}, {3.568271940673606*^9, 3.5682720536490045`*^9}, {
   3.5682722259669075`*^9, 3.5682723223126764`*^9}, {3.568272370813162*^9, 
   3.568272411466833*^9}, {3.5682725893539457`*^9, 3.568272622238803*^9}, {
   3.5682726936713285`*^9, 3.5682726960113325`*^9}, {3.568272741079812*^9, 
   3.56827275143823*^9}, {3.5682729218373294`*^9, 3.568272991070251*^9}, {
   3.5682734285574193`*^9, 3.5682734302890224`*^9}, {3.5682734898811274`*^9, 
   3.5682735111595645`*^9}, {3.5682736115925407`*^9, 
   3.5682736313577757`*^9}, {3.5682738601945777`*^9, 3.568273987272401*^9}, {
   3.5682763944346285`*^9, 3.5682764313910933`*^9}, 3.5682765128232365`*^9, {
   3.5682766527398825`*^9, 3.568276660976697*^9}, {3.5682774496048822`*^9, 
   3.568277461663703*^9}, {3.5683784032731285`*^9, 3.568378422991563*^9}, {
   3.568378459308427*^9, 3.5683784812420654`*^9}, {3.568378513440522*^9, 
   3.568378867046343*^9}, {3.568378915468828*^9, 3.56837896203491*^9}, {
   3.5683790764299107`*^9, 3.568379270041851*^9}, {3.568380233530743*^9, 
   3.568380234388745*^9}, 3.5683804214954734`*^9, 3.568380471867962*^9, 
   3.5683805399464817`*^9, {3.5683812396077104`*^9, 3.568381293537005*^9}, {
   3.568381339541486*^9, 3.5683813935175805`*^9}, 3.5684646396646404`*^9, 
   3.5684646988823442`*^9, 3.568465086387025*^9, {3.5684655859935026`*^9, 
   3.568465605087936*^9}, 3.568821235038273*^9, 3.5688233741452303`*^9, 
   3.568823449742963*^9, 3.5688274593312387`*^9, {3.5688284679978104`*^9, 
   3.5688284689026117`*^9}, {3.5688285131598897`*^9, 3.5688285201487017`*^9}, 
   3.5688285866360188`*^9, {3.5688292397623663`*^9, 3.568829246688778*^9}, {
   3.5688780328231993`*^9, 3.568878075458074*^9}, 3.5688953020257144`*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"wyc", "=", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"If", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"u", "[", 
          RowBox[{"[", 
           RowBox[{"a", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], 
         "\[NotEqual]", "\"\<H\>\""}], ",", 
        RowBox[{
         RowBox[{
          RowBox[{
           RowBox[{"(", 
            RowBox[{
             RowBox[{"a", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "3"}], "]"}], "]"}], "+", 
             RowBox[{"u2", "[", 
              RowBox[{"[", 
               RowBox[{"a", "[", 
                RowBox[{"[", 
                 RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}]}], ")"}],
            "/", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{"u2", "[", 
              RowBox[{"[", 
               RowBox[{"a", "[", 
                RowBox[{"[", 
                 RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "+", 
             RowBox[{"2", "*", 
              RowBox[{"a", "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", "3"}], "]"}], "]"}]}], "+", 
             RowBox[{"u2", "[", 
              RowBox[{"[", 
               RowBox[{"a", "[", 
                RowBox[{"[", 
                 RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], 
            ")"}]}], "*", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"y", "[", 
             RowBox[{"[", 
              RowBox[{"a", "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}], "-", 
            RowBox[{"y", "[", 
             RowBox[{"[", 
              RowBox[{"a", "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}]}], 
           ")"}]}], "+", 
         RowBox[{"y", "[", 
          RowBox[{"[", 
           RowBox[{"a", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}]}], ",", 
        "0"}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "nb"}], "}"}]}], "]"}]}], ";"}], 
  "   "}]], "Input",
 CellChangeTimes->{{3.5681365021632137`*^9, 3.5681365171392403`*^9}, {
   3.568136560351316*^9, 3.5681365863097615`*^9}, {3.5681369135895367`*^9, 
   3.5681369265999594`*^9}, {3.568137821916732*^9, 3.5681378378443604`*^9}, {
   3.5682072492149415`*^9, 3.5682072529901485`*^9}, 3.568209802830227*^9, {
   3.5682216503836985`*^9, 3.5682217583046885`*^9}, {3.5682217912831464`*^9, 
   3.568221843465238*^9}, {3.5682219760654707`*^9, 3.5682221133457117`*^9}, {
   3.568222159615393*^9, 3.5682221693498106`*^9}, {3.5682222565227633`*^9, 
   3.5682222578487663`*^9}, 3.5682224736283445`*^9, 3.5682227227607822`*^9, 
   3.5682227558640404`*^9, {3.568222818810151*^9, 3.5682228486062036`*^9}, {
   3.5682230057764797`*^9, 3.5682230079916835`*^9}, {3.568223041952943*^9, 
   3.5682230548229656`*^9}, 3.5682232619601297`*^9, 3.5682233000709963`*^9, {
   3.568223536645412*^9, 3.5682235381118145`*^9}, {3.568223591151908*^9, 
   3.568223591526308*^9}, {3.5682236491528096`*^9, 3.568223657483224*^9}, {
   3.568223713222122*^9, 3.5682237532361927`*^9}, {3.568223800379475*^9, 
   3.568223827648323*^9}, {3.5682238887380304`*^9, 3.5682238893152313`*^9}, {
   3.5682240287170763`*^9, 3.5682240612431335`*^9}, 3.5682241088076167`*^9, {
   3.568224198772975*^9, 3.5682242396450467`*^9}, {3.568224282170721*^9, 
   3.5682243185031853`*^9}, {3.568224437890195*^9, 3.568224439590598*^9}, {
   3.568224572253231*^9, 3.5682245867144566`*^9}, {3.5682247081762695`*^9, 
   3.568224733557514*^9}, {3.56822480520844*^9, 3.5682248433505073`*^9}, {
   3.568224884206979*^9, 3.5682248959069996`*^9}, {3.568224982830352*^9, 
   3.5682250047639904`*^9}, {3.56822503870965*^9, 3.5682250803305235`*^9}, {
   3.5682251223413973`*^9, 3.568225141451431*^9}, {3.568225234177994*^9, 
   3.5682252929588966`*^9}, {3.5682253775890455`*^9, 3.568225482140429*^9}, {
   3.56822551694409*^9, 3.5682255297205124`*^9}, {3.568225573197789*^9, 
   3.5682255809666023`*^9}, 3.568225660698343*^9, 3.5682257301964645`*^9, {
   3.568225765078126*^9, 3.5682257828933573`*^9}, {3.568225927026843*^9, 
   3.568225927416844*^9}, {3.568267062420238*^9, 3.568267081358671*^9}, {
   3.568267188328059*^9, 3.5682672830358257`*^9}, {3.5682676367040462`*^9, 
   3.568267747979042*^9}, 3.5682678404716043`*^9, 3.568267882030078*^9, {
   3.5682679252577534`*^9, 3.568267996612279*^9}, {3.5682680488879704`*^9, 
   3.568268061133992*^9}, {3.568268100914062*^9, 3.5682681150944867`*^9}, {
   3.5682681672609787`*^9, 3.5682681709113846`*^9}, {3.568268272904364*^9, 
   3.5682683182536435`*^9}, {3.5682685202739983`*^9, 3.568268599194537*^9}, {
   3.5682686348873997`*^9, 3.568268637773405*^9}, {3.5682686944951043`*^9, 
   3.568268872273017*^9}, {3.5682689420207396`*^9, 3.568269016947671*^9}, {
   3.5682695040741262`*^9, 3.568269601153097*^9}, {3.5682696652380095`*^9, 
   3.5682696728352227`*^9}, {3.5682697489633565`*^9, 
   3.5682698308635006`*^9}, {3.568271156491429*^9, 3.568271160063835*^9}, {
   3.568271244927984*^9, 3.568271369010602*^9}, {3.5682716250070515`*^9, 
   3.5682717312276382`*^9}, {3.5682717843301315`*^9, 
   3.5682718576814604`*^9}, {3.568271940673606*^9, 3.5682720536490045`*^9}, {
   3.5682722259669075`*^9, 3.5682723223126764`*^9}, {3.568272370813162*^9, 
   3.568272411466833*^9}, {3.5682725893539457`*^9, 3.568272622238803*^9}, {
   3.5682726936713285`*^9, 3.5682726960113325`*^9}, {3.568272741079812*^9, 
   3.56827275143823*^9}, {3.5682729218373294`*^9, 3.568272991070251*^9}, {
   3.5682734285574193`*^9, 3.5682734302890224`*^9}, {3.5682734898811274`*^9, 
   3.5682735111595645`*^9}, {3.5682736115925407`*^9, 
   3.5682736313577757`*^9}, {3.5682738601945777`*^9, 3.568273987272401*^9}, {
   3.5682763944346285`*^9, 3.5682764313910933`*^9}, 3.5682765128232365`*^9, {
   3.5682766527398825`*^9, 3.568276660976697*^9}, {3.5682774496048822`*^9, 
   3.568277461663703*^9}, {3.5683784032731285`*^9, 3.568378422991563*^9}, {
   3.568378459308427*^9, 3.5683784812420654`*^9}, {3.568378513440522*^9, 
   3.568378867046343*^9}, {3.568378915468828*^9, 3.56837896203491*^9}, {
   3.5683790764299107`*^9, 3.568379270041851*^9}, {3.568380233530743*^9, 
   3.568380234388745*^9}, 3.5683804214954734`*^9, 3.568380471867962*^9, 
   3.5683805399464817`*^9, {3.5683812396077104`*^9, 3.568381293537005*^9}, {
   3.568381339541486*^9, 3.5683813935175805`*^9}, 3.5684646396646404`*^9, 
   3.5684646988823442`*^9, 3.568465086387025*^9, {3.5684655859935026`*^9, 
   3.568465605087936*^9}, 3.568821235038273*^9, 3.5688233741452303`*^9, 
   3.568823449742963*^9, 3.5688274593312387`*^9, {3.5688284679978104`*^9, 
   3.5688284689026117`*^9}, {3.5688285131598897`*^9, 3.5688285201487017`*^9}, 
   3.5688285866360188`*^9, {3.5688292397623663`*^9, 3.568829246688778*^9}, {
   3.5688780328231993`*^9, 3.5688780904809003`*^9}, 3.5688953042565184`*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"wzc", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"If", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"u", "[", 
         RowBox[{"[", 
          RowBox[{"a", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], 
        "\[NotEqual]", "\"\<H\>\""}], ",", 
       RowBox[{
        RowBox[{
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"a", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "3"}], "]"}], "]"}], "+", 
            RowBox[{"u2", "[", 
             RowBox[{"[", 
              RowBox[{"a", "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}]}], ")"}], 
          "/", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"u2", "[", 
             RowBox[{"[", 
              RowBox[{"a", "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "+", 
            RowBox[{"2", "*", 
             RowBox[{"a", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "3"}], "]"}], "]"}]}], "+", 
            RowBox[{"u2", "[", 
             RowBox[{"[", 
              RowBox[{"a", "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], 
           ")"}]}], "*", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"z", "[", 
            RowBox[{"[", 
             RowBox[{"a", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}], "-", 
           RowBox[{"z", "[", 
            RowBox[{"[", 
             RowBox[{"a", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}]}], ")"}]}],
         "+", 
        RowBox[{"z", "[", 
         RowBox[{"[", 
          RowBox[{"a", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}]}], ",", "0"}],
       "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "nb"}], "}"}]}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.5681365021632137`*^9, 3.5681365171392403`*^9}, {
   3.568136560351316*^9, 3.5681365863097615`*^9}, {3.5681369135895367`*^9, 
   3.5681369265999594`*^9}, {3.568137821916732*^9, 3.5681378378443604`*^9}, {
   3.5682072492149415`*^9, 3.5682072529901485`*^9}, 3.568209802830227*^9, {
   3.5682216503836985`*^9, 3.5682217583046885`*^9}, {3.5682217912831464`*^9, 
   3.568221843465238*^9}, {3.5682219760654707`*^9, 3.5682221133457117`*^9}, {
   3.568222159615393*^9, 3.5682221693498106`*^9}, {3.5682222565227633`*^9, 
   3.5682222578487663`*^9}, 3.5682224736283445`*^9, 3.5682227227607822`*^9, 
   3.5682227558640404`*^9, {3.568222818810151*^9, 3.5682228486062036`*^9}, {
   3.5682230057764797`*^9, 3.5682230079916835`*^9}, {3.568223041952943*^9, 
   3.5682230548229656`*^9}, 3.5682232619601297`*^9, 3.5682233000709963`*^9, {
   3.568223536645412*^9, 3.5682235381118145`*^9}, {3.568223591151908*^9, 
   3.568223591526308*^9}, {3.5682236491528096`*^9, 3.568223657483224*^9}, {
   3.568223713222122*^9, 3.5682237532361927`*^9}, {3.568223800379475*^9, 
   3.568223827648323*^9}, {3.5682238887380304`*^9, 3.5682238893152313`*^9}, {
   3.5682240287170763`*^9, 3.5682240612431335`*^9}, 3.5682241088076167`*^9, {
   3.568224198772975*^9, 3.5682242396450467`*^9}, {3.568224282170721*^9, 
   3.5682243185031853`*^9}, {3.568224437890195*^9, 3.568224439590598*^9}, {
   3.568224572253231*^9, 3.5682245867144566`*^9}, {3.5682247081762695`*^9, 
   3.568224733557514*^9}, {3.56822480520844*^9, 3.5682248433505073`*^9}, {
   3.568224884206979*^9, 3.5682248959069996`*^9}, {3.568224982830352*^9, 
   3.5682250047639904`*^9}, {3.56822503870965*^9, 3.5682250803305235`*^9}, {
   3.5682251223413973`*^9, 3.568225141451431*^9}, {3.568225234177994*^9, 
   3.5682252929588966`*^9}, {3.5682253775890455`*^9, 3.568225482140429*^9}, {
   3.56822551694409*^9, 3.5682255297205124`*^9}, {3.568225573197789*^9, 
   3.5682255809666023`*^9}, 3.568225660698343*^9, 3.5682257301964645`*^9, {
   3.568225765078126*^9, 3.5682257828933573`*^9}, {3.568225927026843*^9, 
   3.568225927416844*^9}, {3.568267062420238*^9, 3.568267081358671*^9}, {
   3.568267188328059*^9, 3.5682672830358257`*^9}, {3.5682676367040462`*^9, 
   3.568267747979042*^9}, 3.5682678404716043`*^9, 3.568267882030078*^9, {
   3.5682679252577534`*^9, 3.568267996612279*^9}, {3.5682680488879704`*^9, 
   3.568268061133992*^9}, {3.568268100914062*^9, 3.5682681150944867`*^9}, {
   3.5682681672609787`*^9, 3.5682681709113846`*^9}, {3.568268272904364*^9, 
   3.5682683182536435`*^9}, {3.5682685202739983`*^9, 3.568268599194537*^9}, {
   3.5682686348873997`*^9, 3.568268637773405*^9}, {3.5682686944951043`*^9, 
   3.568268872273017*^9}, {3.5682689420207396`*^9, 3.568269016947671*^9}, {
   3.5682695040741262`*^9, 3.568269601153097*^9}, {3.5682696652380095`*^9, 
   3.5682696728352227`*^9}, {3.5682697489633565`*^9, 
   3.5682698308635006`*^9}, {3.568271156491429*^9, 3.568271160063835*^9}, {
   3.568271244927984*^9, 3.568271369010602*^9}, {3.5682716250070515`*^9, 
   3.5682717312276382`*^9}, {3.5682717843301315`*^9, 
   3.5682718576814604`*^9}, {3.568271940673606*^9, 3.5682720536490045`*^9}, {
   3.5682722259669075`*^9, 3.5682723223126764`*^9}, {3.568272370813162*^9, 
   3.568272411466833*^9}, {3.5682725893539457`*^9, 3.568272622238803*^9}, {
   3.5682726936713285`*^9, 3.5682726960113325`*^9}, {3.568272741079812*^9, 
   3.56827275143823*^9}, {3.5682729218373294`*^9, 3.568272991070251*^9}, {
   3.5682734285574193`*^9, 3.5682734302890224`*^9}, {3.5682734898811274`*^9, 
   3.5682735111595645`*^9}, {3.5682736115925407`*^9, 
   3.5682736313577757`*^9}, {3.5682738601945777`*^9, 3.568273987272401*^9}, {
   3.5682763944346285`*^9, 3.5682764313910933`*^9}, 3.5682765128232365`*^9, {
   3.5682766527398825`*^9, 3.568276660976697*^9}, {3.5682774496048822`*^9, 
   3.568277461663703*^9}, {3.5683784032731285`*^9, 3.568378422991563*^9}, {
   3.568378459308427*^9, 3.5683784812420654`*^9}, {3.568378513440522*^9, 
   3.568378867046343*^9}, {3.568378915468828*^9, 3.56837896203491*^9}, {
   3.5683790764299107`*^9, 3.568379270041851*^9}, {3.568380233530743*^9, 
   3.568380234388745*^9}, 3.5683804214954734`*^9, 3.568380471867962*^9, 
   3.5683805399464817`*^9, {3.5683812396077104`*^9, 3.568381293537005*^9}, {
   3.568381339541486*^9, 3.5683813935175805`*^9}, 3.5684646396646404`*^9, 
   3.5684646988823442`*^9, 3.568465086387025*^9, {3.5684655859935026`*^9, 
   3.568465605087936*^9}, 3.568821235038273*^9, 3.5688233741452303`*^9, 
   3.568823449742963*^9, 3.5688274593312387`*^9, {3.5688284679978104`*^9, 
   3.5688284689026117`*^9}, {3.5688285131598897`*^9, 3.5688285201487017`*^9}, 
   3.5688285866360188`*^9, {3.5688292397623663`*^9, 3.568829246688778*^9}, {
   3.5688780328231993`*^9, 3.568878099435316*^9}, 3.5688953066901226`*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"Plot7", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"Graphics3D", "[", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"Darker", "[", "Gray", "]"}], ",", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"Opacity", "[", "0.5", "]"}], ",", 
          RowBox[{"If", "[", 
           RowBox[{
            RowBox[{
             RowBox[{"wxc", "[", 
              RowBox[{"[", "i", "]"}], "]"}], "\[NotEqual]", "0"}], ",", 
            RowBox[{"Sphere", "[", 
             RowBox[{
              RowBox[{"{", 
               RowBox[{
                RowBox[{"wxc", "[", 
                 RowBox[{"[", "i", "]"}], "]"}], ",", 
                RowBox[{"wyc", "[", 
                 RowBox[{"[", "i", "]"}], "]"}], ",", 
                RowBox[{"wzc", "[", 
                 RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", 
              RowBox[{"a", "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", "3"}], "]"}], "]"}]}], "]"}]}], "]"}]}], 
         "}"}]}], "}"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "nb"}], "}"}]}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.5681365021632137`*^9, 3.5681365171392403`*^9}, {
   3.568136560351316*^9, 3.5681365863097615`*^9}, {3.5681369135895367`*^9, 
   3.5681369265999594`*^9}, {3.568137821916732*^9, 3.5681378378443604`*^9}, {
   3.5682072492149415`*^9, 3.5682072529901485`*^9}, 3.568209802830227*^9, {
   3.5682216503836985`*^9, 3.5682217583046885`*^9}, {3.5682217912831464`*^9, 
   3.568221843465238*^9}, {3.5682219760654707`*^9, 3.5682221133457117`*^9}, {
   3.568222159615393*^9, 3.5682221693498106`*^9}, {3.5682222565227633`*^9, 
   3.5682222578487663`*^9}, 3.5682224736283445`*^9, 3.5682227227607822`*^9, 
   3.5682227558640404`*^9, {3.568222818810151*^9, 3.5682228486062036`*^9}, {
   3.5682230057764797`*^9, 3.5682230079916835`*^9}, {3.568223041952943*^9, 
   3.5682230548229656`*^9}, 3.5682232619601297`*^9, 3.5682233000709963`*^9, {
   3.568223536645412*^9, 3.5682235381118145`*^9}, {3.568223591151908*^9, 
   3.568223591526308*^9}, {3.5682236491528096`*^9, 3.568223657483224*^9}, {
   3.568223713222122*^9, 3.5682237532361927`*^9}, {3.568223800379475*^9, 
   3.568223827648323*^9}, {3.5682238887380304`*^9, 3.5682238893152313`*^9}, {
   3.5682240287170763`*^9, 3.5682240612431335`*^9}, 3.5682241088076167`*^9, {
   3.568224198772975*^9, 3.5682242396450467`*^9}, {3.568224282170721*^9, 
   3.5682243185031853`*^9}, {3.568224437890195*^9, 3.568224439590598*^9}, {
   3.568224572253231*^9, 3.5682245867144566`*^9}, {3.5682247081762695`*^9, 
   3.568224733557514*^9}, {3.56822480520844*^9, 3.5682248433505073`*^9}, {
   3.568224884206979*^9, 3.5682248959069996`*^9}, {3.568224982830352*^9, 
   3.5682250047639904`*^9}, {3.56822503870965*^9, 3.5682250803305235`*^9}, {
   3.5682251223413973`*^9, 3.568225141451431*^9}, {3.568225234177994*^9, 
   3.5682252929588966`*^9}, {3.5682253775890455`*^9, 3.568225482140429*^9}, {
   3.56822551694409*^9, 3.5682255297205124`*^9}, {3.568225573197789*^9, 
   3.5682255809666023`*^9}, 3.568225660698343*^9, 3.5682257301964645`*^9, {
   3.568225765078126*^9, 3.5682257828933573`*^9}, {3.568225927026843*^9, 
   3.568225927416844*^9}, {3.568267062420238*^9, 3.568267081358671*^9}, {
   3.568267188328059*^9, 3.5682672830358257`*^9}, {3.5682676367040462`*^9, 
   3.568267747979042*^9}, 3.5682678404716043`*^9, 3.568267882030078*^9, {
   3.5682679252577534`*^9, 3.568267996612279*^9}, {3.5682680488879704`*^9, 
   3.568268061133992*^9}, {3.568268100914062*^9, 3.5682681150944867`*^9}, {
   3.5682681672609787`*^9, 3.5682681709113846`*^9}, {3.568268272904364*^9, 
   3.5682683182536435`*^9}, {3.5682685202739983`*^9, 3.568268599194537*^9}, {
   3.5682686348873997`*^9, 3.568268637773405*^9}, {3.5682686944951043`*^9, 
   3.568268872273017*^9}, {3.5682689420207396`*^9, 3.568269016947671*^9}, {
   3.5682695040741262`*^9, 3.568269601153097*^9}, {3.5682696652380095`*^9, 
   3.5682696728352227`*^9}, {3.5682697489633565`*^9, 
   3.5682698308635006`*^9}, {3.568271156491429*^9, 3.568271160063835*^9}, {
   3.568271244927984*^9, 3.568271369010602*^9}, {3.5682716250070515`*^9, 
   3.5682717312276382`*^9}, {3.5682717843301315`*^9, 
   3.5682718576814604`*^9}, {3.568271940673606*^9, 3.5682720536490045`*^9}, {
   3.5682722259669075`*^9, 3.5682723223126764`*^9}, {3.568272370813162*^9, 
   3.568272411466833*^9}, {3.5682725893539457`*^9, 3.568272622238803*^9}, {
   3.5682726936713285`*^9, 3.5682726960113325`*^9}, {3.568272741079812*^9, 
   3.56827275143823*^9}, {3.5682729218373294`*^9, 3.568272991070251*^9}, {
   3.5682734285574193`*^9, 3.5682734302890224`*^9}, {3.5682734898811274`*^9, 
   3.5682735111595645`*^9}, {3.5682736115925407`*^9, 
   3.5682736313577757`*^9}, {3.5682738601945777`*^9, 3.568273987272401*^9}, {
   3.5682763944346285`*^9, 3.5682764313910933`*^9}, 3.5682765128232365`*^9, {
   3.5682766527398825`*^9, 3.568276660976697*^9}, {3.5682774496048822`*^9, 
   3.568277461663703*^9}, {3.5683784032731285`*^9, 3.568378422991563*^9}, {
   3.568378459308427*^9, 3.5683784812420654`*^9}, {3.568378513440522*^9, 
   3.568378867046343*^9}, {3.568378915468828*^9, 3.56837896203491*^9}, {
   3.5683790764299107`*^9, 3.568379270041851*^9}, {3.568380233530743*^9, 
   3.568380234388745*^9}, 3.5683804214954734`*^9, 3.568380471867962*^9, 
   3.5683805399464817`*^9, {3.5683812396077104`*^9, 3.568381293537005*^9}, {
   3.568381339541486*^9, 3.5683813935175805`*^9}, 3.5684646396646404`*^9, 
   3.5684646988823442`*^9, 3.568465086387025*^9, {3.5684655859935026`*^9, 
   3.568465605087936*^9}, 3.568821235038273*^9, 3.5688233741452303`*^9, 
   3.568823449742963*^9, 3.5688274593312387`*^9, {3.5688284679978104`*^9, 
   3.5688284689026117`*^9}, {3.5688285131598897`*^9, 3.5688285201487017`*^9}, 
   3.5688285866360188`*^9, {3.5688292397623663`*^9, 3.568829246688778*^9}, {
   3.5688780328231993`*^9, 3.56887811871695*^9}, {3.568895331556566*^9, 
   3.5688953367201753`*^9}, 3.568895995276532*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"wxh", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"If", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"u", "[", 
         RowBox[{"[", 
          RowBox[{"a", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "\[Equal]", 
        "\"\<H\>\""}], ",", 
       RowBox[{
        RowBox[{
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"a", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "3"}], "]"}], "]"}], "+", 
            RowBox[{"u2", "[", 
             RowBox[{"[", 
              RowBox[{"a", "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ")"}], 
          "/", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"u2", "[", 
             RowBox[{"[", 
              RowBox[{"a", "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "+", 
            RowBox[{"a", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "3"}], "]"}], "]"}], "+", 
            RowBox[{"u2", "[", 
             RowBox[{"[", 
              RowBox[{"a", "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], 
           ")"}]}], "*", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"x", "[", 
            RowBox[{"[", 
             RowBox[{"a", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "-", 
           RowBox[{"x", "[", 
            RowBox[{"[", 
             RowBox[{"a", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ")"}]}],
         "+", 
        RowBox[{"x", "[", 
         RowBox[{"[", 
          RowBox[{"a", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ",", "0"}],
       "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "nb"}], "}"}]}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.5681365021632137`*^9, 3.5681365171392403`*^9}, {
   3.568136560351316*^9, 3.5681365863097615`*^9}, {3.5681369135895367`*^9, 
   3.5681369265999594`*^9}, {3.568137821916732*^9, 3.5681378378443604`*^9}, {
   3.5682072492149415`*^9, 3.5682072529901485`*^9}, 3.568209802830227*^9, {
   3.5682216503836985`*^9, 3.5682217583046885`*^9}, {3.5682217912831464`*^9, 
   3.568221843465238*^9}, {3.5682219760654707`*^9, 3.5682221133457117`*^9}, {
   3.568222159615393*^9, 3.5682221693498106`*^9}, {3.5682222565227633`*^9, 
   3.5682222578487663`*^9}, 3.5682224736283445`*^9, 3.5682227227607822`*^9, 
   3.5682227558640404`*^9, {3.568222818810151*^9, 3.5682228486062036`*^9}, {
   3.5682230057764797`*^9, 3.5682230079916835`*^9}, {3.568223041952943*^9, 
   3.5682230548229656`*^9}, 3.5682232619601297`*^9, 3.5682233000709963`*^9, {
   3.568223536645412*^9, 3.5682235381118145`*^9}, {3.568223591151908*^9, 
   3.568223591526308*^9}, {3.5682236491528096`*^9, 3.568223657483224*^9}, {
   3.568223713222122*^9, 3.5682237532361927`*^9}, {3.568223800379475*^9, 
   3.568223827648323*^9}, {3.5682238887380304`*^9, 3.5682238893152313`*^9}, {
   3.5682240287170763`*^9, 3.5682240612431335`*^9}, 3.5682241088076167`*^9, {
   3.568224198772975*^9, 3.5682242396450467`*^9}, {3.568224282170721*^9, 
   3.5682243185031853`*^9}, {3.568224437890195*^9, 3.568224439590598*^9}, {
   3.568224572253231*^9, 3.5682245867144566`*^9}, {3.5682247081762695`*^9, 
   3.568224733557514*^9}, {3.56822480520844*^9, 3.5682248433505073`*^9}, {
   3.568224884206979*^9, 3.5682248959069996`*^9}, {3.568224982830352*^9, 
   3.5682250047639904`*^9}, {3.56822503870965*^9, 3.5682250803305235`*^9}, {
   3.5682251223413973`*^9, 3.568225141451431*^9}, {3.568225234177994*^9, 
   3.5682252929588966`*^9}, {3.5682253775890455`*^9, 3.568225482140429*^9}, {
   3.56822551694409*^9, 3.5682255297205124`*^9}, {3.568225573197789*^9, 
   3.5682255809666023`*^9}, 3.568225660698343*^9, 3.5682257301964645`*^9, {
   3.568225765078126*^9, 3.5682257828933573`*^9}, {3.568225927026843*^9, 
   3.568225927416844*^9}, {3.568267062420238*^9, 3.568267081358671*^9}, {
   3.568267188328059*^9, 3.5682672830358257`*^9}, {3.5682676367040462`*^9, 
   3.568267747979042*^9}, 3.5682678404716043`*^9, 3.568267882030078*^9, {
   3.5682679252577534`*^9, 3.568267996612279*^9}, {3.5682680488879704`*^9, 
   3.568268061133992*^9}, {3.568268100914062*^9, 3.5682681150944867`*^9}, {
   3.5682681672609787`*^9, 3.5682681709113846`*^9}, {3.568268272904364*^9, 
   3.5682683182536435`*^9}, {3.5682685202739983`*^9, 3.568268599194537*^9}, {
   3.5682686348873997`*^9, 3.568268637773405*^9}, {3.5682686944951043`*^9, 
   3.568268872273017*^9}, {3.5682689420207396`*^9, 3.568269016947671*^9}, {
   3.5682695040741262`*^9, 3.568269601153097*^9}, {3.5682696652380095`*^9, 
   3.5682696728352227`*^9}, {3.5682697489633565`*^9, 
   3.5682698308635006`*^9}, {3.568271156491429*^9, 3.568271160063835*^9}, {
   3.568271244927984*^9, 3.568271369010602*^9}, {3.5682716250070515`*^9, 
   3.5682717312276382`*^9}, {3.5682717843301315`*^9, 
   3.5682718576814604`*^9}, {3.568271940673606*^9, 3.5682720536490045`*^9}, {
   3.5682722259669075`*^9, 3.5682723223126764`*^9}, {3.568272370813162*^9, 
   3.568272411466833*^9}, {3.5682725893539457`*^9, 3.568272622238803*^9}, {
   3.5682726936713285`*^9, 3.5682726960113325`*^9}, {3.568272741079812*^9, 
   3.56827275143823*^9}, {3.5682729218373294`*^9, 3.568272991070251*^9}, {
   3.5682734285574193`*^9, 3.5682734302890224`*^9}, {3.5682734898811274`*^9, 
   3.5682735111595645`*^9}, {3.5682736115925407`*^9, 
   3.5682736313577757`*^9}, {3.5682738601945777`*^9, 3.568273987272401*^9}, {
   3.5682763944346285`*^9, 3.5682764313910933`*^9}, 3.5682765128232365`*^9, {
   3.5682766527398825`*^9, 3.568276660976697*^9}, {3.5682774496048822`*^9, 
   3.568277461663703*^9}, {3.5683784032731285`*^9, 3.568378422991563*^9}, {
   3.568378459308427*^9, 3.5683784812420654`*^9}, {3.568378513440522*^9, 
   3.568378867046343*^9}, {3.568378915468828*^9, 3.56837896203491*^9}, {
   3.5683790764299107`*^9, 3.568379270041851*^9}, {3.568380233530743*^9, 
   3.568380234388745*^9}, 3.5683804214954734`*^9, 3.568380471867962*^9, 
   3.5683805399464817`*^9, {3.5683812396077104`*^9, 3.568381293537005*^9}, {
   3.568381339541486*^9, 3.5683813935175805`*^9}, 3.5684646396646404`*^9, 
   3.5684646988823442`*^9, 3.568465086387025*^9, {3.5684655859935026`*^9, 
   3.568465605087936*^9}, 3.568821235038273*^9, 3.5688233741452303`*^9, 
   3.568823449742963*^9, 3.5688274593312387`*^9, {3.5688284679978104`*^9, 
   3.5688284689026117`*^9}, {3.5688285131598897`*^9, 3.5688285201487017`*^9}, 
   3.5688285866360188`*^9, {3.5688292397623663`*^9, 3.568829246688778*^9}, {
   3.5688780328231993`*^9, 3.5688781280925665`*^9}, 3.568895312524533*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"wyh", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"If", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"u", "[", 
         RowBox[{"[", 
          RowBox[{"a", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "\[Equal]", 
        "\"\<H\>\""}], ",", 
       RowBox[{
        RowBox[{
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"a", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "3"}], "]"}], "]"}], "+", 
            RowBox[{"u2", "[", 
             RowBox[{"[", 
              RowBox[{"a", "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ")"}], 
          "/", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"u2", "[", 
             RowBox[{"[", 
              RowBox[{"a", "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "+", 
            RowBox[{"a", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "3"}], "]"}], "]"}], "+", 
            RowBox[{"u2", "[", 
             RowBox[{"[", 
              RowBox[{"a", "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], 
           ")"}]}], "*", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"y", "[", 
            RowBox[{"[", 
             RowBox[{"a", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "-", 
           RowBox[{"y", "[", 
            RowBox[{"[", 
             RowBox[{"a", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ")"}]}],
         "+", 
        RowBox[{"y", "[", 
         RowBox[{"[", 
          RowBox[{"a", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ",", "0"}],
       "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "nb"}], "}"}]}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.5681365021632137`*^9, 3.5681365171392403`*^9}, {
   3.568136560351316*^9, 3.5681365863097615`*^9}, {3.5681369135895367`*^9, 
   3.5681369265999594`*^9}, {3.568137821916732*^9, 3.5681378378443604`*^9}, {
   3.5682072492149415`*^9, 3.5682072529901485`*^9}, 3.568209802830227*^9, {
   3.5682216503836985`*^9, 3.5682217583046885`*^9}, {3.5682217912831464`*^9, 
   3.568221843465238*^9}, {3.5682219760654707`*^9, 3.5682221133457117`*^9}, {
   3.568222159615393*^9, 3.5682221693498106`*^9}, {3.5682222565227633`*^9, 
   3.5682222578487663`*^9}, 3.5682224736283445`*^9, 3.5682227227607822`*^9, 
   3.5682227558640404`*^9, {3.568222818810151*^9, 3.5682228486062036`*^9}, {
   3.5682230057764797`*^9, 3.5682230079916835`*^9}, {3.568223041952943*^9, 
   3.5682230548229656`*^9}, 3.5682232619601297`*^9, 3.5682233000709963`*^9, {
   3.568223536645412*^9, 3.5682235381118145`*^9}, {3.568223591151908*^9, 
   3.568223591526308*^9}, {3.5682236491528096`*^9, 3.568223657483224*^9}, {
   3.568223713222122*^9, 3.5682237532361927`*^9}, {3.568223800379475*^9, 
   3.568223827648323*^9}, {3.5682238887380304`*^9, 3.5682238893152313`*^9}, {
   3.5682240287170763`*^9, 3.5682240612431335`*^9}, 3.5682241088076167`*^9, {
   3.568224198772975*^9, 3.5682242396450467`*^9}, {3.568224282170721*^9, 
   3.5682243185031853`*^9}, {3.568224437890195*^9, 3.568224439590598*^9}, {
   3.568224572253231*^9, 3.5682245867144566`*^9}, {3.5682247081762695`*^9, 
   3.568224733557514*^9}, {3.56822480520844*^9, 3.5682248433505073`*^9}, {
   3.568224884206979*^9, 3.5682248959069996`*^9}, {3.568224982830352*^9, 
   3.5682250047639904`*^9}, {3.56822503870965*^9, 3.5682250803305235`*^9}, {
   3.5682251223413973`*^9, 3.568225141451431*^9}, {3.568225234177994*^9, 
   3.5682252929588966`*^9}, {3.5682253775890455`*^9, 3.568225482140429*^9}, {
   3.56822551694409*^9, 3.5682255297205124`*^9}, {3.568225573197789*^9, 
   3.5682255809666023`*^9}, 3.568225660698343*^9, 3.5682257301964645`*^9, {
   3.568225765078126*^9, 3.5682257828933573`*^9}, {3.568225927026843*^9, 
   3.568225927416844*^9}, {3.568267062420238*^9, 3.568267081358671*^9}, {
   3.568267188328059*^9, 3.5682672830358257`*^9}, {3.5682676367040462`*^9, 
   3.568267747979042*^9}, 3.5682678404716043`*^9, 3.568267882030078*^9, {
   3.5682679252577534`*^9, 3.568267996612279*^9}, {3.5682680488879704`*^9, 
   3.568268061133992*^9}, {3.568268100914062*^9, 3.5682681150944867`*^9}, {
   3.5682681672609787`*^9, 3.5682681709113846`*^9}, {3.568268272904364*^9, 
   3.5682683182536435`*^9}, {3.5682685202739983`*^9, 3.568268599194537*^9}, {
   3.5682686348873997`*^9, 3.568268637773405*^9}, {3.5682686944951043`*^9, 
   3.568268872273017*^9}, {3.5682689420207396`*^9, 3.568269016947671*^9}, {
   3.5682695040741262`*^9, 3.568269601153097*^9}, {3.5682696652380095`*^9, 
   3.5682696728352227`*^9}, {3.5682697489633565`*^9, 
   3.5682698308635006`*^9}, {3.568271156491429*^9, 3.568271160063835*^9}, {
   3.568271244927984*^9, 3.568271369010602*^9}, {3.5682716250070515`*^9, 
   3.5682717312276382`*^9}, {3.5682717843301315`*^9, 
   3.5682718576814604`*^9}, {3.568271940673606*^9, 3.5682720536490045`*^9}, {
   3.5682722259669075`*^9, 3.5682723223126764`*^9}, {3.568272370813162*^9, 
   3.568272411466833*^9}, {3.5682725893539457`*^9, 3.568272622238803*^9}, {
   3.5682726936713285`*^9, 3.5682726960113325`*^9}, {3.568272741079812*^9, 
   3.56827275143823*^9}, {3.5682729218373294`*^9, 3.568272991070251*^9}, {
   3.5682734285574193`*^9, 3.5682734302890224`*^9}, {3.5682734898811274`*^9, 
   3.5682735111595645`*^9}, {3.5682736115925407`*^9, 
   3.5682736313577757`*^9}, {3.5682738601945777`*^9, 3.568273987272401*^9}, {
   3.5682763944346285`*^9, 3.5682764313910933`*^9}, 3.5682765128232365`*^9, {
   3.5682766527398825`*^9, 3.568276660976697*^9}, {3.5682774496048822`*^9, 
   3.568277461663703*^9}, {3.5683784032731285`*^9, 3.568378422991563*^9}, {
   3.568378459308427*^9, 3.5683784812420654`*^9}, {3.568378513440522*^9, 
   3.568378867046343*^9}, {3.568378915468828*^9, 3.56837896203491*^9}, {
   3.5683790764299107`*^9, 3.568379270041851*^9}, {3.568380233530743*^9, 
   3.568380234388745*^9}, 3.5683804214954734`*^9, 3.568380471867962*^9, 
   3.5683805399464817`*^9, {3.5683812396077104`*^9, 3.568381293537005*^9}, {
   3.568381339541486*^9, 3.5683813935175805`*^9}, 3.5684646396646404`*^9, 
   3.5684646988823442`*^9, 3.568465086387025*^9, {3.5684655859935026`*^9, 
   3.568465605087936*^9}, 3.568821235038273*^9, 3.5688233741452303`*^9, 
   3.568823449742963*^9, 3.5688274593312387`*^9, {3.5688284679978104`*^9, 
   3.5688284689026117`*^9}, {3.5688285131598897`*^9, 3.5688285201487017`*^9}, 
   3.5688285866360188`*^9, {3.5688292397623663`*^9, 3.568829246688778*^9}, {
   3.5688780328231993`*^9, 3.5688781383261843`*^9}, 3.5688953148957367`*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"wzh", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"If", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"u", "[", 
         RowBox[{"[", 
          RowBox[{"a", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "\[Equal]", 
        "\"\<H\>\""}], ",", 
       RowBox[{
        RowBox[{
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"a", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "3"}], "]"}], "]"}], "+", 
            RowBox[{"u2", "[", 
             RowBox[{"[", 
              RowBox[{"a", "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ")"}], 
          "/", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"u2", "[", 
             RowBox[{"[", 
              RowBox[{"a", "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "+", 
            RowBox[{"a", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "3"}], "]"}], "]"}], "+", 
            RowBox[{"u2", "[", 
             RowBox[{"[", 
              RowBox[{"a", "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], 
           ")"}]}], "*", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"z", "[", 
            RowBox[{"[", 
             RowBox[{"a", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "-", 
           RowBox[{"z", "[", 
            RowBox[{"[", 
             RowBox[{"a", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ")"}]}],
         "+", 
        RowBox[{"z", "[", 
         RowBox[{"[", 
          RowBox[{"a", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ",", "0"}],
       "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "nb"}], "}"}]}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.5681365021632137`*^9, 3.5681365171392403`*^9}, {
   3.568136560351316*^9, 3.5681365863097615`*^9}, {3.5681369135895367`*^9, 
   3.5681369265999594`*^9}, {3.568137821916732*^9, 3.5681378378443604`*^9}, {
   3.5682072492149415`*^9, 3.5682072529901485`*^9}, 3.568209802830227*^9, {
   3.5682216503836985`*^9, 3.5682217583046885`*^9}, {3.5682217912831464`*^9, 
   3.568221843465238*^9}, {3.5682219760654707`*^9, 3.5682221133457117`*^9}, {
   3.568222159615393*^9, 3.5682221693498106`*^9}, {3.5682222565227633`*^9, 
   3.5682222578487663`*^9}, 3.5682224736283445`*^9, 3.5682227227607822`*^9, 
   3.5682227558640404`*^9, {3.568222818810151*^9, 3.5682228486062036`*^9}, {
   3.5682230057764797`*^9, 3.5682230079916835`*^9}, {3.568223041952943*^9, 
   3.5682230548229656`*^9}, 3.5682232619601297`*^9, 3.5682233000709963`*^9, {
   3.568223536645412*^9, 3.5682235381118145`*^9}, {3.568223591151908*^9, 
   3.568223591526308*^9}, {3.5682236491528096`*^9, 3.568223657483224*^9}, {
   3.568223713222122*^9, 3.5682237532361927`*^9}, {3.568223800379475*^9, 
   3.568223827648323*^9}, {3.5682238887380304`*^9, 3.5682238893152313`*^9}, {
   3.5682240287170763`*^9, 3.5682240612431335`*^9}, 3.5682241088076167`*^9, {
   3.568224198772975*^9, 3.5682242396450467`*^9}, {3.568224282170721*^9, 
   3.5682243185031853`*^9}, {3.568224437890195*^9, 3.568224439590598*^9}, {
   3.568224572253231*^9, 3.5682245867144566`*^9}, {3.5682247081762695`*^9, 
   3.568224733557514*^9}, {3.56822480520844*^9, 3.5682248433505073`*^9}, {
   3.568224884206979*^9, 3.5682248959069996`*^9}, {3.568224982830352*^9, 
   3.5682250047639904`*^9}, {3.56822503870965*^9, 3.5682250803305235`*^9}, {
   3.5682251223413973`*^9, 3.568225141451431*^9}, {3.568225234177994*^9, 
   3.5682252929588966`*^9}, {3.5682253775890455`*^9, 3.568225482140429*^9}, {
   3.56822551694409*^9, 3.5682255297205124`*^9}, {3.568225573197789*^9, 
   3.5682255809666023`*^9}, 3.568225660698343*^9, 3.5682257301964645`*^9, {
   3.568225765078126*^9, 3.5682257828933573`*^9}, {3.568225927026843*^9, 
   3.568225927416844*^9}, {3.568267062420238*^9, 3.568267081358671*^9}, {
   3.568267188328059*^9, 3.5682672830358257`*^9}, {3.5682676367040462`*^9, 
   3.568267747979042*^9}, 3.5682678404716043`*^9, 3.568267882030078*^9, {
   3.5682679252577534`*^9, 3.568267996612279*^9}, {3.5682680488879704`*^9, 
   3.568268061133992*^9}, {3.568268100914062*^9, 3.5682681150944867`*^9}, {
   3.5682681672609787`*^9, 3.5682681709113846`*^9}, {3.568268272904364*^9, 
   3.5682683182536435`*^9}, {3.5682685202739983`*^9, 3.568268599194537*^9}, {
   3.5682686348873997`*^9, 3.568268637773405*^9}, {3.5682686944951043`*^9, 
   3.568268872273017*^9}, {3.5682689420207396`*^9, 3.568269016947671*^9}, {
   3.5682695040741262`*^9, 3.568269601153097*^9}, {3.5682696652380095`*^9, 
   3.5682696728352227`*^9}, {3.5682697489633565`*^9, 
   3.5682698308635006`*^9}, {3.568271156491429*^9, 3.568271160063835*^9}, {
   3.568271244927984*^9, 3.568271369010602*^9}, {3.5682716250070515`*^9, 
   3.5682717312276382`*^9}, {3.5682717843301315`*^9, 
   3.5682718576814604`*^9}, {3.568271940673606*^9, 3.5682720536490045`*^9}, {
   3.5682722259669075`*^9, 3.5682723223126764`*^9}, {3.568272370813162*^9, 
   3.568272411466833*^9}, {3.5682725893539457`*^9, 3.568272622238803*^9}, {
   3.5682726936713285`*^9, 3.5682726960113325`*^9}, {3.568272741079812*^9, 
   3.56827275143823*^9}, {3.5682729218373294`*^9, 3.568272991070251*^9}, {
   3.5682734285574193`*^9, 3.5682734302890224`*^9}, {3.5682734898811274`*^9, 
   3.5682735111595645`*^9}, {3.5682736115925407`*^9, 
   3.5682736313577757`*^9}, {3.5682738601945777`*^9, 3.568273987272401*^9}, {
   3.5682763944346285`*^9, 3.5682764313910933`*^9}, 3.5682765128232365`*^9, {
   3.5682766527398825`*^9, 3.568276660976697*^9}, {3.5682774496048822`*^9, 
   3.568277461663703*^9}, {3.5683784032731285`*^9, 3.568378422991563*^9}, {
   3.568378459308427*^9, 3.5683784812420654`*^9}, {3.568378513440522*^9, 
   3.568378867046343*^9}, {3.568378915468828*^9, 3.56837896203491*^9}, {
   3.5683790764299107`*^9, 3.568379270041851*^9}, {3.568380233530743*^9, 
   3.568380234388745*^9}, 3.5683804214954734`*^9, 3.568380471867962*^9, 
   3.5683805399464817`*^9, {3.5683812396077104`*^9, 3.568381293537005*^9}, {
   3.568381339541486*^9, 3.5683813935175805`*^9}, 3.5684646396646404`*^9, 
   3.5684646988823442`*^9, 3.568465086387025*^9, {3.5684655859935026`*^9, 
   3.568465605087936*^9}, 3.568821235038273*^9, 3.5688233741452303`*^9, 
   3.568823449742963*^9, 3.5688274593312387`*^9, {3.5688284679978104`*^9, 
   3.5688284689026117`*^9}, {3.5688285131598897`*^9, 3.5688285201487017`*^9}, 
   3.5688285866360188`*^9, {3.5688292397623663`*^9, 3.568829246688778*^9}, {
   3.5688780328231993`*^9, 3.5688781532398105`*^9}, 3.568895317251341*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"Plot8", "=", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"Graphics3D", "[", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"Opacity", "[", "0.5", "]"}], ",", 
         RowBox[{"If", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"wxh", "[", 
             RowBox[{"[", "i", "]"}], "]"}], "\[NotEqual]", "0"}], ",", 
           RowBox[{"Sphere", "[", 
            RowBox[{
             RowBox[{"{", 
              RowBox[{
               RowBox[{"wxh", "[", 
                RowBox[{"[", "i", "]"}], "]"}], ",", 
               RowBox[{"wyh", "[", 
                RowBox[{"[", "i", "]"}], "]"}], ",", 
               RowBox[{"wzh", "[", 
                RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", 
             RowBox[{"a", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "3"}], "]"}], "]"}]}], "]"}]}], "]"}]}], 
        "}"}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "nb"}], "}"}]}], "]"}]}], ";"}], 
  "\[IndentingNewLine]"}]], "Input",
 CellChangeTimes->{{3.5681365021632137`*^9, 3.5681365171392403`*^9}, {
   3.568136560351316*^9, 3.5681365863097615`*^9}, {3.5681369135895367`*^9, 
   3.5681369265999594`*^9}, {3.568137821916732*^9, 3.5681378378443604`*^9}, {
   3.5682072492149415`*^9, 3.5682072529901485`*^9}, 3.568209802830227*^9, {
   3.5682216503836985`*^9, 3.5682217583046885`*^9}, {3.5682217912831464`*^9, 
   3.568221843465238*^9}, {3.5682219760654707`*^9, 3.5682221133457117`*^9}, {
   3.568222159615393*^9, 3.5682221693498106`*^9}, {3.5682222565227633`*^9, 
   3.5682222578487663`*^9}, 3.5682224736283445`*^9, 3.5682227227607822`*^9, 
   3.5682227558640404`*^9, {3.568222818810151*^9, 3.5682228486062036`*^9}, {
   3.5682230057764797`*^9, 3.5682230079916835`*^9}, {3.568223041952943*^9, 
   3.5682230548229656`*^9}, 3.5682232619601297`*^9, 3.5682233000709963`*^9, {
   3.568223536645412*^9, 3.5682235381118145`*^9}, {3.568223591151908*^9, 
   3.568223591526308*^9}, {3.5682236491528096`*^9, 3.568223657483224*^9}, {
   3.568223713222122*^9, 3.5682237532361927`*^9}, {3.568223800379475*^9, 
   3.568223827648323*^9}, {3.5682238887380304`*^9, 3.5682238893152313`*^9}, {
   3.5682240287170763`*^9, 3.5682240612431335`*^9}, 3.5682241088076167`*^9, {
   3.568224198772975*^9, 3.5682242396450467`*^9}, {3.568224282170721*^9, 
   3.5682243185031853`*^9}, {3.568224437890195*^9, 3.568224439590598*^9}, {
   3.568224572253231*^9, 3.5682245867144566`*^9}, {3.5682247081762695`*^9, 
   3.568224733557514*^9}, {3.56822480520844*^9, 3.5682248433505073`*^9}, {
   3.568224884206979*^9, 3.5682248959069996`*^9}, {3.568224982830352*^9, 
   3.5682250047639904`*^9}, {3.56822503870965*^9, 3.5682250803305235`*^9}, {
   3.5682251223413973`*^9, 3.568225141451431*^9}, {3.568225234177994*^9, 
   3.5682252929588966`*^9}, {3.5682253775890455`*^9, 3.568225482140429*^9}, {
   3.56822551694409*^9, 3.5682255297205124`*^9}, {3.568225573197789*^9, 
   3.5682255809666023`*^9}, 3.568225660698343*^9, 3.5682257301964645`*^9, {
   3.568225765078126*^9, 3.5682257828933573`*^9}, {3.568225927026843*^9, 
   3.568225927416844*^9}, {3.568267062420238*^9, 3.568267081358671*^9}, {
   3.568267188328059*^9, 3.5682672830358257`*^9}, {3.5682676367040462`*^9, 
   3.568267747979042*^9}, 3.5682678404716043`*^9, 3.568267882030078*^9, {
   3.5682679252577534`*^9, 3.568267996612279*^9}, {3.5682680488879704`*^9, 
   3.568268061133992*^9}, {3.568268100914062*^9, 3.5682681150944867`*^9}, {
   3.5682681672609787`*^9, 3.5682681709113846`*^9}, {3.568268272904364*^9, 
   3.5682683182536435`*^9}, {3.5682685202739983`*^9, 3.568268599194537*^9}, {
   3.5682686348873997`*^9, 3.568268637773405*^9}, {3.5682686944951043`*^9, 
   3.568268872273017*^9}, {3.5682689420207396`*^9, 3.568269016947671*^9}, {
   3.5682695040741262`*^9, 3.568269601153097*^9}, {3.5682696652380095`*^9, 
   3.5682696728352227`*^9}, {3.5682697489633565`*^9, 
   3.5682698308635006`*^9}, {3.568271156491429*^9, 3.568271160063835*^9}, {
   3.568271244927984*^9, 3.568271369010602*^9}, {3.5682716250070515`*^9, 
   3.5682717312276382`*^9}, {3.5682717843301315`*^9, 
   3.5682718576814604`*^9}, {3.568271940673606*^9, 3.5682720536490045`*^9}, {
   3.5682722259669075`*^9, 3.5682723223126764`*^9}, {3.568272370813162*^9, 
   3.568272411466833*^9}, {3.5682725893539457`*^9, 3.568272622238803*^9}, {
   3.5682726936713285`*^9, 3.5682726960113325`*^9}, {3.568272741079812*^9, 
   3.56827275143823*^9}, {3.5682729218373294`*^9, 3.568272991070251*^9}, {
   3.5682734285574193`*^9, 3.5682734302890224`*^9}, {3.5682734898811274`*^9, 
   3.5682735111595645`*^9}, {3.5682736115925407`*^9, 
   3.5682736313577757`*^9}, {3.5682738601945777`*^9, 3.568273987272401*^9}, {
   3.5682763944346285`*^9, 3.5682764313910933`*^9}, 3.5682765128232365`*^9, {
   3.5682766527398825`*^9, 3.568276660976697*^9}, {3.5682774496048822`*^9, 
   3.568277461663703*^9}, {3.5683784032731285`*^9, 3.568378422991563*^9}, {
   3.568378459308427*^9, 3.5683784812420654`*^9}, {3.568378513440522*^9, 
   3.568378867046343*^9}, {3.568378915468828*^9, 3.56837896203491*^9}, {
   3.5683790764299107`*^9, 3.568379270041851*^9}, {3.568380233530743*^9, 
   3.568380234388745*^9}, 3.5683804214954734`*^9, 3.568380471867962*^9, 
   3.5683805399464817`*^9, {3.5683812396077104`*^9, 3.568381293537005*^9}, {
   3.568381339541486*^9, 3.5683813935175805`*^9}, 3.5684646396646404`*^9, 
   3.5684646988823442`*^9, 3.568465086387025*^9, {3.5684655859935026`*^9, 
   3.568465605087936*^9}, 3.568821235038273*^9, 3.5688233741452303`*^9, 
   3.568823449742963*^9, 3.5688274593312387`*^9, {3.5688284679978104`*^9, 
   3.5688284689026117`*^9}, {3.5688285131598897`*^9, 3.5688285201487017`*^9}, 
   3.5688285866360188`*^9, {3.5688292397623663`*^9, 3.568829246688778*^9}, {
   3.5688780328231993`*^9, 3.568878147967001*^9}, {3.56889532230575*^9, 
   3.5688953242557535`*^9}, 3.56889600529175*^9, 3.5688961595292206`*^9}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"Plot9", "=", " ", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"Graphics3D", "[", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"Darker", "[", 
         RowBox[{"Red", ",", "0.8"}], "]"}], ",", 
        RowBox[{"Cylinder", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{"x", "[", 
               RowBox[{"[", 
                RowBox[{"a", "[", 
                 RowBox[{"[", 
                  RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], ",", 
              RowBox[{"y", "[", 
               RowBox[{"[", 
                RowBox[{"a", "[", 
                 RowBox[{"[", 
                  RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], ",", 
              RowBox[{"z", "[", 
               RowBox[{"[", 
                RowBox[{"a", "[", 
                 RowBox[{"[", 
                  RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}]}], 
             "}"}], ",", 
            RowBox[{"{", 
             RowBox[{
              RowBox[{"x", "[", 
               RowBox[{"[", 
                RowBox[{"a", "[", 
                 RowBox[{"[", 
                  RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}], ",", 
              RowBox[{"y", "[", 
               RowBox[{"[", 
                RowBox[{"a", "[", 
                 RowBox[{"[", 
                  RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}], ",", 
              RowBox[{"z", "[", 
               RowBox[{"[", 
                RowBox[{"a", "[", 
                 RowBox[{"[", 
                  RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], 
             "}"}]}], "}"}], ",", "0.05"}], "]"}]}], "}"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "nb"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"Show", "[", 
  RowBox[{"Plot1", ",", "Plot7", ",", "Plot8", ",", "Plot9", ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"AspectRatio", "\[Rule]", "Automatic"}], ",", 
     RowBox[{"Boxed", "->", "False"}]}], "}"}]}], 
  "]"}], "\[IndentingNewLine]", 
 RowBox[{"Show", "[", 
  RowBox[{"Plot1", ",", "Plot2", ",", "Plot9", ",", 
   RowBox[{"{", 
    RowBox[{"AspectRatio", "\[Rule]", "Automatic"}], "}"}]}], "]"}]}], "Input",
 CellChangeTimes->{{3.5681365021632137`*^9, 3.5681365171392403`*^9}, {
   3.568136560351316*^9, 3.5681365863097615`*^9}, {3.5681369135895367`*^9, 
   3.5681369265999594`*^9}, {3.568137821916732*^9, 3.5681378378443604`*^9}, {
   3.5682072492149415`*^9, 3.5682072529901485`*^9}, 3.568209802830227*^9, {
   3.5682216503836985`*^9, 3.5682217583046885`*^9}, {3.5682217912831464`*^9, 
   3.568221843465238*^9}, {3.5682219760654707`*^9, 3.5682221133457117`*^9}, {
   3.568222159615393*^9, 3.5682221693498106`*^9}, {3.5682222565227633`*^9, 
   3.5682222578487663`*^9}, 3.5682224736283445`*^9, 3.5682227227607822`*^9, 
   3.5682227558640404`*^9, {3.568222818810151*^9, 3.5682228486062036`*^9}, {
   3.5682230057764797`*^9, 3.5682230079916835`*^9}, {3.568223041952943*^9, 
   3.5682230548229656`*^9}, 3.5682232619601297`*^9, 3.5682233000709963`*^9, {
   3.568223536645412*^9, 3.5682235381118145`*^9}, {3.568223591151908*^9, 
   3.568223591526308*^9}, {3.5682236491528096`*^9, 3.568223657483224*^9}, {
   3.568223713222122*^9, 3.5682237532361927`*^9}, {3.568223800379475*^9, 
   3.568223827648323*^9}, {3.5682238887380304`*^9, 3.5682238893152313`*^9}, {
   3.5682240287170763`*^9, 3.5682240612431335`*^9}, 3.5682241088076167`*^9, {
   3.568224198772975*^9, 3.5682242396450467`*^9}, {3.568224282170721*^9, 
   3.5682243185031853`*^9}, {3.568224437890195*^9, 3.568224439590598*^9}, {
   3.568224572253231*^9, 3.5682245867144566`*^9}, {3.5682247081762695`*^9, 
   3.568224733557514*^9}, {3.56822480520844*^9, 3.5682248433505073`*^9}, {
   3.568224884206979*^9, 3.5682248959069996`*^9}, {3.568224982830352*^9, 
   3.5682250047639904`*^9}, {3.56822503870965*^9, 3.5682250803305235`*^9}, {
   3.5682251223413973`*^9, 3.568225141451431*^9}, {3.568225234177994*^9, 
   3.5682252929588966`*^9}, {3.5682253775890455`*^9, 3.568225482140429*^9}, {
   3.56822551694409*^9, 3.5682255297205124`*^9}, {3.568225573197789*^9, 
   3.5682255809666023`*^9}, 3.568225660698343*^9, 3.5682257301964645`*^9, {
   3.568225765078126*^9, 3.5682257828933573`*^9}, {3.568225927026843*^9, 
   3.568225927416844*^9}, {3.568267062420238*^9, 3.568267081358671*^9}, {
   3.568267188328059*^9, 3.5682672830358257`*^9}, {3.5682676367040462`*^9, 
   3.568267747979042*^9}, 3.5682678404716043`*^9, 3.568267882030078*^9, {
   3.5682679252577534`*^9, 3.568267996612279*^9}, {3.5682680488879704`*^9, 
   3.568268061133992*^9}, {3.568268100914062*^9, 3.5682681150944867`*^9}, {
   3.5682681672609787`*^9, 3.5682681709113846`*^9}, {3.568268272904364*^9, 
   3.5682683182536435`*^9}, {3.5682685202739983`*^9, 3.568268599194537*^9}, {
   3.5682686348873997`*^9, 3.568268637773405*^9}, {3.5682686944951043`*^9, 
   3.568268872273017*^9}, {3.5682689420207396`*^9, 3.568269016947671*^9}, {
   3.5682695040741262`*^9, 3.568269601153097*^9}, {3.5682696652380095`*^9, 
   3.5682696728352227`*^9}, {3.5682697489633565`*^9, 
   3.5682698308635006`*^9}, {3.568271156491429*^9, 3.568271160063835*^9}, {
   3.568271244927984*^9, 3.568271369010602*^9}, {3.5682716250070515`*^9, 
   3.5682717312276382`*^9}, {3.5682717843301315`*^9, 
   3.5682718576814604`*^9}, {3.568271940673606*^9, 3.5682720536490045`*^9}, {
   3.5682722259669075`*^9, 3.5682723223126764`*^9}, {3.568272370813162*^9, 
   3.568272411466833*^9}, {3.5682725893539457`*^9, 3.568272622238803*^9}, {
   3.5682726936713285`*^9, 3.5682726960113325`*^9}, {3.568272741079812*^9, 
   3.56827275143823*^9}, {3.5682729218373294`*^9, 3.568272991070251*^9}, {
   3.5682734285574193`*^9, 3.5682734302890224`*^9}, {3.5682734898811274`*^9, 
   3.5682735111595645`*^9}, {3.5682736115925407`*^9, 
   3.5682736313577757`*^9}, {3.5682738601945777`*^9, 3.568273987272401*^9}, {
   3.5682763944346285`*^9, 3.5682764313910933`*^9}, 3.5682765128232365`*^9, {
   3.5682766527398825`*^9, 3.568276660976697*^9}, {3.5682774496048822`*^9, 
   3.568277461663703*^9}, {3.5683784032731285`*^9, 3.568378422991563*^9}, {
   3.568378459308427*^9, 3.5683784812420654`*^9}, {3.568378513440522*^9, 
   3.568378867046343*^9}, {3.568378915468828*^9, 3.56837896203491*^9}, {
   3.5683790764299107`*^9, 3.568379270041851*^9}, {3.568380233530743*^9, 
   3.568380234388745*^9}, 3.5683804214954734`*^9, 3.568380471867962*^9, 
   3.5683805399464817`*^9, {3.5683812396077104`*^9, 3.568381293537005*^9}, {
   3.568381339541486*^9, 3.568381420224828*^9}, {3.5684645223056345`*^9, 
   3.5684645572340956`*^9}, {3.5684651558539467`*^9, 
   3.5684651618911576`*^9}, {3.5688209923798466`*^9, 3.568820997683856*^9}, {
   3.568821158754139*^9, 3.568821160750942*^9}, 3.56882222404881*^9, {
   3.568822649274357*^9, 3.5688226508031597`*^9}, 3.568822723109287*^9}],

Cell[BoxData[
 Graphics3DBox[{
   {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], 
    SphereBox[{1.305415161896, 0.246471293002, -0.553079336502}, 0.255]}, 
   {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], 
    SphereBox[{-1.305405713266, -0.24642971903, 0.5530264241739999}, 0.255]}, 
   {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], 
    SphereBox[{2.7502127341, -0.782858679746, 0.576912560814}, 0.1]}, 
   {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], 
    SphereBox[{1.360699096026, -0.427807510236, -2.54459164804}, 0.1]}, 
   {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], 
    SphereBox[{1.71326338612, 2.30909399392, -0.49999315370999997`}, 0.1]}, 
   {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], 
    SphereBox[{-2.7517623094199997`, 0.7821254660579999, -0.575640775216}, 
     0.1]}, 
   {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], 
    SphereBox[{-1.3610033419119998`, 0.42886386707, 2.54425149736}, 0.1]}, 
   {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], 
    SphereBox[{-1.711419013544, -2.30945304186, 0.49911443112}, 0.1]}, 
   {RGBColor[
    0.33333333333333337`, 0.33333333333333337`, 0.33333333333333337`], 
    {Opacity[0.5], 
     SphereBox[{4.724315000093071*^-6, 
      0.00002078698599999873, -0.00002645616399998829}, 
      1.1849933116283475`]}}, 
   {RGBColor[
    0.33333333333333337`, 0.33333333333333337`, 0.33333333333333337`], 
    {Opacity[0.5]}}, 
   {RGBColor[
    0.33333333333333337`, 0.33333333333333337`, 0.33333333333333337`], 
    {Opacity[0.5]}}, 
   {RGBColor[
    0.33333333333333337`, 0.33333333333333337`, 0.33333333333333337`], 
    {Opacity[0.5]}}, 
   {RGBColor[
    0.33333333333333337`, 0.33333333333333337`, 0.33333333333333337`], 
    {Opacity[0.5]}}, 
   {RGBColor[
    0.33333333333333337`, 0.33333333333333337`, 0.33333333333333337`], 
    {Opacity[0.5]}}, 
   {RGBColor[
    0.33333333333333337`, 0.33333333333333337`, 0.33333333333333337`], 
    {Opacity[0.5]}}, 
   {Opacity[0.5]}, 
   {Opacity[0.5], 
    SphereBox[{2.3591482390706986`, -0.5042491265771325, 0.2710567653637229}, 
     1.2799931097314476`]}, 
   {Opacity[0.5], 
    SphereBox[{
     1.3457353291281724`, -0.24529968052776002`, -2.0055466890844267`}, 
     1.2799901618360665`]}, 
   {Opacity[0.5], 
    SphereBox[{1.6028674028467824`, 1.7507851729787047`, -0.5143624729799201},
      1.2799286426226604`]}, 
   {Opacity[0.5], 
    SphereBox[{-2.360272332739936, 0.5037231367188904, -0.27014080602189083`},
      1.2799743090211648`]}, 
   {Opacity[0.5], 
    SphereBox[{-1.3459551096881213`, 0.2460867423902468, 2.0053001420404883`},
      1.2800520427696944`]}, 
   {Opacity[0.5], 
    SphereBox[{-1.601522783976179, -1.7510514216809103`, 0.5137068712224904}, 
     1.2799875536973593`]}, 
   {RGBColor[0.19999999999999996`, 0., 0.], 
    CylinderBox[{{-1.305405713266, -0.24642971903, 0.5530264241739999}, {
     1.305415161896, 0.246471293002, -0.553079336502}}, 0.05]}, 
   {RGBColor[0.19999999999999996`, 0., 0.], 
    CylinderBox[{{2.7502127341, -0.782858679746, 0.576912560814}, {
     1.305415161896, 0.246471293002, -0.553079336502}}, 0.05]}, 
   {RGBColor[0.19999999999999996`, 0., 0.], 
    CylinderBox[{{1.360699096026, -0.427807510236, -2.54459164804}, {
     1.305415161896, 0.246471293002, -0.553079336502}}, 0.05]}, 
   {RGBColor[0.19999999999999996`, 0., 0.], 
    CylinderBox[{{1.71326338612, 2.30909399392, -0.49999315370999997`}, {
     1.305415161896, 0.246471293002, -0.553079336502}}, 0.05]}, 
   {RGBColor[0.19999999999999996`, 0., 0.], 
    CylinderBox[{{-2.7517623094199997`, 
     0.7821254660579999, -0.575640775216}, {-1.305405713266, -0.24642971903, 
     0.5530264241739999}}, 0.05]}, 
   {RGBColor[0.19999999999999996`, 0., 0.], 
    CylinderBox[{{-1.3610033419119998`, 0.42886386707, 
     2.54425149736}, {-1.305405713266, -0.24642971903, 0.5530264241739999}}, 
     0.05]}, 
   {RGBColor[0.19999999999999996`, 0., 0.], 
    CylinderBox[{{-1.711419013544, -2.30945304186, 
     0.49911443112}, {-1.305405713266, -0.24642971903, 0.5530264241739999}}, 
     0.05]}},
  AspectRatio->Automatic,
  Boxed->False]], "Output",
 CellChangeTimes->{{3.5690806271651382`*^9, 3.569080654012786*^9}, {
   3.5690807608729734`*^9, 3.5690808419619155`*^9}, 3.5690808953608093`*^9, {
   3.569080926155264*^9, 3.5690809603817234`*^9}, {3.5690809941401825`*^9, 
   3.56908100939701*^9}, 3.569082172815853*^9, 3.569088595674734*^9, 
   3.5691725087367716`*^9, 3.5691760356713657`*^9, 3.5691760999278784`*^9, 
   3.5691762427617292`*^9, 3.5691764704909296`*^9}],

Cell[BoxData[
 Graphics3DBox[{
   {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], 
    SphereBox[{1.305415161896, 0.246471293002, -0.553079336502}, 0.255]}, 
   {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], 
    SphereBox[{-1.305405713266, -0.24642971903, 0.5530264241739999}, 0.255]}, 
   {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], 
    SphereBox[{2.7502127341, -0.782858679746, 0.576912560814}, 0.1]}, 
   {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], 
    SphereBox[{1.360699096026, -0.427807510236, -2.54459164804}, 0.1]}, 
   {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], 
    SphereBox[{1.71326338612, 2.30909399392, -0.49999315370999997`}, 0.1]}, 
   {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], 
    SphereBox[{-2.7517623094199997`, 0.7821254660579999, -0.575640775216}, 
     0.1]}, 
   {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], 
    SphereBox[{-1.3610033419119998`, 0.42886386707, 2.54425149736}, 0.1]}, 
   {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], 
    SphereBox[{-1.711419013544, -2.30945304186, 0.49911443112}, 0.1]}, 
   Text3DBox["1", {1.0054151618959999`, 0.246471293002, -0.553079336502}, {1, 
    0}], Text3DBox["2", {-1.605405713266, -0.24642971903, 
    0.5530264241739999}, {1, 0}], 
   Text3DBox["3", {2.4502127341, -0.782858679746, 0.576912560814}, {1, 0}], 
   Text3DBox["4", {1.0606990960259999`, -0.427807510236, -2.54459164804}, {1, 
    0}], Text3DBox["5", {1.41326338612, 
    2.30909399392, -0.49999315370999997`}, {1, 0}], 
   Text3DBox["6", {-3.0517623094199995`, 
    0.7821254660579999, -0.575640775216}, {1, 0}], 
   Text3DBox["7", {-1.6610033419119998`, 0.42886386707, 2.54425149736}, {1, 
    0}], Text3DBox["8", {-2.0114190135439998`, -2.30945304186, 
    0.49911443112}, {1, 0}], 
   {RGBColor[0.19999999999999996`, 0., 0.], 
    CylinderBox[{{-1.305405713266, -0.24642971903, 0.5530264241739999}, {
     1.305415161896, 0.246471293002, -0.553079336502}}, 0.05]}, 
   {RGBColor[0.19999999999999996`, 0., 0.], 
    CylinderBox[{{2.7502127341, -0.782858679746, 0.576912560814}, {
     1.305415161896, 0.246471293002, -0.553079336502}}, 0.05]}, 
   {RGBColor[0.19999999999999996`, 0., 0.], 
    CylinderBox[{{1.360699096026, -0.427807510236, -2.54459164804}, {
     1.305415161896, 0.246471293002, -0.553079336502}}, 0.05]}, 
   {RGBColor[0.19999999999999996`, 0., 0.], 
    CylinderBox[{{1.71326338612, 2.30909399392, -0.49999315370999997`}, {
     1.305415161896, 0.246471293002, -0.553079336502}}, 0.05]}, 
   {RGBColor[0.19999999999999996`, 0., 0.], 
    CylinderBox[{{-2.7517623094199997`, 
     0.7821254660579999, -0.575640775216}, {-1.305405713266, -0.24642971903, 
     0.5530264241739999}}, 0.05]}, 
   {RGBColor[0.19999999999999996`, 0., 0.], 
    CylinderBox[{{-1.3610033419119998`, 0.42886386707, 
     2.54425149736}, {-1.305405713266, -0.24642971903, 0.5530264241739999}}, 
     0.05]}, 
   {RGBColor[0.19999999999999996`, 0., 0.], 
    CylinderBox[{{-1.711419013544, -2.30945304186, 
     0.49911443112}, {-1.305405713266, -0.24642971903, 0.5530264241739999}}, 
     0.05]}},
  AspectRatio->Automatic]], "Output",
 CellChangeTimes->{{3.5690806271651382`*^9, 3.569080654012786*^9}, {
   3.5690807608729734`*^9, 3.5690808419619155`*^9}, 3.5690808953608093`*^9, {
   3.569080926155264*^9, 3.5690809603817234`*^9}, {3.5690809941401825`*^9, 
   3.56908100939701*^9}, 3.569082172815853*^9, 3.569088595674734*^9, 
   3.5691725087367716`*^9, 3.5691760356713657`*^9, 3.5691760999278784`*^9, 
   3.5691762427617292`*^9, 3.5691764705377293`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"twc", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"If", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"wxc", "[", 
         RowBox[{"[", "i", "]"}], "]"}], "\[NotEqual]", "0"}], ",", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"wxc", "[", 
          RowBox[{"[", "i", "]"}], "]"}], ",", 
         RowBox[{"wyc", "[", 
          RowBox[{"[", "i", "]"}], "]"}], ",", 
         RowBox[{"wzc", "[", 
          RowBox[{"[", "i", "]"}], "]"}]}], "}"}]}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "cw"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"twh", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"If", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"wxh", "[", 
         RowBox[{"[", 
          RowBox[{"i", "+", "cw"}], "]"}], "]"}], "\[NotEqual]", "0"}], ",", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"wxh", "[", 
          RowBox[{"[", 
           RowBox[{"i", "+", "cw"}], "]"}], "]"}], ",", 
         RowBox[{"wyh", "[", 
          RowBox[{"[", 
           RowBox[{"i", "+", "cw"}], "]"}], "]"}], ",", 
         RowBox[{"wzh", "[", 
          RowBox[{"[", 
           RowBox[{"i", "+", "cw"}], "]"}], "]"}]}], "}"}]}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "bh"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"tw", "=", 
   RowBox[{"Join", "[", 
    RowBox[{"twc", ",", "twh"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"MatrixForm", "[", "tw", "]"}], "\[IndentingNewLine]"}], "Input",
 CellChangeTimes->{{3.568895397419882*^9, 3.5688955157148895`*^9}, {
   3.568895549629349*^9, 3.5688955855718126`*^9}, {3.568895617317868*^9, 
   3.568895845234268*^9}, {3.5688960797194805`*^9, 3.5688961015127187`*^9}, {
   3.568896219074525*^9, 3.5688962276233397`*^9}, {3.568896262910602*^9, 
   3.5688962903198504`*^9}, 3.5688964250728865`*^9, 3.568896456912543*^9, {
   3.5688981540863237`*^9, 3.5688981869399815`*^9}, {3.5689021233572955`*^9, 
   3.5689021341525145`*^9}, {3.5689111169043026`*^9, 
   3.5689111218495116`*^9}, {3.568911235261711*^9, 3.568911242437723*^9}, 
   3.5689112948850155`*^9, {3.5689116510180407`*^9, 3.568911800419503*^9}, {
   3.568911871820829*^9, 3.5689119012892804`*^9}, {3.5689119993042526`*^9, 
   3.5689120006770554`*^9}, {3.568912144306507*^9, 3.568912165085744*^9}, {
   3.5689528310632005`*^9, 3.5689528569592457`*^9}, {3.5689534065794115`*^9, 
   3.5689534297298517`*^9}, {3.568954496214126*^9, 3.5689545054337425`*^9}, 
   3.5689569073699613`*^9, {3.56917645947731*^9, 3.569176464453719*^9}}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"4.724315000093071`*^-6", "0.00002078698599999873`", 
      RowBox[{"-", "0.00002645616399998829`"}]},
     {"2.3591482390706986`", 
      RowBox[{"-", "0.5042491265771325`"}], "0.2710567653637229`"},
     {"1.3457353291281724`", 
      RowBox[{"-", "0.24529968052776002`"}], 
      RowBox[{"-", "2.0055466890844267`"}]},
     {"1.6028674028467824`", "1.7507851729787047`", 
      RowBox[{"-", "0.5143624729799201`"}]},
     {
      RowBox[{"-", "2.360272332739936`"}], "0.5037231367188904`", 
      RowBox[{"-", "0.27014080602189083`"}]},
     {
      RowBox[{"-", "1.3459551096881213`"}], "0.2460867423902468`", 
      "2.0053001420404883`"},
     {
      RowBox[{"-", "1.601522783976179`"}], 
      RowBox[{"-", "1.7510514216809103`"}], "0.5137068712224904`"}
    },
    GridBoxAlignment->{
     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{3.5691764705689297`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"Vee", "=", "0.0"}], ";", " ", 
  RowBox[{"Vne", "=", "0.0"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"bw", "=", 
  RowBox[{"Length", "[", "tw", "]"}]}]}], "Input",
 CellChangeTimes->{{3.5688982251756487`*^9, 3.5688986440675845`*^9}, {
   3.5688986873108606`*^9, 3.56889871560931*^9}, {3.5688987970414534`*^9, 
   3.5688988526711507`*^9}, {3.568898887147211*^9, 3.568898888644814*^9}, {
   3.5688989206248703`*^9, 3.568899041759083*^9}, {3.5688990876855636`*^9, 
   3.568899258349863*^9}, 3.568899300267137*^9, {3.568899351326027*^9, 
   3.5688993831032825`*^9}, {3.5688994370481772`*^9, 3.568899468217032*^9}, {
   3.5688995813016305`*^9, 3.5688997165538683`*^9}, {3.5688998495129013`*^9, 
   3.568899932364647*^9}, {3.568899967667509*^9, 3.568900024529609*^9}, {
   3.5689000794261055`*^9, 3.568900136257005*^9}, {3.5689001767858763`*^9, 
   3.56890020149632*^9}, {3.568902068726*^9, 3.568902074139209*^9}, 
   3.56890216040736*^9, 3.568902344113283*^9, {3.5689024237982235`*^9, 
   3.56890243907065*^9}, 3.5689025770684924`*^9, 3.568910195129884*^9, 
   3.568910620995032*^9, {3.568912393532545*^9, 3.5689124022841606`*^9}}],

Cell[BoxData["7"], "Output",
 CellChangeTimes->{{3.5690806271807384`*^9, 3.5690806540907855`*^9}, {
   3.5690807609509735`*^9, 3.569080842039916*^9}, 3.5690808954388094`*^9, {
   3.5690809262332635`*^9, 3.5690809604597235`*^9}, {3.569080994218183*^9, 
   3.56908100947501*^9}, 3.569082172893853*^9, 3.5690885957371345`*^9, 
   3.5691725088147717`*^9, 3.569176035749366*^9, 3.5691761000058784`*^9, 
   3.5691762428397293`*^9, 3.5691764705689297`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"For", "[", 
  RowBox[{
   RowBox[{"i", "=", "1"}], ",", 
   RowBox[{"i", "<", "bw"}], ",", 
   RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
   RowBox[{"For", "[", 
    RowBox[{
     RowBox[{"j", "=", 
      RowBox[{"i", "+", "1"}]}], ",", 
     RowBox[{"j", "<", 
      RowBox[{"bw", "+", "1"}]}], ",", 
     RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", 
     RowBox[{"Vee", "=", 
      RowBox[{"Vee", "+", 
       RowBox[{"4", "/", 
        RowBox[{"Sqrt", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"(", 
            RowBox[{
             RowBox[{"tw", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "1"}], "]"}], "]"}], "-", 
             RowBox[{"tw", "[", 
              RowBox[{"[", 
               RowBox[{"j", ",", "1"}], "]"}], "]"}]}], ")"}], "^", "2"}], 
          "+", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{
             RowBox[{"tw", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "2"}], "]"}], "]"}], "-", 
             RowBox[{"tw", "[", 
              RowBox[{"[", 
               RowBox[{"j", ",", "2"}], "]"}], "]"}]}], ")"}], "^", "2"}], 
          "+", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{
             RowBox[{"tw", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "3"}], "]"}], "]"}], "-", 
             RowBox[{"tw", "[", 
              RowBox[{"[", 
               RowBox[{"j", ",", "3"}], "]"}], "]"}]}], ")"}], "^", "2"}]}], 
         "]"}]}]}]}]}], "]"}]}], 
  "]"}], "\[IndentingNewLine]", "Vee"}], "Input",
 CellChangeTimes->{{3.5688982251756487`*^9, 3.5688986440675845`*^9}, {
   3.5688986873108606`*^9, 3.56889871560931*^9}, {3.5688987970414534`*^9, 
   3.5688988526711507`*^9}, {3.568898887147211*^9, 3.568898888644814*^9}, {
   3.5688989206248703`*^9, 3.568899041759083*^9}, {3.5688990876855636`*^9, 
   3.568899258349863*^9}, 3.568899300267137*^9, {3.568899351326027*^9, 
   3.5688993831032825`*^9}, {3.5688994370481772`*^9, 3.568899468217032*^9}, {
   3.5688995813016305`*^9, 3.5688997165538683`*^9}, {3.5688998495129013`*^9, 
   3.568899932364647*^9}, {3.568899967667509*^9, 3.568900024529609*^9}, {
   3.5689000794261055`*^9, 3.568900136257005*^9}, {3.5689001767858763`*^9, 
   3.56890020149632*^9}, {3.568902068726*^9, 3.568902074139209*^9}, 
   3.56890216040736*^9, 3.568902344113283*^9, {3.5689024237982235`*^9, 
   3.56890243907065*^9}, 3.5689025770684924`*^9, {3.568910195129884*^9, 
   3.568910199825492*^9}, 3.568910617609826*^9, {3.56891263539537*^9, 
   3.568912638842976*^9}, {3.568952611212014*^9, 3.5689526227248344`*^9}, {
   3.5689527476498537`*^9, 3.5689527480086546`*^9}, {3.5689530637688093`*^9, 
   3.5689530642368097`*^9}}],

Cell[BoxData["27.70567990471211`"], "Output",
 CellChangeTimes->{{3.5690806271963387`*^9, 3.5690806540907855`*^9}, {
   3.5690807609665737`*^9, 3.569080842055516*^9}, 3.5690808954544096`*^9, {
   3.5690809262488637`*^9, 3.5690809604753237`*^9}, {3.5690809942337832`*^9, 
   3.56908100949061*^9}, 3.5690821729094534`*^9, 3.569088595752734*^9, 
   3.569172508830372*^9, 3.569176035764966*^9, 3.5691761000214787`*^9, 
   3.5691762428553295`*^9, 3.56917647058453*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"For", "[", 
  RowBox[{
   RowBox[{"i", "=", "1"}], ",", 
   RowBox[{"i", "<", "bc"}], ",", 
   RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
   RowBox[{"For", "[", 
    RowBox[{
     RowBox[{"j", "=", 
      RowBox[{"i", "+", "1"}]}], ",", 
     RowBox[{"j", "<", 
      RowBox[{"bc", "+", "1"}]}], ",", 
     RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", 
     RowBox[{"Vee", "=", 
      RowBox[{"Vee", "+", 
       RowBox[{"4", "/", 
        RowBox[{"Sqrt", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"(", 
            RowBox[{
             RowBox[{"t", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "1"}], "]"}], "]"}], "-", 
             RowBox[{"t", "[", 
              RowBox[{"[", 
               RowBox[{"j", ",", "1"}], "]"}], "]"}]}], ")"}], "^", "2"}], 
          "+", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{
             RowBox[{"t", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "2"}], "]"}], "]"}], "-", 
             RowBox[{"t", "[", 
              RowBox[{"[", 
               RowBox[{"j", ",", "2"}], "]"}], "]"}]}], ")"}], "^", "2"}], 
          "+", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{
             RowBox[{"t", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "3"}], "]"}], "]"}], "-", 
             RowBox[{"t", "[", 
              RowBox[{"[", 
               RowBox[{"j", ",", "3"}], "]"}], "]"}]}], ")"}], "^", "2"}]}], 
         "]"}]}]}]}]}], "]"}]}], 
  "]"}], "\[IndentingNewLine]", "Vee"}], "Input",
 CellChangeTimes->{{3.5688982251756487`*^9, 3.5688986440675845`*^9}, {
   3.5688986873108606`*^9, 3.56889871560931*^9}, {3.5688987970414534`*^9, 
   3.5688988526711507`*^9}, {3.568898887147211*^9, 3.568898888644814*^9}, {
   3.5688989206248703`*^9, 3.568899041759083*^9}, {3.5688990876855636`*^9, 
   3.568899258349863*^9}, 3.568899300267137*^9, {3.568899351326027*^9, 
   3.5688993831032825`*^9}, {3.5688994370481772`*^9, 3.568899468217032*^9}, {
   3.5688995813016305`*^9, 3.5688997165538683`*^9}, {3.5688998495129013`*^9, 
   3.568899932364647*^9}, {3.568899967667509*^9, 3.568900024529609*^9}, {
   3.5689000794261055`*^9, 3.568900136257005*^9}, {3.5689001767858763`*^9, 
   3.56890020149632*^9}, {3.568902068726*^9, 3.568902074139209*^9}, 
   3.56890216040736*^9, 3.568902344113283*^9, {3.5689024237982235`*^9, 
   3.56890243907065*^9}, 3.5689025770684924`*^9, {3.568910195129884*^9, 
   3.5689102048487005`*^9}, 3.5689106144742203`*^9, {3.568912644739786*^9, 
   3.56891264654939*^9}, {3.568952467364841*^9, 3.5689524768652573`*^9}, {
   3.5689525940994635`*^9, 3.568952604472802*^9}, {3.568954607114721*^9, 
   3.5689546170051384`*^9}}],

Cell[BoxData["29.09554043001005`"], "Output",
 CellChangeTimes->{{3.5690806271963387`*^9, 3.5690806541063857`*^9}, {
   3.5690807609665737`*^9, 3.5690808420711155`*^9}, 3.56908089547001*^9, {
   3.5690809262644634`*^9, 3.569080960490924*^9}, {3.569080994249383*^9, 
   3.5690810095062103`*^9}, 3.569082172925053*^9, 3.5690885957683344`*^9, 
   3.569172508845972*^9, 3.569176035780566*^9, 3.569176100037079*^9, 
   3.5691762428709297`*^9, 3.5691764706001296`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"For", "[", 
  RowBox[{
   RowBox[{"i", "=", "1"}], ",", 
   RowBox[{"i", "<", 
    RowBox[{"bw", "+", "1"}]}], ",", 
   RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
   RowBox[{"For", "[", 
    RowBox[{
     RowBox[{"j", "=", "1"}], ",", 
     RowBox[{"j", "<", 
      RowBox[{"bc", "+", "1"}]}], ",", 
     RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", 
     RowBox[{"Vee", "=", 
      RowBox[{"Vee", "+", 
       RowBox[{"4", "/", 
        RowBox[{"Sqrt", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"(", 
            RowBox[{
             RowBox[{"tw", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "1"}], "]"}], "]"}], "-", 
             RowBox[{"t", "[", 
              RowBox[{"[", 
               RowBox[{"j", ",", "1"}], "]"}], "]"}]}], ")"}], "^", "2"}], 
          "+", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{
             RowBox[{"tw", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "2"}], "]"}], "]"}], "-", 
             RowBox[{"t", "[", 
              RowBox[{"[", 
               RowBox[{"j", ",", "2"}], "]"}], "]"}]}], ")"}], "^", "2"}], 
          "+", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{
             RowBox[{"tw", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "3"}], "]"}], "]"}], "-", 
             RowBox[{"t", "[", 
              RowBox[{"[", 
               RowBox[{"j", ",", "3"}], "]"}], "]"}]}], ")"}], "^", "2"}]}], 
         "]"}]}]}]}]}], "]"}]}], 
  "]"}], "\[IndentingNewLine]", "Vee"}], "Input",
 CellChangeTimes->{{3.5688982251756487`*^9, 3.5688986440675845`*^9}, {
   3.5688986873108606`*^9, 3.56889871560931*^9}, {3.5688987970414534`*^9, 
   3.5688988526711507`*^9}, {3.568898887147211*^9, 3.568898888644814*^9}, {
   3.5688989206248703`*^9, 3.568899041759083*^9}, {3.5688990876855636`*^9, 
   3.568899258349863*^9}, 3.568899300267137*^9, {3.568899351326027*^9, 
   3.5688993831032825`*^9}, {3.5688994370481772`*^9, 3.568899468217032*^9}, {
   3.5688995813016305`*^9, 3.5688997165538683`*^9}, {3.5688998495129013`*^9, 
   3.568899932364647*^9}, {3.568899967667509*^9, 3.568900024529609*^9}, {
   3.5689000794261055`*^9, 3.568900136257005*^9}, {3.5689001767858763`*^9, 
   3.56890020149632*^9}, {3.568902068726*^9, 3.568902074139209*^9}, 
   3.56890216040736*^9, 3.568902344113283*^9, {3.5689024237982235`*^9, 
   3.56890243907065*^9}, 3.5689025770684924`*^9, {3.568910195129884*^9, 
   3.5689102112915125`*^9}, {3.5689106097162123`*^9, 3.568910610371413*^9}, {
   3.568912654084203*^9, 3.5689126558314056`*^9}, 3.568952500358899*^9, 
   3.568952653737689*^9, {3.568953090039255*^9, 3.5689530906476564`*^9}, {
   3.568954639047977*^9, 3.5689546392507772`*^9}}],

Cell[BoxData["56.81346433123603`"], "Output",
 CellChangeTimes->{{3.5690806272119384`*^9, 3.5690806541219854`*^9}, {
   3.569080760997773*^9, 3.5690808420867157`*^9}, 3.5690808954856095`*^9, {
   3.5690809262800636`*^9, 3.569080960506524*^9}, {3.569080994264983*^9, 
   3.56908100952181*^9}, 3.569082172925053*^9, 3.5690885957683344`*^9, 
   3.5691725088615723`*^9, 3.569176035780566*^9, 3.569176100052679*^9, 
   3.56917624288653*^9, 3.5691764706001296`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"sc", "=", "0.2577"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"For", "[", 
  RowBox[{
   RowBox[{"i", "=", "1"}], " ", ",", 
   RowBox[{"i", "<", 
    RowBox[{"bc", "+", "1"}]}], ",", 
   RowBox[{"i", "++"}], ",", 
   RowBox[{"Vee", "=", 
    RowBox[{"Vee", "+", 
     RowBox[{"3", "*", 
      RowBox[{"sc", "/", 
       RowBox[{"u2", "[", 
        RowBox[{"[", "i", "]"}], "]"}]}]}]}]}]}], 
  "]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"sw", "=", "0.258"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"For", "[", 
  RowBox[{
   RowBox[{"i", "=", "1"}], " ", ",", 
   RowBox[{"i", "<", 
    RowBox[{"nb", "+", "1"}]}], ",", 
   RowBox[{"i", "++"}], ",", 
   RowBox[{"Vee", "=", 
    RowBox[{"Vee", "+", 
     RowBox[{"3", "*", 
      RowBox[{"sw", "/", 
       RowBox[{"a", "[", 
        RowBox[{"[", 
         RowBox[{"i", ",", "3"}], "]"}], "]"}]}]}]}]}]}], 
  "]"}], "\[IndentingNewLine]", "Vee", "\[IndentingNewLine]"}], "Input",
 CellChangeTimes->{{3.5689548890696163`*^9, 3.568955116299615*^9}, {
   3.5689552553894596`*^9, 3.5689552575266633`*^9}, {3.5689552936875267`*^9, 
   3.5689553203323736`*^9}, {3.568955404229321*^9, 3.568955433245372*^9}, 
   3.5689554882822685`*^9, 3.5689562157115464`*^9, {3.568956283961666*^9, 
   3.5689562858180695`*^9}, {3.5689566111878405`*^9, 
   3.5689566291278725`*^9}, {3.5689566695319433`*^9, 
   3.5689567220104356`*^9}, {3.569078870680053*^9, 3.569078873987259*^9}, {
   3.5690806436387672`*^9, 3.5690806493171773`*^9}, {3.569080750654955*^9, 
   3.569080836080705*^9}, {3.569080885922793*^9, 3.5690809556861153`*^9}, {
   3.5690809867613697`*^9, 3.569081004342601*^9}}],

Cell[BoxData["67.18219419652564`"], "Output",
 CellChangeTimes->{{3.5690806272119384`*^9, 3.5690806541219854`*^9}, {
   3.5690807610133734`*^9, 3.569080842102316*^9}, 3.5690808954856095`*^9, {
   3.569080926295664*^9, 3.5690809605221243`*^9}, {3.5690809942805834`*^9, 
   3.56908100952181*^9}, 3.5690821729406533`*^9, 3.5690885957839346`*^9, 
   3.5691725088615723`*^9, 3.569176035796166*^9, 3.569176100068279*^9, 
   3.5691762429021297`*^9, 3.5691764706157293`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"For", "[", 
  RowBox[{
   RowBox[{"i", "=", "1"}], ",", 
   RowBox[{"i", "<", 
    RowBox[{"bw", "+", "1"}]}], ",", 
   RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
   RowBox[{"For", "[", 
    RowBox[{
     RowBox[{"j", "=", "1"}], ",", 
     RowBox[{"j", "<", 
      RowBox[{"n", "+", "1"}]}], ",", 
     RowBox[{"j", "++"}], ",", 
     RowBox[{"If", "[", " ", 
      RowBox[{
       RowBox[{"i", "\[NotEqual]", "j"}], ",", 
       RowBox[{"Vne", "=", 
        RowBox[{"Vne", "-", 
         RowBox[{"2", "*", 
          RowBox[{
           RowBox[{"cn", "[", 
            RowBox[{"[", "j", "]"}], "]"}], "/", 
           RowBox[{"Sqrt", "[", 
            RowBox[{
             RowBox[{
              RowBox[{"(", 
               RowBox[{
                RowBox[{"tw", "[", 
                 RowBox[{"[", 
                  RowBox[{"i", ",", "1"}], "]"}], "]"}], "-", 
                RowBox[{"t", "[", 
                 RowBox[{"[", 
                  RowBox[{"j", ",", "1"}], "]"}], "]"}]}], ")"}], "^", "2"}], 
             "+", 
             RowBox[{
              RowBox[{"(", 
               RowBox[{
                RowBox[{"tw", "[", 
                 RowBox[{"[", 
                  RowBox[{"i", ",", "2"}], "]"}], "]"}], "-", 
                RowBox[{"t", "[", 
                 RowBox[{"[", 
                  RowBox[{"j", ",", "2"}], "]"}], "]"}]}], ")"}], "^", "2"}], 
             "+", 
             RowBox[{
              RowBox[{"(", 
               RowBox[{
                RowBox[{"tw", "[", 
                 RowBox[{"[", 
                  RowBox[{"i", ",", "3"}], "]"}], "]"}], "-", 
                RowBox[{"t", "[", 
                 RowBox[{"[", 
                  RowBox[{"j", ",", "3"}], "]"}], "]"}]}], ")"}], "^", 
              "2"}]}], "]"}]}]}]}]}], ",", "0"}], "]"}]}], "]"}]}], 
  "]"}], "\[IndentingNewLine]", "Vne"}], "Input",
 CellChangeTimes->{{3.5688982251756487`*^9, 3.5688986440675845`*^9}, {
   3.5688986873108606`*^9, 3.56889871560931*^9}, {3.5688987970414534`*^9, 
   3.5688988526711507`*^9}, {3.568898887147211*^9, 3.568898888644814*^9}, {
   3.5688989206248703`*^9, 3.568899041759083*^9}, {3.5688990876855636`*^9, 
   3.568899258349863*^9}, 3.568899300267137*^9, {3.568899351326027*^9, 
   3.5688993831032825`*^9}, {3.5688994370481772`*^9, 3.568899468217032*^9}, {
   3.5688995813016305`*^9, 3.5688997165538683`*^9}, {3.5688998495129013`*^9, 
   3.568899932364647*^9}, {3.568899967667509*^9, 3.568900024529609*^9}, {
   3.5689000794261055`*^9, 3.568900136257005*^9}, {3.5689001767858763`*^9, 
   3.56890020149632*^9}, {3.568902068726*^9, 3.568902074139209*^9}, 
   3.56890216040736*^9, 3.568902344113283*^9, {3.5689024237982235`*^9, 
   3.56890243907065*^9}, 3.5689025770684924`*^9, {3.568910195129884*^9, 
   3.568910216923122*^9}, {3.568910624395838*^9, 3.5689106253786397`*^9}, {
   3.568912441143829*^9, 3.568912441424629*^9}, 3.5689124745746875`*^9, {
   3.5689126628670177`*^9, 3.568912665066622*^9}, {3.568953170722597*^9, 
   3.568953194637439*^9}, {3.569079934960723*^9, 3.5690799595307655`*^9}, 
   3.5690800004340377`*^9, 3.569080091116997*^9, {3.569084566967658*^9, 
   3.5690845922397027`*^9}, {3.5691723518316965`*^9, 
   3.5691724380530477`*^9}, {3.5691724907811403`*^9, 
   3.5691724921227427`*^9}, {3.56917596963645*^9, 3.569175990337686*^9}, 
   3.5691760280741525`*^9, 3.5691760941870685`*^9, {3.569176194869645*^9, 
   3.5691762330273123`*^9}}],

Cell[BoxData[
 RowBox[{"-", "109.28584159078939`"}]], "Output",
 CellChangeTimes->{3.569176100068279*^9, 3.5691762429021297`*^9, 
  3.5691764706313295`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"For", "[", 
  RowBox[{
   RowBox[{"i", "=", "1"}], ",", 
   RowBox[{"i", "<", 
    RowBox[{"bc", "+", "1"}]}], ",", 
   RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
   RowBox[{"For", "[", 
    RowBox[{
     RowBox[{"j", "=", 
      RowBox[{"i", "+", "1"}]}], ",", 
     RowBox[{"j", "<", 
      RowBox[{"bc", "+", "1"}]}], ",", 
     RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", 
     RowBox[{"Vne", "=", 
      RowBox[{"Vne", "-", 
       RowBox[{"2", "*", 
        RowBox[{
         RowBox[{"cn", "[", 
          RowBox[{"[", "j", "]"}], "]"}], "/", 
         RowBox[{"Sqrt", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"(", 
             RowBox[{
              RowBox[{"t", "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", "1"}], "]"}], "]"}], "-", 
              RowBox[{"t", "[", 
               RowBox[{"[", 
                RowBox[{"j", ",", "1"}], "]"}], "]"}]}], ")"}], "^", "2"}], 
           "+", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{
              RowBox[{"t", "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", "2"}], "]"}], "]"}], "-", 
              RowBox[{"t", "[", 
               RowBox[{"[", 
                RowBox[{"j", ",", "2"}], "]"}], "]"}]}], ")"}], "^", "2"}], 
           "+", 
           RowBox[{
            RowBox[{"(", 
             RowBox[{
              RowBox[{"t", "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", "3"}], "]"}], "]"}], "-", 
              RowBox[{"t", "[", 
               RowBox[{"[", 
                RowBox[{"j", ",", "3"}], "]"}], "]"}]}], ")"}], "^", "2"}]}], 
          "]"}]}]}]}]}]}], "]"}]}], 
  "]"}], "\[IndentingNewLine]", "Vne"}], "Input",
 CellChangeTimes->{{3.5688982251756487`*^9, 3.5688986440675845`*^9}, {
   3.5688986873108606`*^9, 3.56889871560931*^9}, {3.5688987970414534`*^9, 
   3.5688988526711507`*^9}, {3.568898887147211*^9, 3.568898888644814*^9}, {
   3.5688989206248703`*^9, 3.568899041759083*^9}, {3.5688990876855636`*^9, 
   3.568899258349863*^9}, 3.568899300267137*^9, {3.568899351326027*^9, 
   3.5688993831032825`*^9}, {3.5688994370481772`*^9, 3.568899468217032*^9}, {
   3.5688995813016305`*^9, 3.5688997165538683`*^9}, {3.5688998495129013`*^9, 
   3.568899932364647*^9}, {3.568899967667509*^9, 3.568900024529609*^9}, {
   3.5689000794261055`*^9, 3.568900136257005*^9}, {3.5689001767858763`*^9, 
   3.56890020149632*^9}, {3.568902068726*^9, 3.568902074139209*^9}, 
   3.56890216040736*^9, 3.568902344113283*^9, {3.5689024237982235`*^9, 
   3.56890243907065*^9}, 3.5689025770684924`*^9, {3.568910195129884*^9, 
   3.5689102216343307`*^9}, 3.568910629060246*^9, {3.5689126699494305`*^9, 
   3.5689126747230387`*^9}, {3.5690803094237804`*^9, 3.569080333557023*^9}, {
   3.5690821082629395`*^9, 3.569082126140571*^9}}],

Cell[BoxData[
 RowBox[{"-", "113.45542316668322`"}]], "Output",
 CellChangeTimes->{{3.569080627227538*^9, 3.569080654153186*^9}, {
   3.569080761044574*^9, 3.5690808421335163`*^9}, 3.5690808955168095`*^9, {
   3.569080926326864*^9, 3.5690809605533237`*^9}, {3.5690809943117833`*^9, 
   3.56908100955301*^9}, 3.5690821729718533`*^9, 3.569088595815134*^9, 
   3.5691725090799723`*^9, 3.5691760359677663`*^9, 3.5691761000838785`*^9, 
   3.5691762429177294`*^9, 3.5691764706469297`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"For", "[", 
  RowBox[{
   RowBox[{"i", "=", "1"}], ",", 
   RowBox[{"i", "<", 
    RowBox[{"bc", "+", "1"}]}], ",", 
   RowBox[{"i", "++"}], ",", 
   RowBox[{"Vne", "=", 
    RowBox[{"Vne", "-", 
     RowBox[{"3", "*", 
      RowBox[{
       RowBox[{"cn", "[", 
        RowBox[{"[", "i", "]"}], "]"}], "/", 
       RowBox[{"u2", "[", 
        RowBox[{"[", "i", "]"}], "]"}]}]}]}]}]}], 
  "]"}], "\[IndentingNewLine]", "Vne"}], "Input",
 CellChangeTimes->{{3.5688982251756487`*^9, 3.5688986440675845`*^9}, {
   3.5688986873108606`*^9, 3.56889871560931*^9}, {3.5688987970414534`*^9, 
   3.5688988526711507`*^9}, {3.568898887147211*^9, 3.568898888644814*^9}, {
   3.5688989206248703`*^9, 3.568899041759083*^9}, {3.5688990876855636`*^9, 
   3.568899258349863*^9}, 3.568899300267137*^9, {3.568899351326027*^9, 
   3.5688993831032825`*^9}, {3.5688994370481772`*^9, 3.568899468217032*^9}, {
   3.5688995813016305`*^9, 3.5688997165538683`*^9}, {3.5688998495129013`*^9, 
   3.568899932364647*^9}, {3.568899967667509*^9, 3.568900024529609*^9}, {
   3.5689000794261055`*^9, 3.568900136257005*^9}, {3.5689001767858763`*^9, 
   3.56890020149632*^9}, {3.568902068726*^9, 3.568902074139209*^9}, 
   3.56890216040736*^9, 3.568902344113283*^9, {3.5689024237982235`*^9, 
   3.56890243907065*^9}, 3.5689025770684924`*^9, {3.568910195129884*^9, 
   3.5689102262831383`*^9}, {3.5689126851126575`*^9, 3.5689126875306616`*^9}, 
   3.5690795661448746`*^9}],

Cell[BoxData[
 RowBox[{"-", "255.18770663125017`"}]], "Output",
 CellChangeTimes->{{3.569080627227538*^9, 3.569080654168786*^9}, {
   3.5690807610601735`*^9, 3.569080842149116*^9}, 3.5690808955324097`*^9, {
   3.5690809263424635`*^9, 3.569080960568924*^9}, {3.569080994327383*^9, 
   3.56908100956861*^9}, 3.5690821729874535`*^9, 3.5690885958307343`*^9, 
   3.5691725090955725`*^9, 3.5691760359677663`*^9, 3.5691761000994787`*^9, 
   3.5691762429333296`*^9, 3.56917647066253*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"For", "[", 
  RowBox[{
   RowBox[{"i", "=", 
    RowBox[{"cw", "+", "1"}]}], ",", 
   RowBox[{"i", "<", " ", 
    RowBox[{"nb", "+", "1"}]}], ",", 
   RowBox[{"i", "++"}], ",", 
   RowBox[{"Vne", "=", 
    RowBox[{"Vne", "-", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"3", "-", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"u2", "[", 
            RowBox[{"[", 
             RowBox[{"a", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "/", 
           RowBox[{"a", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "3"}], "]"}], "]"}]}], ")"}], "^", "2"}]}], 
       ")"}], "/", 
      RowBox[{"a", "[", 
       RowBox[{"[", 
        RowBox[{"i", ",", "3"}], "]"}], "]"}]}]}]}]}], 
  "]"}], "\[IndentingNewLine]", "Vne"}], "Input",
 CellChangeTimes->{{3.5688982251756487`*^9, 3.5688986440675845`*^9}, {
   3.5688986873108606`*^9, 3.56889871560931*^9}, {3.5688987970414534`*^9, 
   3.5688988526711507`*^9}, {3.568898887147211*^9, 3.568898888644814*^9}, {
   3.5688989206248703`*^9, 3.568899041759083*^9}, {3.5688990876855636`*^9, 
   3.568899258349863*^9}, 3.568899300267137*^9, {3.568899351326027*^9, 
   3.5688993831032825`*^9}, {3.5688994370481772`*^9, 3.568899468217032*^9}, {
   3.5688995813016305`*^9, 3.5688997165538683`*^9}, {3.5688998495129013`*^9, 
   3.568899932364647*^9}, {3.568899967667509*^9, 3.568900024529609*^9}, {
   3.5689000794261055`*^9, 3.568900136257005*^9}, {3.5689001767858763`*^9, 
   3.56890020149632*^9}, {3.568902068726*^9, 3.568902074139209*^9}, 
   3.56890216040736*^9, 3.568902344113283*^9, {3.5689024237982235`*^9, 
   3.56890243907065*^9}, 3.5689025770684924`*^9, {3.568910195129884*^9, 
   3.5689102484819775`*^9}, {3.568910416821873*^9, 3.568910421314681*^9}, {
   3.5689104658527594`*^9, 3.568910491858005*^9}, {3.5689105421368933`*^9, 
   3.5689105444144974`*^9}, {3.568912692943871*^9, 3.56891270358309*^9}, 
   3.569074021101745*^9, {3.5690789214581423`*^9, 3.569078932050561*^9}, {
   3.5690789673222227`*^9, 3.569078983749052*^9}, {3.569079021095517*^9, 
   3.569079050376769*^9}, {3.5690791845214043`*^9, 3.5690792071882443`*^9}, {
   3.5690792779031687`*^9, 3.5690793371832724`*^9}, 3.5690793688825283`*^9, 
   3.5690795751460905`*^9, {3.569080188960369*^9, 3.5690801963079815`*^9}, {
   3.569080598476688*^9, 3.569080609318707*^9}}],

Cell[BoxData[
 RowBox[{"-", "268.32305093620107`"}]], "Output",
 CellChangeTimes->{{3.5690806272431383`*^9, 3.569080654184386*^9}, {
   3.5690807610601735`*^9, 3.569080842149116*^9}, 3.56908089554801*^9, {
   3.5690809263580637`*^9, 3.569080960568924*^9}, {3.569080994327383*^9, 
   3.56908100958421*^9}, 3.5690821730030537`*^9, 3.5690885958307343`*^9, 
   3.5691725091111727`*^9, 3.569176035983366*^9, 3.569176100115079*^9, 
   3.56917624294893*^9, 3.5691764706781297`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{"Vee", "\[IndentingNewLine]", "Vnn", "\[IndentingNewLine]", \
"Vne", "\[IndentingNewLine]", 
 RowBox[{"Vtot", "=", 
  RowBox[{
  "Vne", "+", "Vee", "+", 
   "Vnn"}]}], "\[IndentingNewLine]", "Ttot", "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"-", "Vtot"}], "/", "Ttot"}], "\[IndentingNewLine]", 
 RowBox[{"Etot", "=", 
  RowBox[{"Ttot", "+", "Vtot"}]}], "\[IndentingNewLine]"}], "Input",
 CellChangeTimes->{{3.5689074950860367`*^9, 3.5689075109980645`*^9}, {
  3.568910298074465*^9, 3.5689103015064707`*^9}, {3.568910724064413*^9, 
  3.5689107693044925`*^9}, {3.568955368910859*^9, 3.568955378629676*^9}, {
  3.568955601366867*^9, 3.568955616218093*^9}}],

Cell[BoxData["67.18219419652564`"], "Output",
 CellChangeTimes->{{3.5690806272431383`*^9, 3.5690806541999855`*^9}, {
   3.5690807610913734`*^9, 3.569080842180316*^9}, 3.56908089556361*^9, {
   3.569080926373664*^9, 3.569080960600124*^9}, {3.5690809943585835`*^9, 
   3.5690810095998096`*^9}, 3.569082173018654*^9, 3.5690885958463345`*^9, 
   3.5691725091267724`*^9, 3.569176035998966*^9, 3.569176100130679*^9, 
   3.56917624296453*^9, 3.5691764706937294`*^9}],

Cell[BoxData["41.956156943012274`"], "Output",
 CellChangeTimes->{{3.5690806272431383`*^9, 3.5690806541999855`*^9}, {
   3.5690807610913734`*^9, 3.569080842180316*^9}, 3.56908089556361*^9, {
   3.569080926373664*^9, 3.569080960600124*^9}, {3.5690809943585835`*^9, 
   3.5690810095998096`*^9}, 3.569082173018654*^9, 3.5690885958463345`*^9, 
   3.5691725091267724`*^9, 3.569176035998966*^9, 3.569176100130679*^9, 
   3.56917624296453*^9, 3.5691764706937294`*^9}],

Cell[BoxData[
 RowBox[{"-", "268.32305093620107`"}]], "Output",
 CellChangeTimes->{{3.5690806272431383`*^9, 3.5690806541999855`*^9}, {
   3.5690807610913734`*^9, 3.569080842180316*^9}, 3.56908089556361*^9, {
   3.569080926373664*^9, 3.569080960600124*^9}, {3.5690809943585835`*^9, 
   3.5690810095998096`*^9}, 3.569082173018654*^9, 3.5690885958463345`*^9, 
   3.5691725091267724`*^9, 3.569176035998966*^9, 3.569176100130679*^9, 
   3.56917624296453*^9, 3.5691764706937294`*^9}],

Cell[BoxData[
 RowBox[{"-", "159.18469979666315`"}]], "Output",
 CellChangeTimes->{{3.5690806272431383`*^9, 3.5690806541999855`*^9}, {
   3.5690807610913734`*^9, 3.569080842180316*^9}, 3.56908089556361*^9, {
   3.569080926373664*^9, 3.569080960600124*^9}, {3.5690809943585835`*^9, 
   3.5690810095998096`*^9}, 3.569082173018654*^9, 3.5690885958463345`*^9, 
   3.5691725091267724`*^9, 3.569176035998966*^9, 3.569176100130679*^9, 
   3.56917624296453*^9, 3.5691764706937294`*^9}],

Cell[BoxData["79.59237019992709`"], "Output",
 CellChangeTimes->{{3.5690806272431383`*^9, 3.5690806541999855`*^9}, {
   3.5690807610913734`*^9, 3.569080842180316*^9}, 3.56908089556361*^9, {
   3.569080926373664*^9, 3.569080960600124*^9}, {3.5690809943585835`*^9, 
   3.5690810095998096`*^9}, 3.569082173018654*^9, 3.5690885958463345`*^9, 
   3.5691725091267724`*^9, 3.569176035998966*^9, 3.569176100130679*^9, 
   3.56917624296453*^9, 3.5691764706937294`*^9}],

Cell[BoxData["1.999999489860763`"], "Output",
 CellChangeTimes->{{3.5690806272431383`*^9, 3.5690806541999855`*^9}, {
   3.5690807610913734`*^9, 3.569080842180316*^9}, 3.56908089556361*^9, {
   3.569080926373664*^9, 3.569080960600124*^9}, {3.5690809943585835`*^9, 
   3.5690810095998096`*^9}, 3.569082173018654*^9, 3.5690885958463345`*^9, 
   3.5691725091267724`*^9, 3.569176035998966*^9, 3.569176100130679*^9, 
   3.56917624296453*^9, 3.5691764706937294`*^9}],

Cell[BoxData[
 RowBox[{"-", "79.59232959673606`"}]], "Output",
 CellChangeTimes->{{3.5690806272431383`*^9, 3.5690806541999855`*^9}, {
   3.5690807610913734`*^9, 3.569080842180316*^9}, 3.56908089556361*^9, {
   3.569080926373664*^9, 3.569080960600124*^9}, {3.5690809943585835`*^9, 
   3.5690810095998096`*^9}, 3.569082173018654*^9, 3.5690885958463345`*^9, 
   3.5691725091267724`*^9, 3.569176035998966*^9, 3.569176100130679*^9, 
   3.56917624296453*^9, 3.5691764706937294`*^9}]
}, Open  ]]
},
WindowSize->{908, 939},
WindowMargins->{{Automatic, 38}, {-4, Automatic}},
PrintingCopies->1,
PrintingPageRange->{32000, 32000},
PrintingOptions->{"Magnification"->1.,
"PaperOrientation"->"Portrait",
"PaperSize"->{595.2, 841.8}},
Magnification->0.8999999761581421,
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, 1294, 28, 70, "Input"],
Cell[CellGroupData[{
Cell[1876, 52, 442, 7, 29, "Input"],
Cell[2321, 61, 1868, 35, 146, "Output"]
}, Open  ]],
Cell[4204, 99, 785, 16, 50, "Input"],
Cell[4992, 117, 789, 13, 29, "Input"],
Cell[CellGroupData[{
Cell[5806, 134, 875, 13, 29, "Input"],
Cell[6684, 149, 622, 10, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[7343, 164, 251, 5, 29, "Input"],
Cell[7597, 171, 446, 6, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[8080, 182, 1560, 29, 29, "Input"],
Cell[9643, 213, 622, 10, 28, "Output"]
}, Open  ]],
Cell[10280, 226, 1390, 22, 29, "Input"],
Cell[11673, 250, 1566, 28, 29, "Input"],
Cell[CellGroupData[{
Cell[13264, 282, 1234, 18, 29, "Input"],
Cell[14501, 302, 1009, 24, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[15547, 331, 697, 14, 29, "Input"],
Cell[16247, 347, 553, 10, 28, "Output"]
}, Open  ]],
Cell[16815, 360, 908, 24, 29, "Input"],
Cell[CellGroupData[{
Cell[17748, 388, 3671, 89, 130, "Input"],
Cell[21422, 479, 584, 10, 28, "Output"],
Cell[22009, 491, 584, 10, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[22630, 506, 519, 12, 29, "Input"],
Cell[23152, 520, 614, 10, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[23803, 535, 1070, 31, 29, "Input"],
Cell[24876, 568, 572, 10, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[25485, 583, 701, 21, 90, "Input"],
Cell[26189, 606, 444, 6, 28, "Output"],
Cell[26636, 614, 444, 6, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[27117, 625, 1428, 32, 110, "Input"],
Cell[28548, 659, 457, 6, 28, "Output"]
}, Open  ]],
Cell[29020, 668, 935, 19, 29, "Input"],
Cell[29958, 689, 3535, 66, 70, "Input"],
Cell[33496, 757, 3216, 56, 29, "Input"],
Cell[CellGroupData[{
Cell[36737, 817, 722, 11, 29, "Input"],
Cell[37462, 830, 4589, 111, 66, "Output"]
}, Open  ]],
Cell[42066, 944, 1254, 35, 70, "Input"],
Cell[43323, 981, 3847, 78, 90, "Input"],
Cell[CellGroupData[{
Cell[47195, 1063, 3292, 75, 130, "Input"],
Cell[50490, 1140, 598, 10, 28, "Output"]
}, Open  ]],
Cell[51103, 1153, 1708, 43, 90, "Input"],
Cell[52814, 1198, 511, 11, 50, "Input"],
Cell[53328, 1211, 2048, 48, 90, "Input"],
Cell[55379, 1261, 2180, 51, 70, "Input"],
Cell[57562, 1314, 1345, 34, 50, "Input"],
Cell[CellGroupData[{
Cell[58932, 1352, 923, 20, 29, "Input"],
Cell[59858, 1374, 467, 6, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[60362, 1385, 611, 17, 29, "Input"],
Cell[60976, 1404, 458, 6, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[61471, 1415, 507, 10, 29, "Input"],
Cell[61981, 1427, 461, 6, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[62479, 1438, 823, 23, 29, "Input"],
Cell[63305, 1463, 467, 6, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[63809, 1474, 204, 4, 29, "Input"],
Cell[64016, 1480, 461, 6, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[64514, 1491, 1679, 48, 110, "Input"],
Cell[66196, 1541, 466, 6, 28, "Output"]
}, Open  ]],
Cell[66677, 1550, 1971, 35, 90, "Input"],
Cell[68651, 1587, 1638, 24, 29, "Input"],
Cell[70292, 1613, 1655, 24, 29, "Input"],
Cell[71950, 1639, 2877, 57, 70, "Input"],
Cell[74830, 1698, 2696, 53, 90, "Input"],
Cell[CellGroupData[{
Cell[77551, 1755, 939, 17, 70, "Input"],
Cell[78493, 1774, 1087, 23, 28, "Output"],
Cell[79583, 1799, 450, 6, 28, "Output"],
Cell[80036, 1807, 450, 6, 28, "Output"]
}, Open  ]],
Cell[80501, 1816, 6921, 129, 90, "Input"],
Cell[87425, 1947, 7013, 131, 90, "Input"],
Cell[94441, 2080, 6921, 129, 90, "Input"],
Cell[101365, 2211, 5982, 97, 90, "Input"],
Cell[107350, 2310, 6882, 128, 90, "Input"],
Cell[114235, 2440, 6884, 128, 90, "Input"],
Cell[121122, 2570, 6882, 128, 90, "Input"],
Cell[128007, 2700, 5931, 96, 70, "Input"],
Cell[CellGroupData[{
Cell[133963, 2800, 6916, 125, 150, "Input"],
Cell[140882, 2927, 4582, 94, 351, "Output"],
Cell[145467, 3023, 3577, 61, 350, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[149081, 3089, 2643, 63, 110, "Input"],
Cell[151727, 3154, 1365, 34, 117, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[153129, 3193, 1163, 18, 50, "Input"],
Cell[154295, 3213, 449, 6, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[154781, 3224, 2746, 64, 90, "Input"],
Cell[157530, 3290, 463, 6, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[158030, 3301, 2740, 64, 90, "Input"],
Cell[160773, 3367, 462, 6, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[161272, 3378, 2764, 64, 90, "Input"],
Cell[164039, 3444, 459, 6, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[164535, 3455, 1658, 41, 130, "Input"],
Cell[166196, 3498, 466, 6, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[166699, 3509, 3439, 76, 110, "Input"],
Cell[170141, 3587, 156, 3, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[170334, 3595, 2844, 67, 90, "Input"],
Cell[173181, 3664, 482, 7, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[173700, 3676, 1457, 29, 50, "Input"],
Cell[175160, 3707, 479, 7, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[175676, 3719, 2411, 49, 50, "Input"],
Cell[178090, 3770, 474, 7, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[178601, 3782, 669, 13, 170, "Input"],
Cell[179273, 3797, 459, 6, 28, "Output"],
Cell[179735, 3805, 460, 6, 28, "Output"],
Cell[180198, 3813, 477, 7, 28, "Output"],
Cell[180678, 3822, 477, 7, 28, "Output"],
Cell[181158, 3831, 459, 6, 28, "Output"],
Cell[181620, 3839, 459, 6, 28, "Output"],
Cell[182082, 3847, 476, 7, 28, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)