(* 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[    126817,       2830]
NotebookOptionsPosition[    123275,       2716]
NotebookOutlinePosition[    123812,       2737]
CellTagsIndexPosition[    123769,       2734]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{
    "Kimball", " ", "model", " ", "and", " ", "energy", " ", "derived", " ", 
     "from", " ", "experimental", " ", "molecule", " ", "data", " ", 
     "source"}], ",", " ", 
    RowBox[{"ES", " ", "end", " ", "january", " ", "2013"}]}], " ", "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"Clear", "[", "Plot9", "]"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{"ChemicalData", "[", 
    RowBox[{"\"\<Butane\>\"", ",", "\"\<MoleculePlot\>\""}], 
    "]"}]}]}]], "Input",
 CellChangeTimes->{{3.568278170270548*^9, 3.568278174310955*^9}, {
   3.5682782109554195`*^9, 3.5682782183342323`*^9}, 3.568284752749406*^9, 
   3.5682847957898817`*^9, {3.5682872948362713`*^9, 3.568287369232802*^9}, {
   3.5683911898567343`*^9, 3.568391205612762*^9}, 3.5685665220166903`*^9, 
   3.568615530503583*^9}],

Cell[BoxData[
 Graphics3DBox[
  {AbsoluteThickness[3], EdgeForm[None], Specularity[
    GrayLevel[1], 100], GraphicsComplex3DBox[CompressedData["
1:eJxTTMoPSmViYGCQAGJmIH56Qen2zzZnhyMb9fIWC9oc2Bs0TbGv1fdAGgi0
OR84nLIqIUTIxuF5j8Zb3jZfhy9/r1S8fJXqkMmQ38iy19mh9sLmyK9LAw5A
xQ9kKOVUVO11PrDtM1BgaYCDrqL8l5w4L4dls48obDBLPKB+qG15+CW3A0vu
7+Obo2x4oObqqZb4VYoHZs0EAtXkA0Ii9sfuxHodEIiw3HLCLNFhS+TXnbcu
uTkc/qoR069k6HDufosTyxpFB8UNRRkTVZMdMia+rbHnzXVY9I55wq0jP/aL
rHN/WCWTfAAk/dYkxaFd7Oa574uTHJbuPSiu9j/4wL75Uvp3P+Q4sB7tNyy3
NXE4E+t9gr3XGmbOgT3r4jdIb/th75cEdIFMssNkCZYwPtOUA7JRKdb3Fycd
COidnif0P9jBDeTBDzkHJvJXma22MzmgFsy6eFKX9QHOnwvSN39yc7hX2NX3
ZEfIgSIg9Unb64AayN+vVB0c+bxmZBYbHfDS+6/4ty3qgA0owD+5HVixxM9B
eGeIg9/nvuASHS8HjuuLC2xfqR4IzfrmoF1i5BB1ec9jkfYoB5D3RGoyHb48
eLfm7jcjB1B0leZHHwjaIdf6mj/VYdE1k/cWk8MdBI/v2tEbFHQA6Pvgx7mZ
Dm8emUkduG3rIPLJ8Xxaq/EBqDkHPqm8nHXyu9EBzZj+Q1/zoh1cunOe/+ZP
PbA9qf6m7eTwA80HTi10DQpyAPs/N/OAgdZK4Qt3bA/csK2MWNFq7AAAVG0d
6A==
    "], {
     {RGBColor[0.65, 0.7, 0.7], SphereBox[5, 24.], SphereBox[6, 24.], 
      SphereBox[7, 24.], SphereBox[8, 24.], SphereBox[9, 24.], 
      SphereBox[10, 24.], SphereBox[11, 24.], SphereBox[12, 24.], 
      SphereBox[13, 24.], SphereBox[14, 24.]}, 
     {RGBColor[0.4, 0.4, 0.4], SphereBox[1, 34.], SphereBox[2, 34.], 
      SphereBox[3, 34.], SphereBox[4, 34.]}, 
     {RGBColor[0.65, 0.7, 0.7], CylinderBox[{15, 5}, 15.], 
      CylinderBox[{16, 6}, 15.], CylinderBox[{17, 7}, 15.], 
      CylinderBox[{18, 8}, 15.], CylinderBox[{19, 9}, 15.], 
      CylinderBox[{20, 10}, 15.], CylinderBox[{21, 11}, 15.], 
      CylinderBox[{22, 12}, 15.], CylinderBox[{23, 13}, 15.], 
      CylinderBox[{24, 14}, 15.]}, 
     {RGBColor[0.4, 0.4, 0.4], CylinderBox[{1, 2}, 15.], 
      CylinderBox[{1, 3}, 15.], CylinderBox[{1, 15}, 15.], 
      CylinderBox[{1, 16}, 15.], CylinderBox[{2, 4}, 15.], 
      CylinderBox[{2, 17}, 15.], CylinderBox[{2, 18}, 15.], 
      CylinderBox[{3, 19}, 15.], CylinderBox[{3, 20}, 15.], 
      CylinderBox[{3, 21}, 15.], CylinderBox[{4, 22}, 15.], 
      CylinderBox[{4, 23}, 15.], CylinderBox[{4, 24}, 15.]}}]},
  Boxed->False,
  ImageSize->175.42451257407296`,
  Lighting->"Neutral",
  ViewPoint->{1.3, -2.4, -2.}]], "Output",
 CellChangeTimes->{3.568628146262823*^9, 3.5687219049551935`*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"t", "=", 
    RowBox[{"0.01889726", "*", 
     RowBox[{"ChemicalData", "[", 
      RowBox[{"\"\<Butane\>\"", ",", "\"\<AtomPositions\>\""}], "]"}]}]}], 
   ";"}], " "}]], "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.5682847420789876`*^9, 3.568284795821082*^9, {
   3.568287382898426*^9, 3.5682873928356433`*^9}, 3.568566526821499*^9, 
   3.568615530503583*^9, {3.568623871740513*^9, 3.5686239073709755`*^9}, {
   3.5686239374634285`*^9, 3.5686239652158775`*^9}, {3.5686240254007826`*^9, 
   3.5686240291135893`*^9}, {3.568624068550459*^9, 3.5686241298273664`*^9}, 
   3.5686241856286645`*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"bt", "=", 
  RowBox[{"ChemicalData", "\[NoBreak]", "[", 
   RowBox[{"\"\<Butane\>\"", ",", "\"\<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.568284752702606*^9, 3.568284795821082*^9, 
   3.568566526821499*^9, 3.568615530519183*^9, {3.5686196140291557`*^9, 
   3.568619615745158*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"\<\"C\"\>", ",", "\<\"C\"\>"}], "}"}], ",", "1"}], "}"}], 
     ",", "3"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"\<\"C\"\>", ",", "\<\"H\"\>"}], "}"}], ",", "1"}], "}"}], 
     ",", "10"}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.568628146309623*^9, 3.568721905001993*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"\"\<C\>\"", ",", "\"\<C\>\""}], "}"}], ",", "1"}], "}"}], 
     ",", "3"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"\"\<C\>\"", ",", "\"\<H\>\""}], "}"}], ",", "1"}], "}"}], 
     ",", "10"}], "}"}]}], "}"}], "\[IndentingNewLine]", 
 RowBox[{"bc", "=", 
  RowBox[{"bt", "[", 
   RowBox[{"[", 
    RowBox[{"1", ",", "2"}], "]"}], "]"}]}], "\[IndentingNewLine]", 
 RowBox[{"bh", "=", 
  RowBox[{"bt", "[", 
   RowBox[{"[", 
    RowBox[{"2", ",", "2"}], "]"}], "]"}]}]}], "Input",
 CellChangeTimes->{{3.568619649566018*^9, 3.5686197006249075`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"\<\"C\"\>", ",", "\<\"C\"\>"}], "}"}], ",", "1"}], "}"}], 
     ",", "3"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"\<\"C\"\>", ",", "\<\"H\"\>"}], "}"}], ",", "1"}], "}"}], 
     ",", "10"}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.568628146325223*^9, 3.5687219050175934`*^9}],

Cell[BoxData["3"], "Output",
 CellChangeTimes->{3.568628146325223*^9, 3.5687219050175934`*^9}],

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

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"u", "=", 
  RowBox[{"ChemicalData", "\[NoBreak]", "[", 
   RowBox[{"\"\<Butane\>\"", ",", "\"\<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.568284752702606*^9, 3.568284795852282*^9, 
   3.568566526837099*^9, 3.568615530519183*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"\<\"C\"\>", ",", "\<\"C\"\>", ",", "\<\"C\"\>", ",", "\<\"C\"\>", 
   ",", "\<\"H\"\>", ",", "\<\"H\"\>", ",", "\<\"H\"\>", ",", "\<\"H\"\>", 
   ",", "\<\"H\"\>", ",", "\<\"H\"\>", ",", "\<\"H\"\>", ",", "\<\"H\"\>", 
   ",", "\<\"H\"\>", ",", "\<\"H\"\>"}], "}"}]], "Output",
 CellChangeTimes->{3.568628146356423*^9, 3.568721905033193*^9}]
}, Open  ]],

