(* 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[    142570,       3167]
NotebookOptionsPosition[    139080,       3059]
NotebookOutlinePosition[    139616,       3080]
CellTagsIndexPosition[    139573,       3077]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell[BoxData[
 RowBox[{"ChemicalData", "[", 
  RowBox[{"\"\<Aspirin\>\"", ",", "\"\<MoleculePlot\>\""}], "]"}]], "Input"],

Cell[BoxData[
 Graphics3DBox[
  {AbsoluteThickness[3], EdgeForm[None], Specularity[
    GrayLevel[1], 100], GraphicsComplex3DBox[CompressedData["
1:eJwtUm1Ik1EUnl9hlmGiYhGmBjYZw0Ltg1BPaNmI0Fmsptva3rm9tqWbS01E
6EdIpcXQhD7E0lTsx5DQCNmPdWhDJWwKkhYhJFEgWRKWsR9R997Tgcu9nM/n
ee7JklxVtmiFQpHNTgw7+dwmavBt62rOwZsa9L76pTTMHUADf/issLfjq3Yy
IGHKWPlK289GMPk126OdpTiZwQJKO1bxh0WGhw+YNcr4vb2k0NSpA+HfacPm
u+mxunA2jPevj0517QfdDvV4f14uJiwNuYscevg00hfKdOdj1jNPXXeiDsSc
UBmI/jPV0FjEMhU1uKZlHtVl+Dyf/SESOAmqG2nvw9P/89trQeC+aAM7ty96
yt9wkf+aBQYXC9aPxCUTT6WW8OaaiO+0DG9mmQ1fQr9za2RgdB8KXYJ6uuea
UOBI0uGmkinTYUShS3ERnlppSxlzmUHMOyaTv8cDgce785aHC4HDS9lio/63
ZeC0E0YbQPQttVP9RjMMM/dSUAY2XZ5wRKHIn21FpqZfs1ZKOP/KlL+tFRmL
eF/seXw6dfrwreOHUMHNcAWDLM07LiOHkbEs0f/8MAMvK49yUjynBUKZTHC9
Buu6v7WXjOQj/1aVrCa+L60YWlDZdnmKcZHj/X0GxRgwgOC/aqb9eFdPfCok
EPM/2kn3xCa435B8HV+rgcP2PpKQ0xvqlaGA9++qw6sVGWmO1DLqFzFjEheq
N510WDDS/hgl2pNzbpCSLhx9PgPw5MSLjT8FZhD9Yu2kY4yD7j4L6b/pof8d
rIWcs3FDPf4E0ifcgvG+PdW1MZWkj8WKxXxAfQuq+YJ26rDyzr2G5C4tivrU
JoobLVQ/b8Z//FB4nw==
    "], {
     {RGBColor[0.65, 0.7, 0.7], SphereBox[14, 24.], SphereBox[15, 24.], 
      SphereBox[16, 24.], SphereBox[17, 24.], SphereBox[18, 24.], 
      SphereBox[19, 24.], SphereBox[20, 24.], SphereBox[21, 24.]}, 
     {RGBColor[0.4, 0.4, 0.4], SphereBox[5, 34.], SphereBox[6, 34.], 
      SphereBox[7, 34.], SphereBox[8, 34.], SphereBox[9, 34.], 
      SphereBox[10, 34.], SphereBox[11, 34.], SphereBox[12, 34.], 
      SphereBox[13, 34.]}, 
     {RGBColor[0.800498, 0.201504, 0.192061], SphereBox[1, 31.], 
      SphereBox[2, 31.], SphereBox[3, 31.], SphereBox[4, 31.]}, 
     {RGBColor[0.65, 0.7, 0.7], CylinderBox[{25, 21}, 15.], 
      CylinderBox[{28, 14}, 15.], CylinderBox[{29, 15}, 15.], 
      CylinderBox[{30, 16}, 15.], CylinderBox[{31, 17}, 15.], 
      CylinderBox[{32, 18}, 15.], CylinderBox[{33, 19}, 15.], 
      CylinderBox[{34, 20}, 15.]}, 
     {RGBColor[0.4, 0.4, 0.4], CylinderBox[{22, 5}, 15.], 
      CylinderBox[{23, 12}, 15.], CylinderBox[{24, 11}, 15.], 
      CylinderBox[{26, 11}, 15.], CylinderBox[{27, 12}, 15.], 
      CylinderBox[{5, 6}, 15.], CylinderBox[{5, 7}, 15.], 
      CylinderBox[{6, 8}, 15.], CylinderBox[{6, 11}, 15.], 
      CylinderBox[{7, 9}, 15.], CylinderBox[{7, 28}, 15.], 
      CylinderBox[{8, 10}, 15.], CylinderBox[{8, 29}, 15.], 
      CylinderBox[{9, 10}, 15.], CylinderBox[{9, 30}, 15.], 
      CylinderBox[{10, 31}, 15.], CylinderBox[{12, 13}, 15.], 
      CylinderBox[{13, 32}, 15.], CylinderBox[{13, 33}, 15.], 
      CylinderBox[{13, 34}, 15.]}, 
     {RGBColor[0.800498, 0.201504, 0.192061], CylinderBox[{1, 22}, 15.], 
      CylinderBox[{1, 23}, 15.], CylinderBox[{2, 24}, 15.], 
      CylinderBox[{2, 25}, 15.], CylinderBox[{3, 26}, 15.], 
      CylinderBox[{4, 27}, 15.]}}]},
  Boxed->False,
  ImageSize->282.44217559121745`,
  Lighting->"Neutral",
  ViewPoint->{1.3, -2.4, -2.}]], "Output",
 CellChangeTimes->{
  3.5682249012266088`*^9, 3.5682250183360147`*^9, 3.5682250935437465`*^9, 
   3.568225796293781*^9, 3.5682259554500933`*^9, 3.568266119351782*^9, {
   3.5682661961351166`*^9, 3.568266218552356*^9}, 3.568266325599744*^9, 
   3.5682664258611197`*^9, 3.5682665066068616`*^9, 3.568267092637491*^9, 
   3.568267862545643*^9, {3.568267931310564*^9, 3.5682679541802044`*^9}, 
   3.5682680040378914`*^9, 3.568268075439217*^9, 3.5682681858250113`*^9, 
   3.568268325726057*^9, 3.5682688789810286`*^9, 3.5682690287880917`*^9, 
   3.5682695424033937`*^9, 3.5682696911496553`*^9, 3.568271398229453*^9, 
   3.568272429219664*^9, 3.5682726342040243`*^9, 3.5682730124422884`*^9, 
   3.5682734354994316`*^9, 3.5682735174307756`*^9, 3.5682736370985856`*^9, 
   3.568273995368815*^9, 3.5682764646347523`*^9, 3.568276518657647*^9, 
   3.568276667949909*^9, {3.568277700484523*^9, 3.568277718237354*^9}}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"t", "=", 
   RowBox[{"0.01889726", "*", 
    RowBox[{"ChemicalData", "[", 
     RowBox[{"\"\<Aspirin\>\"", ",", "\"\<AtomPositions\>\""}], "]"}]}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.567862912170579*^9, 3.567862949329844*^9}, {
   3.5678630095303497`*^9, 3.5678630107471523`*^9}, {3.5679573013548217`*^9, 
   3.5679573132576427`*^9}, 3.568114367533684*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"ChemicalData", "\[NoBreak]", "[", 
  RowBox[{"\"\<Aspirin\>\"", ",", "\"\<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}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"\<\"C\"\>", ",", "\<\"C\"\>"}], "}"}], ",", "2"}], "}"}], 
     ",", "3"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"\<\"C\"\>", ",", "\<\"C\"\>"}], "}"}], ",", "1"}], "}"}], 
     ",", "5"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"\<\"C\"\>", ",", "\<\"H\"\>"}], "}"}], ",", "1"}], "}"}], 
     ",", "7"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"\<\"C\"\>", ",", "\<\"O\"\>"}], "}"}], ",", "2"}], "}"}], 
     ",", "2"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"\<\"C\"\>", ",", "\<\"O\"\>"}], "}"}], ",", "1"}], "}"}], 
     ",", "3"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"\<\"H\"\>", ",", "\<\"O\"\>"}], "}"}], ",", "1"}], "}"}], 
     ",", "1"}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{
  3.5682249012578087`*^9, 3.568225018367214*^9, 3.5682250935905466`*^9, 
   3.5682257963249807`*^9, 3.5682259554656935`*^9, 3.5682661193985815`*^9, {
   3.568266196181916*^9, 3.5682662185679555`*^9}, 3.5682663256621437`*^9, 
   3.5682664258923197`*^9, 3.5682665066536617`*^9, 3.568267092684291*^9, 
   3.5682678625768433`*^9, {3.568267931341764*^9, 3.568267954211404*^9}, 
   3.5682680040534916`*^9, 3.568268075486017*^9, 3.568268185856211*^9, 
   3.5682683257572565`*^9, 3.5682688790122285`*^9, 3.5682690288348913`*^9, 
   3.568269542418994*^9, 3.568269691165255*^9, 3.5682713982606535`*^9, 
   3.5682724292820644`*^9, 3.5682726342352247`*^9, 3.5682730124578886`*^9, 
   3.568273435515032*^9, 3.5682735174619756`*^9, 3.5682736371453857`*^9, 
   3.568273995415615*^9, 3.568276464681552*^9, 3.5682765186888466`*^9, 
   3.568276667996709*^9, {3.5682777005313225`*^9, 3.568277718299754*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"u", "=", 
  RowBox[{"ChemicalData", "\[NoBreak]", "[", 
   RowBox[{"\"\<Aspirin\>\"", ",", "\"\<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}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"\<\"O\"\>", ",", "\<\"O\"\>", ",", "\<\"O\"\>", ",", "\<\"O\"\>", 
   ",", "\<\"C\"\>", ",", "\<\"C\"\>", ",", "\<\"C\"\>", ",", "\<\"C\"\>", 
   ",", "\<\"C\"\>", ",", "\<\"C\"\>", ",", "\<\"C\"\>", ",", "\<\"C\"\>", 
   ",", "\<\"C\"\>", ",", "\<\"H\"\>", ",", "\<\"H\"\>", ",", "\<\"H\"\>", 
   ",", "\<\"H\"\>", ",", "\<\"H\"\>", ",", "\<\"H\"\>", ",", "\<\"H\"\>", 
   ",", "\<\"H\"\>"}], "}"}]], "Output",
 CellChangeTimes->{{3.5682777005469227`*^9, 3.5682777183153543`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

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

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"\<\"O\"\>", ",", "\<\"O\"\>", ",", "\<\"O\"\>", ",", "\<\"O\"\>", 
   ",", "\<\"C\"\>", ",", "\<\"C\"\>", ",", "\<\"C\"\>", ",", "\<\"C\"\>", 
   ",", "\<\"C\"\>", ",", "\<\"C\"\>", ",", "\<\"C\"\>", ",", "\<\"C\"\>", 
   ",", "\<\"C\"\>", ",", "\<\"H\"\>", ",", "\<\"H\"\>", ",", "\<\"H\"\>", 
   ",", "\<\"H\"\>", ",", "\<\"H\"\>", ",", "\<\"H\"\>", ",", "\<\"H\"\>", 
   ",", "\<\"H\"\>"}], "}"}]], "Output",
 CellChangeTimes->{{3.568277700562523*^9, 3.568277718330954*^9}}],

Cell[BoxData["21"], "Output",
 CellChangeTimes->{{3.568277700562523*^9, 3.568277718330954*^9}}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"u1", "=", 
   RowBox[{"{", 
    RowBox[{
    "0.1887", ",", "0.1887", ",", "0.1887", ",", "0.1887", ",", "0.2555", ",",
      "0.2555", ",", "0.2555", ",", "0.2555", ",", "0.2555", ",", "0.2555", 
     ",", "0.2555", ",", "0.2555", ",", "0.2555", ",", "0.1", ",", "0.1", ",",
      "0.1", ",", "0.1", ",", "0.1", ",", "0.1", ",", "0.1", ",", "0.1"}], 
    "}"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.5679480369450207`*^9, 3.567948069112277*^9}, {
   3.5679486805553513`*^9, 3.5679487948723516`*^9}, 3.567950206581231*^9, {
   3.5681216077716007`*^9, 3.5681216839309344`*^9}, {3.568266380090639*^9, 
   3.5682664108850937`*^9}, {3.568266481350417*^9, 3.56826649971165*^9}}],

Cell[BoxData[
 RowBox[{"(*", 
  RowBox[{
   RowBox[{"For", "[", 
    RowBox[{
     RowBox[{"i", "=", "1"}], ",", 
     RowBox[{"i", "<", "15"}], ",", 
     RowBox[{"i", "++"}], ",", 
     RowBox[{"If", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"u", "[", 
         RowBox[{"[", "i", "]"}], "]"}], "=", "\"\<C\>\""}], ",", 
       RowBox[{
        RowBox[{"u1", "[", 
         RowBox[{"[", "i", "]"}], "]"}], "=", "0.26"}], ",", 
       RowBox[{
        RowBox[{"u1", "[", 
         RowBox[{"[", "i", "]"}], "]"}], "=", "0.18"}]}], "]"}]}], "]"}], 
   "\[IndentingNewLine]", "u1"}], "*)"}]], "Input",
 CellChangeTimes->{{3.5679472243083935`*^9, 3.567947394894693*^9}, {
  3.567947453425996*^9, 3.567947471241227*^9}, {3.5679475415037503`*^9, 
  3.5679475448733563`*^9}, {3.5679488138107853`*^9, 3.567948820877598*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"ChemicalData", "\[NoBreak]", "[", 
  RowBox[{"\"\<Aspirin\>\"", ",", "\"\<MolecularFormulaString\>\""}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.567869377256353*^9, 3.5678694101568108`*^9}, {
   3.5678695244894114`*^9, 3.56786953511303*^9}, 3.56786988243204*^9, 
   3.5678699252853155`*^9, {3.567870069179968*^9, 3.5678700766835814`*^9}, {
   3.5678703216664114`*^9, 3.567870325738019*^9}, 3.5679271251072817`*^9, {
   3.567942840466694*^9, 3.567942854381918*^9}, {3.567942932272855*^9, 
   3.567942950509287*^9}, 3.5679430681802936`*^9}],

Cell[BoxData["\<\"(CH3CO2)C6H4CO2H\"\>"], "Output",
 CellChangeTimes->{
  3.568224901289009*^9, 3.5682250183984146`*^9, 3.5682250936217465`*^9, 
   3.568225796356181*^9, 3.568225955512493*^9, 3.568266119429782*^9, {
   3.5682661962131166`*^9, 3.568266218599156*^9}, 3.568266325693344*^9, 
   3.56826642592352*^9, 3.568266506700462*^9, 3.5682670927154913`*^9, 
   3.568267862608043*^9, {3.568267931372964*^9, 3.5682679542426043`*^9}, 
   3.568268004084692*^9, 3.568268075517217*^9, 3.568268185887411*^9, 
   3.568268325788457*^9, 3.5682688790434284`*^9, 3.5682690288660917`*^9, 
   3.568269542450194*^9, 3.568269691196455*^9, 3.5682713982918534`*^9, 
   3.5682724293444643`*^9, 3.568272634266424*^9, 3.5682730124890885`*^9, 
   3.5682734355462313`*^9, 3.5682735174931755`*^9, 3.568273637176586*^9, 
   3.5682739954468155`*^9, 3.568276464697152*^9, 3.568276518720047*^9, 
   3.568276668012309*^9, {3.568277700578123*^9, 3.568277718346554*^9}}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"v", "=", 
   RowBox[{"0.01889726", "*", 
    RowBox[{"ChemicalData", "\[NoBreak]", "[", 
     RowBox[{"\"\<Aspirin\>\"", ",", "\"\<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}],

Cell[BoxData[
 RowBox[{"(*", 
  RowBox[{
   RowBox[{"lp", "=", 
    RowBox[{"ListPlot", "[", 
     RowBox[{"v", ",", " ", 
      RowBox[{"Axes", "\[Rule]", "False"}]}], "]"}]}], ";"}], "*)"}]], "Input",
 CellChangeTimes->{{3.567946177515355*^9, 3.5679462297286463`*^9}, 
   3.5679463060439806`*^9, {3.5679467119254932`*^9, 3.5679467232199135`*^9}, {
   3.5680116987964935`*^9, 3.568011700100568*^9}, 3.5681220290347404`*^9, {
   3.5682660959829407`*^9, 3.5682661026441517`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"Plot1", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"Graphics", "[", 
        RowBox[{"{", 
         RowBox[{"Circle", "[", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{
             RowBox[{
              RowBox[{"v", "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", "1"}], "]"}], "]"}], "/", "0.5292"}], ",", 
             RowBox[{
              RowBox[{"v", "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", "2"}], "]"}], "]"}], "/", "0.5292"}]}], 
            "}"}], ",", 
           RowBox[{"If", "[", 
            RowBox[{
             RowBox[{"i", "<", "5"}], ",", "0.18", ",", "0.26"}], "]"}]}], 
          "]"}], "}"}], "]"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "13"}], "}"}]}], "]"}]}], ";"}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"Plot1", "=", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"Graphics3D", "[", 
       RowBox[{"{", 
        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", ",", "21"}], "}"}]}], "]"}]}], ";"}]}]], "Input",
 CellChangeTimes->{{3.5679440895764875`*^9, 3.5679441000441055`*^9}, {
   3.5679441319773617`*^9, 3.5679441523041973`*^9}, {3.567944240210352*^9, 
   3.5679442813476243`*^9}, {3.567944323966899*^9, 3.567944345931738*^9}, {
   3.567944397567828*^9, 3.567944456941533*^9}, {3.567944491557994*^9, 
   3.5679445468444905`*^9}, {3.5679446505690727`*^9, 3.567944683875131*^9}, {
   3.567944736244423*^9, 3.5679448177233667`*^9}, {3.567945549707852*^9, 
   3.5679455499730525`*^9}, {3.5679456039959474`*^9, 3.567945606195551*^9}, {
   3.567945704460124*^9, 3.5679457164877453`*^9}, {3.5679457857362666`*^9, 
   3.5679458050023003`*^9}, {3.567946353670864*^9, 3.5679464127169676`*^9}, {
   3.5679466013368993`*^9, 3.5679466311173515`*^9}, {3.5679467590999765`*^9, 
   3.567946915739851*^9}, {3.5679469466123056`*^9, 3.567946946815106*^9}, {
   3.5679471385082426`*^9, 3.5679471497870626`*^9}, {3.5679483907224417`*^9, 
   3.5679484036080647`*^9}, {3.567948441937332*^9, 3.5679484421869326`*^9}, {
   3.567948496740228*^9, 3.56794849772303*^9}, {3.567948847912445*^9, 
   3.5679488949933276`*^9}, {3.568011728877214*^9, 3.5680117308213253`*^9}, {
   3.568011766325356*^9, 3.568011773436763*^9}, {3.5681173314452896`*^9, 
   3.5681173370612993`*^9}, {3.568117919800323*^9, 3.5681179678952074`*^9}, {
   3.56811809459863*^9, 3.568118105409449*^9}, {3.5681181818963833`*^9, 
   3.5681182389144835`*^9}, {3.568118274888147*^9, 3.5681182883977704`*^9}, {
   3.568118518513775*^9, 3.56811852739019*^9}, {3.5681188250543127`*^9, 
   3.568118862634779*^9}, {3.568118972396572*^9, 3.5681189862493963`*^9}, {
   3.568119020694257*^9, 3.568119047136303*^9}, {3.5681190903015785`*^9, 
   3.568119136930061*^9}, {3.5681191671005135`*^9, 3.568119202808976*^9}, {
   3.5681194788046613`*^9, 3.568119487603077*^9}, {3.5681195353235607`*^9, 
   3.5681195410331707`*^9}, 3.56811959201406*^9, {3.5681196426985493`*^9, 
   3.568119659702579*^9}, {3.568119726813897*^9, 3.56811974027672*^9}, {
   3.5682659184078283`*^9, 3.5682659267694435`*^9}, {3.568265972477524*^9, 
   3.568266002725977*^9}}],

Cell[BoxData[
 RowBox[{
  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", ",", "21"}], "}"}]}], "]"}]}], ";"}], 
  "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"Show", "[", 
     RowBox[{"Plot2", ",", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"AspectRatio", "\[Rule]", "Automatic"}], ",", 
        RowBox[{"Axes", " ", "\[Rule]", " ", "False"}], ",", 
        RowBox[{"Frame", " ", "\[Rule]", " ", "Automatic"}], ",", 
        RowBox[{"PlotRange", " ", "\[Rule]", " ", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{"0", ",", "14"}], "}"}], ",", 
           RowBox[{"{", 
            RowBox[{
             RowBox[{"-", "6"}], ",", "6"}], "}"}]}], "}"}]}]}], "}"}]}], 
     "]"}], ";"}], "*)"}], "\[IndentingNewLine]"}]], "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}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Show", "[", 
  RowBox[{"Plot1", ",", 
   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}],