Cell[CellGroupData[{

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

Cell[BoxData["14"], "Output",
 CellChangeTimes->{3.5686281463720226`*^9, 3.5687219050487933`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"DeleteFile", "[", " ", "\"\<C:\\\\Kimb_Mathem\\\\tmp3\>\"", "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"For", "[", 
  RowBox[{
   RowBox[{"i", "=", "1"}], ",", 
   RowBox[{"i", "<", 
    RowBox[{"n", "+", "1"}]}], ",", 
   RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
      RowBox[{"u", "[", 
       RowBox[{"[", "i", "]"}], "]"}], ",", 
      RowBox[{"Extract", "[", 
       RowBox[{"t", ",", 
        RowBox[{"{", "i", "}"}]}], "]"}]}], "}"}], ">>>", " ", 
    "\"\<C:\\\\Kimb_Mathem\\\\tmp3\>\""}]}], "]"}], "\[IndentingNewLine]", 
 RowBox[{"u1", "=", 
  RowBox[{"Table", "[", 
   RowBox[{
    RowBox[{"If", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"u", "[", 
        RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\<C\>\""}], ",", 
      "0.2556", ",", 
      RowBox[{"If", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"u", "[", 
          RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\<H\>\""}], ",", 
        "0.1", ",", 
        RowBox[{"If", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"u", "[", 
            RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\<N\>\""}], ",", 
          "0.213", ",", 
          RowBox[{"If", "[", 
           RowBox[{
            RowBox[{
             RowBox[{"u", "[", 
              RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\<O\>\""}], ",", 
            "0.178"}], "]"}]}], "]"}]}], "]"}]}], "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], 
   "]"}]}], "\[IndentingNewLine]", 
 RowBox[{"u2", "=", 
  RowBox[{"Table", "[", 
   RowBox[{
    RowBox[{"If", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"u", "[", 
        RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\<C\>\""}], ",", 
      "0.2556", ",", 
      RowBox[{"If", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"u", "[", 
          RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\<H\>\""}], ",", 
        "0.564", ",", 
        RowBox[{"If", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"u", "[", 
            RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\<N\>\""}], ",", 
          "0.213", ",", 
          RowBox[{"If", "[", 
           RowBox[{
            RowBox[{
             RowBox[{"u", "[", 
              RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\<O\>\""}], ",", 
            "0.178"}], "]"}]}], "]"}]}], "]"}]}], "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}]}], "Input",
 CellChangeTimes->{{3.5686154384478216`*^9, 3.568615482112298*^9}, {
   3.5686157311667356`*^9, 3.568615736252344*^9}, {3.56861579063404*^9, 
   3.568615829727709*^9}, {3.56861597290476*^9, 3.568615973981162*^9}, {
   3.56861739044805*^9, 3.5686174033960724`*^9}, {3.5686241445225925`*^9, 
   3.568624159108618*^9}, {3.568624862622655*^9, 3.568624863137456*^9}, {
   3.5686270815141525`*^9, 3.5686273366682005`*^9}, {3.5686274562268105`*^9, 
   3.5686274982844844`*^9}, {3.5686278155110416`*^9, 3.568627853465908*^9}, 
   3.568628035393428*^9, {3.5687219594772887`*^9, 3.568721964718898*^9}}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"DeleteFile", "::", "nffil"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"File not found during \
\[NoBreak]\\!\\(DeleteFile[\\\"C:\\\\\\\\Kimb_Mathem\\\\\\\\tmp3\\\"]\\)\
\[NoBreak]. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \
ButtonStyle->\\\"Link\\\", ButtonFrame->None, \
ButtonData:>\\\"paclet:ref/message/General/nffil\\\", ButtonNote -> \
\\\"DeleteFile::nffil\\\"]\\)\"\>"}]], "Message", "MSG",
 CellChangeTimes->{3.5687219051579933`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "0.2556`", ",", "0.2556`", ",", "0.2556`", ",", "0.2556`", ",", "0.1`", ",",
    "0.1`", ",", "0.1`", ",", "0.1`", ",", "0.1`", ",", "0.1`", ",", "0.1`", 
   ",", "0.1`", ",", "0.1`", ",", "0.1`"}], "}"}]], "Output",
 CellChangeTimes->{3.568628146387623*^9, 3.5687219051735935`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "0.2556`", ",", "0.2556`", ",", "0.2556`", ",", "0.2556`", ",", "0.564`", 
   ",", "0.564`", ",", "0.564`", ",", "0.564`", ",", "0.564`", ",", "0.564`", 
   ",", "0.564`", ",", "0.564`", ",", "0.564`", ",", "0.564`"}], 
  "}"}]], "Output",
 CellChangeTimes->{3.568628146387623*^9, 3.5687219051735935`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"ChemicalData", "[", 
  RowBox[{"\"\<Butane\>\"", ",", "\"\<MolecularFormulaString\>\""}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.568284795883482*^9, 3.568284795883482*^9}, 
   3.568566526837099*^9, 3.568615530519183*^9}],

Cell[BoxData["\<\"CH3CH2CH2CH3\"\>"], "Output",
 CellChangeTimes->{3.568628146403223*^9, 3.5687219051891937`*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"v", "=", 
   RowBox[{"0.01889726", "*", 
    RowBox[{"ChemicalData", "\[NoBreak]", "[", 
     RowBox[{"\"\<Butane\>\"", ",", "\"\<VertexCoordinates\>\""}], "]"}]}]}], 
  ";"}]], "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.568284752733806*^9, 3.568284795899082*^9, 3.568566526852699*^9, 
   3.5686155305347834`*^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.5"}], ",", 
          RowBox[{"t", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "2"}], "]"}], "]"}], ",", 
          RowBox[{"t", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "3"}], "]"}], "]"}]}], "}"}], " ", ",", 
        RowBox[{"{", 
         RowBox[{"1", ",", "0"}], "}"}]}], "]"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 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.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.568615628331355*^9, {
   3.568628116965971*^9, 3.5686281181983733`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Show", "[", 
  RowBox[{"Plot1", ",", "Plot2", ",", 
   RowBox[{"{", 
    RowBox[{"AspectRatio", "\[Rule]", "Automatic"}], "}"}]}], "]"}]], "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.568120119762987*^9, 3.5681201279374013`*^9}, {3.5681212588145876`*^9, 
   3.5681212809978266`*^9}, 3.568121331261115*^9, {3.5682660696812944`*^9, 
   3.5682660793689113`*^9}, {3.5682661730782757`*^9, 3.5682661751998796`*^9}, 
   3.568266213092346*^9, {3.568266308392914*^9, 3.5682663102025166`*^9}, {
   3.568266350887388*^9, 3.5682663517297897`*^9}, 3.568266419574309*^9, {
   3.5683892670285573`*^9, 3.5683892701017623`*^9}}],

Cell[BoxData[
 Graphics3DBox[{
   {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], 
    SphereBox[{0.7380230406699999, -0.530425301214, -1.115758481084}, 
     0.2556]}, 
   {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], 
    SphereBox[{-0.7379380029999999, 0.530476323816, 1.115832180398}, 0.2556]}, 
   {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], 
    SphereBox[{3.31346446566, 0.746001463842, -1.2584271249060002`}, 0.2556]}, 
   {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], 
    SphereBox[{-3.31346446566, -0.7460694939779999, 1.25828350573}, 0.2556]}, 
   {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], 
    SphereBox[{0.9965980287019999, -2.60232277734, -0.862533307358}, 0.1]}, 
   {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], 
    SphereBox[{-0.323882028866, -0.1669149622376, -2.8945878005}, 0.1]}, 
   {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], 
    SphereBox[{-0.996514880758, 2.6023794691199997`, 0.862608896398}, 0.1]}, 
   {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], 
    SphereBox[{0.32380455010000003`, 0.1669756224422, 2.89456890324}, 0.1]}, 
   {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], 
    SphereBox[{4.39198778564, -0.0292540923156, -2.8892020813999997`}, 0.1]}, 
   {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], 
    SphereBox[{3.0548176680399997`, 
     2.8180350002399996`, -1.5115294664419998`}, 0.1]}, 
   {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], 
    SphereBox[{4.37522591602, 0.382482432126, 0.520636520534}, 0.1]}, 
   {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], 
    SphereBox[{-4.39198778564, 0.029186062179599997`, 2.88905090332}, 0.1]}, 
   {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], 
    SphereBox[{-3.05489325708, -2.81797830846, 1.5115067897299999`}, 0.1]}, 
   {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], 
    SphereBox[{-4.37520701876, -0.382571249248, -0.520455106838}, 0.1]}, 
   Text3DBox["1", {0.2380230406699999, -0.530425301214, -1.115758481084}, {1, 
    0}], Text3DBox["2", {-1.237938003, 0.530476323816, 1.115832180398}, {1, 
    0}], Text3DBox["3", {2.81346446566, 
    0.746001463842, -1.2584271249060002`}, {1, 0}], 
   Text3DBox["4", {-3.81346446566, -0.7460694939779999, 1.25828350573}, {1, 
    0}], Text3DBox["5", {
    0.49659802870199987`, -2.60232277734, -0.862533307358}, {1, 0}], 
   Text3DBox["6", {-0.823882028866, -0.1669149622376, -2.8945878005}, {1, 0}],
    Text3DBox["7", {-1.496514880758, 2.6023794691199997`, 0.862608896398}, {1,
     0}], Text3DBox["8", {-0.17619544989999997`, 0.1669756224422, 
    2.89456890324}, {1, 0}], 
   Text3DBox["9", {
    3.8919877856399996`, -0.0292540923156, -2.8892020813999997`}, {1, 0}], 
   Text3DBox["10", {2.5548176680399997`, 
    2.8180350002399996`, -1.5115294664419998`}, {1, 0}], 
   Text3DBox["11", {3.8752259160199998`, 0.382482432126, 0.520636520534}, {1, 
    0}], Text3DBox["12", {-4.89198778564, 0.029186062179599997`, 
    2.88905090332}, {1, 0}], 
   Text3DBox["13", {-3.55489325708, -2.81797830846, 1.5115067897299999`}, {1, 
    0}], Text3DBox["14", {-4.87520701876, -0.382571249248, -0.520455106838}, {
    1, 0}]},
  AspectRatio->Automatic]], "Output",
 CellChangeTimes->{3.568628146434423*^9, 3.5687219052359934`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"ChemicalData", "\[NoBreak]", "[", 
  RowBox[{"\"\<Butane\>\"", ",", "\"\<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.568284752733806*^9, 3.568284795914682*^9, 
   3.568566526852699*^9, 3.5686155305347834`*^9}],