Cell[BoxData[
 Graphics3DBox[{
   SphereBox[{-2.093816408, -0.96496968464, -0.31766294059999994`}, 0.1887], 
   SphereBox[{2.9740507787999997`, -2.8338331096, 5.428804852799999}, 0.1887],
    SphereBox[{-0.7338939893599999, -3.3497283075999995`, 
    3.5349214555999997`}, 0.1887], 
   SphereBox[{-3.544181113, 1.70381475612, -3.1853221456}, 0.1887], 
   SphereBox[{0.17769093578, 0.22247744198, -0.23797319517999999`}, 0.2555], 
   SphereBox[{1.9210954515999998`, -0.36732493988, 1.66696509912}, 0.2555], 
   SphereBox[{0.7848976940999999, 2.0244634637999996`, -2.0410930526}, 
    0.2555], SphereBox[{4.2586865136, 0.8596741519199999, 1.72195612572}, 
    0.2555], SphereBox[{3.0970719413999994`, 
    3.2310535147999997`, -1.9577561359999998`}, 0.2555], 
   SphereBox[{4.8303286286, 2.6452384548, -0.08513404602599998}, 0.2555], 
   SphereBox[{1.29688115928, -2.279954419, 3.5948257697999995`}, 0.2555], 
   SphereBox[{-3.8767728889999997`, -0.185710932924, -1.947362643}, 0.2555], 
   SphereBox[{-6.291842717, -1.6678721676000001`, -2.2315774334}, 0.2555], 
   SphereBox[{-0.5526503687, 2.4848007174, -3.5122447436}, 0.1], 
   SphereBox[{5.6156987542, 0.40901229543999995`, 3.1783301593999997`}, 0.1], 
   SphereBox[{3.560243784, 4.633986097199999, -3.3656020059999996`}, 0.1], 
   SphereBox[{6.6372846298, 3.5927470712, -0.041141224745999995`}, 0.1], 
   SphereBox[{-7.498432768, -0.7529802219599999, -3.6278959747999995`}, 0.1], 
   SphereBox[{-7.2703428398, -1.73992741998, -0.42057741855999997`}, 0.1], 
   SphereBox[{-5.853426284999999, -3.5825425508, -2.8517855065999997`}, 0.1], 
   SphereBox[{2.5617125655999997`, -4.082564050399999, 6.6973779166}, 0.1]},
  AspectRatio->Automatic,
  ImageSize->{651., Automatic},
  ViewPoint->{0.5384478532598622, -3.27921562820733, 0.6378234654202138},
  ViewVertical->{0.039479586127532953`, -0.6882474948061027, 
   1.229747061122278}]], "Output",
 CellChangeTimes->{
  3.5682249013202085`*^9, 3.5682250184296145`*^9, 3.5682250936529465`*^9, 
   3.568225796387381*^9, 3.5682259555436935`*^9, 3.568266119663782*^9, {
   3.568266196244316*^9, 3.568266218630356*^9}, 3.5682663259429445`*^9, 
   3.5682664259703197`*^9, 3.568266506731662*^9, 3.568267092746691*^9, 
   3.5682678626236434`*^9, {3.5682679314041643`*^9, 3.5682679542582045`*^9}, 
   3.568268004115892*^9, 3.568268075532817*^9, 3.5682681859030113`*^9, 
   3.568268325804057*^9, 3.568268879074629*^9, 3.568269028881692*^9, 
   3.568269542481394*^9, 3.5682696912276554`*^9, 3.5682713983074536`*^9, 
   3.568272429406864*^9, 3.5682726342820244`*^9, 3.5682730125202885`*^9, 
   3.5682734355618315`*^9, 3.5682735175087757`*^9, 3.5682736371921864`*^9, 
   3.568273995478015*^9, 3.5682764647283525`*^9, 3.568276518735647*^9, 
   3.5682766680435095`*^9, {3.5682777006093225`*^9, 3.568277718377754*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"ChemicalData", "\[NoBreak]", "[", 
  RowBox[{"\"\<Aspirin\>\"", ",", "\"\<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}}],

Cell[BoxData[
 GraphicsBox[{{
    {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[
      StyleBox["\<\"O\"\>",
       StripOnInput->False,
       ScriptSizeMultipliers->0.8,
       ScriptMinSize->8,
       FontSize->Scaled[0.05129653480986163],
       FontWeight->Plain], {373.21, -6.}], 
     LineBox[{{391.66318322568156`, -16.654262832379658`}, {
      416.51, -31.000000000000004`}}], 
     {GrayLevel[0], InsetBox[
       StyleBox["\<\"C\"\>",
        StripOnInput->False,
        ScriptSizeMultipliers->0.8,
        ScriptMinSize->8,
        FontSize->Scaled[0.05129653480986163],
        FontWeight->Plain], {459.80999999999995`, -56.00000000000001}], 
      LineBox[{{416.51, -31.000000000000004`}, {
       442.7564020679032, -46.15381181749611}}]}}, 
    {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[
      StyleBox["\<\"O\"\>",
       StripOnInput->False,
       ScriptSizeMultipliers->0.8,
       ScriptMinSize->8,
       FontSize->Scaled[0.05129653480986163],
       FontWeight->Plain], {373.21, -6.}], 
     LineBox[{{354.7562841069346, -16.65334019920644}, {
      329.905, -31.000000000000004`}}], 
     {GrayLevel[0], InsetBox[
       StyleBox["\<\"C\"\>",
        StripOnInput->False,
        ScriptSizeMultipliers->0.8,
        ScriptMinSize->8,
        FontSize->Scaled[0.05129653480986163],
        FontWeight->Plain], {286.6, -56.00000000000001}], 
      LineBox[{{329.905, -31.000000000000004`}, {
       303.6540901992186, -46.15466447337572}}]}}, 
    {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[
      StyleBox["\<\"O\"\>",
       StripOnInput->False,
       ScriptSizeMultipliers->0.8,
       ScriptMinSize->8,
       FontSize->Scaled[0.05129653480986163],
       FontWeight->Plain], {633.01, 144.}], 
     LineBox[{{614.5568167743185, 133.34573716762034`}, {589.71, 119.}}], 
     {GrayLevel[0], InsetBox[
       StyleBox["\<\"C\"\>",
        StripOnInput->False,
        ScriptSizeMultipliers->0.8,
        ScriptMinSize->8,
        FontSize->Scaled[0.05129653480986163],
        FontWeight->Plain], {546.41, 94.}], 
      LineBox[{{589.71, 119.}, {563.4635979320967, 103.84618818250388`}}]}}, 
    {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[
      StyleBox["\<\"O\"\>",
       StripOnInput->False,
       ScriptSizeMultipliers->0.8,
       ScriptMinSize->8,
       FontSize->Scaled[0.05129653480986163],
       FontWeight->Plain], {633.01, 144.}], 
     LineBox[{{633.01, 165.30805687203792`}, {633.01, 175.}}], 
     {GrayLevel[0], InsetBox[
       StyleBox["\<\"H\"\>",
        StripOnInput->False,
        ScriptSizeMultipliers->0.8,
        ScriptMinSize->8,
        FontSize->Scaled[0.05129653480986163],
        FontWeight->Plain], {633.01, 206.}], 
      LineBox[{{633.01, 175.}, {633.01, 186.30805687203792`}}]}}, 
    {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[
      StyleBox["\<\"O\"\>",
       StripOnInput->False,
       ScriptSizeMultipliers->0.8,
       ScriptMinSize->8,
       FontSize->Scaled[0.05129653480986163],
       FontWeight->Plain], {459.80999999999995`, 144.}], 
     LineBox[{{{480.76318322568153`, 137.67573716762035`}, {
      505.60999999999996`, 123.33}}, {{475.76318322568153`, 
      129.01573716762033`}, {500.60999999999996`, 114.67}}}], 
     {GrayLevel[0], InsetBox[
       StyleBox["\<\"C\"\>",
        StripOnInput->False,
        ScriptSizeMultipliers->0.8,
        ScriptMinSize->8,
        FontSize->Scaled[0.05129653480986163],
        FontWeight->Plain], {546.41, 94.}], 
      LineBox[{{{505.60999999999996`, 123.33}, {531.8564020679032, 
       108.17618818250388`}}, {{500.60999999999996`, 114.67}, {
       526.8564020679032, 99.51618818250388}}}]}}, 
    {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[
      StyleBox["\<\"O\"\>",
       StripOnInput->False,
       ScriptSizeMultipliers->0.8,
       ScriptMinSize->8,
       FontSize->Scaled[0.05129653480986163],
       FontWeight->Plain], {286.6, -156.}], 
     LineBox[{{{281.6, -134.69194312796208`}, {281.6, -106.}}, {{
      291.6, -134.69194312796208`}, {291.6, -106.}}}], 
     {GrayLevel[0], InsetBox[
       StyleBox["\<\"C\"\>",
        StripOnInput->False,
        ScriptSizeMultipliers->0.8,
        ScriptMinSize->8,
        FontSize->Scaled[0.05129653480986163],
        FontWeight->Plain], {286.6, -56.00000000000001}], 
      LineBox[{{{281.6, -106.}, {281.6, -75.69194312796209}}, {{
       291.6, -106.}, {291.6, -75.69194312796209}}}]}}, 
    {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[
      StyleBox["\<\"C\"\>",
       StripOnInput->False,
       ScriptSizeMultipliers->0.8,
       ScriptMinSize->8,
       FontSize->Scaled[0.05129653480986163],
       FontWeight->Plain], {459.80999999999995`, -56.00000000000001}], 
     LineBox[{{476.8635979320967, -46.15381181749612}, {
      503.10999999999996`, -31.000000000000004`}}], 
     {GrayLevel[0], InsetBox[
       StyleBox["\<\"C\"\>",
        StripOnInput->False,
        ScriptSizeMultipliers->0.8,
        ScriptMinSize->8,
        FontSize->Scaled[0.05129653480986163],
        FontWeight->Plain], {546.41, -6.}], 
      LineBox[{{503.10999999999996`, -31.000000000000004`}, {
       529.3564020679032, -15.846188182503885`}}]}}, 
    {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[
      StyleBox["\<\"C\"\>",
       StripOnInput->False,
       ScriptSizeMultipliers->0.8,
       ScriptMinSize->8,
       FontSize->Scaled[0.05129653480986163],
       FontWeight->Plain], {459.80999999999995`, -56.00000000000001}], 
     LineBox[{{{464.80999999999995`, -75.69194312796209}, {
      464.80999999999995`, -106.}}, {{
      454.80999999999995`, -75.69194312796209}, {
      454.80999999999995`, -106.}}}], 
     {GrayLevel[0], InsetBox[
       StyleBox["\<\"C\"\>",
        StripOnInput->False,
        ScriptSizeMultipliers->0.8,
        ScriptMinSize->8,
        FontSize->Scaled[0.05129653480986163],
        FontWeight->Plain], {459.80999999999995`, -156.}], 
      LineBox[{{{464.80999999999995`, -106.}, {
       464.80999999999995`, -136.30805687203792`}}, {{
       454.80999999999995`, -106.}, {
       454.80999999999995`, -136.30805687203792`}}}]}}, 
    {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[
      StyleBox["\<\"C\"\>",
       StripOnInput->False,
       ScriptSizeMultipliers->0.8,
       ScriptMinSize->8,
       FontSize->Scaled[0.05129653480986163],
       FontWeight->Plain], {546.41, -6.}], 
     LineBox[{{{565.9635979320967, -11.516188182503875`}, {
      592.21, -26.67}}, {{560.9635979320967, -20.17618818250388}, {
      587.21, -35.330000000000005`}}}], 
     {GrayLevel[0], InsetBox[
       StyleBox["\<\"C\"\>",
        StripOnInput->False,
        ScriptSizeMultipliers->0.8,
        ScriptMinSize->8,
        FontSize->Scaled[0.05129653480986163],
        FontWeight->Plain], {633.01, -56.00000000000001}], 
      LineBox[{{{592.21, -26.67}, {618.4564020679032, -41.82381181749611}}, {{
       587.21, -35.330000000000005`}, {
       613.4564020679032, -50.48381181749611}}}]}}, 
    {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[
      StyleBox["\<\"C\"\>",
       StripOnInput->False,
       ScriptSizeMultipliers->0.8,
       ScriptMinSize->8,
       FontSize->Scaled[0.05129653480986163],
       FontWeight->Plain], {546.41, -6.}], 
     LineBox[{{546.41, 13.691943127962084`}, {546.41, 44.}}], 
     {GrayLevel[0], InsetBox[
       StyleBox["\<\"C\"\>",
        StripOnInput->False,
        ScriptSizeMultipliers->0.8,
        ScriptMinSize->8,
        FontSize->Scaled[0.05129653480986163],
        FontWeight->Plain], {546.41, 94.}], 
      LineBox[{{546.41, 44.}, {546.41, 74.30805687203792}}]}}, 
    {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[
      StyleBox["\<\"C\"\>",
       StripOnInput->False,
       ScriptSizeMultipliers->0.8,
       ScriptMinSize->8,
       FontSize->Scaled[0.05129653480986163],
       FontWeight->Plain], {459.80999999999995`, -156.}], 
     LineBox[{{476.8635979320967, -165.84618818250388`}, {
      503.10999999999996`, -181.}}], 
     {GrayLevel[0], InsetBox[
       StyleBox["\<\"C\"\>",
        StripOnInput->False,
        ScriptSizeMultipliers->0.8,
        ScriptMinSize->8,
        FontSize->Scaled[0.05129653480986163],
        FontWeight->Plain], {546.41, -206.}], 
      LineBox[{{503.10999999999996`, -181.}, {
       529.3564020679032, -196.15381181749612`}}]}}, 
    {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[
      StyleBox["\<\"C\"\>",
       StripOnInput->False,
       ScriptSizeMultipliers->0.8,
       ScriptMinSize->8,
       FontSize->Scaled[0.05129653480986163],
       FontWeight->Plain], {459.80999999999995`, -156.}], 
     LineBox[{{442.75576690847316`, -165.84508800441955`}, {
      432.9599999999999, -171.5}}], 
     {GrayLevel[0], InsetBox[
       StyleBox["\<\"H\"\>",
        StripOnInput->False,
        ScriptSizeMultipliers->0.8,
        ScriptMinSize->8,
        FontSize->Scaled[0.05129653480986163],
        FontWeight->Plain], {406.10999999999996`, -187.}], 
      LineBox[{{432.9599999999999, -171.5}, {
       423.16423309152674`, -177.15491199558042`}}]}}, 
    {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[
      StyleBox["\<\"C\"\>",
       StripOnInput->False,
       ScriptSizeMultipliers->0.8,
       ScriptMinSize->8,
       FontSize->Scaled[0.05129653480986163],
       FontWeight->Plain], {633.01, -56.00000000000001}], 
     LineBox[{{633.01, -75.69194312796209}, {633.01, -106.}}], 
     {GrayLevel[0], InsetBox[
       StyleBox["\<\"C\"\>",
        StripOnInput->False,
        ScriptSizeMultipliers->0.8,
        ScriptMinSize->8,
        FontSize->Scaled[0.05129653480986163],
        FontWeight->Plain], {633.01, -156.}], 
      LineBox[{{633.01, -106.}, {633.01, -136.30805687203792`}}]}}, 
    {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[
      StyleBox["\<\"C\"\>",
       StripOnInput->False,
       ScriptSizeMultipliers->0.8,
       ScriptMinSize->8,
       FontSize->Scaled[0.05129653480986163],
       FontWeight->Plain], {633.01, -56.00000000000001}], 
     LineBox[{{650.0642330915267, -46.15491199558042}, {
      659.8599999999999, -40.5}}], 
     {GrayLevel[0], InsetBox[
       StyleBox["\<\"H\"\>",
        StripOnInput->False,
        ScriptSizeMultipliers->0.8,
        ScriptMinSize->8,
        FontSize->Scaled[0.05129653480986163],
        FontWeight->Plain], {686.7099999999999, -25.}], 
      LineBox[{{659.8599999999999, -40.5}, {
       669.6557669084731, -34.84508800441955}}]}}, 
    {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[
      StyleBox["\<\"C\"\>",
       StripOnInput->False,
       ScriptSizeMultipliers->0.8,
       ScriptMinSize->8,
       FontSize->Scaled[0.05129653480986163],
       FontWeight->Plain], {546.41, -206.}], 
     LineBox[{{{560.9635979320967, -191.8238118174961}, {587.21, -176.67}}, {{
      565.9635979320967, -200.48381181749613`}, {592.21, -185.33}}}], 
     {GrayLevel[0], InsetBox[
       StyleBox["\<\"C\"\>",
        StripOnInput->False,
        ScriptSizeMultipliers->0.8,
        ScriptMinSize->8,
        FontSize->Scaled[0.05129653480986163],
        FontWeight->Plain], {633.01, -156.}], 
      LineBox[{{{587.21, -176.67}, {
       613.4564020679032, -161.51618818250387`}}, {{592.21, -185.33}, {
       618.4564020679032, -170.1761881825039}}}]}}, 
    {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[
      StyleBox["\<\"C\"\>",
       StripOnInput->False,
       ScriptSizeMultipliers->0.8,
       ScriptMinSize->8,
       FontSize->Scaled[0.05129653480986163],
       FontWeight->Plain], {546.41, -206.}], 
     LineBox[{{546.41, -225.69194312796208`}, {546.41, -237.}}], 
     {GrayLevel[0], InsetBox[
       StyleBox["\<\"H\"\>",
        StripOnInput->False,
        ScriptSizeMultipliers->0.8,
        ScriptMinSize->8,
        FontSize->Scaled[0.05129653480986163],
        FontWeight->Plain], {546.41, -268.}], 
      LineBox[{{546.41, -237.}, {546.41, -248.30805687203792`}}]}}, 
    {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[
      StyleBox["\<\"C\"\>",
       StripOnInput->False,
       ScriptSizeMultipliers->0.8,
       ScriptMinSize->8,
       FontSize->Scaled[0.05129653480986163],
       FontWeight->Plain], {633.01, -156.}], 
     LineBox[{{650.0642330915267, -165.84508800441958`}, {
      659.8599999999999, -171.5}}], 
     {GrayLevel[0], InsetBox[
       StyleBox["\<\"H\"\>",
        StripOnInput->False,
        ScriptSizeMultipliers->0.8,
        ScriptMinSize->8,
        FontSize->Scaled[0.05129653480986163],
        FontWeight->Plain], {686.7099999999999, -187.}], 
      LineBox[{{659.8599999999999, -171.5}, {
       669.6557669084731, -177.15491199558045`}}]}}, 
    {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[
      StyleBox["\<\"C\"\>",
       StripOnInput->False,
       ScriptSizeMultipliers->0.8,
       ScriptMinSize->8,
       FontSize->Scaled[0.05129653480986163],
       FontWeight->Plain], {286.6, -56.00000000000001}], 
     LineBox[{{269.5464020679033, -46.15381181749612}, {
      243.3, -31.000000000000004`}}], 
     {GrayLevel[0], InsetBox[
       StyleBox["\<\"C\"\>",
        StripOnInput->False,
        ScriptSizeMultipliers->0.8,
        ScriptMinSize->8,
        FontSize->Scaled[0.05129653480986163],
        FontWeight->Plain], {200., -6.}], 
      LineBox[{{243.3, -31.000000000000004`}, {
       217.05359793209672`, -15.846188182503885`}}]}}, 
    {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[
      StyleBox["\<\"C\"\>",
       StripOnInput->False,
       ScriptSizeMultipliers->0.8,
       ScriptMinSize->8,
       FontSize->Scaled[0.05129653480986163],
       FontWeight->Plain], {200., -6.}], 
     LineBox[{{209.84646325904197`, 11.05343910896655}, {215.5, 20.845}}], 
     {GrayLevel[0], InsetBox[
       StyleBox["\<\"H\"\>",
        StripOnInput->False,
        ScriptSizeMultipliers->0.8,
        ScriptMinSize->8,
        FontSize->Scaled[0.05129653480986163],
        FontWeight->Plain], {231., 47.69}], 
      LineBox[{{215.5, 20.845}, {221.15353674095803`, 
       30.636560891033447`}}]}}, 
    {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[
      StyleBox["\<\"C\"\>",
       StripOnInput->False,
       ScriptSizeMultipliers->0.8,
       ScriptMinSize->8,
       FontSize->Scaled[0.05129653480986163],
       FontWeight->Plain], {200., -6.}], 
     LineBox[{{182.94656089103344`, 3.846463259041963}, {173.155, 9.5}}], 
     {GrayLevel[0], InsetBox[
       StyleBox["\<\"H\"\>",
        StripOnInput->False,
        ScriptSizeMultipliers->0.8,
        ScriptMinSize->8,
        FontSize->Scaled[0.05129653480986163],
        FontWeight->Plain], {146.31, 25.}], 
      LineBox[{{173.155, 9.5}, {163.36343910896656`, 15.153536740958037`}}]}}, 
    {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[
      StyleBox["\<\"C\"\>",
       StripOnInput->False,
       ScriptSizeMultipliers->0.8,
       ScriptMinSize->8,
       FontSize->Scaled[0.05129653480986163],
       FontWeight->Plain], {200., -6.}], 
     LineBox[{{190.15353674095803`, -23.05343910896655}, {184.5, -32.845}}], 
     {GrayLevel[0], InsetBox[
       StyleBox["\<\"H\"\>",
        StripOnInput->False,
        ScriptSizeMultipliers->0.8,
        ScriptMinSize->8,
        FontSize->Scaled[0.05129653480986163],
        FontWeight->Plain], {169., -59.69}], 
      LineBox[{{184.5, -32.845}, {
       178.84646325904197`, -42.63656089103345}}]}}}, {}, {}},
  ImageSize->{{0., 194.94494193548383`}, {0., 178.81290322580645`}},
  PlotRange->{{130.098, 702.9219999999999}, {-297.38388625592415`, 
   228.03791469194312`}}]], "Output",
 CellChangeTimes->{
  3.568224901382609*^9, 3.568225018492015*^9, 3.568225093715347*^9, 
   3.5682257964497814`*^9, 3.5682259555904937`*^9, 3.568266119772982*^9, {
   3.568266196306717*^9, 3.5682662186927557`*^9}, 3.5682663260209446`*^9, 
   3.56826642603272*^9, 3.568266506778462*^9, 3.5682670927934914`*^9, 
   3.5682678626860437`*^9, {3.5682679314509645`*^9, 3.5682679543206043`*^9}, 
   3.568268004162692*^9, 3.5682680755952177`*^9, 3.568268185965411*^9, 
   3.568268325866457*^9, 3.5682688791214285`*^9, 3.568269028944092*^9, 
   3.5682695425437937`*^9, 3.568269691274455*^9, 3.5682713983698535`*^9, 
   3.5682724294692645`*^9, 3.5682726343444242`*^9, 3.568273012567089*^9, 
   3.568273435624232*^9, 3.5682735175711756`*^9, 3.568273637254586*^9, 
   3.568273995524815*^9, 3.5682764647907524`*^9, 3.568276518798047*^9, 
   3.568276668090309*^9, {3.568277700656123*^9, 3.568277718424554*^9}}]
}, Open  ]],

Cell[BoxData[{
 RowBox[{
  RowBox[{"x", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"t", "[", 
      RowBox[{"[", 
       RowBox[{"i", ",", "1"}], "]"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "21"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"y", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"t", "[", 
      RowBox[{"[", 
       RowBox[{"i", ",", "2"}], "]"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "21"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"z", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"t", "[", 
      RowBox[{"[", 
       RowBox[{"i", ",", "3"}], "]"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "21"}], "}"}]}], "]"}]}], ";"}]}], "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}}],

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", ",", "21"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"j", ",", "1", ",", "21"}], "}"}]}], "]"}]}], ";"}], 
   "\[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", ",", "21"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"j", ",", "1", ",", "21"}], "}"}]}], "]"}]}], ";"}], "*)"}], 
   "\[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}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"0.`", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", \
"\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", \
"\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", \
"\<\" \"\>", "\<\" \"\>", "\<\" \"\>"},
     {"0", 
      "0.`", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", \
"\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", \
"\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", \
"\<\" \"\>", "\<\" \"\>"},
     {"0", "0", 
      "0.`", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", \
"\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", \
"\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", \
"\<\" \"\>"},
     {"0", "0", "0", 
      "0.`", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", \
"\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", \
"\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>"},
     {"2.5643959809065504`", "0", "0", "0", 
      "0.`", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", \
"\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", \
"\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>"},
     {"0", "0", "0", "0", "2.6487952092451184`", 
      "0.`", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", \
"\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", \
"\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>"},
     {"0", "0", "0", "0", "2.620514241595812`", "0", 
      "0.`", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", \
"\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", \
"\<\" \"\>", "\<\" \"\>", "\<\" \"\>"},
     {"0", "0", "0", "0", "0", "2.6406216611674673`", "0", 
      "0.`", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", \
"\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", \
"\<\" \"\>", "\<\" \"\>"},
     {"0", "0", "0", "0", "0", "0", "2.6093973140780804`", "0", 
      "0.`", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", \
"\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", \
"\<\" \"\>"},
     {"0", "0", "0", "0", "0", "0", "0", "2.6039565042799775`", 
      "2.6180281740472`", 
      "0.`", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", \
"\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>"},
     {"0", "2.5462047862305086`", "2.296095033730954`", "0", "0", 
      "2.78647479311946`", "0", "0", "0", "0", 
      "0.`", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", \
"\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>"},
     {"2.5381290619802077`", "0", "0", "2.283302071364557`", "0", "0", "0", 
      "0", "0", "0", "0", 
      "0.`", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", \
"\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>"},
     {"0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", 
      "2.847831147869261`", 
      "0.`", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", \
"\<\" \"\>", "\<\" \"\>", "\<\" \"\>"},
     {"0", "0", "0", "0", "0", "0", "2.0408901257255825`", "0", "0", "0", "0",
       "0", "0", 
      "0.`", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", \
"\<\" \"\>", "\<\" \"\>"},
     {"0", "0", "0", "0", "0", "0", "0", "2.04098105233772`", "0", "0", "0", 
      "0", "0", "0", 
      "0.`", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", \
"\<\" \"\>"},
     {"0", "0", "0", "0", "0", "0", "0", "0", "2.040778768071901`", "0", "0", 
      "0", "0", "0", "0", 
      "0.`", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>"},
     {"0", "0", "0", "0", "0", "0", "0", "0", "0", "2.040783657507358`", "0", 
      "0", "0", "0", "0", "0", 
      "0.`", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>"},
     {"0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", 
      "2.0597553962610586`", "0", "0", "0", "0", 
      "0.`", "\<\" \"\>", "\<\" \"\>", "\<\" \"\>"},
     {"0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", 
      "2.0597027706602815`", "0", "0", "0", "0", "0", 
      "0.`", "\<\" \"\>", "\<\" \"\>"},
     {"0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", "0", 
      "2.059813025061571`", "0", "0", "0", "0", "0", "0", "0.`", "\<\" \"\>"},
     {"0", "1.827191665582771`", "0", "0", "0", "0", "0", "0", "0", "0", "0", 
      "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.568224901398209*^9, 3.568225018507615*^9, 3.5682250937309465`*^9, 
   3.568225796480981*^9, 3.5682259556216936`*^9, 3.568266119788582*^9, {
   3.568266196337917*^9, 3.568266218723956*^9}, 3.5682663260365443`*^9, 
   3.56826642606392*^9, 3.5682665067940617`*^9, 3.5682670928246913`*^9, 
   3.5682678627172437`*^9, {3.5682679314821644`*^9, 3.5682679543362045`*^9}, 
   3.568268004193892*^9, 3.568268075610818*^9, 3.568268185996611*^9, 
   3.568268325882057*^9, 3.568268879152629*^9, 3.568269028975292*^9, 
   3.568269542559394*^9, 3.5682696913056555`*^9, 3.5682713984010534`*^9, 
   3.5682724295160646`*^9, 3.5682726343600245`*^9, 3.5682730125982885`*^9, 
   3.5682734356554317`*^9, 3.568273517586776*^9, 3.568273637270186*^9, 
   3.568273995556015*^9, 3.5682764648063526`*^9, 3.5682765188136473`*^9, 
   3.5682766681215096`*^9, {3.568277700671723*^9, 3.568277718440154*^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", ",", "21"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"j", ",", "1", ",", "21"}], "}"}]}], "]"}]}], ";", 
   "\[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}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"u2", "=", 
   RowBox[{"{", 
    RowBox[{
    "0.1887", ",", "0.1887", ",", "0.1887", ",", "0.1887", ",", "0.2555", ",",
      "0.2555", ",", "0.2555", ",", "0.2555", ",", "0.2555", ",", "0.2555", 
     ",", "0.2555", ",", "0.2555", ",", "0.2555", ",", "0.563", ",", "0.563", 
     ",", "0.563", ",", "0.563", ",", "0.563", ",", "0.563", ",", "0.563", 
     ",", "0.50"}], "}"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.5679585418846006`*^9, 3.5679585738334565`*^9}, {
  3.567959075187137*^9, 3.5679590758579383`*^9}, {3.5681145227851562`*^9, 
  3.5681146476945763`*^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", ",", "21"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"j", ",", "1", ",", "21"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"MatrixForm", "[", "Gij", "]"}], ";"}], "*)"}]}]}], "Input",
 CellChangeTimes->{{3.567958596890297*^9, 3.567958603270708*^9}, {
   3.5679587604565845`*^9, 3.5679587707682023`*^9}, {3.5679588272871017`*^9, 
   3.567958828254304*^9}, {3.567959580066424*^9, 3.5679595814704266`*^9}, {
   3.567959624994503*^9, 3.5679596300177116`*^9}, {3.5679597079086485`*^9, 
   3.567959902877791*^9}, {3.5679599961659546`*^9, 3.5679600446196404`*^9}, {
   3.5679609611992497`*^9, 3.567960974006872*^9}, 3.567961975154231*^9, {
   3.5679622637079377`*^9, 3.5679622945179915`*^9}, 3.5681142413606625`*^9, 
   3.568121878697276*^9, {3.5681219332193723`*^9, 3.56812197189184*^9}}],

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", ",", "21"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"j", ",", "1", ",", "21"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"MatrixForm", "[", "Lij", "]"}], ";"}], "*)"}]}]}], "Input",
 CellChangeTimes->{{3.5679602339571724`*^9, 3.567960257419614*^9}, {
   3.567960343406965*^9, 3.567960353500182*^9}, {3.5679604155258913`*^9, 
   3.567960421422702*^9}, {3.567960854323462*^9, 3.567960874197897*^9}, {
   3.5679611596939983`*^9, 3.567961170270817*^9}, {3.5679614000124207`*^9, 
   3.567961407016833*^9}, {3.5679615066074076`*^9, 3.567961509649413*^9}, 
   3.5679616124379935`*^9, 3.5679616766477065`*^9, {3.5679617424486217`*^9, 
   3.5679617456310277`*^9}, 3.567961801276325*^9, {3.5679618582788253`*^9, 
   3.5679618970760937`*^9}, {3.5679620279135237`*^9, 3.567962030113127*^9}, 
   3.567962107894864*^9, {3.567962149468937*^9, 3.5679621572533503`*^9}, 
   3.5679623389156694`*^9, {3.567962645815009*^9, 3.567962648888214*^9}, {
   3.567962691647889*^9, 3.567962693020692*^9}, {3.568048161486014*^9, 
   3.568048162827616*^9}, 3.5680482096900983`*^9, 3.5681142588638926`*^9, 
   3.568121899632513*^9, 3.5681219958534822`*^9}],

Cell[BoxData[{
 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[{"j", "<", "14"}], ",", 
           RowBox[{"0.5", "*", 
            RowBox[{"Lij", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "j"}], "]"}], "]"}]}], ",", 
           RowBox[{"Lij", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "j"}], "]"}], "]"}]}], "]"}], ",", "0"}], 
        "]"}], ",", "\"\< \>\""}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "21"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"j", ",", "1", ",", "21"}], "}"}]}], "]"}]}], ";", " ", 
  RowBox[{"Qij", "=", "Lij"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"MatrixForm", "[", "Lij", "]"}], ";"}], "*)"}]}]}], "Input",
 CellChangeTimes->{{3.5679627617388124`*^9, 3.567962837461345*^9}, {
   3.567963118121438*^9, 3.567963118698639*^9}, {3.5679631520670977`*^9, 
   3.5679631541731014`*^9}, {3.5679631898191643`*^9, 
   3.5679632154656096`*^9}, {3.5679635823158536`*^9, 3.56796361468591*^9}, {
   3.567963691516045*^9, 3.567963693590849*^9}, {3.567963737957327*^9, 
   3.567963807299449*^9}, {3.568048114904332*^9, 3.5680481161835346`*^9}, 
   3.5680482221389203`*^9, {3.5680485003498087`*^9, 3.5680485123618298`*^9}, 
   3.5681142659931054`*^9, 3.5681219084621286`*^9}],

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", ",", "21"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"j", ",", "1", ",", "21"}], "}"}]}], "]"}]}], 
  ";"}], "\[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}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Tsigma", "=", 
  RowBox[{
   RowBox[{"Sum", "[", 
    RowBox[{
     RowBox[{"If", "[", 
      RowBox[{
       RowBox[{"i", "\[LessEqual]", "j"}], ",", 
       RowBox[{"Nij", "[", 
        RowBox[{"[", 
         RowBox[{"i", ",", "j"}], "]"}], "]"}], ",", "0"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "21"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"j", ",", "1", ",", "21"}], "}"}]}], "]"}], "+", 
   RowBox[{"Sum", "[", 
    RowBox[{
     RowBox[{"2.25", "/", 
      RowBox[{
       RowBox[{"u2", "[", 
        RowBox[{"[", "i", "]"}], "]"}], "^", "2"}]}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "13"}], "}"}]}], "]"}]}]}]], "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}}],

Cell[BoxData["601.9609558806221`"], "Output",
 CellChangeTimes->{
  3.568224901445009*^9, 3.5682250185544147`*^9, 3.568225093777747*^9, 
   3.5682257965277815`*^9, 3.5682259556840935`*^9, 3.5682661198197823`*^9, {
   3.5682661964003167`*^9, 3.568266218770756*^9}, 3.5682663260833445`*^9, 
   3.5682664261107206`*^9, 3.5682665068408623`*^9, 3.5682670928714914`*^9, 
   3.568267862764044*^9, {3.568267931528964*^9, 3.5682679543986044`*^9}, 
   3.5682680042562923`*^9, 3.568268075673218*^9, 3.5682681860434113`*^9, 
   3.568268325928857*^9, 3.5682688791994286`*^9, 3.568269029022092*^9, 
   3.5682695426217937`*^9, 3.568269691352455*^9, 3.5682713984478536`*^9, 
   3.5682724295784645`*^9, 3.5682726343912244`*^9, 3.568273012629489*^9, 
   3.568273435686632*^9, 3.568273517633576*^9, 3.568273637316986*^9, 
   3.5682739955872154`*^9, 3.568276464837552*^9, 3.568276518844847*^9, 
   3.5682766681527095`*^9, {3.568277700718523*^9, 3.5682777184869547`*^9}}]
}, Open  ]],

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

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

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

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"For", "[", 
  RowBox[{
   RowBox[{"i", "=", "2"}], ",", 
   RowBox[{"i", "<", "23"}], ",", 
   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}}],

Cell[CellGroupData[{

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"C\"\>", "\[InvisibleSpace]", "5", 
   "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"O\"\>", 
   "\[InvisibleSpace]", "1", "\[InvisibleSpace]", "\<\"  \"\>", 
   "\[InvisibleSpace]", "1.0600979904532752`"}],
  SequenceForm["C", 5, " ", "O", 1, "  ", 1.0600979904532752`],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.5682249014606094`*^9, 3.5682250185856147`*^9, 3.568225093808947*^9, 
   3.568225796558981*^9, 3.5682259556996937`*^9, 3.5682661198353825`*^9, {
   3.5682661964159164`*^9, 3.568266218801956*^9}, 3.5682663260989447`*^9, 
   3.56826642614192*^9, 3.568266506856462*^9, 3.5682670929026914`*^9, 
   3.5682678627952433`*^9, {3.5682679315601645`*^9, 3.5682679544142046`*^9}, 
   3.5682680042718925`*^9, 3.5682680756888175`*^9, 3.5682681860746117`*^9, 
   3.568268325944457*^9, 3.568268879230629*^9, 3.5682690290532923`*^9, 
   3.568269542637394*^9, 3.5682696913836555`*^9, 3.5682713984790535`*^9, 
   3.568272429609665*^9, 3.5682726344224243`*^9, 3.568273012660689*^9, 
   3.568273435702232*^9, 3.5682735176491756`*^9, 3.5682736373325863`*^9, 
   3.5682739956184154`*^9, 3.5682764648687525`*^9, 3.568276518876047*^9, 
   3.5682766681683097`*^9, {3.568277700734123*^9, 3.5682777185025544`*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"C\"\>", "\[InvisibleSpace]", "6", 
   "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", 
   "\[InvisibleSpace]", "5", "\[InvisibleSpace]", "\<\"  \"\>", 
   "\[InvisibleSpace]", "1.0688976046225591`"}],
  SequenceForm["C", 6, " ", "C", 5, "  ", 1.0688976046225591`],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.5682249014606094`*^9, 3.5682250185856147`*^9, 3.568225093808947*^9, 
   3.568225796558981*^9, 3.5682259556996937`*^9, 3.5682661198353825`*^9, {
   3.5682661964159164`*^9, 3.568266218801956*^9}, 3.5682663260989447`*^9, 
   3.56826642614192*^9, 3.568266506856462*^9, 3.5682670929026914`*^9, 
   3.5682678627952433`*^9, {3.5682679315601645`*^9, 3.5682679544142046`*^9}, 
   3.5682680042718925`*^9, 3.5682680756888175`*^9, 3.5682681860746117`*^9, 
   3.568268325944457*^9, 3.568268879230629*^9, 3.5682690290532923`*^9, 
   3.568269542637394*^9, 3.5682696913836555`*^9, 3.5682713984790535`*^9, 
   3.568272429609665*^9, 3.5682726344224243`*^9, 3.568273012660689*^9, 
   3.568273435702232*^9, 3.5682735176491756`*^9, 3.5682736373325863`*^9, 
   3.5682739956184154`*^9, 3.5682764648687525`*^9, 3.568276518876047*^9, 
   3.5682766681683097`*^9, {3.568277700734123*^9, 3.5682777185025544`*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"C\"\>", "\[InvisibleSpace]", "7", 
   "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", 
   "\[InvisibleSpace]", "5", "\[InvisibleSpace]", "\<\"  \"\>", 
   "\[InvisibleSpace]", "1.0547571207979058`"}],
  SequenceForm["C", 7, " ", "C", 5, "  ", 1.0547571207979058`],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.5682249014606094`*^9, 3.5682250185856147`*^9, 3.568225093808947*^9, 
   3.568225796558981*^9, 3.5682259556996937`*^9, 3.5682661198353825`*^9, {
   3.5682661964159164`*^9, 3.568266218801956*^9}, 3.5682663260989447`*^9, 
   3.56826642614192*^9, 3.568266506856462*^9, 3.5682670929026914`*^9, 
   3.5682678627952433`*^9, {3.5682679315601645`*^9, 3.5682679544142046`*^9}, 
   3.5682680042718925`*^9, 3.5682680756888175`*^9, 3.5682681860746117`*^9, 
   3.568268325944457*^9, 3.568268879230629*^9, 3.5682690290532923`*^9, 
   3.568269542637394*^9, 3.5682696913836555`*^9, 3.5682713984790535`*^9, 
   3.568272429609665*^9, 3.5682726344224243`*^9, 3.568273012660689*^9, 
   3.568273435702232*^9, 3.5682735176491756`*^9, 3.5682736373325863`*^9, 
   3.5682739956184154`*^9, 3.5682764648687525`*^9, 3.568276518876047*^9, 
   3.5682766681683097`*^9, {3.568277700734123*^9, 3.5682777185025544`*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"C\"\>", "\[InvisibleSpace]", "8", 
   "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", 
   "\[InvisibleSpace]", "6", "\[InvisibleSpace]", "\<\"  \"\>", 
   "\[InvisibleSpace]", "1.0648108305837336`"}],
  SequenceForm["C", 8, " ", "C", 6, "  ", 1.0648108305837336`],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.5682249014606094`*^9, 3.5682250185856147`*^9, 3.568225093808947*^9, 
   3.568225796558981*^9, 3.5682259556996937`*^9, 3.5682661198353825`*^9, {
   3.5682661964159164`*^9, 3.568266218801956*^9}, 3.5682663260989447`*^9, 
   3.56826642614192*^9, 3.568266506856462*^9, 3.5682670929026914`*^9, 
   3.5682678627952433`*^9, {3.5682679315601645`*^9, 3.5682679544142046`*^9}, 
   3.5682680042718925`*^9, 3.5682680756888175`*^9, 3.5682681860746117`*^9, 
   3.568268325944457*^9, 3.568268879230629*^9, 3.5682690290532923`*^9, 
   3.568269542637394*^9, 3.5682696913836555`*^9, 3.5682713984790535`*^9, 
   3.568272429609665*^9, 3.5682726344224243`*^9, 3.568273012660689*^9, 
   3.568273435702232*^9, 3.5682735176491756`*^9, 3.5682736373325863`*^9, 
   3.5682739956184154`*^9, 3.5682764648687525`*^9, 3.568276518876047*^9, 
   3.5682766681683097`*^9, {3.568277700734123*^9, 3.5682777185025544`*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"C\"\>", "\[InvisibleSpace]", "9", 
   "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", 
   "\[InvisibleSpace]", "7", "\[InvisibleSpace]", "\<\"  \"\>", 
   "\[InvisibleSpace]", "1.0491986570390401`"}],
  SequenceForm["C", 9, " ", "C", 7, "  ", 1.0491986570390401`],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.5682249014606094`*^9, 3.5682250185856147`*^9, 3.568225093808947*^9, 
   3.568225796558981*^9, 3.5682259556996937`*^9, 3.5682661198353825`*^9, {
   3.5682661964159164`*^9, 3.568266218801956*^9}, 3.5682663260989447`*^9, 
   3.56826642614192*^9, 3.568266506856462*^9, 3.5682670929026914`*^9, 
   3.5682678627952433`*^9, {3.5682679315601645`*^9, 3.5682679544142046`*^9}, 
   3.5682680042718925`*^9, 3.5682680756888175`*^9, 3.5682681860746117`*^9, 
   3.568268325944457*^9, 3.568268879230629*^9, 3.5682690290532923`*^9, 
   3.568269542637394*^9, 3.5682696913836555`*^9, 3.5682713984790535`*^9, 
   3.568272429609665*^9, 3.5682726344224243`*^9, 3.568273012660689*^9, 
   3.568273435702232*^9, 3.5682735176491756`*^9, 3.5682736373325863`*^9, 
   3.5682739956184154`*^9, 3.5682764648687525`*^9, 3.568276518876047*^9, 
   3.5682766681683097`*^9, {3.568277700734123*^9, 3.568277718518154*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"C\"\>", "\[InvisibleSpace]", "10", 
   "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", 
   "\[InvisibleSpace]", "8", "\[InvisibleSpace]", "\<\"  \"\>", 
   "\[InvisibleSpace]", "1.0464782521399887`"}],
  SequenceForm["C", 10, " ", "C", 8, "  ", 1.0464782521399887`],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.5682249014606094`*^9, 3.5682250185856147`*^9, 3.568225093808947*^9, 
   3.568225796558981*^9, 3.5682259556996937`*^9, 3.5682661198353825`*^9, {
   3.5682661964159164`*^9, 3.568266218801956*^9}, 3.5682663260989447`*^9, 
   3.56826642614192*^9, 3.568266506856462*^9, 3.5682670929026914`*^9, 
   3.5682678627952433`*^9, {3.5682679315601645`*^9, 3.5682679544142046`*^9}, 
   3.5682680042718925`*^9, 3.5682680756888175`*^9, 3.5682681860746117`*^9, 
   3.568268325944457*^9, 3.568268879230629*^9, 3.5682690290532923`*^9, 
   3.568269542637394*^9, 3.5682696913836555`*^9, 3.5682713984790535`*^9, 
   3.568272429609665*^9, 3.5682726344224243`*^9, 3.568273012660689*^9, 
   3.568273435702232*^9, 3.5682735176491756`*^9, 3.5682736373325863`*^9, 
   3.5682739956184154`*^9, 3.5682764648687525`*^9, 3.568276518876047*^9, 
   3.5682766681683097`*^9, {3.568277700734123*^9, 3.568277718518154*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"C\"\>", "\[InvisibleSpace]", "10", 
   "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", 
   "\[InvisibleSpace]", "9", "\[InvisibleSpace]", "\<\"  \"\>", 
   "\[InvisibleSpace]", "1.0535140870236`"}],
  SequenceForm["C", 10, " ", "C", 9, "  ", 1.0535140870236],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.5682249014606094`*^9, 3.5682250185856147`*^9, 3.568225093808947*^9, 
   3.568225796558981*^9, 3.5682259556996937`*^9, 3.5682661198353825`*^9, {
   3.5682661964159164`*^9, 3.568266218801956*^9}, 3.5682663260989447`*^9, 
   3.56826642614192*^9, 3.568266506856462*^9, 3.5682670929026914`*^9, 
   3.5682678627952433`*^9, {3.5682679315601645`*^9, 3.5682679544142046`*^9}, 
   3.5682680042718925`*^9, 3.5682680756888175`*^9, 3.5682681860746117`*^9, 
   3.568268325944457*^9, 3.568268879230629*^9, 3.5682690290532923`*^9, 
   3.568269542637394*^9, 3.5682696913836555`*^9, 3.5682713984790535`*^9, 
   3.568272429609665*^9, 3.5682726344224243`*^9, 3.568273012660689*^9, 
   3.568273435702232*^9, 3.5682735176491756`*^9, 3.5682736373325863`*^9, 
   3.5682739956184154`*^9, 3.5682764648687525`*^9, 3.568276518876047*^9, 
   3.5682766681683097`*^9, {3.568277700734123*^9, 3.568277718518154*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"C\"\>", "\[InvisibleSpace]", "11", 
   "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"O\"\>", 
   "\[InvisibleSpace]", "2", "\[InvisibleSpace]", "\<\"  \"\>", 
   "\[InvisibleSpace]", "1.0510023931152543`"}],
  SequenceForm["C", 11, " ", "O", 2, "  ", 1.0510023931152543`],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.5682249014606094`*^9, 3.5682250185856147`*^9, 3.568225093808947*^9, 
   3.568225796558981*^9, 3.5682259556996937`*^9, 3.5682661198353825`*^9, {
   3.5682661964159164`*^9, 3.568266218801956*^9}, 3.5682663260989447`*^9, 
   3.56826642614192*^9, 3.568266506856462*^9, 3.5682670929026914`*^9, 
   3.5682678627952433`*^9, {3.5682679315601645`*^9, 3.5682679544142046`*^9}, 
   3.5682680042718925`*^9, 3.5682680756888175`*^9, 3.5682681860746117`*^9, 
   3.568268325944457*^9, 3.568268879230629*^9, 3.5682690290532923`*^9, 
   3.568269542637394*^9, 3.5682696913836555`*^9, 3.5682713984790535`*^9, 
   3.568272429609665*^9, 3.5682726344224243`*^9, 3.568273012660689*^9, 
   3.568273435702232*^9, 3.5682735176491756`*^9, 3.5682736373325863`*^9, 
   3.5682739956184154`*^9, 3.5682764648687525`*^9, 3.568276518876047*^9, 
   3.5682766681683097`*^9, {3.568277700734123*^9, 3.568277718518154*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"C\"\>", "\[InvisibleSpace]", "11", 
   "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"O\"\>", 
   "\[InvisibleSpace]", "3", "\[InvisibleSpace]", "\<\"  \"\>", 
   "\[InvisibleSpace]", "0.9259475168654768`"}],
  SequenceForm["C", 11, " ", "O", 3, "  ", 0.9259475168654768],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.5682249014606094`*^9, 3.5682250185856147`*^9, 3.568225093808947*^9, 
   3.568225796558981*^9, 3.5682259556996937`*^9, 3.5682661198353825`*^9, {
   3.5682661964159164`*^9, 3.568266218801956*^9}, 3.5682663260989447`*^9, 
   3.56826642614192*^9, 3.568266506856462*^9, 3.5682670929026914`*^9, 
   3.5682678627952433`*^9, {3.5682679315601645`*^9, 3.5682679544142046`*^9}, 
   3.5682680042718925`*^9, 3.5682680756888175`*^9, 3.5682681860746117`*^9, 
   3.568268325944457*^9, 3.568268879230629*^9, 3.5682690290532923`*^9, 
   3.568269542637394*^9, 3.5682696913836555`*^9, 3.5682713984790535`*^9, 
   3.568272429609665*^9, 3.5682726344224243`*^9, 3.568273012660689*^9, 
   3.568273435702232*^9, 3.5682735176491756`*^9, 3.5682736373325863`*^9, 
   3.5682739956184154`*^9, 3.5682764648687525`*^9, 3.568276518876047*^9, 
   3.5682766681683097`*^9, {3.568277700734123*^9, 3.568277718518154*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"C\"\>", "\[InvisibleSpace]", "11", 
   "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", 
   "\[InvisibleSpace]", "6", "\[InvisibleSpace]", "\<\"  \"\>", 
   "\[InvisibleSpace]", "1.13773739655973`"}],
  SequenceForm["C", 11, " ", "C", 6, "  ", 1.13773739655973],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.5682249014606094`*^9, 3.5682250185856147`*^9, 3.568225093808947*^9, 
   3.568225796558981*^9, 3.5682259556996937`*^9, 3.5682661198353825`*^9, {
   3.5682661964159164`*^9, 3.568266218801956*^9}, 3.5682663260989447`*^9, 
   3.56826642614192*^9, 3.568266506856462*^9, 3.5682670929026914`*^9, 
   3.5682678627952433`*^9, {3.5682679315601645`*^9, 3.5682679544142046`*^9}, 
   3.5682680042718925`*^9, 3.5682680756888175`*^9, 3.5682681860746117`*^9, 
   3.568268325944457*^9, 3.568268879230629*^9, 3.5682690290532923`*^9, 
   3.568269542637394*^9, 3.5682696913836555`*^9, 3.5682713984790535`*^9, 
   3.568272429609665*^9, 3.5682726344224243`*^9, 3.568273012660689*^9, 
   3.568273435702232*^9, 3.5682735176491756`*^9, 3.5682736373325863`*^9, 
   3.5682739956184154`*^9, 3.5682764648687525`*^9, 3.568276518876047*^9, 
   3.5682766681683097`*^9, {3.568277700734123*^9, 3.568277718518154*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"C\"\>", "\[InvisibleSpace]", "12", 
   "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"O\"\>", 
   "\[InvisibleSpace]", "1", "\[InvisibleSpace]", "\<\"  \"\>", 
   "\[InvisibleSpace]", "1.046964530990104`"}],
  SequenceForm["C", 12, " ", "O", 1, "  ", 1.046964530990104],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.5682249014606094`*^9, 3.5682250185856147`*^9, 3.568225093808947*^9, 
   3.568225796558981*^9, 3.5682259556996937`*^9, 3.5682661198353825`*^9, {
   3.5682661964159164`*^9, 3.568266218801956*^9}, 3.5682663260989447`*^9, 
   3.56826642614192*^9, 3.568266506856462*^9, 3.5682670929026914`*^9, 
   3.5682678627952433`*^9, {3.5682679315601645`*^9, 3.5682679544142046`*^9}, 
   3.5682680042718925`*^9, 3.5682680756888175`*^9, 3.5682681860746117`*^9, 
   3.568268325944457*^9, 3.568268879230629*^9, 3.5682690290532923`*^9, 
   3.568269542637394*^9, 3.5682696913836555`*^9, 3.5682713984790535`*^9, 
   3.568272429609665*^9, 3.5682726344224243`*^9, 3.568273012660689*^9, 
   3.568273435702232*^9, 3.5682735176491756`*^9, 3.5682736373325863`*^9, 
   3.5682739956184154`*^9, 3.5682764648687525`*^9, 3.568276518876047*^9, 
   3.5682766681683097`*^9, {3.568277700734123*^9, 3.568277718518154*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"C\"\>", "\[InvisibleSpace]", "12", 
   "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"O\"\>", 
   "\[InvisibleSpace]", "4", "\[InvisibleSpace]", "\<\"  \"\>", 
   "\[InvisibleSpace]", "0.9195510356822784`"}],
  SequenceForm["C", 12, " ", "O", 4, "  ", 0.9195510356822784],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.5682249014606094`*^9, 3.5682250185856147`*^9, 3.568225093808947*^9, 
   3.568225796558981*^9, 3.5682259556996937`*^9, 3.5682661198353825`*^9, {
   3.5682661964159164`*^9, 3.568266218801956*^9}, 3.5682663260989447`*^9, 
   3.56826642614192*^9, 3.568266506856462*^9, 3.5682670929026914`*^9, 
   3.5682678627952433`*^9, {3.5682679315601645`*^9, 3.5682679544142046`*^9}, 
   3.5682680042718925`*^9, 3.5682680756888175`*^9, 3.5682681860746117`*^9, 
   3.568268325944457*^9, 3.568268879230629*^9, 3.5682690290532923`*^9, 
   3.568269542637394*^9, 3.5682696913836555`*^9, 3.5682713984790535`*^9, 
   3.568272429609665*^9, 3.5682726344224243`*^9, 3.568273012660689*^9, 
   3.568273435702232*^9, 3.5682735176491756`*^9, 3.5682736373325863`*^9, 
   3.5682739956184154`*^9, 3.5682764648687525`*^9, 3.568276518876047*^9, 
   3.5682766681683097`*^9, {3.568277700734123*^9, 3.5682777185337543`*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"C\"\>", "\[InvisibleSpace]", "13", 
   "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", 
   "\[InvisibleSpace]", "12", "\[InvisibleSpace]", "\<\"  \"\>", 
   "\[InvisibleSpace]", "1.1684155739346305`"}],
  SequenceForm["C", 13, " ", "C", 12, "  ", 1.1684155739346305`],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.5682249014606094`*^9, 3.5682250185856147`*^9, 3.568225093808947*^9, 
   3.568225796558981*^9, 3.5682259556996937`*^9, 3.5682661198353825`*^9, {
   3.5682661964159164`*^9, 3.568266218801956*^9}, 3.5682663260989447`*^9, 
   3.56826642614192*^9, 3.568266506856462*^9, 3.5682670929026914`*^9, 
   3.5682678627952433`*^9, {3.5682679315601645`*^9, 3.5682679544142046`*^9}, 
   3.5682680042718925`*^9, 3.5682680756888175`*^9, 3.5682681860746117`*^9, 
   3.568268325944457*^9, 3.568268879230629*^9, 3.5682690290532923`*^9, 
   3.568269542637394*^9, 3.5682696913836555`*^9, 3.5682713984790535`*^9, 
   3.568272429609665*^9, 3.5682726344224243`*^9, 3.568273012660689*^9, 
   3.568273435702232*^9, 3.5682735176491756`*^9, 3.5682736373325863`*^9, 
   3.5682739956184154`*^9, 3.5682764648687525`*^9, 3.568276518876047*^9, 
   3.5682766681683097`*^9, {3.568277700734123*^9, 3.5682777185337543`*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"H\"\>", "\[InvisibleSpace]", "14", 
   "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", 
   "\[InvisibleSpace]", "7", "\[InvisibleSpace]", "\<\"  \"\>", 
   "\[InvisibleSpace]", "1.2223901257255825`"}],
  SequenceForm["H", 14, " ", "C", 7, "  ", 1.2223901257255825`],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.5682249014606094`*^9, 3.5682250185856147`*^9, 3.568225093808947*^9, 
   3.568225796558981*^9, 3.5682259556996937`*^9, 3.5682661198353825`*^9, {
   3.5682661964159164`*^9, 3.568266218801956*^9}, 3.5682663260989447`*^9, 
   3.56826642614192*^9, 3.568266506856462*^9, 3.5682670929026914`*^9, 
   3.5682678627952433`*^9, {3.5682679315601645`*^9, 3.5682679544142046`*^9}, 
   3.5682680042718925`*^9, 3.5682680756888175`*^9, 3.5682681860746117`*^9, 
   3.568268325944457*^9, 3.568268879230629*^9, 3.5682690290532923`*^9, 
   3.568269542637394*^9, 3.5682696913836555`*^9, 3.5682713984790535`*^9, 
   3.568272429609665*^9, 3.5682726344224243`*^9, 3.568273012660689*^9, 
   3.568273435702232*^9, 3.5682735176491756`*^9, 3.5682736373325863`*^9, 
   3.5682739956184154`*^9, 3.5682764648687525`*^9, 3.568276518876047*^9, 
   3.5682766681683097`*^9, {3.568277700734123*^9, 3.5682777185337543`*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"H\"\>", "\[InvisibleSpace]", "15", 
   "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", 
   "\[InvisibleSpace]", "8", "\[InvisibleSpace]", "\<\"  \"\>", 
   "\[InvisibleSpace]", "1.2224810523377199`"}],
  SequenceForm["H", 15, " ", "C", 8, "  ", 1.2224810523377199`],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.5682249014606094`*^9, 3.5682250185856147`*^9, 3.568225093808947*^9, 
   3.568225796558981*^9, 3.5682259556996937`*^9, 3.5682661198353825`*^9, {
   3.5682661964159164`*^9, 3.568266218801956*^9}, 3.5682663260989447`*^9, 
   3.56826642614192*^9, 3.568266506856462*^9, 3.5682670929026914`*^9, 
   3.5682678627952433`*^9, {3.5682679315601645`*^9, 3.5682679544142046`*^9}, 
   3.5682680042718925`*^9, 3.5682680756888175`*^9, 3.5682681860746117`*^9, 
   3.568268325944457*^9, 3.568268879230629*^9, 3.5682690290532923`*^9, 
   3.568269542637394*^9, 3.5682696913836555`*^9, 3.5682713984790535`*^9, 
   3.568272429609665*^9, 3.5682726344224243`*^9, 3.568273012660689*^9, 
   3.568273435702232*^9, 3.5682735176491756`*^9, 3.5682736373325863`*^9, 
   3.5682739956184154`*^9, 3.5682764648687525`*^9, 3.568276518876047*^9, 
   3.5682766681683097`*^9, {3.568277700734123*^9, 3.5682777185337543`*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"H\"\>", "\[InvisibleSpace]", "16", 
   "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", 
   "\[InvisibleSpace]", "9", "\[InvisibleSpace]", "\<\"  \"\>", 
   "\[InvisibleSpace]", "1.2222787680719012`"}],
  SequenceForm["H", 16, " ", "C", 9, "  ", 1.2222787680719012`],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.5682249014606094`*^9, 3.5682250185856147`*^9, 3.568225093808947*^9, 
   3.568225796558981*^9, 3.5682259556996937`*^9, 3.5682661198353825`*^9, {
   3.5682661964159164`*^9, 3.568266218801956*^9}, 3.5682663260989447`*^9, 
   3.56826642614192*^9, 3.568266506856462*^9, 3.5682670929026914`*^9, 
   3.5682678627952433`*^9, {3.5682679315601645`*^9, 3.5682679544142046`*^9}, 
   3.5682680042718925`*^9, 3.5682680756888175`*^9, 3.5682681860746117`*^9, 
   3.568268325944457*^9, 3.568268879230629*^9, 3.5682690290532923`*^9, 
   3.568269542637394*^9, 3.5682696913836555`*^9, 3.5682713984790535`*^9, 
   3.568272429609665*^9, 3.5682726344224243`*^9, 3.568273012660689*^9, 
   3.568273435702232*^9, 3.5682735176491756`*^9, 3.5682736373325863`*^9, 
   3.5682739956184154`*^9, 3.5682764648687525`*^9, 3.568276518876047*^9, 
   3.5682766681683097`*^9, {3.568277700734123*^9, 3.5682777185337543`*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"H\"\>", "\[InvisibleSpace]", "17", 
   "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", 
   "\[InvisibleSpace]", "10", "\[InvisibleSpace]", "\<\"  \"\>", 
   "\[InvisibleSpace]", "1.2222836575073581`"}],
  SequenceForm["H", 17, " ", "C", 10, "  ", 1.2222836575073581`],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.5682249014606094`*^9, 3.5682250185856147`*^9, 3.568225093808947*^9, 
   3.568225796558981*^9, 3.5682259556996937`*^9, 3.5682661198353825`*^9, {
   3.5682661964159164`*^9, 3.568266218801956*^9}, 3.5682663260989447`*^9, 
   3.56826642614192*^9, 3.568266506856462*^9, 3.5682670929026914`*^9, 
   3.5682678627952433`*^9, {3.5682679315601645`*^9, 3.5682679544142046`*^9}, 
   3.5682680042718925`*^9, 3.5682680756888175`*^9, 3.5682681860746117`*^9, 
   3.568268325944457*^9, 3.568268879230629*^9, 3.5682690290532923`*^9, 
   3.568269542637394*^9, 3.5682696913836555`*^9, 3.5682713984790535`*^9, 
   3.568272429609665*^9, 3.5682726344224243`*^9, 3.568273012660689*^9, 
   3.568273435702232*^9, 3.5682735176491756`*^9, 3.5682736373325863`*^9, 
   3.5682739956184154`*^9, 3.5682764648687525`*^9, 3.568276518876047*^9, 
   3.5682766681683097`*^9, {3.568277700734123*^9, 3.5682777185337543`*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"H\"\>", "\[InvisibleSpace]", "18", 
   "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", 
   "\[InvisibleSpace]", "13", "\[InvisibleSpace]", "\<\"  \"\>", 
   "\[InvisibleSpace]", "1.2412553962610586`"}],
  SequenceForm["H", 18, " ", "C", 13, "  ", 1.2412553962610586`],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.5682249014606094`*^9, 3.5682250185856147`*^9, 3.568225093808947*^9, 
   3.568225796558981*^9, 3.5682259556996937`*^9, 3.5682661198353825`*^9, {
   3.5682661964159164`*^9, 3.568266218801956*^9}, 3.5682663260989447`*^9, 
   3.56826642614192*^9, 3.568266506856462*^9, 3.5682670929026914`*^9, 
   3.5682678627952433`*^9, {3.5682679315601645`*^9, 3.5682679544142046`*^9}, 
   3.5682680042718925`*^9, 3.5682680756888175`*^9, 3.5682681860746117`*^9, 
   3.568268325944457*^9, 3.568268879230629*^9, 3.5682690290532923`*^9, 
   3.568269542637394*^9, 3.5682696913836555`*^9, 3.5682713984790535`*^9, 
   3.568272429609665*^9, 3.5682726344224243`*^9, 3.568273012660689*^9, 
   3.568273435702232*^9, 3.5682735176491756`*^9, 3.5682736373325863`*^9, 
   3.5682739956184154`*^9, 3.5682764648687525`*^9, 3.568276518876047*^9, 
   3.5682766681683097`*^9, {3.568277700734123*^9, 3.5682777185337543`*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"H\"\>", "\[InvisibleSpace]", "19", 
   "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", 
   "\[InvisibleSpace]", "13", "\[InvisibleSpace]", "\<\"  \"\>", 
   "\[InvisibleSpace]", "1.2412027706602815`"}],
  SequenceForm["H", 19, " ", "C", 13, "  ", 1.2412027706602815`],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.5682249014606094`*^9, 3.5682250185856147`*^9, 3.568225093808947*^9, 
   3.568225796558981*^9, 3.5682259556996937`*^9, 3.5682661198353825`*^9, {
   3.5682661964159164`*^9, 3.568266218801956*^9}, 3.5682663260989447`*^9, 
   3.56826642614192*^9, 3.568266506856462*^9, 3.5682670929026914`*^9, 
   3.5682678627952433`*^9, {3.5682679315601645`*^9, 3.5682679544142046`*^9}, 
   3.5682680042718925`*^9, 3.5682680756888175`*^9, 3.5682681860746117`*^9, 
   3.568268325944457*^9, 3.568268879230629*^9, 3.5682690290532923`*^9, 
   3.568269542637394*^9, 3.5682696913836555`*^9, 3.5682713984790535`*^9, 
   3.568272429609665*^9, 3.5682726344224243`*^9, 3.568273012660689*^9, 
   3.568273435702232*^9, 3.5682735176491756`*^9, 3.5682736373325863`*^9, 
   3.5682739956184154`*^9, 3.5682764648687525`*^9, 3.568276518876047*^9, 
   3.5682766681683097`*^9, {3.568277700734123*^9, 3.5682777185337543`*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"H\"\>", "\[InvisibleSpace]", "20", 
   "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", 
   "\[InvisibleSpace]", "13", "\[InvisibleSpace]", "\<\"  \"\>", 
   "\[InvisibleSpace]", "1.2413130250615712`"}],
  SequenceForm["H", 20, " ", "C", 13, "  ", 1.2413130250615712`],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.5682249014606094`*^9, 3.5682250185856147`*^9, 3.568225093808947*^9, 
   3.568225796558981*^9, 3.5682259556996937`*^9, 3.5682661198353825`*^9, {
   3.5682661964159164`*^9, 3.568266218801956*^9}, 3.5682663260989447`*^9, 
   3.56826642614192*^9, 3.568266506856462*^9, 3.5682670929026914`*^9, 
   3.5682678627952433`*^9, {3.5682679315601645`*^9, 3.5682679544142046`*^9}, 
   3.5682680042718925`*^9, 3.5682680756888175`*^9, 3.5682681860746117`*^9, 
   3.568268325944457*^9, 3.568268879230629*^9, 3.5682690290532923`*^9, 
   3.568269542637394*^9, 3.5682696913836555`*^9, 3.5682713984790535`*^9, 
   3.568272429609665*^9, 3.5682726344224243`*^9, 3.568273012660689*^9, 
   3.568273435702232*^9, 3.5682735176491756`*^9, 3.5682736373325863`*^9, 
   3.5682739956184154`*^9, 3.5682764648687525`*^9, 3.568276518876047*^9, 
   3.5682766681683097`*^9, {3.568277700734123*^9, 3.5682777185493546`*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"H\"\>", "\[InvisibleSpace]", "21", 
   "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"O\"\>", 
   "\[InvisibleSpace]", "2", "\[InvisibleSpace]", "\<\"  \"\>", 
   "\[InvisibleSpace]", "1.1384916655827708`"}],
  SequenceForm["H", 21, " ", "O", 2, "  ", 1.1384916655827708`],
  Editable->False]], "Print",
 CellChangeTimes->{
  3.5682249014606094`*^9, 3.5682250185856147`*^9, 3.568225093808947*^9, 
   3.568225796558981*^9, 3.5682259556996937`*^9, 3.5682661198353825`*^9, {
   3.5682661964159164`*^9, 3.568266218801956*^9}, 3.5682663260989447`*^9, 
   3.56826642614192*^9, 3.568266506856462*^9, 3.5682670929026914`*^9, 
   3.5682678627952433`*^9, {3.5682679315601645`*^9, 3.5682679544142046`*^9}, 
   3.5682680042718925`*^9, 3.5682680756888175`*^9, 3.5682681860746117`*^9, 
   3.568268325944457*^9, 3.568268879230629*^9, 3.5682690290532923`*^9, 
   3.568269542637394*^9, 3.5682696913836555`*^9, 3.5682713984790535`*^9, 
   3.568272429609665*^9, 3.5682726344224243`*^9, 3.568273012660689*^9, 
   3.568273435702232*^9, 3.5682735176491756`*^9, 3.5682736373325863`*^9, 
   3.5682739956184154`*^9, 3.5682764648687525`*^9, 3.568276518876047*^9, 
   3.5682766681683097`*^9, {3.568277700734123*^9, 3.5682777185493546`*^9}}]
}, Open  ]]
}, Open  ]],

Cell[BoxData[{
 RowBox[{
  RowBox[{"DeleteFile", "[", " ", "\"\<C:\\\\Kimb_Mathem\\\\tmp4\>\"", "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"For", "[", 
  RowBox[{
   RowBox[{"i", "=", "2"}], ",", 
   RowBox[{"i", "<", "23"}], ",", 
   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}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"a", "=", 
  RowBox[{
  "ReadList", "[", "\"\<C:\\\\Kimb_Mathem\\\\tmp4\>\"", "]"}]}]], "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}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"5", ",", "1", ",", "1.0600979904532752`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"6", ",", "5", ",", "1.0688976046225591`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"7", ",", "5", ",", "1.0547571207979058`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"8", ",", "6", ",", "1.0648108305837336`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"9", ",", "7", ",", "1.0491986570390401`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"10", ",", "8", ",", "1.0464782521399887`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"10", ",", "9", ",", "1.0535140870236`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"11", ",", "2", ",", "1.0510023931152543`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"11", ",", "3", ",", "0.9259475168654768`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"11", ",", "6", ",", "1.13773739655973`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"12", ",", "1", ",", "1.046964530990104`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"12", ",", "4", ",", "0.9195510356822784`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"13", ",", "12", ",", "1.1684155739346305`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"14", ",", "7", ",", "1.2223901257255825`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"15", ",", "8", ",", "1.2224810523377199`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"16", ",", "9", ",", "1.2222787680719012`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"17", ",", "10", ",", "1.2222836575073581`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"18", ",", "13", ",", "1.2412553962610586`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"19", ",", "13", ",", "1.2412027706602815`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"20", ",", "13", ",", "1.2413130250615712`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"21", ",", "2", ",", "1.1384916655827708`"}], "}"}]}], 
  "}"}]], "Output",
 CellChangeTimes->{
  3.5682249015230093`*^9, 3.568225018663615*^9, 3.568225093886947*^9, 
   3.5682257966057816`*^9, 3.5682259557620935`*^9, 3.5682661198977823`*^9, {
   3.568266196462717*^9, 3.5682662188643565`*^9}, 3.5682663261613445`*^9, 
   3.5682664261887207`*^9, 3.5682665069344625`*^9, 3.5682670929650917`*^9, 
   3.568267862842044*^9, {3.568267931606964*^9, 3.5682679544766045`*^9}, 
   3.5682680043342924`*^9, 3.568268075751218*^9, 3.5682681861214113`*^9, 
   3.568268326006857*^9, 3.5682688792930293`*^9, 3.568269029100092*^9, 
   3.5682695426997943`*^9, 3.5682696914304557`*^9, 3.568271398525854*^9, 
   3.568272429672065*^9, 3.568272634484825*^9, 3.568273012723089*^9, 
   3.568273435749032*^9, 3.568273517711576*^9, 3.568273637394986*^9, 
   3.5682739956652155`*^9, 3.568276464915552*^9, 3.5682765189228477`*^9, 
   3.5682766682307096`*^9, {3.568277700796523*^9, 3.5682777185649543`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{" ", 
  RowBox[{
   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", ",", "13"}], "}"}]}], "]"}]}], ";"}], "\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", ",", "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", ",", "13"}], "}"}]}], "]"}]}], ";"}], "\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", ",", "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", ",", "13"}], "}"}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", "\n", 
   RowBox[{
    RowBox[{"Plot7", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"Graphics3D", "[", 
        RowBox[{"{", 
         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", ",", "13"}], "}"}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{"q", "=", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{
         RowBox[{"a", "[", 
          RowBox[{"[", 
           RowBox[{"i", ",", "3"}], "]"}], "]"}], "+", 
         RowBox[{"u2", "[", 
          RowBox[{"[", "a2", "]"}], "]"}]}], ")"}], "/", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"u2", "[", 
          RowBox[{"[", "a1", "]"}], "]"}], "+", 
         RowBox[{"a", "[", 
          RowBox[{"[", 
           RowBox[{"i", ",", "3"}], "]"}], "]"}], "+", 
         RowBox[{"u2", "[", 
          RowBox[{"[", "a2", "]"}], "]"}]}], ")"}]}]}], ";"}], "*)"}], 
   "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{"q", "=", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"a3", "+", 
         RowBox[{"u2", "[", 
          RowBox[{"[", "a2", "]"}], "]"}]}], ")"}], "/", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"u2", "[", 
          RowBox[{"[", "a1", "]"}], "]"}], "+", "a3", "+", 
         RowBox[{"u2", "[", 
          RowBox[{"[", "a2", "]"}], "]"}]}], ")"}]}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"wx", "=", 
      RowBox[{
       RowBox[{"q", "*", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"v", "[", 
           RowBox[{"[", 
            RowBox[{"a1", ",", "1"}], "]"}], "]"}], "-", 
          RowBox[{"v", "[", 
           RowBox[{"[", 
            RowBox[{"a2", ",", "1"}], "]"}], "]"}]}], ")"}]}], "+", 
       RowBox[{"v", "[", 
        RowBox[{"[", 
         RowBox[{"a2", ",", "1"}], "]"}], "]"}]}]}], ";", "\n", 
     RowBox[{"wy", "=", 
      RowBox[{
       RowBox[{"q", "*", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"v", "[", 
           RowBox[{"[", 
            RowBox[{"a1", ",", "2"}], "]"}], "]"}], "-", 
          RowBox[{"v", "[", 
           RowBox[{"[", 
            RowBox[{"a2", ",", "2"}], "]"}], "]"}]}], ")"}]}], "+", 
       RowBox[{"v", "[", 
        RowBox[{"[", 
         RowBox[{"a2", ",", "2"}], "]"}], "]"}]}]}], ";"}], "*)"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"wx", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{
        RowBox[{
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"a", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "3"}], "]"}], "]"}], "-", "0.2555", "+", 
            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", ",", "14", ",", "20"}], "}"}]}], "]"}]}], ";"}], "\n", 
   "      ", 
   RowBox[{
    RowBox[{"wy", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{
        RowBox[{
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"a", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "3"}], "]"}], "]"}], "-", "0.2555", "+", 
            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", ",", "14", ",", "20"}], "}"}]}], "]"}]}], ";"}], "\n", 
   "      ", 
   RowBox[{
    RowBox[{"wz", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{
        RowBox[{
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"a", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "3"}], "]"}], "]"}], "-", "0.2555", "+", 
            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", ",", "14", ",", "20"}], "}"}]}], "]"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"Plot8", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"Graphics3D", "[", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"Opacity", "[", "0.5", "]"}], ",", 
          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", ",", "7"}], "}"}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"wx", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{
        RowBox[{
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"a", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "3"}], "]"}], "]"}], "-", "0.1887", "+", 
            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", ",", "21", ",", "21"}], "}"}]}], "]"}]}], ";"}], "\n", 
   "      ", 
   RowBox[{
    RowBox[{"wy", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{
        RowBox[{
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"a", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "3"}], "]"}], "]"}], "-", "0.1887", "+", 
            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", ",", "21", ",", "21"}], "}"}]}], "]"}]}], ";"}], "\n", 
   "      ", 
   RowBox[{
    RowBox[{"wz", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{
        RowBox[{
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"a", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "3"}], "]"}], "]"}], "-", "0.1887", "+", 
            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", ",", "21", ",", "21"}], "}"}]}], "]"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"Plot9", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"Graphics3D", "[", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"Opacity", "[", "0.5", "]"}], ",", 
          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", ",", "1"}], "}"}]}], "]"}]}], ";"}], "\n", 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"Show", "[", 
    RowBox[{"Plot1", ",", "Plot7", ",", "Plot8", ",", "Plot9", ",", 
     RowBox[{"{", 
      RowBox[{"AspectRatio", "\[Rule]", "Automatic"}], "}"}]}], 
    "]"}]}]}]], "Input",
 CellChangeTimes->{{3.5681365021632137`*^9, 3.5681365171392403`*^9}, {
   3.568136560351316*^9, 3.5681365863097615`*^9}, {3.5681369135895367`*^9, 
   3.5681369265999594`*^9}, {3.568137821916732*^9, 3.5681378378443604`*^9}, {
   3.5682072492149415`*^9, 3.5682072529901485`*^9}, 3.568209802830227*^9, {
   3.5682216503836985`*^9, 3.5682217583046885`*^9}, {3.5682217912831464`*^9, 
   3.568221843465238*^9}, {3.5682219760654707`*^9, 3.5682221133457117`*^9}, {
   3.568222159615393*^9, 3.5682221693498106`*^9}, {3.5682222565227633`*^9, 
   3.5682222578487663`*^9}, 3.5682224736283445`*^9, 3.5682227227607822`*^9, 
   3.5682227558640404`*^9, {3.568222818810151*^9, 3.5682228486062036`*^9}, {
   3.5682230057764797`*^9, 3.5682230079916835`*^9}, {3.568223041952943*^9, 
   3.5682230548229656`*^9}, 3.5682232619601297`*^9, 3.5682233000709963`*^9, {
   3.568223536645412*^9, 3.5682235381118145`*^9}, {3.568223591151908*^9, 
   3.568223591526308*^9}, {3.5682236491528096`*^9, 3.568223657483224*^9}, {
   3.568223713222122*^9, 3.5682237532361927`*^9}, {3.568223800379475*^9, 
   3.568223827648323*^9}, {3.5682238887380304`*^9, 3.5682238893152313`*^9}, {
   3.5682240287170763`*^9, 3.5682240612431335`*^9}, 3.5682241088076167`*^9, {
   3.568224198772975*^9, 3.5682242396450467`*^9}, {3.568224282170721*^9, 
   3.5682243185031853`*^9}, {3.568224437890195*^9, 3.568224439590598*^9}, {
   3.568224572253231*^9, 3.5682245867144566`*^9}, {3.5682247081762695`*^9, 
   3.568224733557514*^9}, {3.56822480520844*^9, 3.5682248433505073`*^9}, {
   3.568224884206979*^9, 3.5682248959069996`*^9}, {3.568224982830352*^9, 
   3.5682250047639904`*^9}, {3.56822503870965*^9, 3.5682250803305235`*^9}, {
   3.5682251223413973`*^9, 3.568225141451431*^9}, {3.568225234177994*^9, 
   3.5682252929588966`*^9}, {3.5682253775890455`*^9, 3.568225482140429*^9}, {
   3.56822551694409*^9, 3.5682255297205124`*^9}, {3.568225573197789*^9, 
   3.5682255809666023`*^9}, 3.568225660698343*^9, 3.5682257301964645`*^9, {
   3.568225765078126*^9, 3.5682257828933573`*^9}, {3.568225927026843*^9, 
   3.568225927416844*^9}, {3.568267062420238*^9, 3.568267081358671*^9}, {
   3.568267188328059*^9, 3.5682672830358257`*^9}, {3.5682676367040462`*^9, 
   3.568267747979042*^9}, 3.5682678404716043`*^9, 3.568267882030078*^9, {
   3.5682679252577534`*^9, 3.568267996612279*^9}, {3.5682680488879704`*^9, 
   3.568268061133992*^9}, {3.568268100914062*^9, 3.5682681150944867`*^9}, {
   3.5682681672609787`*^9, 3.5682681709113846`*^9}, {3.568268272904364*^9, 
   3.5682683182536435`*^9}, {3.5682685202739983`*^9, 3.568268599194537*^9}, {
   3.5682686348873997`*^9, 3.568268637773405*^9}, {3.5682686944951043`*^9, 
   3.568268872273017*^9}, {3.5682689420207396`*^9, 3.568269016947671*^9}, {
   3.5682695040741262`*^9, 3.568269601153097*^9}, {3.5682696652380095`*^9, 
   3.5682696728352227`*^9}, {3.5682697489633565`*^9, 
   3.5682698308635006`*^9}, {3.568271156491429*^9, 3.568271160063835*^9}, {
   3.568271244927984*^9, 3.568271369010602*^9}, {3.5682716250070515`*^9, 
   3.5682717312276382`*^9}, {3.5682717843301315`*^9, 
   3.5682718576814604`*^9}, {3.568271940673606*^9, 3.5682720536490045`*^9}, {
   3.5682722259669075`*^9, 3.5682723223126764`*^9}, {3.568272370813162*^9, 
   3.568272411466833*^9}, {3.5682725893539457`*^9, 3.568272622238803*^9}, {
   3.5682726936713285`*^9, 3.5682726960113325`*^9}, {3.568272741079812*^9, 
   3.56827275143823*^9}, {3.5682729218373294`*^9, 3.568272991070251*^9}, {
   3.5682734285574193`*^9, 3.5682734302890224`*^9}, {3.5682734898811274`*^9, 
   3.5682735111595645`*^9}, {3.5682736115925407`*^9, 
   3.5682736313577757`*^9}, {3.5682738601945777`*^9, 3.568273987272401*^9}, {
   3.5682763944346285`*^9, 3.5682764313910933`*^9}, 3.5682765128232365`*^9, {
   3.5682766527398825`*^9, 3.568276660976697*^9}, {3.5682774496048822`*^9, 
   3.568277461663703*^9}}],

Cell[BoxData[
 Graphics3DBox[{
   SphereBox[{-2.093816408, -0.96496968464, -0.31766294059999994`}, 0.1887], 
   SphereBox[{2.9740507787999997`, -2.8338331096, 5.428804852799999}, 0.1887],
    SphereBox[{-0.7338939893599999, -3.3497283075999995`, 
    3.5349214555999997`}, 0.1887], 
   SphereBox[{-3.544181113, 1.70381475612, -3.1853221456}, 0.1887], 
   SphereBox[{0.17769093578, 0.22247744198, -0.23797319517999999`}, 0.2555], 
   SphereBox[{1.9210954515999998`, -0.36732493988, 1.66696509912}, 0.2555], 
   SphereBox[{0.7848976940999999, 2.0244634637999996`, -2.0410930526}, 
    0.2555], SphereBox[{4.2586865136, 0.8596741519199999, 1.72195612572}, 
    0.2555], SphereBox[{3.0970719413999994`, 
    3.2310535147999997`, -1.9577561359999998`}, 0.2555], 
   SphereBox[{4.8303286286, 2.6452384548, -0.08513404602599998}, 0.2555], 
   SphereBox[{1.29688115928, -2.279954419, 3.5948257697999995`}, 0.2555], 
   SphereBox[{-3.8767728889999997`, -0.185710932924, -1.947362643}, 0.2555], 
   SphereBox[{-6.291842717, -1.6678721676000001`, -2.2315774334}, 0.2555], 
   SphereBox[{-0.5526503687, 2.4848007174, -3.5122447436}, 0.1], 
   SphereBox[{5.6156987542, 0.40901229543999995`, 3.1783301593999997`}, 0.1], 
   SphereBox[{3.560243784, 4.633986097199999, -3.3656020059999996`}, 0.1], 
   SphereBox[{6.6372846298, 3.5927470712, -0.041141224745999995`}, 0.1], 
   SphereBox[{-7.498432768, -0.7529802219599999, -3.6278959747999995`}, 0.1], 
   SphereBox[{-7.2703428398, -1.73992741998, -0.42057741855999997`}, 0.1], 
   SphereBox[{-5.853426284999999, -3.5825425508, -2.8517855065999997`}, 0.1], 
   SphereBox[{2.5617125655999997`, -4.082564050399999, 6.6973779166}, 0.1], 
   SphereBox[{-0.9876480052521825, -0.3867120377970378, -0.2788559877420218}, 
    1.0600979904532752`], 
   SphereBox[{1.04939319369, -0.07242374895000003, 0.7144959519700002}, 
    1.0688976046225591`], 
   SphereBox[{0.48129431493999997`, 1.12347045289, -1.13953312389}, 
    1.0547571207979058`], 
   SphereBox[{3.0898909826, 0.24617460601999996`, 1.69446061242}, 
    1.0648108305837336`], 
   SphereBox[{1.9409848177499995`, 2.6277584892999997`, -1.9994245942999997`},
     1.0491986570390401`], 
   SphereBox[{4.5445075711000005`, 1.7524563033599998`, 0.818411039847}, 
    1.0464782521399887`], 
   SphereBox[{3.963700285, 2.9381459847999998`, -1.021445091013}, 
    1.0535140870236], 
   SphereBox[{2.1574663460697274`, -2.5641593025497444`, 4.535872646241579}, 
    1.0510023931152543`], 
   SphereBox[{0.25195304375062255`, -2.8304027609574742`, 3.564002218413227}, 
    0.9259475168654768], 
   SphereBox[{1.6089883054399998`, -1.3236396794400003`, 2.6308954344599997`},
     1.13773739655973], 
   SphereBox[{-2.961832190519656, -0.5855948079314673, -1.1110670856396943`}, 
    1.046964530990104], 
   SphereBox[{-3.705611869286717, 0.7866917752873251, -2.5844511885350654`}, 
    0.9195510356822784], 
   SphereBox[{-5.084307803, -0.9267915502620002, -2.0894700382}, 
    1.1684155739346305`], 
   {Opacity[0.5], 
    SphereBox[{-0.016226052857627127`, 
     2.3001822328857275`, -2.9222386262175655`}, 1.0600979904532752`]}, 
   {Opacity[0.5], 
    SphereBox[{5.071492565877527, 0.5897424051616867, 2.5942766498381515`}, 
     1.0688976046225591`]}, 
   {Opacity[0.5], 
    SphereBox[{3.3744783501865765`, 4.071308585722703, -2.8009539105654824`}, 
     1.0547571207979058`]}, 
   {Opacity[0.5], 
    SphereBox[{5.912566220015856, 3.212728444693234, -0.05878548805887894}, 
     1.0648108305837336`]}, 
   {Opacity[0.5], 
    SphereBox[{-7.018961293602498, -1.1165374962196812`, \
-3.0730307086569866`}, 1.0491986570390401`]}, 
   {Opacity[0.5], 
    SphereBox[{-6.881499186278548, -1.7112935680897472`, \
-1.1402460684552256`}, 1.0464782521399887`]}, 
   {Opacity[0.5], 
    SphereBox[{-6.02763814074903, -2.8217173256063903`, -2.6053358037688286`},
      1.0535140870236]}, 
   {Opacity[0.5], 
    SphereBox[{2.7171299379424174`, -3.6118958577149876`, 6.219230886107493}, 
     1.0600979904532752`]}},
  AspectRatio->Automatic,
  ImageSize->{802.3527740428565, 643.},
  ViewPoint->{1.397094377009007, -1.250935673201176, -2.816609174742311},
  ViewVertical->{
   0.12364663907566328`, -1.4853556447128842`, -0.29221139691535325`}]], \
"Output",
 CellChangeTimes->{
  3.568267862857644*^9, {3.5682679316225643`*^9, 3.5682679544922047`*^9}, 
   3.5682680043498926`*^9, 3.5682680757668176`*^9, 3.568268186324212*^9, 
   3.5682683262096577`*^9, 3.56826887966743*^9, 3.5682690294900928`*^9, 
   3.568269543152195*^9, 3.568269691742456*^9, 3.568271398557054*^9, 
   3.5682724298748655`*^9, 3.568272634703225*^9, 3.568273012941489*^9, 
   3.5682734359674325`*^9, 3.568273517742776*^9, 3.5682736374261866`*^9, 
   3.5682739957120156`*^9, 3.5682764649467525`*^9, 3.568276518954047*^9, 
   3.5682766682619095`*^9, {3.5682777008277235`*^9, 3.568277718596154*^9}}]
}, Open  ]]
},
WindowSize->{1013, 789},
WindowMargins->{{409, Automatic}, {Automatic, 5}},
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, 121, 2, 29, "Input"],
Cell[703, 26, 3788, 65, 258, "Output"]
}, Open  ]],
Cell[4506, 94, 400, 9, 29, "Input"],
Cell[CellGroupData[{
Cell[4931, 107, 557, 8, 29, "Input"],
Cell[5491, 117, 2108, 57, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[7636, 179, 724, 11, 29, "Input"],
Cell[8363, 192, 521, 8, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[8921, 205, 604, 11, 50, InheritFromParent],
Cell[9528, 218, 517, 8, 28, "Output"],
Cell[10048, 228, 95, 1, 28, "Output"]
}, Open  ]],
Cell[10158, 232, 711, 13, 50, "Input"],
Cell[10872, 247, 827, 22, 50, "Input"],
Cell[CellGroupData[{
Cell[11724, 273, 564, 9, 29, "Input"],
Cell[12291, 284, 941, 13, 28, "Output"]
}, Open  ]],
Cell[13247, 300, 616, 12, 29, "Input"],
Cell[13866, 314, 478, 10, 29, "Input"],
Cell[14347, 326, 3767, 79, 50, "Input"],
Cell[18117, 407, 3784, 74, 70, "Input"],
Cell[CellGroupData[{
Cell[21926, 485, 2500, 35, 29, "Input"],
Cell[24429, 522, 2801, 42, 488, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[27267, 569, 597, 9, 29, "Input"],
Cell[27867, 580, 16362, 390, 174, "Output"]
}, Open  ]],
Cell[44244, 973, 1204, 34, 70, "Input"],
Cell[CellGroupData[{
Cell[45473, 1011, 4738, 110, 110, "Input"],
Cell[50214, 1123, 6150, 106, 337, "Output"]
}, Open  ]],
Cell[56379, 1232, 1935, 49, 70, "Input"],
Cell[58317, 1283, 606, 12, 50, "Input"],
Cell[58926, 1297, 1656, 42, 90, "Input"],
Cell[60585, 1341, 511, 11, 50, "Input"],
Cell[61099, 1354, 1998, 47, 90, "Input"],
Cell[63100, 1403, 1675, 41, 50, "Input"],
Cell[64778, 1446, 1217, 33, 50, "Input"],
Cell[CellGroupData[{
Cell[66020, 1483, 1027, 27, 29, "Input"],
Cell[67050, 1512, 950, 13, 28, "Output"]
}, Open  ]],
Cell[68015, 1528, 1971, 35, 90, "Input"],
Cell[69989, 1565, 1638, 24, 29, "Input"],
Cell[71630, 1591, 1655, 24, 29, "Input"],
Cell[CellGroupData[{
Cell[73310, 1619, 2699, 53, 70, "Input"],
Cell[CellGroupData[{
Cell[76034, 1676, 1273, 20, 21, "Print"],
Cell[77310, 1698, 1273, 20, 21, "Print"],
Cell[78586, 1720, 1273, 20, 21, "Print"],
Cell[79862, 1742, 1273, 20, 21, "Print"],
Cell[81138, 1764, 1271, 20, 21, "Print"],
Cell[82412, 1786, 1273, 20, 21, "Print"],
Cell[83688, 1808, 1266, 20, 21, "Print"],
Cell[84957, 1830, 1273, 20, 21, "Print"],
Cell[86233, 1852, 1272, 20, 21, "Print"],
Cell[87508, 1874, 1268, 20, 21, "Print"],
Cell[88779, 1896, 1270, 20, 21, "Print"],
Cell[90052, 1918, 1274, 20, 21, "Print"],
Cell[91329, 1940, 1277, 20, 21, "Print"],
Cell[92609, 1962, 1275, 20, 21, "Print"],
Cell[93887, 1984, 1275, 20, 21, "Print"],
Cell[95165, 2006, 1275, 20, 21, "Print"],
Cell[96443, 2028, 1277, 20, 21, "Print"],
Cell[97723, 2050, 1277, 20, 21, "Print"],
Cell[99003, 2072, 1277, 20, 21, "Print"],
Cell[100283, 2094, 1277, 20, 21, "Print"],
Cell[101563, 2116, 1275, 20, 21, "Print"]
}, Open  ]]
}, Open  ]],
Cell[102865, 2140, 2648, 52, 90, "Input"],
Cell[CellGroupData[{
Cell[105538, 2196, 675, 11, 29, "Input"],
Cell[106216, 2209, 2801, 58, 67, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[109054, 2272, 25170, 693, 750, "Input"],
Cell[134227, 2967, 4837, 89, 592, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)