Cell[BoxData[
 GraphicsBox[{{
    {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[
      StyleBox["\<\"C\"\>",
       StripOnInput->False,
       ScriptSizeMultipliers->0.8,
       ScriptMinSize->8,
       FontSize->Scaled[0.0770769541448646],
       FontWeight->Plain], {286.6, -25.}], 
     LineBox[{{303.9208810206263, -15.000645987399642`}, {329.905, 0.}}], 
     {GrayLevel[0], InsetBox[
       StyleBox["\<\"C\"\>",
        StripOnInput->False,
        ScriptSizeMultipliers->0.8,
        ScriptMinSize->8,
        FontSize->Scaled[0.0770769541448646],
        FontWeight->Plain], {373.21, 25.}], 
      LineBox[{{329.905, 0.}, {355.8891189793737, 15.000645987399652`}}]}}, 
    {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[
      StyleBox["\<\"C\"\>",
       StripOnInput->False,
       ScriptSizeMultipliers->0.8,
       ScriptMinSize->8,
       FontSize->Scaled[0.0770769541448646],
       FontWeight->Plain], {286.6, -25.}], 
     LineBox[{{269.27961894742526`, -14.999779992739736`}, {243.3, 0.}}], 
     {GrayLevel[0], InsetBox[
       StyleBox["\<\"C\"\>",
        StripOnInput->False,
        ScriptSizeMultipliers->0.8,
        ScriptMinSize->8,
        FontSize->Scaled[0.0770769541448646],
        FontWeight->Plain], {200., 25.}], 
      LineBox[{{243.3, 0.}, {217.3203810525748, 14.999779992739736`}}]}}, 
    {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[
      StyleBox["\<\"C\"\>",
       StripOnInput->False,
       ScriptSizeMultipliers->0.8,
       ScriptMinSize->8,
       FontSize->Scaled[0.0770769541448646],
       FontWeight->Plain], {286.6, -25.}], 
     LineBox[{{273.74402848350627`, -40.32070482605489}, {266.675, -48.745}}], 
     {GrayLevel[0], InsetBox[
       StyleBox["\<\"H\"\>",
        StripOnInput->False,
        ScriptSizeMultipliers->0.8,
        ScriptMinSize->8,
        FontSize->Scaled[0.0770769541448646],
        FontWeight->Plain], {246.74999999999997`, -72.49}], 
      LineBox[{{266.675, -48.745}, {
       259.6059715164937, -57.169295173945116`}}]}}, 
    {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[
      StyleBox["\<\"C\"\>",
       StripOnInput->False,
       ScriptSizeMultipliers->0.8,
       ScriptMinSize->8,
       FontSize->Scaled[0.0770769541448646],
       FontWeight->Plain], {286.6, -25.}], 
     LineBox[{{299.4578643243525, -40.31911632622926}, {
      306.53000000000003`, -48.745}}], 
     {GrayLevel[0], InsetBox[
       StyleBox["\<\"H\"\>",
        StripOnInput->False,
        ScriptSizeMultipliers->0.8,
        ScriptMinSize->8,
        FontSize->Scaled[0.0770769541448646],
        FontWeight->Plain], {326.46000000000004`, -72.49}], 
      LineBox[{{306.53000000000003`, -48.745}, {
       313.60213567564756`, -57.17088367377073}}]}}, 
    {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[
      StyleBox["\<\"C\"\>",
       StripOnInput->False,
       ScriptSizeMultipliers->0.8,
       ScriptMinSize->8,
       FontSize->Scaled[0.0770769541448646],
       FontWeight->Plain], {373.21, 25.}], 
     LineBox[{{390.53038105257474`, 14.999779992739736`}, {416.51, 0.}}], 
     {GrayLevel[0], InsetBox[
       StyleBox["\<\"C\"\>",
        StripOnInput->False,
        ScriptSizeMultipliers->0.8,
        ScriptMinSize->8,
        FontSize->Scaled[0.0770769541448646],
        FontWeight->Plain], {459.80999999999995`, -25.}], 
      LineBox[{{416.51, 0.}, {442.4896189474252, -14.999779992739725`}}]}}, 
    {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[
      StyleBox["\<\"C\"\>",
       StripOnInput->False,
       ScriptSizeMultipliers->0.8,
       ScriptMinSize->8,
       FontSize->Scaled[0.0770769541448646],
       FontWeight->Plain], {373.21, 25.}], 
     LineBox[{{386.06597151649373`, 40.32070482605489}, {393.135, 48.745}}], 
     {GrayLevel[0], InsetBox[
       StyleBox["\<\"H\"\>",
        StripOnInput->False,
        ScriptSizeMultipliers->0.8,
        ScriptMinSize->8,
        FontSize->Scaled[0.0770769541448646],
        FontWeight->Plain], {413.06, 72.49}], 
      LineBox[{{393.135, 48.745}, {400.20402848350625`, 57.1692951739451}}]}}, 
    {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[
      StyleBox["\<\"C\"\>",
       StripOnInput->False,
       ScriptSizeMultipliers->0.8,
       ScriptMinSize->8,
       FontSize->Scaled[0.0770769541448646],
       FontWeight->Plain], {373.21, 25.}], 
     LineBox[{{360.3521356756475, 40.31911632622926}, {353.28, 48.745}}], 
     {GrayLevel[0], InsetBox[
       StyleBox["\<\"H\"\>",
        StripOnInput->False,
        ScriptSizeMultipliers->0.8,
        ScriptMinSize->8,
        FontSize->Scaled[0.0770769541448646],
        FontWeight->Plain], {333.34999999999997`, 72.49}], 
      LineBox[{{353.28, 48.745}, {346.20786432435244`, 57.17088367377073}}]}}, 
    {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[
      StyleBox["\<\"C\"\>",
       StripOnInput->False,
       ScriptSizeMultipliers->0.8,
       ScriptMinSize->8,
       FontSize->Scaled[0.0770769541448646],
       FontWeight->Plain], {200., 25.}], 
     LineBox[{{210.00049938704143`, 42.320219744846895`}, {215.5, 51.845}}], 
     {GrayLevel[0], InsetBox[
       StyleBox["\<\"H\"\>",
        StripOnInput->False,
        ScriptSizeMultipliers->0.8,
        ScriptMinSize->8,
        FontSize->Scaled[0.0770769541448646],
        FontWeight->Plain], {231., 78.69}], 
      LineBox[{{215.5, 51.845}, {220.99950061295857`, 61.3697802551531}}]}}, 
    {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[
      StyleBox["\<\"C\"\>",
       StripOnInput->False,
       ScriptSizeMultipliers->0.8,
       ScriptMinSize->8,
       FontSize->Scaled[0.0770769541448646],
       FontWeight->Plain], {200., 25.}], 
     LineBox[{{182.6797802551531, 35.000499387041415`}, {173.155, 40.5}}], 
     {GrayLevel[0], InsetBox[
       StyleBox["\<\"H\"\>",
        StripOnInput->False,
        ScriptSizeMultipliers->0.8,
        ScriptMinSize->8,
        FontSize->Scaled[0.0770769541448646],
        FontWeight->Plain], {146.31, 56.00000000000001}], 
      LineBox[{{173.155, 40.5}, {163.6302197448469, 45.999500612958585`}}]}}, 
    {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[
      StyleBox["\<\"C\"\>",
       StripOnInput->False,
       ScriptSizeMultipliers->0.8,
       ScriptMinSize->8,
       FontSize->Scaled[0.0770769541448646],
       FontWeight->Plain], {200., 25.}], 
     LineBox[{{189.99950061295857`, 7.679780255153105}, {
      184.5, -1.8449999999999989`}}], 
     {GrayLevel[0], InsetBox[
       StyleBox["\<\"H\"\>",
        StripOnInput->False,
        ScriptSizeMultipliers->0.8,
        ScriptMinSize->8,
        FontSize->Scaled[0.0770769541448646],
        FontWeight->Plain], {169., -28.689999999999998`}], 
      LineBox[{{184.5, -1.8449999999999989`}, {
       179.00049938704143`, -11.369780255153103`}}]}}, 
    {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[
      StyleBox["\<\"C\"\>",
       StripOnInput->False,
       ScriptSizeMultipliers->0.8,
       ScriptMinSize->8,
       FontSize->Scaled[0.0770769541448646],
       FontWeight->Plain], {459.80999999999995`, -25.}], 
     LineBox[{{449.8095006129585, -42.320219744846895`}, {
      444.30999999999995`, -51.845}}], 
     {GrayLevel[0], InsetBox[
       StyleBox["\<\"H\"\>",
        StripOnInput->False,
        ScriptSizeMultipliers->0.8,
        ScriptMinSize->8,
        FontSize->Scaled[0.0770769541448646],
        FontWeight->Plain], {428.81, -78.69}], 
      LineBox[{{444.30999999999995`, -51.845}, {
       438.8104993870414, -61.36978025515309}}]}}, 
    {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[
      StyleBox["\<\"C\"\>",
       StripOnInput->False,
       ScriptSizeMultipliers->0.8,
       ScriptMinSize->8,
       FontSize->Scaled[0.0770769541448646],
       FontWeight->Plain], {459.80999999999995`, -25.}], 
     LineBox[{{477.13021974484684`, -35.00049938704141}, {486.655, -40.5}}], 
     {GrayLevel[0], InsetBox[
       StyleBox["\<\"H\"\>",
        StripOnInput->False,
        ScriptSizeMultipliers->0.8,
        ScriptMinSize->8,
        FontSize->Scaled[0.0770769541448646],
        FontWeight->Plain], {513.5, -56.00000000000001}], 
      LineBox[{{486.655, -40.5}, {496.1797802551531, -45.99950061295859}}]}}, 
    {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[
      StyleBox["\<\"C\"\>",
       StripOnInput->False,
       ScriptSizeMultipliers->0.8,
       ScriptMinSize->8,
       FontSize->Scaled[0.0770769541448646],
       FontWeight->Plain], {459.80999999999995`, -25.}], 
     LineBox[{{469.8104993870414, -7.679780255153105}, {475.30999999999995`, 
      1.8449999999999989`}}], 
     {GrayLevel[0], InsetBox[
       StyleBox["\<\"H\"\>",
        StripOnInput->False,
        ScriptSizeMultipliers->0.8,
        ScriptMinSize->8,
        FontSize->Scaled[0.0770769541448646],
        FontWeight->Plain], {490.81, 28.689999999999998`}], 
      LineBox[{{475.30999999999995`, 1.8449999999999989`}, {480.8095006129586,
        11.36978025515312}}]}}}, {}, {}},
  ImageSize->{{0., 142.0302016249526}, {0., 73.11756166951619}},
  PlotRange->{{123.81, 531.5}, {-108.69, 101.19}}]], "Output",
 CellChangeTimes->{3.568628146465623*^9, 3.568721905267194*^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[CellGroupData[{

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"}], "]"}], ",", 
         "\"\< \>\""}], "]"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "n"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"j", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{"Dij", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        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", ",", "\"\< \>\""}], "]"}],
         ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "1", ",", "n"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"j", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}], "*)"}], 
   "\[IndentingNewLine]", 
   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.568286391407484*^9, 3.568286398053096*^9}, {
   3.56856731902209*^9, 3.568567323936099*^9}, {3.5686259672509956`*^9, 
   3.568625969107399*^9}, {3.568626068916374*^9, 3.5686260698991756`*^9}}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"0.`", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", \
"\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", \
"\<\" \"\>", "\<\" \"\>"},
     {"2.8781886909715286`", 
      "0.`", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", \
"\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", \
"\<\" \"\>"},
     {"2.8779364416254674`", "0", 
      "0.`", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", \
"\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>"},
     {"0", "2.87805457989909`", "0", 
      "0.`", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", \
"\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>"},
     {"2.103269636689885`", "0", "0", "0", 
      "0.`", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", \
"\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>"},
     {"2.1033344695656426`", "0", "0", "0", "0", 
      "0.`", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", \
"\<\" \"\>", "\<\" \"\>", "\<\" \"\>"},
     {"0", "2.1032752261155525`", "0", "0", "0", "0", 
      "0.`", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", \
"\<\" \"\>", "\<\" \"\>"},
     {"0", "2.103172446127338`", "0", "0", "0", "0", "0", 
      "0.`", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", \
"\<\" \"\>"},
     {"0", "0", "2.103250029792444`", "0", "0", "0", "0", "0", 
      "0.`", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>"},
     {"0", "0", "2.1033977125517866`", "0", "0", "0", "0", "0", "0", 
      "0.`", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>"},
     {"0", "0", "2.1034616512826116`", "0", "0", "0", "0", "0", "0", "0", 
      "0.`", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>"},
     {"0", "0", "0", "2.1032441689296015`", "0", "0", "0", "0", "0", "0", "0",
       "0.`", "\<\" \"\>", "\<\" \"\>"},
     {"0", "0", "0", "2.1032801137733586`", "0", "0", "0", "0", "0", "0", "0",
       "0", "0.`", "\<\" \"\>"},
     {"0", "0", "0", "2.1031736197531177`", "0", "0", "0", "0", "0", "0", "0",
       "0", "0", "0.`"}
    },
    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.568628146481223*^9, 3.5687219052827935`*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"(*", 
  RowBox[{
   RowBox[{"Eij", "=", 
    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"}]}], ",", 
          RowBox[{"0.529177", "*", "sq"}], ",", "0"}], "]"}], ",", 
        "\"\< \>\""}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "n"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"j", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";", 
   "\[IndentingNewLine]", 
   RowBox[{"MatrixForm", "[", "Eij", "]"}], ";"}], "*)"}]], "Input",
 CellChangeTimes->{{3.5679574512554846`*^9, 3.5679574566530943`*^9}, {
   3.5679575983637433`*^9, 3.5679576252893906`*^9}, {3.5679577071427345`*^9, 
   3.5679577424611964`*^9}, 3.567958973147358*^9, {3.5679591782097178`*^9, 
   3.56795917902092*^9}, {3.567959225290601*^9, 3.567959254322252*^9}, {
   3.5682778318055534`*^9, 3.568277834176758*^9}}],

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[{"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}, {3.5685673432021327`*^9, 
   3.5685673517353477`*^9}, {3.568625843839179*^9, 3.568625849782789*^9}, {
   3.568626005627063*^9, 3.568626007389866*^9}, {3.568626080553994*^9, 
   3.5686260810375953`*^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[{"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}, {3.568567050623619*^9, 3.56856705698843*^9}, {
   3.568625871295227*^9, 3.568625877160837*^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[{"(*", 
   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.568285867917365*^9, 3.568285868634966*^9}, {
   3.5685662802006655`*^9, 3.5685663279991493`*^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[{
    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}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{"Total", " ", "kinetic", " ", "energy", " ", "in", " ", "Eh"}], 
    ",", " ", 
    RowBox[{"=", " ", 
     RowBox[{"|", "Etot", "|"}]}]}], " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{"Tsigma", "=", 
   RowBox[{
    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"}], "}"}]}], "]"}], "+", 
    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.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.568286113290196*^9, 3.5682861137113967`*^9}, {3.5682870401502237`*^9, 
   3.5682870902263117`*^9}, 3.568566812692001*^9, {3.56856697153148*^9, 
   3.568566972342681*^9}, {3.5686252771057835`*^9, 3.5686252869962006`*^9}, {
   3.5686253620323324`*^9, 3.568625362874734*^9}, {3.5686254979241714`*^9, 
   3.568625499718174*^9}, {3.5686260279195023`*^9, 3.568626034424713*^9}, {
   3.568626091162013*^9, 3.5686260974956245`*^9}, {3.5687218450510883`*^9, 
   3.5687218755647416`*^9}}],

Cell[BoxData["156.23310767412`"], "Output",
 CellChangeTimes->{3.5686281465124235`*^9, 3.5687219053139935`*^9}]
}, Open  ]],

Cell[BoxData[{
 RowBox[{
  RowBox[{"Qij", "=", "Lij"}], ";"}], "\[IndentingNewLine]", 
 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, 3.5685663358303633`*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 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, {3.568567024103572*^9, 3.5685670312327847`*^9}}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"\<\" \"\>", "\<\"C\"\>", "\<\"C\"\>", "\<\"C\"\>", "\<\"C\"\>", "\<\"H\
\"\>", "\<\"H\"\>", "\<\"H\"\>", "\<\"H\"\>", "\<\"H\"\>", "\<\"H\"\>", \
"\<\"H\"\>", "\<\"H\"\>", "\<\"H\"\>", "\<\"H\"\>"},
     {"\<\"C\"\>", 
      "0", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\
\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", \
"\<\" \"\>", "\<\" \"\>"},
     {"\<\"C\"\>", "1.1834943454857645`", 
      "0", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\
\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", \
"\<\" \"\>"},
     {"\<\"C\"\>", "1.1833682208127339`", "0", 
      "0", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\
\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>"},
     {"\<\"C\"\>", "0", "1.1834272899495453`", "0", 
      "0", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\
\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>"},
     {"\<\"H\"\>", "1.2836696366898848`", "0", "0", "0", 
      "0", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\
\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>"},
     {"\<\"H\"\>", "1.2837344695656425`", "0", "0", "0", "0", 
      "0", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\
\<\" \"\>", "\<\" \"\>", "\<\" \"\>"},
     {"\<\"H\"\>", "0", "1.2836752261155524`", "0", "0", "0", "0", 
      "0", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\
\<\" \"\>", "\<\" \"\>"},
     {"\<\"H\"\>", "0", "1.283572446127338`", "0", "0", "0", "0", "0", 
      "0", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\
\<\" \"\>"},
     {"\<\"H\"\>", "0", "0", "1.283650029792444`", "0", "0", "0", "0", "0", 
      "0", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>"},
     {"\<\"H\"\>", "0", "0", "1.2837977125517865`", "0", "0", "0", "0", "0", 
      "0", "0", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>"},
     {"\<\"H\"\>", "0", "0", "1.2838616512826115`", "0", "0", "0", "0", "0", 
      "0", "0", "0", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>"},
     {"\<\"H\"\>", "0", "0", "0", "1.2836441689296014`", "0", "0", "0", "0", 
      "0", "0", "0", "0", "\<\" \"\>", "\<\" \"\>"},
     {"\<\"H\"\>", "0", "0", "0", "1.2836801137733584`", "0", "0", "0", "0", 
      "0", "0", "0", "0", "0", "\<\" \"\>"},
     {"\<\"H\"\>", "0", "0", "0", "1.2835736197531176`", "0", "0", "0", "0", 
      "0", "0", "0", "0", "0", "0"}
    },
    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.5686281465280232`*^9, 3.5687219053295937`*^9}]
}, Open  ]],

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[CellGroupData[{

Cell[BoxData[
 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}],

Cell[CellGroupData[{

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"C\"\>", "\[InvisibleSpace]", "2", 
   "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", 
   "\[InvisibleSpace]", "1", "\[InvisibleSpace]", "\<\"  \"\>", 
   "\[InvisibleSpace]", "1.1834943454857645`"}],
  SequenceForm["C", 2, " ", "C", 1, "  ", 1.1834943454857645`],
  Editable->False]], "Print",
 CellChangeTimes->{3.568628146559223*^9, 3.568721905345194*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"C\"\>", "\[InvisibleSpace]", "3", 
   "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", 
   "\[InvisibleSpace]", "1", "\[InvisibleSpace]", "\<\"  \"\>", 
   "\[InvisibleSpace]", "1.1833682208127339`"}],
  SequenceForm["C", 3, " ", "C", 1, "  ", 1.1833682208127339`],
  Editable->False]], "Print",
 CellChangeTimes->{3.568628146559223*^9, 3.568721905345194*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"C\"\>", "\[InvisibleSpace]", "4", 
   "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", 
   "\[InvisibleSpace]", "2", "\[InvisibleSpace]", "\<\"  \"\>", 
   "\[InvisibleSpace]", "1.1834272899495453`"}],
  SequenceForm["C", 4, " ", "C", 2, "  ", 1.1834272899495453`],
  Editable->False]], "Print",
 CellChangeTimes->{3.568628146559223*^9, 3.568721905345194*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"H\"\>", "\[InvisibleSpace]", "5", 
   "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", 
   "\[InvisibleSpace]", "1", "\[InvisibleSpace]", "\<\"  \"\>", 
   "\[InvisibleSpace]", "1.2836696366898848`"}],
  SequenceForm["H", 5, " ", "C", 1, "  ", 1.2836696366898848`],
  Editable->False]], "Print",
 CellChangeTimes->{3.568628146559223*^9, 3.568721905345194*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"H\"\>", "\[InvisibleSpace]", "6", 
   "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", 
   "\[InvisibleSpace]", "1", "\[InvisibleSpace]", "\<\"  \"\>", 
   "\[InvisibleSpace]", "1.2837344695656425`"}],
  SequenceForm["H", 6, " ", "C", 1, "  ", 1.2837344695656425`],
  Editable->False]], "Print",
 CellChangeTimes->{3.568628146559223*^9, 3.568721905345194*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"H\"\>", "\[InvisibleSpace]", "7", 
   "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", 
   "\[InvisibleSpace]", "2", "\[InvisibleSpace]", "\<\"  \"\>", 
   "\[InvisibleSpace]", "1.2836752261155524`"}],
  SequenceForm["H", 7, " ", "C", 2, "  ", 1.2836752261155524`],
  Editable->False]], "Print",
 CellChangeTimes->{3.568628146559223*^9, 3.568721905360794*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"H\"\>", "\[InvisibleSpace]", "8", 
   "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", 
   "\[InvisibleSpace]", "2", "\[InvisibleSpace]", "\<\"  \"\>", 
   "\[InvisibleSpace]", "1.283572446127338`"}],
  SequenceForm["H", 8, " ", "C", 2, "  ", 1.283572446127338],
  Editable->False]], "Print",
 CellChangeTimes->{3.568628146559223*^9, 3.568721905360794*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"H\"\>", "\[InvisibleSpace]", "9", 
   "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", 
   "\[InvisibleSpace]", "3", "\[InvisibleSpace]", "\<\"  \"\>", 
   "\[InvisibleSpace]", "1.283650029792444`"}],
  SequenceForm["H", 9, " ", "C", 3, "  ", 1.283650029792444],
  Editable->False]], "Print",
 CellChangeTimes->{3.568628146559223*^9, 3.568721905360794*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"H\"\>", "\[InvisibleSpace]", "10", 
   "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", 
   "\[InvisibleSpace]", "3", "\[InvisibleSpace]", "\<\"  \"\>", 
   "\[InvisibleSpace]", "1.2837977125517865`"}],
  SequenceForm["H", 10, " ", "C", 3, "  ", 1.2837977125517865`],
  Editable->False]], "Print",
 CellChangeTimes->{3.568628146559223*^9, 3.568721905360794*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"H\"\>", "\[InvisibleSpace]", "11", 
   "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", 
   "\[InvisibleSpace]", "3", "\[InvisibleSpace]", "\<\"  \"\>", 
   "\[InvisibleSpace]", "1.2838616512826115`"}],
  SequenceForm["H", 11, " ", "C", 3, "  ", 1.2838616512826115`],
  Editable->False]], "Print",
 CellChangeTimes->{3.568628146559223*^9, 3.568721905360794*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"H\"\>", "\[InvisibleSpace]", "12", 
   "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", 
   "\[InvisibleSpace]", "4", "\[InvisibleSpace]", "\<\"  \"\>", 
   "\[InvisibleSpace]", "1.2836441689296014`"}],
  SequenceForm["H", 12, " ", "C", 4, "  ", 1.2836441689296014`],
  Editable->False]], "Print",
 CellChangeTimes->{3.568628146559223*^9, 3.568721905360794*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"H\"\>", "\[InvisibleSpace]", "13", 
   "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", 
   "\[InvisibleSpace]", "4", "\[InvisibleSpace]", "\<\"  \"\>", 
   "\[InvisibleSpace]", "1.2836801137733584`"}],
  SequenceForm["H", 13, " ", "C", 4, "  ", 1.2836801137733584`],
  Editable->False]], "Print",
 CellChangeTimes->{3.568628146559223*^9, 3.568721905360794*^9}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"H\"\>", "\[InvisibleSpace]", "14", 
   "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", 
   "\[InvisibleSpace]", "4", "\[InvisibleSpace]", "\<\"  \"\>", 
   "\[InvisibleSpace]", "1.2835736197531176`"}],
  SequenceForm["H", 14, " ", "C", 4, "  ", 1.2835736197531176`],
  Editable->False]], "Print",
 CellChangeTimes->{3.568628146559223*^9, 3.5687219053763943`*^9}]
}, Open  ]]
}, Open  ]],

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", "]"}]}]}], "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.568389047208571*^9, 
   3.5683890542753835`*^9}, {3.5683902124526176`*^9, 3.5683902131078186`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"2", ",", "1", ",", "1.1834943454857645`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"3", ",", "1", ",", "1.1833682208127339`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"4", ",", "2", ",", "1.1834272899495453`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"5", ",", "1", ",", "1.2836696366898848`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"6", ",", "1", ",", "1.2837344695656425`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"7", ",", "2", ",", "1.2836752261155524`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"8", ",", "2", ",", "1.283572446127338`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"9", ",", "3", ",", "1.283650029792444`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"10", ",", "3", ",", "1.2837977125517865`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"11", ",", "3", ",", "1.2838616512826115`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"12", ",", "4", ",", "1.2836441689296014`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"13", ",", "4", ",", "1.2836801137733584`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"14", ",", "4", ",", "1.2835736197531176`"}], "}"}]}], 
  "}"}]], "Output",
 CellChangeTimes->{3.5686281466060233`*^9, 3.568721905391994*^9}],

Cell[BoxData["13"], "Output",
 CellChangeTimes->{3.5686281466060233`*^9, 3.568721905391994*^9}]
}, Open  ]],

Cell[BoxData[" "], "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.5682852978767633`*^9, 3.568285382460112*^9}, {
   3.568285675366227*^9, 3.568285707861084*^9}, {3.568285985822372*^9, 
   3.568285991485182*^9}, {3.568286358600627*^9, 3.568286369505046*^9}, {
   3.5682869898869357`*^9, 3.568287001134555*^9}, 3.56838891559114*^9}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"wx", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     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"}], "]"}], "]"}], "]"}], "]"}]}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "bc"}], "}"}]}], "]"}]}], ";", "     ", 
  RowBox[{"wy", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     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"}], "]"}], "]"}], "]"}], "]"}]}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "bc"}], "}"}]}], "]"}]}], ";", "      ", 
  RowBox[{"wz", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     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"}], "]"}], "]"}], "]"}], "]"}]}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "bc"}], "}"}]}], "]"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"Plot7", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"Graphics3D", "[", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"Darker", "[", "Gray", "]"}], ",", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"Opacity", "[", "0.3", "]"}], ",", 
          RowBox[{"Sphere", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{"wx", "[", 
               RowBox[{"[", "i", "]"}], "]"}], ",", 
              RowBox[{"wy", "[", 
               RowBox[{"[", "i", "]"}], "]"}], ",", 
              RowBox[{"wz", "[", 
               RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", 
            RowBox[{"a", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "3"}], "]"}], "]"}]}], "]"}]}], "}"}]}], 
       "}"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "bc"}], "}"}]}], "]"}]}], 
  ";"}], "\[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.5682852978767633`*^9, 3.568285382460112*^9}, {
   3.568285675366227*^9, 3.568285707861084*^9}, {3.568285985822372*^9, 
   3.568285991485182*^9}, {3.568286358600627*^9, 3.568286369505046*^9}, {
   3.5682869898869357`*^9, 3.568287001134555*^9}, {3.56838891559114*^9, 
   3.568388937930379*^9}, {3.568391616907484*^9, 3.568391624754298*^9}, {
   3.568391667342373*^9, 3.5683916742531853`*^9}, {3.568391794825797*^9, 
   3.5683917959489985`*^9}, {3.5683918276482544`*^9, 3.56839185905111*^9}, {
   3.5683919288612328`*^9, 3.5683919590316854`*^9}, {3.568524054858098*^9, 
   3.5685240555756993`*^9}, {3.5685666362556906`*^9, 3.568566666394944*^9}, {
   3.5685667055978127`*^9, 3.5685667164398317`*^9}, {3.5686174374353323`*^9, 
   3.568617488572222*^9}, {3.568618796213319*^9, 3.568618805526535*^9}, {
   3.568619728205756*^9, 3.5686197473469896`*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"wx", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     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"}], "]"}], "]"}], "]"}], "]"}]}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", 
       RowBox[{"bc", "+", "1"}], ",", "nb"}], "}"}]}], "]"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"wy", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     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"}], "]"}], "]"}], "]"}], "]"}]}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", 
       RowBox[{"bc", "+", "1"}], ",", "nb"}], "}"}]}], "]"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"wz", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     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"}], "]"}], "]"}], "]"}], "]"}]}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", 
       RowBox[{"bc", "+", "1"}], ",", "nb"}], "}"}]}], "]"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"Plot8", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"Graphics3D", "[", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"Opacity", "[", "0.4", "]"}], ",", 
        RowBox[{"Sphere", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"wx", "[", 
             RowBox[{"[", "i", "]"}], "]"}], ",", 
            RowBox[{"wy", "[", 
             RowBox[{"[", "i", "]"}], "]"}], ",", 
            RowBox[{"wz", "[", 
             RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", 
          RowBox[{"a", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "3"}], "]"}], "]"}]}], "]"}]}], "}"}], "]"}], 
     ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "bh"}], "}"}]}], "]"}]}], 
  ";"}], "\[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.5682852978767633`*^9, 3.568285382460112*^9}, {
   3.568285675366227*^9, 3.568285707861084*^9}, {3.568285985822372*^9, 
   3.568285991485182*^9}, {3.568286358600627*^9, 3.568286369505046*^9}, {
   3.5682869898869357`*^9, 3.568287001134555*^9}, {3.56838891559114*^9, 
   3.5683889530312057`*^9}, {3.5683893397246847`*^9, 3.568389482371335*^9}, {
   3.568389512432588*^9, 3.568389840313964*^9}, {3.568390052801937*^9, 
   3.568390055422742*^9}, 3.568390114578046*^9, 3.568390152064912*^9, {
   3.568390244027073*^9, 3.568390244635474*^9}, 3.5683904737686768`*^9, {
   3.56852316213093*^9, 3.5685231796497607`*^9}, {3.568524015343229*^9, 
   3.5685240160296297`*^9}, 3.568524093889367*^9, 3.568524191046337*^9, {
   3.5685666461617084`*^9, 3.568566671137352*^9}, {3.568566721712641*^9, 
   3.568566754238698*^9}, {3.5686174565765657`*^9, 3.5686174617401752`*^9}, {
   3.5686174951554337`*^9, 3.568617503891449*^9}, {3.5686185965173683`*^9, 
   3.5686185984517717`*^9}, {3.5686188095825424`*^9, 
   3.5686188240749683`*^9}, {3.5686197519957976`*^9, 3.5686197778606434`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"PLot9", "=", " ", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"Graphics3D", "[", 
        RowBox[{"{", 
         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.1"}], "]"}], "}"}], "]"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "nb"}], "}"}]}], "]"}]}], ";"}], "*)"}], 
  "\[IndentingNewLine]", 
  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"}], "}"}]}], "]"}]}], ";"}]}]], "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.5682852978767633`*^9, 3.568285382460112*^9}, {
   3.568285675366227*^9, 3.568285707861084*^9}, {3.568285985822372*^9, 
   3.568285991485182*^9}, {3.568286358600627*^9, 3.568286369505046*^9}, {
   3.5682869898869357`*^9, 3.568287001134555*^9}, {3.56838891559114*^9, 
   3.5683889530312057`*^9}, {3.5683893397246847`*^9, 3.568389482371335*^9}, {
   3.568389512432588*^9, 3.568389840313964*^9}, {3.568390052801937*^9, 
   3.568390055422742*^9}, 3.568390114578046*^9, 3.568390152064912*^9, {
   3.568390244027073*^9, 3.568390244635474*^9}, 3.5683904737686768`*^9, {
   3.568390648613784*^9, 3.568390694009864*^9}, {3.5683907253191185`*^9, 
   3.5683908630517607`*^9}, {3.5683910922941628`*^9, 
   3.5683911307638307`*^9}, {3.5683913309433823`*^9, 3.568391333236586*^9}, {
   3.56839141488713*^9, 3.568391447397587*^9}, {3.568391486709656*^9, 
   3.568391491342864*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Show", "[", 
  RowBox[{"Plot1", ",", "Plot7", ",", "Plot8", ",", "Plot9", ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"AspectRatio", "\[Rule]", "Automatic"}], ",", 
     RowBox[{"Boxed", "\[Rule]", "False"}]}], "}"}]}], "]"}]], "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.5682852978767633`*^9, 3.568285382460112*^9}, {
   3.568285675366227*^9, 3.568285707861084*^9}, {3.568285985822372*^9, 
   3.568285991485182*^9}, {3.568286358600627*^9, 3.568286369505046*^9}, {
   3.5682869898869357`*^9, 3.568287001134555*^9}, {3.56838891559114*^9, 
   3.568388924030755*^9}, {3.568390166448137*^9, 3.568390171252945*^9}, {
   3.5683908961862183`*^9, 3.5683909036430316`*^9}, 3.5683909701615486`*^9, 
   3.5683913684302483`*^9, {3.5683915370821447`*^9, 3.5683915698422017`*^9}, {
   3.5685232430326724`*^9, 3.5685232511134863`*^9}, {3.568524281323696*^9, 
   3.5685242837417*^9}, 3.5685243967170987`*^9}],

Cell[BoxData[
 Graphics3DBox[{
   {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], 
    SphereBox[{0.7380230406699999, -0.530425301214, -1.115758481084}, 
     0.2556]}, 
   {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], 
    SphereBox[{-0.7379380029999999, 0.530476323816, 1.115832180398}, 0.2556]}, 
   {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], 
    SphereBox[{3.31346446566, 0.746001463842, -1.2584271249060002`}, 0.2556]}, 
   {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], 
    SphereBox[{-3.31346446566, -0.7460694939779999, 1.25828350573}, 0.2556]}, 
   {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], 
    SphereBox[{0.9965980287019999, -2.60232277734, -0.862533307358}, 0.1]}, 
   {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], 
    SphereBox[{-0.323882028866, -0.1669149622376, -2.8945878005}, 0.1]}, 
   {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], 
    SphereBox[{-0.996514880758, 2.6023794691199997`, 0.862608896398}, 0.1]}, 
   {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], 
    SphereBox[{0.32380455010000003`, 0.1669756224422, 2.89456890324}, 0.1]}, 
   {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], 
    SphereBox[{4.39198778564, -0.0292540923156, -2.8892020813999997`}, 0.1]}, 
   {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], 
    SphereBox[{3.0548176680399997`, 
     2.8180350002399996`, -1.5115294664419998`}, 0.1]}, 
   {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], 
    SphereBox[{4.37522591602, 0.382482432126, 0.520636520534}, 0.1]}, 
   {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], 
    SphereBox[{-4.39198778564, 0.029186062179599997`, 2.88905090332}, 0.1]}, 
   {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], 
    SphereBox[{-3.05489325708, -2.81797830846, 1.5115067897299999`}, 0.1]}, 
   {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], 
    SphereBox[{-4.37520701876, -0.382571249248, -0.520455106838}, 0.1]}, 
   {RGBColor[
    0.33333333333333337`, 0.33333333333333337`, 0.33333333333333337`], 
    {Opacity[0.3], 
     SphereBox[{0.000042518835000171507`, 0.000025511300999814246`, 
      0.00003684965699979337}, 1.1834943454857645`]}}, 
   {RGBColor[
    0.33333333333333337`, 0.33333333333333337`, 0.33333333333333337`], 
    {Opacity[0.3], 
     SphereBox[{2.025743753165, 0.10778808131399986`, -1.187092802995}, 
      1.1833682208127339`]}}, 
   {RGBColor[
    0.33333333333333337`, 0.33333333333333337`, 0.33333333333333337`], 
    {Opacity[0.3], 
     SphereBox[{-2.0257012343299996`, -0.10779658508099987`, 1.187057843064}, 
      1.1834272899495453`]}}, 
   {Opacity[0.4], 
    SphereBox[{0.9272601317885236, -2.0467353452054677`, -0.9304366305886971},
      1.1834943454857645`]}, 
   {Opacity[0.4], 
    SphereBox[{-0.039136797955585334`, -0.26438868035007307`, \
-2.417602446688057}, 1.1833682208127339`]}, 
   {Opacity[0.4], 
    SphereBox[{-0.9271766613729637, 2.0467919932433376`, 0.9305115324404376}, 
     1.1834272899495453`]}, 
   {Opacity[0.4], 
    SphereBox[{0.03908096452286969, 0.264454265219365, 2.417571634884694}, 
     1.2836696366898848`]}, 
   {Opacity[0.4], 
    SphereBox[{4.102774832423616, 0.17863566281543153`, -2.451899305964699}, 
     1.2837344695656425`]}, 
   {Opacity[0.4], 
    SphereBox[{3.1241705978198397`, 
     2.2624449149438033`, -1.4436632137716365`}, 1.2836752261155524`]}, 
   {Opacity[0.4], 
    SphereBox[{4.0905364102056705`, 0.47995258745914215`, 
     0.043617177007423624`}, 1.283572446127338]}, 
   {Opacity[0.4], 
    SphereBox[{-4.102774026507937, -0.17870427225333063`, 
     2.4517489362788423`}, 1.283650029792444]}, 
   {Opacity[0.4], 
    SphereBox[{-3.124229795130316, -2.2623906034674315`, 1.443604311481499}, 
     1.2837977125517865`]}, 
   {Opacity[0.4], 
    SphereBox[{-4.090483592066063, -0.4800491788407958, -0.04345759789905013},
      1.2838616512826115`]}, 
   {RGBColor[0.19999999999999996`, 0., 0.], 
    CylinderBox[{{-0.7379380029999999, 0.530476323816, 1.115832180398}, {
     0.7380230406699999, -0.530425301214, -1.115758481084}}, 0.05]}, 
   {RGBColor[0.19999999999999996`, 0., 0.], 
    CylinderBox[{{3.31346446566, 0.746001463842, -1.2584271249060002`}, {
     0.7380230406699999, -0.530425301214, -1.115758481084}}, 0.05]}, 
   {RGBColor[0.19999999999999996`, 0., 0.], 
    CylinderBox[{{-3.31346446566, -0.7460694939779999, 
     1.25828350573}, {-0.7379380029999999, 0.530476323816, 1.115832180398}}, 
     0.05]}, 
   {RGBColor[0.19999999999999996`, 0., 0.], 
    CylinderBox[{{0.9965980287019999, -2.60232277734, -0.862533307358}, {
     0.7380230406699999, -0.530425301214, -1.115758481084}}, 0.05]}, 
   {RGBColor[0.19999999999999996`, 0., 0.], 
    CylinderBox[{{-0.323882028866, -0.1669149622376, -2.8945878005}, {
     0.7380230406699999, -0.530425301214, -1.115758481084}}, 0.05]}, 
   {RGBColor[0.19999999999999996`, 0., 0.], 
    CylinderBox[{{-0.996514880758, 2.6023794691199997`, 
     0.862608896398}, {-0.7379380029999999, 0.530476323816, 1.115832180398}}, 
     0.05]}, 
   {RGBColor[0.19999999999999996`, 0., 0.], 
    CylinderBox[{{0.32380455010000003`, 0.1669756224422, 
     2.89456890324}, {-0.7379380029999999, 0.530476323816, 1.115832180398}}, 
     0.05]}, 
   {RGBColor[0.19999999999999996`, 0., 0.], 
    CylinderBox[{{4.39198778564, -0.0292540923156, -2.8892020813999997`}, {
     3.31346446566, 0.746001463842, -1.2584271249060002`}}, 0.05]}, 
   {RGBColor[0.19999999999999996`, 0., 0.], 
    CylinderBox[{{3.0548176680399997`, 
     2.8180350002399996`, -1.5115294664419998`}, {3.31346446566, 
     0.746001463842, -1.2584271249060002`}}, 0.05]}, 
   {RGBColor[0.19999999999999996`, 0., 0.], 
    CylinderBox[{{4.37522591602, 0.382482432126, 0.520636520534}, {
     3.31346446566, 0.746001463842, -1.2584271249060002`}}, 0.05]}, 
   {RGBColor[0.19999999999999996`, 0., 0.], 
    CylinderBox[{{-4.39198778564, 0.029186062179599997`, 
     2.88905090332}, {-3.31346446566, -0.7460694939779999, 1.25828350573}}, 
     0.05]}, 
   {RGBColor[0.19999999999999996`, 0., 0.], 
    CylinderBox[{{-3.05489325708, -2.81797830846, 
     1.5115067897299999`}, {-3.31346446566, -0.7460694939779999, 
     1.25828350573}}, 0.05]}, 
   {RGBColor[0.19999999999999996`, 0., 0.], 
    CylinderBox[{{-4.37520701876, -0.382571249248, -0.520455106838}, \
{-3.31346446566, -0.7460694939779999, 1.25828350573}}, 0.05]}},
  AspectRatio->Automatic,
  Boxed->False]], "Output",
 CellChangeTimes->{3.5686281466528234`*^9, 3.568721905423194*^9}]
}, Open  ]]
},
WindowSize->{1245, 910},
WindowMargins->{{Automatic, 226}, {Automatic, 16}},
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) (November 20, 2012)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[579, 22, 885, 20, 70, "Input"],
Cell[1467, 44, 2272, 42, 161, "Output"]
}, Open  ]],
Cell[3754, 89, 854, 16, 29, "Input"],
Cell[CellGroupData[{
Cell[4633, 109, 779, 12, 29, "Input"],
Cell[5415, 123, 500, 17, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[5952, 145, 763, 25, 70, "Input"],
Cell[6718, 172, 502, 17, 28, "Output"],
Cell[7223, 191, 94, 1, 28, "Output"],
Cell[7320, 194, 95, 1, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[7452, 200, 867, 13, 29, "Input"],
Cell[8322, 215, 381, 6, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[8740, 226, 213, 5, 50, "Input"],
Cell[8956, 233, 97, 1, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[9090, 239, 3126, 85, 110, "Input"],
Cell[12219, 326, 500, 11, 19, "Message"],
Cell[12722, 339, 324, 6, 28, "Output"],
Cell[13049, 347, 347, 7, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[13433, 359, 247, 5, 29, "Input"],
Cell[13683, 366, 113, 1, 28, "Output"]
}, Open  ]],
Cell[13811, 370, 763, 14, 29, "Input"],
Cell[14577, 386, 4522, 93, 70, "Input"],
Cell[CellGroupData[{
Cell[19124, 483, 2568, 36, 29, "Input"],
Cell[21695, 521, 3300, 54, 306, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[25032, 580, 718, 11, 29, "Input"],
Cell[25753, 593, 9117, 221, 79, "Output"]
}, Open  ]],
Cell[34885, 817, 1254, 35, 70, "Input"],
Cell[CellGroupData[{
Cell[36164, 856, 4981, 113, 90, "Input"],
Cell[41148, 971, 2791, 53, 232, "Output"]
}, Open  ]],
Cell[43954, 1027, 1985, 50, 50, "Input"],
Cell[45942, 1079, 1890, 45, 90, "Input"],
Cell[47835, 1126, 511, 11, 50, "Input"],
Cell[48349, 1139, 2123, 48, 90, "Input"],
Cell[50475, 1189, 1935, 47, 50, "Input"],
Cell[52413, 1238, 1261, 33, 50, "Input"],
Cell[CellGroupData[{
Cell[53699, 1275, 2015, 47, 50, "Input"],
Cell[55717, 1324, 111, 1, 28, "Output"]
}, Open  ]],
Cell[55843, 1328, 1995, 35, 90, "Input"],
Cell[CellGroupData[{
Cell[57863, 1367, 1642, 22, 29, "Input"],
Cell[59508, 1391, 3163, 57, 247, "Output"]
}, Open  ]],
Cell[62686, 1451, 1655, 24, 29, "Input"],
Cell[CellGroupData[{
Cell[64366, 1479, 2749, 55, 70, "Input"],
Cell[CellGroupData[{
Cell[67140, 1538, 427, 8, 21, "Print"],
Cell[67570, 1548, 427, 8, 21, "Print"],
Cell[68000, 1558, 427, 8, 21, "Print"],
Cell[68430, 1568, 427, 8, 21, "Print"],
Cell[68860, 1578, 427, 8, 21, "Print"],
Cell[69290, 1588, 427, 8, 21, "Print"],
Cell[69720, 1598, 424, 8, 21, "Print"],
Cell[70147, 1608, 424, 8, 21, "Print"],
Cell[70574, 1618, 429, 8, 21, "Print"],
Cell[71006, 1628, 429, 8, 21, "Print"],
Cell[71438, 1638, 429, 8, 21, "Print"],
Cell[71870, 1648, 429, 8, 21, "Print"],
Cell[72302, 1658, 431, 8, 21, "Print"]
}, Open  ]]
}, Open  ]],
Cell[72760, 1670, 2696, 53, 90, "Input"],
Cell[CellGroupData[{
Cell[75481, 1727, 867, 15, 50, "Input"],
Cell[76351, 1744, 1256, 30, 48, "Output"],
Cell[77610, 1776, 95, 1, 28, "Output"]
}, Open  ]],
Cell[77720, 1780, 4130, 55, 29, "Input"],
Cell[81853, 1837, 10795, 254, 110, "Input"],
Cell[92651, 2093, 10888, 256, 110, "Input"],
Cell[103542, 2351, 8460, 166, 70, "Input"],
Cell[CellGroupData[{
Cell[112027, 2521, 4724, 66, 29, "Input"],
Cell[116754, 2589, 6505, 124, 348, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)