(* 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[ 114937, 2693] NotebookOptionsPosition[ 111342, 2582] NotebookOutlinePosition[ 111878, 2603] CellTagsIndexPosition[ 111835, 2600] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ RowBox[{ "Kimball", " ", "model", " ", "and", " ", "energy", " ", "derived", " ", "from", " ", "experimental", " ", "molecule", " ", "data", " ", "source"}], ",", " ", RowBox[{"ES", " ", "end", " ", "january", " ", "2013"}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"ChemicalData", "[", RowBox[{"\"\\"", ",", "\"\\""}], "]"}]}]], "Input", CellChangeTimes->{{3.568278170270548*^9, 3.568278174310955*^9}, { 3.5682782109554195`*^9, 3.5682782183342323`*^9}, 3.568284752749406*^9, 3.5682847957898817`*^9, {3.5682872948362713`*^9, 3.568287369232802*^9}}], Cell[BoxData[ Graphics3DBox[ {AbsoluteThickness[3], EdgeForm[None], Specularity[ GrayLevel[1], 100], GraphicsComplex3DBox[CompressedData[" 1:eJw1Uw1MjGEcv1Vi1kKMskQfWqQdktZa7l9XXSrdR+e4Ps51d2/dW1cdopwK axixfLXaRIQsSz4WWrP2X8VsIp3M1FoqreUjLGqkzfO8D8/27N37fPz+v4// 423IV3EOIpFoLZmOZBa21S4X35ZD3JBtSVNWKt76bt9msGqALg80a+Dy15tP E6bDoMXr6CflTh26OAQeX/pOjJl0bNHDZyXZ8VVDRfvPgDRJBqpHblzsBBlG l1vGZrRx7FySESi8jFdj9uOC8+4fJPj8aszDydsyEM7rN4GIjjgO/jQQAh9V GExHmRH3ms9+Ka6VwmrnI5XSSRXjecgIuyPmv712RooHxv03zKuRoW//r9Ls 4k1IqgzZojlUUcJbjThRLAnRGdQQeiJSPyoCWJNGmFanwEoip2VzFuM3YQC6 HOBnBArXWBMDAk4vjx5OGtcgezgKdatygF7zmtWz+7G5SOU4NSVjPDWmXs14 B+fDN2Jb6AEZxl6Yct7lPofVfW9iPhVy2LqVXFigQQFnIgNyx2YavtfFM/wS Hojacsvc9XiwTJxYf9IRiFprhD0fKK0leRl4kBLs43AqgAA/0UHzs/1ttT4K FPAf5TF/V2xEhsfBiy4yCrIgqjvzzVj2IubLLWC+L8xl++u0WEKWW7VR2EFg KupSmT86C+pHX/n039Pgu5fTySMPEhiPKQvj4ZPD8OLlcN1KgvFNB4HXgBmF c5eN8HBytreIj8D7l0hDhW6HYdoP+wwYPig/XeWvA+FbqmM5JEayfjtswjNU z2AIUvnzP3Eg1HNNBeFf+i+/JDnzgVcgldEFPFion24KHB6vvO8SGcTyW2Zi eht3YW8RTToJz30hxPt2wCl6gY8CWq7CkMXyV4uh43Ug52H3Z7mczoZGzxRT +EgK43PHhNV5bmUYr4adYSSBY0oU+nwFj57iAT9n79Ug7E/pWH42DrTLHvWI 2kKwpnPV3T1xCcwvmRkEv21K5ufvaIylhukVSNvc2mNG1yDinESJ1M6mwg3w v46gn/Sx4KujCq6QZ+t3TsVwF2ei8G660wEpfvsPyV9RsNyT "], { {RGBColor[0.65, 0.7, 0.7], SphereBox[11, 24.], SphereBox[12, 24.], SphereBox[13, 24.], SphereBox[14, 24.], SphereBox[15, 24.], SphereBox[16, 24.], SphereBox[17, 24.], SphereBox[18, 24.], SphereBox[19, 24.], SphereBox[20, 24.], SphereBox[21, 24.], SphereBox[22, 24.], SphereBox[23, 24.], SphereBox[24, 24.], SphereBox[25, 24.], SphereBox[26, 24.]}, {RGBColor[0.4, 0.4, 0.4], SphereBox[1, 34.], SphereBox[2, 34.], SphereBox[3, 34.], SphereBox[4, 34.], SphereBox[5, 34.], SphereBox[6, 34.], SphereBox[7, 34.], SphereBox[8, 34.], SphereBox[9, 34.], SphereBox[10, 34.]}, {RGBColor[0.65, 0.7, 0.7], CylinderBox[{27, 11}, 15.], CylinderBox[{28, 12}, 15.], CylinderBox[{29, 13}, 15.], CylinderBox[{30, 14}, 15.], CylinderBox[{31, 15}, 15.], CylinderBox[{32, 16}, 15.], CylinderBox[{33, 17}, 15.], CylinderBox[{34, 18}, 15.], CylinderBox[{35, 19}, 15.], CylinderBox[{36, 20}, 15.], CylinderBox[{37, 21}, 15.], CylinderBox[{38, 22}, 15.], CylinderBox[{39, 23}, 15.], CylinderBox[{40, 24}, 15.], CylinderBox[{41, 25}, 15.], CylinderBox[{42, 26}, 15.]}, {RGBColor[0.4, 0.4, 0.4], CylinderBox[{1, 5}, 15.], CylinderBox[{1, 6}, 15.], CylinderBox[{1, 7}, 15.], CylinderBox[{1, 27}, 15.], CylinderBox[{2, 5}, 15.], CylinderBox[{2, 8}, 15.], CylinderBox[{2, 9}, 15.], CylinderBox[{2, 28}, 15.], CylinderBox[{3, 6}, 15.], CylinderBox[{3, 8}, 15.], CylinderBox[{3, 10}, 15.], CylinderBox[{3, 29}, 15.], CylinderBox[{4, 7}, 15.], CylinderBox[{4, 9}, 15.], CylinderBox[{4, 10}, 15.], CylinderBox[{4, 30}, 15.], CylinderBox[{5, 31}, 15.], CylinderBox[{5, 32}, 15.], CylinderBox[{6, 33}, 15.], CylinderBox[{6, 34}, 15.], CylinderBox[{7, 35}, 15.], CylinderBox[{7, 36}, 15.], CylinderBox[{8, 37}, 15.], CylinderBox[{8, 38}, 15.], CylinderBox[{9, 39}, 15.], CylinderBox[{9, 40}, 15.], CylinderBox[{10, 41}, 15.], CylinderBox[{10, 42}, 15.]}}]}, Boxed->False, ImageSize->217.75316128581244`, Lighting->"Neutral", ViewPoint->{1.3, -2.4, -2.}]], "Output", CellChangeTimes->{3.568286549435762*^9, 3.5682870083729677`*^9}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"t", "=", RowBox[{"0.01889726", "*", RowBox[{"ChemicalData", "[", RowBox[{"\"\\"", ",", "\"\\""}], "]"}]}]}], ";", " ", RowBox[{"(*", " ", RowBox[{"convert", " ", "to", " ", "bohr"}], " ", "*)"}]}]], "Input", CellChangeTimes->{{3.567862912170579*^9, 3.567862949329844*^9}, { 3.5678630095303497`*^9, 3.5678630107471523`*^9}, {3.5679573013548217`*^9, 3.5679573132576427`*^9}, 3.568114367533684*^9, {3.5682782383334675`*^9, 3.5682782435438766`*^9}, 3.5682847420789876`*^9, 3.568284795821082*^9, { 3.568287382898426*^9, 3.5682873928356433`*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"ChemicalData", "\[NoBreak]", "[", RowBox[{"\"\\"", ",", "\"\\""}], "]"}]], "Input", CellChangeTimes->{{3.5678630319943895`*^9, 3.5678630956269007`*^9}, { 3.5678632280087337`*^9, 3.567863228336334*^9}, {3.5678632621727934`*^9, 3.5678633509213495`*^9}, {3.5678644293824434`*^9, 3.567864451908883*^9}, { 3.567864923404111*^9, 3.56786494544695*^9}, {3.5678650143210707`*^9, 3.5678650845367947`*^9}, 3.567869231084096*^9, 3.5678693210962543`*^9, { 3.567869354027912*^9, 3.5678693684423375`*^9}, {3.5682782542610955`*^9, 3.5682782564138994`*^9}, 3.568284752702606*^9, 3.568284795821082*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"\<\"C\"\>", ",", "\<\"C\"\>"}], "}"}], ",", "1"}], "}"}], ",", "12"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"\<\"C\"\>", ",", "\<\"H\"\>"}], "}"}], ",", "1"}], "}"}], ",", "16"}], "}"}]}], "}"}]], "Output", CellChangeTimes->{3.5682865494981623`*^9, 3.5682870084197683`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"u", "=", RowBox[{"ChemicalData", "\[NoBreak]", "[", RowBox[{"\"\\"", ",", "\"\\""}], "]"}]}]], "Input", CellChangeTimes->{{3.567869377256353*^9, 3.5678694101568108`*^9}, { 3.5678695244894114`*^9, 3.56786953511303*^9}, 3.56786988243204*^9, 3.5678699252853155`*^9, {3.567870069179968*^9, 3.5678700766835814`*^9}, { 3.5678703216664114`*^9, 3.567870325738019*^9}, 3.5679271251072817`*^9, { 3.567942840466694*^9, 3.567942854381918*^9}, {3.567942932272855*^9, 3.567942950509287*^9}, {3.567943087867528*^9, 3.567943093733138*^9}, 3.567943185211699*^9, 3.5679432253505697`*^9, 3.5679433277647495`*^9, { 3.5679434430645523`*^9, 3.5679434440941534`*^9}, {3.568278260766307*^9, 3.5682782624355097`*^9}, 3.568284752702606*^9, 3.568284795852282*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{"\<\"C\"\>", ",", "\<\"C\"\>", ",", "\<\"C\"\>", ",", "\<\"C\"\>", ",", "\<\"C\"\>", ",", "\<\"C\"\>", ",", "\<\"C\"\>", ",", "\<\"C\"\>", ",", "\<\"C\"\>", ",", "\<\"C\"\>", ",", "\<\"H\"\>", ",", "\<\"H\"\>", ",", "\<\"H\"\>", ",", "\<\"H\"\>", ",", "\<\"H\"\>", ",", "\<\"H\"\>", ",", "\<\"H\"\>", ",", "\<\"H\"\>", ",", "\<\"H\"\>", ",", "\<\"H\"\>", ",", "\<\"H\"\>", ",", "\<\"H\"\>", ",", "\<\"H\"\>", ",", "\<\"H\"\>", ",", "\<\"H\"\>", ",", "\<\"H\"\>"}], "}"}]], "Output", CellChangeTimes->{3.568286549513762*^9, 3.568287008435368*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"\[IndentingNewLine]", RowBox[{"n", "=", RowBox[{"Length", "[", "u", "]"}]}]}]], "Input", CellChangeTimes->{{3.568277663590458*^9, 3.568277712652544*^9}, 3.5682785838808746`*^9}], Cell[BoxData["26"], "Output", CellChangeTimes->{3.5682865495293617`*^9, 3.568287008435368*^9}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"u1", "=", RowBox[{"{", RowBox[{ "0.2555", ",", "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", ",", "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}, 3.5682795170433135`*^9, {3.5682795636873956`*^9, 3.5682797843809834`*^9}, { 3.568279837187076*^9, 3.5682798459854913`*^9}, {3.568279877201146*^9, 3.568279880711152*^9}, {3.5682799570264864`*^9, 3.5682799599124913`*^9}, { 3.5682849151924915`*^9, 3.568284959231369*^9}, {3.5682849938634295`*^9, 3.568284994877432*^9}}], Cell[BoxData[ RowBox[{"(*", RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "<", "15"}], ",", RowBox[{"i", "++"}], ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"[", "i", "]"}], "]"}], "=", "\"\\""}], ",", 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", "[", RowBox[{"\"\\"", ",", "\"\\""}], "]"}]], "Input", CellChangeTimes->{{3.568284795883482*^9, 3.568284795883482*^9}}], Cell[BoxData["\<\"C10H16\"\>"], "Output", CellChangeTimes->{3.568286549544962*^9, 3.568287008450968*^9}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"v", "=", RowBox[{"0.01889726", "*", RowBox[{"ChemicalData", "\[NoBreak]", "[", RowBox[{"\"\\"", ",", "\"\\""}], "]"}]}]}], ";"}]], "Input", CellChangeTimes->{{3.567869377256353*^9, 3.567869394073182*^9}, { 3.5678694248832364`*^9, 3.5678694526200852`*^9}, {3.5678694891241493`*^9, 3.5678694955045605`*^9}, {3.567869596046737*^9, 3.567869602645549*^9}, 3.56786971175214*^9, {3.5678697908910794`*^9, 3.567869798800293*^9}, { 3.5679434548425727`*^9, 3.5679434621121855`*^9}, 3.568011875875622*^9, 3.56811435956207*^9, {3.5682782821071444`*^9, 3.5682782839479475`*^9}, 3.568284752733806*^9, 3.568284795899082*^9}], 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", ",", "n"}], "}"}]}], "]"}]}], ";"}]}]], "Input", CellChangeTimes->{{3.5679440895764875`*^9, 3.5679441000441055`*^9}, { 3.5679441319773617`*^9, 3.5679441523041973`*^9}, {3.567944240210352*^9, 3.5679442813476243`*^9}, {3.567944323966899*^9, 3.567944345931738*^9}, { 3.567944397567828*^9, 3.567944456941533*^9}, {3.567944491557994*^9, 3.5679445468444905`*^9}, {3.5679446505690727`*^9, 3.567944683875131*^9}, { 3.567944736244423*^9, 3.5679448177233667`*^9}, {3.567945549707852*^9, 3.5679455499730525`*^9}, {3.5679456039959474`*^9, 3.567945606195551*^9}, { 3.567945704460124*^9, 3.5679457164877453`*^9}, {3.5679457857362666`*^9, 3.5679458050023003`*^9}, {3.567946353670864*^9, 3.5679464127169676`*^9}, { 3.5679466013368993`*^9, 3.5679466311173515`*^9}, {3.5679467590999765`*^9, 3.567946915739851*^9}, {3.5679469466123056`*^9, 3.567946946815106*^9}, { 3.5679471385082426`*^9, 3.5679471497870626`*^9}, {3.5679483907224417`*^9, 3.5679484036080647`*^9}, {3.567948441937332*^9, 3.5679484421869326`*^9}, { 3.567948496740228*^9, 3.56794849772303*^9}, {3.567948847912445*^9, 3.5679488949933276`*^9}, {3.568011728877214*^9, 3.5680117308213253`*^9}, { 3.568011766325356*^9, 3.568011773436763*^9}, {3.5681173314452896`*^9, 3.5681173370612993`*^9}, {3.568117919800323*^9, 3.5681179678952074`*^9}, { 3.56811809459863*^9, 3.568118105409449*^9}, {3.5681181818963833`*^9, 3.5681182389144835`*^9}, {3.568118274888147*^9, 3.5681182883977704`*^9}, { 3.568118518513775*^9, 3.56811852739019*^9}, {3.5681188250543127`*^9, 3.568118862634779*^9}, {3.568118972396572*^9, 3.5681189862493963`*^9}, { 3.568119020694257*^9, 3.568119047136303*^9}, {3.5681190903015785`*^9, 3.568119136930061*^9}, {3.5681191671005135`*^9, 3.568119202808976*^9}, { 3.5681194788046613`*^9, 3.568119487603077*^9}, {3.5681195353235607`*^9, 3.5681195410331707`*^9}, 3.56811959201406*^9, {3.5681196426985493`*^9, 3.568119659702579*^9}, {3.568119726813897*^9, 3.56811974027672*^9}, { 3.5682659184078283`*^9, 3.5682659267694435`*^9}, {3.568265972477524*^9, 3.568266002725977*^9}, 3.568277802649102*^9}], 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", ",", "n"}], "}"}]}], "]"}]}], ";"}], "\[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}, 3.568277806081108*^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[{1.19719811278, -2.0715176412, 1.69614246856}, 0.2555], SphereBox[{1.7163625367599997`, 0.49068625316, -2.2937494187999996`}, 0.2555], SphereBox[{-0.30091996824, 2.4556989369999997`, 1.5222687793}, 0.2555], SphereBox[{-2.6074439347999996`, -0.87883597356, -0.92035325378}, 0.2555], SphereBox[{2.9186818069999996`, -1.5428667927, -0.6026525186599999}, 0.2555], SphereBox[{0.89495533634, 0.38508836428, 3.2172585149999997`}, 0.2555], SphereBox[{-1.42978558886, -2.9498622859999997`, 0.74164186596}, 0.2555], SphereBox[{1.43243120526, 2.9477835874, -0.7394497838}, 0.2555], SphereBox[{-0.89230971994, -0.38712926836, -3.2151798163999996`}, 0.2555], SphereBox[{-2.9160361906, 1.5408258886200001`, 0.6048446008199999}, 0.2555], SphereBox[{2.0104794914, -3.5094101546, 2.8619900269999996`}, 0.1], SphereBox[{2.8937374237999998`, 0.81681516624, -3.9049298063999998`}, 0.1], SphereBox[{-0.52574067046, 4.1713811724, 2.5683266066}, 0.1], SphereBox[{-4.4259272646, -1.4875734099399998`, -1.56095147052}, 0.1], SphereBox[{4.716756095999999, -0.88716966522, 0.04983018489399999}, 0.1], SphereBox[{3.1564093378, -3.2452264597999996`, -1.6674375306199998`}, 0.1], SphereBox[{2.7110009196, 1.00671373198, 3.8533402866}, 0.1], SphereBox[{-0.30290418053999996`, 0.04616033700199999, 4.810675478199999}, 0.1], SphereBox[{-2.6350339344, -3.3039969384, 2.3260637334}, 0.1], SphereBox[{-1.21964805766, -4.6436236998, -0.34236165942}, 0.1], SphereBox[{3.2486279666, 3.546070839, -0.08228622894399999}, 0.1], SphereBox[{0.62986457306, 4.39361295, -1.902954082}, 0.1], SphereBox[{-0.6935483392599999, -2.0889031204, -4.288733156999999}, 0.1], SphereBox[{-1.7154743655399998`, 1.04730504646, -4.3783061694}, 0.1], SphereBox[{-4.1024061733999995`, 1.18260942806, 2.2026646255999998`}, 0.1], SphereBox[{-3.7592319318, 2.9657359843999997`, -0.55580621112}, 0.1]}, AspectRatio->Automatic]], "Output", CellChangeTimes->{3.5682865495917625`*^9, 3.568287008466568*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"ChemicalData", "\[NoBreak]", "[", RowBox[{"\"\\"", ",", "\"\\""}], "]"}]], "Input", CellChangeTimes->{{3.5678630319943895`*^9, 3.5678630956269007`*^9}, { 3.5678632280087337`*^9, 3.567863228336334*^9}, {3.5678632621727934`*^9, 3.5678633509213495`*^9}, {3.5678644293824434`*^9, 3.567864451908883*^9}, { 3.567864923404111*^9, 3.56786494544695*^9}, {3.5678650143210707`*^9, 3.5678650845367947`*^9}, {3.567869231084096*^9, 3.5678692412397137`*^9}, 3.5678693168686466`*^9, {3.5678704450470285`*^9, 3.5678704621290584`*^9}, 3.5682783281584253`*^9, 3.568284752733806*^9, 3.568284795914682*^9}], Cell[BoxData[ GraphicsBox[{{ {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[ StyleBox["\<\"C\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {367.3, 118.00999999999999`}], LineBox[{{377.3934230875848, 112.18232161362599`}, {409.11, 93.87}}], {GrayLevel[0], InsetBox[ StyleBox["\<\"C\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {450.91999999999996`, 69.73}], LineBox[{{409.11, 93.87}, {440.8265769124152, 75.55767838637401}}]}}, {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[ StyleBox["\<\"C\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {367.3, 118.00999999999999`}], LineBox[{{357.20657691241524`, 112.18232161362599`}, {325.49, 93.87}}], {GrayLevel[0], InsetBox[ StyleBox["\<\"C\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {283.68, 69.73}], LineBox[{{325.49, 93.87}, {293.7734230875848, 75.55767838637401}}]}}, {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[ StyleBox["\<\"C\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {367.3, 118.00999999999999`}], LineBox[{{359.88507955790544`, 109.01789347052602`}, {332.745, 76.10499999999999}}], {GrayLevel[0], InsetBox[ StyleBox["\<\"C\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {298.19, 34.2}], LineBox[{{332.745, 76.10499999999999}, {305.60492044209457`, 43.19210652947397}}]}}, {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[ StyleBox["\<\"C\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {367.3, 118.00999999999999`}], LineBox[{{367.3, 129.66500000000002`}, {367.3, 149.01}}], {GrayLevel[0], InsetBox[ StyleBox["\<\"H\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {367.3, 180.01}], LineBox[{{367.3, 149.01}, {367.3, 168.35499999999996`}}]}}, {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[ StyleBox["\<\"C\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {450.91999999999996`, -26.82}], LineBox[{{450.91999999999996`, -15.164999999999985`}, { 450.91999999999996`, 21.455000000000002`}}], {GrayLevel[0], InsetBox[ StyleBox["\<\"C\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {450.91999999999996`, 69.73}], LineBox[{{450.91999999999996`, 21.455000000000002`}, { 450.91999999999996`, 58.07499999999999}}]}}, {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[ StyleBox["\<\"C\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {450.91999999999996`, -26.82}], LineBox[{{439.9613606508754, -30.788280309632768`}, { 409.11, -41.959999999999994`}}], {GrayLevel[0], InsetBox[ StyleBox["\<\"C\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {367.3, -57.099999999999994`}], LineBox[{{409.11, -41.959999999999994`}, { 378.25863934912456`, -53.13171969036723}}]}}, {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[ StyleBox["\<\"C\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {450.91999999999996`, -26.82}], LineBox[{{444.35539075492085`, -36.45041694109981}, { 419.84, -72.41499999999999}}], {GrayLevel[0], InsetBox[ StyleBox["\<\"C\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {388.76, -118.00999999999999`}], LineBox[{{419.84, -72.41499999999999}, { 395.3246092450791, -108.37958305890018`}}]}}, {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[ StyleBox["\<\"C\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {450.91999999999996`, -26.82}], LineBox[{{462.05736625178884`, -30.254835945647912`}, { 480.53999999999996`, -35.955}}], {GrayLevel[0], InsetBox[ StyleBox["\<\"H\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {510.16, -45.09}], LineBox[{{480.53999999999996`, -35.955}, { 499.02263374821115`, -41.65516405435209}}]}}, {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[ StyleBox["\<\"C\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {283.68, -26.82}], LineBox[{{283.68, -15.164999999999985`}, {283.68, 21.455000000000002`}}], {GrayLevel[0], InsetBox[ StyleBox["\<\"C\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {283.68, 69.73}], LineBox[{{283.68, 21.455000000000002`}, {283.68, 58.07499999999999}}]}}, {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[ StyleBox["\<\"C\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {283.68, -26.82}], LineBox[{{294.63863934912456`, -30.788280309632764`}, { 325.49, -41.959999999999994`}}], {GrayLevel[0], InsetBox[ StyleBox["\<\"C\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {367.3, -57.099999999999994`}], LineBox[{{325.49, -41.959999999999994`}, { 356.34136065087546`, -53.13171969036723}}]}}, {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[ StyleBox["\<\"C\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {283.68, -26.82}], LineBox[{{273.5847676744785, -32.64454369832911}, { 241.84, -50.959999999999994`}}], {GrayLevel[0], InsetBox[ StyleBox["\<\"C\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {200., -75.1}], LineBox[{{241.84, -50.959999999999994`}, { 210.09523232552155`, -69.27545630167089}}]}}, {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[ StyleBox["\<\"C\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {283.68, -26.82}], LineBox[{{284.7026467305868, -38.43004817666234}, {286.4, -57.7}}], {GrayLevel[0], InsetBox[ StyleBox["\<\"H\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {289.12, -88.58}], LineBox[{{286.4, -57.7}, {288.09735326941313`, -76.96995182333767}}]}}, {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[ StyleBox["\<\"C\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {298.19, -70.33}], LineBox[{{298.19, -58.67499999999998}, {298.19, -18.064999999999998`}}], {GrayLevel[0], InsetBox[ StyleBox["\<\"C\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {298.19, 34.2}], LineBox[{{298.19, -18.064999999999998`}, {298.19, 22.544999999999987`}}]}}, {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[ StyleBox["\<\"C\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {298.19, -70.33}], LineBox[{{308.5031784173135, -75.75930713191462}, { 343.475, -94.16999999999999}}], {GrayLevel[0], InsetBox[ StyleBox["\<\"C\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {388.76, -118.00999999999999`}], LineBox[{{343.475, -94.16999999999999}, { 378.4468215826865, -112.58069286808535`}}]}}, {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[ StyleBox["\<\"C\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {298.19, -70.33}], LineBox[{{286.5487282970084, -70.89552465651563}, {249.095, -72.715}}], {GrayLevel[0], InsetBox[ StyleBox["\<\"C\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {200., -75.1}], LineBox[{{249.095, -72.715}, { 211.64127170299162`, -74.53447534348436}}]}}, {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[ StyleBox["\<\"C\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {298.19, -70.33}], LineBox[{{295.66904086719666`, -81.70909443017044}, { 291.485, -100.595}}], {GrayLevel[0], InsetBox[ StyleBox["\<\"H\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {284.78, -130.85999999999999`}], LineBox[{{291.485, -100.595}, { 287.3009591328033, -119.48090556982955`}}]}}, {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[ StyleBox["\<\"C\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {450.91999999999996`, 69.73}], LineBox[{{454.9054303831555, 80.68241387370936}, {461.52, 98.86000000000001}}], {GrayLevel[0], InsetBox[ StyleBox["\<\"H\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {472.11999999999995`, 127.99000000000001`}], LineBox[{{461.52, 98.86000000000001}, {468.1345696168445, 117.03758612629063`}}]}}, {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[ StyleBox["\<\"C\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {450.91999999999996`, 69.73}], LineBox[{{462.3981440712103, 67.70732017349782}, {481.44999999999993`, 64.35}}], {GrayLevel[0], InsetBox[ StyleBox["\<\"H\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {511.97999999999996`, 58.97}], LineBox[{{481.44999999999993`, 64.35}, {500.5018559287896, 60.99267982650218}}]}}, {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[ StyleBox["\<\"C\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {283.68, 69.73}], LineBox[{{279.69290965322637`, 80.68180969368373}, {273.075, 98.86000000000001}}], {GrayLevel[0], InsetBox[ StyleBox["\<\"H\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {262.46999999999997`, 127.99000000000001`}], LineBox[{{273.075, 98.86000000000001}, {266.4570903467736, 117.03819030631628`}}]}}, {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[ StyleBox["\<\"C\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {283.68, 69.73}], LineBox[{{272.20185592878966`, 67.70732017349782}, {253.15, 64.35}}], {GrayLevel[0], InsetBox[ StyleBox["\<\"H\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {222.62, 58.97}], LineBox[{{253.15, 64.35}, {234.09814407121036`, 60.99267982650218}}]}}, {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[ StyleBox["\<\"C\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {298.19, 34.2}], LineBox[{{291.7301764142509, 43.9010156293555}, {281.01, 60.}}], {GrayLevel[0], InsetBox[ StyleBox["\<\"H\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {263.83, 85.8}], LineBox[{{281.01, 60.}, {270.2898235857491, 76.0989843706445}}]}}, {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[ StyleBox["\<\"C\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {298.19, 34.2}], LineBox[{{287.0342217142418, 30.82544286772229}, {268.52, 25.225}}], {GrayLevel[0], InsetBox[ StyleBox["\<\"H\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {238.85, 16.25}], LineBox[{{268.52, 25.225}, {250.0057782857582, 19.624557132277715`}}]}}, {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[ StyleBox["\<\"C\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {367.3, -57.099999999999994`}], LineBox[{{360.3446533178201, -66.4521215523906}, {348.8, -81.975}}], {GrayLevel[0], InsetBox[ StyleBox["\<\"H\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {330.3, -106.85}], LineBox[{{348.8, -81.975}, {337.25534668217995`, -97.4978784476094}}]}}, {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[ StyleBox["\<\"C\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {367.3, -57.099999999999994`}], LineBox[{{374.25534668217995`, -66.4521215523906}, {385.8, -81.975}}], {GrayLevel[0], InsetBox[ StyleBox["\<\"H\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {404.3, -106.85}], LineBox[{{385.8, -81.975}, {397.3446533178201, -97.4978784476094}}]}}, {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[ StyleBox["\<\"C\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {388.76, -118.00999999999999`}], LineBox[{{382.9439276851723, -128.1101152383869}, { 373.28999999999996`, -144.875}}], {GrayLevel[0], InsetBox[ StyleBox["\<\"H\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {357.82, -171.74}], LineBox[{{373.28999999999996`, -144.875}, { 363.63607231482763`, -161.6398847616131}}]}}, {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[ StyleBox["\<\"C\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {388.76, -118.00999999999999`}], LineBox[{{398.631467454594, -124.20622104938903`}, {415.015, -134.49}}], {GrayLevel[0], InsetBox[ StyleBox["\<\"H\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {441.27, -150.97}], LineBox[{{415.015, -134.49}, { 431.398532545406, -144.77377895061096`}}]}}, {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[ StyleBox["\<\"C\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {200., -75.1}], LineBox[{{190.9537013178601, -67.75129228003222}, { 175.94, -55.55499999999999}}], {GrayLevel[0], InsetBox[ StyleBox["\<\"H\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {151.88, -36.01}], LineBox[{{175.94, -55.55499999999999}, { 160.9262986821399, -43.35870771996777}}]}}, {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[ StyleBox["\<\"C\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {200., -75.1}], LineBox[{{196.36798649615318`, -86.1746333080547}, {190.34, -104.555}}], {GrayLevel[0], InsetBox[ StyleBox["\<\"H\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.04071196156203995], FontWeight->Plain], {180.68, -134.01000000000002`}], LineBox[{{190.34, -104.555}, { 184.31201350384683`, -122.93536669194532`}}]}}}, {}, {}}, ImageSize->{{0., 245.6280566280563}, {0., 243.85135135135104`}}, PlotRange->{{141.077, 522.783}, {-187.28000000000003`, 191.66500000000002`}}]], "Output", CellChangeTimes->{3.5682865496853623`*^9, 3.568287008544568*^9}] }, Open ]], Cell[BoxData[{ RowBox[{ RowBox[{"x", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"y", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"z", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{"i", ",", "3"}], "]"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}]}], "Input", CellChangeTimes->{{3.5678630319943895`*^9, 3.5678630956269007`*^9}, { 3.5678632280087337`*^9, 3.567863228336334*^9}, {3.5678632621727934`*^9, 3.5678633509213495`*^9}, {3.5678644293824434`*^9, 3.567864451908883*^9}, { 3.567864923404111*^9, 3.56786494544695*^9}, {3.5678650143210707`*^9, 3.5678650845367947`*^9}, {3.567869231084096*^9, 3.567869278118179*^9}, { 3.5682778102931156`*^9, 3.5682778156907253`*^9}}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"(*", " ", RowBox[{"Distance", " ", "matrix", " ", "Bohr"}], " ", "*)"}], "Text"], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Dij", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{"i", "\[GreaterEqual]", "j"}], ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"sq", "=", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"x", "[", RowBox[{"[", "i", "]"}], "]"}], "-", RowBox[{"x", "[", RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"y", "[", RowBox[{"[", "i", "]"}], "]"}], "-", RowBox[{"y", "[", RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"z", "[", RowBox[{"[", "i", "]"}], "]"}], "-", RowBox[{"z", "[", RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}]}], "]"}]}], ";", " ", RowBox[{"sq", "<", "3.0"}]}], ",", "sq", ",", "0"}], "]"}], ",", "\"\< \>\""}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{"Dij", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"sq", "=", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"x", "[", RowBox[{"[", "i", "]"}], "]"}], "-", RowBox[{"x", "[", RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"y", "[", RowBox[{"[", "i", "]"}], "]"}], "-", RowBox[{"y", "[", RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"z", "[", RowBox[{"[", "i", "]"}], "]"}], "-", RowBox[{"z", "[", RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}]}], "]"}]}], ";", " ", RowBox[{"sq", "<", "3.0"}]}], ",", "sq", ",", "\"\< \>\""}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}], "*)"}], "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"MatrixForm", "[", "Dij", "]"}], "*)"}]}]}]], "Input", CellChangeTimes->{{3.567863535157673*^9, 3.5678635412104836`*^9}, { 3.5678638182825704`*^9, 3.567863827283786*^9}, {3.5679503910983553`*^9, 3.5679504266352177`*^9}, {3.5679504653700857`*^9, 3.5679505048069553`*^9}, 3.567950534915008*^9, {3.567951107623214*^9, 3.567951136046464*^9}, 3.567951302545556*^9, {3.5679562986474605`*^9, 3.567956305355472*^9}, { 3.5679563590975666`*^9, 3.567956370688387*^9}, {3.567956414930065*^9, 3.567956430467692*^9}, {3.567956484740187*^9, 3.567956500667815*^9}, { 3.5679565456270943`*^9, 3.5679565868267665`*^9}, {3.5679569267045636`*^9, 3.5679569448941956`*^9}, {3.567956976562251*^9, 3.567957000711094*^9}, { 3.5679570348595533`*^9, 3.5679570632672033`*^9}, {3.5679573263928657`*^9, 3.5679573286548696`*^9}, {3.5679573611029263`*^9, 3.5679573637393312`*^9}, {3.5679592820903006`*^9, 3.567959317127962*^9}, { 3.567959397655304*^9, 3.5679594004945087`*^9}, {3.5679594704762316`*^9, 3.5679594724574347`*^9}, {3.567959655898157*^9, 3.567959673073787*^9}, 3.5681143121223865`*^9, 3.568121866388855*^9, 3.568220668221574*^9, { 3.568220718968463*^9, 3.568220723742071*^9}, {3.568220862051914*^9, 3.5682208641111183`*^9}, {3.568220906761593*^9, 3.5682209083527956`*^9}, { 3.5682209451844606`*^9, 3.5682209651680956`*^9}, {3.568221147782016*^9, 3.5682211695908546`*^9}, {3.568221237279373*^9, 3.568221253612602*^9}, { 3.568221305420293*^9, 3.568221306247094*^9}, {3.5682215161142635`*^9, 3.568221546206716*^9}, 3.56822158247678*^9, {3.568277819185131*^9, 3.5682778284359474`*^9}, {3.568286391407484*^9, 3.568286398053096*^9}}], Cell[BoxData[ RowBox[{"(*", RowBox[{ RowBox[{"Eij", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{"i", "\[GreaterEqual]", "j"}], ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"sq", "=", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"x", "[", RowBox[{"[", "i", "]"}], "]"}], "-", RowBox[{"x", "[", RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"y", "[", RowBox[{"[", "i", "]"}], "]"}], "-", RowBox[{"y", "[", RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"z", "[", RowBox[{"[", "i", "]"}], "]"}], "-", RowBox[{"z", "[", RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}]}], "]"}]}], ";", " ", RowBox[{"sq", "<", "3.0"}]}], ",", RowBox[{"0.529177", "*", "sq"}], ",", "0"}], "]"}], ",", "\"\< \>\""}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"MatrixForm", "[", "Eij", "]"}], ";"}], "*)"}]], "Input", CellChangeTimes->{{3.5679574512554846`*^9, 3.5679574566530943`*^9}, { 3.5679575983637433`*^9, 3.5679576252893906`*^9}, {3.5679577071427345`*^9, 3.5679577424611964`*^9}, 3.567958973147358*^9, {3.5679591782097178`*^9, 3.56795917902092*^9}, {3.567959225290601*^9, 3.567959254322252*^9}, { 3.5682778318055534`*^9, 3.568277834176758*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"u2", "=", RowBox[{"{", RowBox[{ "0.2555", ",", "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.563", ",", "0.563", ",", "0.563", ",", "0.563", ",", "0.563", ",", "0.563", ",", "0.563", ",", "0.563", ",", "0.563"}], "}"}]}], ";"}]], "Input", CellChangeTimes->{{3.5679585418846006`*^9, 3.5679585738334565`*^9}, { 3.567959075187137*^9, 3.5679590758579383`*^9}, {3.5681145227851562`*^9, 3.5681146476945763`*^9}, {3.568280038115429*^9, 3.5682801254131823`*^9}, { 3.5682851622189255`*^9, 3.5682851860713673`*^9}, 3.5682852194866257`*^9}], Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"Fij", "=", RowBox[{"Dij", "-", "u2"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{"MatrixForm", "[", "Fij", "]"}], ";"}], "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Gij", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{"i", "\[GreaterEqual]", "j"}], ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"Fij", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], ">", "0"}], ",", RowBox[{"Fij", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], ",", "0"}], "]"}], ",", "\"\< \>\""}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"(*", RowBox[{ RowBox[{"MatrixForm", "[", "Gij", "]"}], ";"}], "*)"}]}]}], "Input", CellChangeTimes->{{3.567958596890297*^9, 3.567958603270708*^9}, { 3.5679587604565845`*^9, 3.5679587707682023`*^9}, {3.5679588272871017`*^9, 3.567958828254304*^9}, {3.567959580066424*^9, 3.5679595814704266`*^9}, { 3.567959624994503*^9, 3.5679596300177116`*^9}, {3.5679597079086485`*^9, 3.567959902877791*^9}, {3.5679599961659546`*^9, 3.5679600446196404`*^9}, { 3.5679609611992497`*^9, 3.567960974006872*^9}, 3.567961975154231*^9, { 3.5679622637079377`*^9, 3.5679622945179915`*^9}, 3.5681142413606625`*^9, 3.568121878697276*^9, {3.5681219332193723`*^9, 3.56812197189184*^9}, { 3.5682778370003624`*^9, 3.5682778397927675`*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"Hij", "=", RowBox[{"Transpose", "[", "Gij", "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"(*", RowBox[{ RowBox[{"MatrixForm", "[", "Hij", "]"}], ";"}], "*)"}]}]}], "Input", CellChangeTimes->{{3.5679601758938704`*^9, 3.5679601965327067`*^9}, { 3.567960427147912*^9, 3.567960432217921*^9}, {3.5679605764089737`*^9, 3.5679605968294096`*^9}, {3.5679607389456596`*^9, 3.5679607643581038`*^9}, 3.5681142501278777`*^9, 3.5681218877140923`*^9}], Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"Kij", "=", RowBox[{"Hij", "-", "u2"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{"MatrixForm", "[", "Kij", "]"}], ";"}], "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Lij", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{"i", "\[LessEqual]", "j"}], ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"Kij", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], ">", "0"}], ",", RowBox[{"Kij", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], ",", "0"}], "]"}], ",", "\"\< \>\""}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"(*", RowBox[{ RowBox[{"MatrixForm", "[", "Lij", "]"}], ";"}], "*)"}]}]}], "Input", CellChangeTimes->{{3.5679602339571724`*^9, 3.567960257419614*^9}, { 3.567960343406965*^9, 3.567960353500182*^9}, {3.5679604155258913`*^9, 3.567960421422702*^9}, {3.567960854323462*^9, 3.567960874197897*^9}, { 3.5679611596939983`*^9, 3.567961170270817*^9}, {3.5679614000124207`*^9, 3.567961407016833*^9}, {3.5679615066074076`*^9, 3.567961509649413*^9}, 3.5679616124379935`*^9, 3.5679616766477065`*^9, {3.5679617424486217`*^9, 3.5679617456310277`*^9}, 3.567961801276325*^9, {3.5679618582788253`*^9, 3.5679618970760937`*^9}, {3.5679620279135237`*^9, 3.567962030113127*^9}, 3.567962107894864*^9, {3.567962149468937*^9, 3.5679621572533503`*^9}, 3.5679623389156694`*^9, {3.567962645815009*^9, 3.567962648888214*^9}, { 3.567962691647889*^9, 3.567962693020692*^9}, {3.568048161486014*^9, 3.568048162827616*^9}, 3.5680482096900983`*^9, 3.5681142588638926`*^9, 3.568121899632513*^9, 3.5681219958534822`*^9, {3.568277844675576*^9, 3.5682778463603787`*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"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", "<", "11"}], ",", RowBox[{"0.5", "*", RowBox[{"Lij", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}]}], ",", RowBox[{"Lij", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}]}], "]"}], ",", "0"}], "]"}], ",", "\"\< \>\""}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";", " ", RowBox[{"Qij", "=", "Lij"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"(*", RowBox[{ RowBox[{"MatrixForm", "[", "Lij", "]"}], ";"}], "*)"}]}]}], "Input", CellChangeTimes->{{3.5679627617388124`*^9, 3.567962837461345*^9}, { 3.567963118121438*^9, 3.567963118698639*^9}, {3.5679631520670977`*^9, 3.5679631541731014`*^9}, {3.5679631898191643`*^9, 3.5679632154656096`*^9}, {3.5679635823158536`*^9, 3.56796361468591*^9}, { 3.567963691516045*^9, 3.567963693590849*^9}, {3.567963737957327*^9, 3.567963807299449*^9}, {3.568048114904332*^9, 3.5680481161835346`*^9}, 3.5680482221389203`*^9, {3.5680485003498087`*^9, 3.5680485123618298`*^9}, 3.5681142659931054`*^9, 3.5681219084621286`*^9, {3.568277848107582*^9, 3.5682778496519847`*^9}, {3.568285867917365*^9, 3.568285868634966*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"Nij", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{"i", "\[LessEqual]", "j"}], ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"Lij", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], ">", "0"}], ",", RowBox[{"2.25", "/", RowBox[{ RowBox[{"Lij", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], "^", "2"}]}], ",", "0"}], "]"}], ",", "\"\< \>\""}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"(*", RowBox[{ RowBox[{"MatrixForm", "[", "Nij", "]"}], ";"}], "*)"}]}]}], "Input", CellChangeTimes->{{3.5679628944794455`*^9, 3.5679629710755796`*^9}, { 3.567963889963994*^9, 3.567963933800071*^9}, {3.5680478391738477`*^9, 3.5680478414358516`*^9}, 3.568048231171336*^9, {3.5680482823550262`*^9, 3.56804828449223*^9}, 3.5680483349271183`*^9, 3.568114273059918*^9, 3.5681219175101447`*^9, {3.568277851305588*^9, 3.568277857295998*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ RowBox[{"Total", " ", "kinetic", " ", "energy", " ", "in", " ", "Eh"}], ",", " ", RowBox[{"=", " ", RowBox[{"|", "Etot", "|"}]}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"Tsigma", "=", RowBox[{ RowBox[{"Sum", "[", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{"i", "\[LessEqual]", "j"}], ",", RowBox[{"Nij", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], ",", "0"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", "n"}], "}"}]}], "]"}], "+", RowBox[{"Sum", "[", RowBox[{ RowBox[{"2.25", "/", RowBox[{ RowBox[{"u2", "[", RowBox[{"[", "i", "]"}], "]"}], "^", "2"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "10"}], "}"}]}], "]"}]}]}]}]], "Input", CellChangeTimes->{{3.567963992456174*^9, 3.567964003890994*^9}, { 3.5679640974287586`*^9, 3.5679642038209453`*^9}, {3.56796433180357*^9, 3.567964353487608*^9}, {3.567964403641696*^9, 3.567964453046983*^9}, { 3.5679646160984693`*^9, 3.5679646168316708`*^9}, {3.567964688747797*^9, 3.567964762255126*^9}, {3.568277859480002*^9, 3.568277863598409*^9}, { 3.568286113290196*^9, 3.5682861137113967`*^9}, {3.5682870401502237`*^9, 3.5682870902263117`*^9}}], Cell[BoxData["388.2993832472005`"], "Output", CellChangeTimes->{3.5682865497165623`*^9, 3.5682870085913687`*^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", "<", RowBox[{"n", "+", "2"}]}], ",", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"j", "=", "2"}], ",", RowBox[{"j", "\[LessEqual]", "i"}], ",", RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"Pij", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], ">", "0"}], ",", RowBox[{"Print", "[", RowBox[{ RowBox[{"u", "[", RowBox[{"[", RowBox[{"i", "-", "1"}], "]"}], "]"}], ",", RowBox[{"i", "-", "1"}], ",", "\"\< \>\"", ",", RowBox[{"u", "[", RowBox[{"[", RowBox[{"j", "-", "1"}], "]"}], "]"}], ",", RowBox[{"j", "-", "1"}], ",", "\"\< \>\"", ",", RowBox[{"Extract", "[", RowBox[{"Pij", ",", RowBox[{"{", RowBox[{"i", ",", "j"}], "}"}]}], "]"}]}], "]"}]}], "]"}]}], "]"}]}], "]"}]], "Input", CellChangeTimes->{{3.568108824079547*^9, 3.568109012824279*^9}, { 3.568109558076436*^9, 3.568109588168889*^9}, {3.5681097523435774`*^9, 3.5681097637783976`*^9}, {3.5681100381828794`*^9, 3.5681101616726966`*^9}, {3.5681102015307665`*^9, 3.5681102288464146`*^9}, {3.568110263743676*^9, 3.568110267331682*^9}, { 3.5681103163469687`*^9, 3.5681103268925867`*^9}, {3.5681106443375444`*^9, 3.5681106916368275`*^9}, 3.568110873002746*^9, {3.5681109388660617`*^9, 3.5681110268970165`*^9}, {3.568111092526332*^9, 3.568111123461186*^9}, { 3.5681111884196997`*^9, 3.568111192756508*^9}, {3.5681112268113675`*^9, 3.568111366213212*^9}, {3.568111434946933*^9, 3.568111463354583*^9}, { 3.568111544396725*^9, 3.568111547625931*^9}, {3.5681116856081734`*^9, 3.5681117131422215`*^9}, {3.568111864696488*^9, 3.5681119835998964`*^9}, 3.5681120218355637`*^9, {3.56811206543764*^9, 3.568112066389242*^9}, { 3.568112332369709*^9, 3.5681123404505234`*^9}, {3.568112505935614*^9, 3.56811251481203*^9}, {3.5681125623609133`*^9, 3.568112569006525*^9}, { 3.568112648067464*^9, 3.568112734210815*^9}, {3.568112764443668*^9, 3.5681127697320776`*^9}, {3.5681128018681335`*^9, 3.568112843426607*^9}, { 3.568112884641879*^9, 3.5681128887602863`*^9}, {3.5681129663080225`*^9, 3.5681129699584293`*^9}, {3.568113133789917*^9, 3.568113227936082*^9}, { 3.5681132946885996`*^9, 3.5681133044698167`*^9}, {3.568113502200164*^9, 3.568113507348173*^9}, 3.56811354556824*^9, {3.5681222886191964`*^9, 3.5681222896019983`*^9}, {3.568207424138049*^9, 3.5682074415164795`*^9}, 3.568277923134914*^9}], Cell[CellGroupData[{ Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"C\"\>", "\[InvisibleSpace]", "5", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", "\[InvisibleSpace]", "1", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "1.2045887149718768`"}], SequenceForm["C", 5, " ", "C", 1, " ", 1.2045887149718768`], Editable->False]], "Print", CellChangeTimes->{3.5682865497477627`*^9, 3.5682870086069684`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"C\"\>", "\[InvisibleSpace]", "5", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", "\[InvisibleSpace]", "2", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "1.1971457193324913`"}], SequenceForm["C", 5, " ", "C", 2, " ", 1.1971457193324913`], Editable->False]], "Print", CellChangeTimes->{3.5682865497477627`*^9, 3.5682870086069684`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"C\"\>", "\[InvisibleSpace]", "6", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", "\[InvisibleSpace]", "1", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "1.197088877556879`"}], SequenceForm["C", 6, " ", "C", 1, " ", 1.197088877556879], Editable->False]], "Print", CellChangeTimes->{3.5682865497477627`*^9, 3.5682870086069684`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"C\"\>", "\[InvisibleSpace]", "6", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", "\[InvisibleSpace]", "3", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "1.2099808163272308`"}], SequenceForm["C", 6, " ", "C", 3, " ", 1.2099808163272308`], Editable->False]], "Print", CellChangeTimes->{3.5682865497477627`*^9, 3.5682870086069684`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"C\"\>", "\[InvisibleSpace]", "7", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", "\[InvisibleSpace]", "1", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "1.2093894227973998`"}], SequenceForm["C", 7, " ", "C", 1, " ", 1.2093894227973998`], Editable->False]], "Print", CellChangeTimes->{3.5682865497477627`*^9, 3.568287008622568*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"C\"\>", "\[InvisibleSpace]", "7", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", "\[InvisibleSpace]", "4", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "1.1969339007869289`"}], SequenceForm["C", 7, " ", "C", 4, " ", 1.1969339007869289`], Editable->False]], "Print", CellChangeTimes->{3.5682865497477627`*^9, 3.568287008622568*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"C\"\>", "\[InvisibleSpace]", "8", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", "\[InvisibleSpace]", "2", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "1.2051327110449268`"}], SequenceForm["C", 8, " ", "C", 2, " ", 1.2051327110449268`], Editable->False]], "Print", CellChangeTimes->{3.5682865497477627`*^9, 3.568287008622568*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"C\"\>", "\[InvisibleSpace]", "8", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", "\[InvisibleSpace]", "3", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "1.190358261769134`"}], SequenceForm["C", 8, " ", "C", 3, " ", 1.190358261769134], Editable->False]], "Print", CellChangeTimes->{3.5682865497477627`*^9, 3.568287008622568*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"C\"\>", "\[InvisibleSpace]", "9", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", "\[InvisibleSpace]", "2", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "1.1957722875393675`"}], SequenceForm["C", 9, " ", "C", 2, " ", 1.1957722875393675`], Editable->False]], "Print", CellChangeTimes->{3.5682865497477627`*^9, 3.568287008622568*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"C\"\>", "\[InvisibleSpace]", "9", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", "\[InvisibleSpace]", "4", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "1.1979175081018865`"}], SequenceForm["C", 9, " ", "C", 4, " ", 1.1979175081018865`], Editable->False]], "Print", CellChangeTimes->{3.5682865497477627`*^9, 3.568287008622568*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"C\"\>", "\[InvisibleSpace]", "10", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", "\[InvisibleSpace]", "3", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "1.2037371872285063`"}], SequenceForm["C", 10, " ", "C", 3, " ", 1.2037371872285063`], Editable->False]], "Print", CellChangeTimes->{3.5682865497477627`*^9, 3.568287008622568*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"C\"\>", "\[InvisibleSpace]", "10", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", "\[InvisibleSpace]", "4", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "1.1829211138696971`"}], SequenceForm["C", 10, " ", "C", 4, " ", 1.1829211138696971`], Editable->False]], "Print", CellChangeTimes->{3.5682865497477627`*^9, 3.568287008622568*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"H\"\>", "\[InvisibleSpace]", "11", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", "\[InvisibleSpace]", "1", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "1.2034203768749958`"}], SequenceForm["H", 11, " ", "C", 1, " ", 1.2034203768749958`], Editable->False]], "Print", CellChangeTimes->{3.5682865497477627`*^9, 3.568287008622568*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"H\"\>", "\[InvisibleSpace]", "12", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", "\[InvisibleSpace]", "2", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "1.2034975108703974`"}], SequenceForm["H", 12, " ", "C", 2, " ", 1.2034975108703974`], Editable->False]], "Print", CellChangeTimes->{3.5682865497477627`*^9, 3.568287008622568*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"H\"\>", "\[InvisibleSpace]", "13", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", "\[InvisibleSpace]", "3", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "1.2034660875171552`"}], SequenceForm["H", 13, " ", "C", 3, " ", 1.2034660875171552`], Editable->False]], "Print", CellChangeTimes->{3.5682865497477627`*^9, 3.568287008622568*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"H\"\>", "\[InvisibleSpace]", "14", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", "\[InvisibleSpace]", "4", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "1.203333069890158`"}], SequenceForm["H", 14, " ", "C", 4, " ", 1.203333069890158], Editable->False]], "Print", CellChangeTimes->{3.5682865497477627`*^9, 3.568287008622568*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"H\"\>", "\[InvisibleSpace]", "15", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", "\[InvisibleSpace]", "5", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "1.203564180530755`"}], SequenceForm["H", 15, " ", "C", 5, " ", 1.203564180530755], Editable->False]], "Print", CellChangeTimes->{3.5682865497477627`*^9, 3.568287008622568*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"H\"\>", "\[InvisibleSpace]", "16", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", "\[InvisibleSpace]", "5", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "1.2034569571985807`"}], SequenceForm["H", 16, " ", "C", 5, " ", 1.2034569571985807`], Editable->False]], "Print", CellChangeTimes->{3.5682865497477627`*^9, 3.5682870086381683`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"H\"\>", "\[InvisibleSpace]", "17", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", "\[InvisibleSpace]", "6", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "1.2036374034442143`"}], SequenceForm["H", 17, " ", "C", 6, " ", 1.2036374034442143`], Editable->False]], "Print", CellChangeTimes->{3.5682865497477627`*^9, 3.5682870086381683`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"H\"\>", "\[InvisibleSpace]", "18", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", "\[InvisibleSpace]", "6", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "1.2035576768401175`"}], SequenceForm["H", 18, " ", "C", 6, " ", 1.2035576768401175`], Editable->False]], "Print", CellChangeTimes->{3.5682865497477627`*^9, 3.5682870086381683`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"H\"\>", "\[InvisibleSpace]", "19", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", "\[InvisibleSpace]", "7", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "1.2034860485372871`"}], SequenceForm["H", 19, " ", "C", 7, " ", 1.2034860485372871`], Editable->False]], "Print", CellChangeTimes->{3.5682865497477627`*^9, 3.5682870086381683`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"H\"\>", "\[InvisibleSpace]", "20", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", "\[InvisibleSpace]", "7", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "1.2033924679450376`"}], SequenceForm["H", 20, " ", "C", 7, " ", 1.2033924679450376`], Editable->False]], "Print", CellChangeTimes->{3.5682865497477627`*^9, 3.5682870086381683`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"H\"\>", "\[InvisibleSpace]", "21", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", "\[InvisibleSpace]", "8", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "1.2034748389037786`"}], SequenceForm["H", 21, " ", "C", 8, " ", 1.2034748389037786`], Editable->False]], "Print", CellChangeTimes->{3.5682865497477627`*^9, 3.5682870086381683`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"H\"\>", "\[InvisibleSpace]", "22", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", "\[InvisibleSpace]", "8", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "1.203449058905057`"}], SequenceForm["H", 22, " ", "C", 8, " ", 1.203449058905057], Editable->False]], "Print", CellChangeTimes->{3.5682865497477627`*^9, 3.5682870086381683`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"H\"\>", "\[InvisibleSpace]", "23", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", "\[InvisibleSpace]", "9", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "1.203394434694922`"}], SequenceForm["H", 23, " ", "C", 9, " ", 1.203394434694922], Editable->False]], "Print", CellChangeTimes->{3.5682865497477627`*^9, 3.5682870086381683`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"H\"\>", "\[InvisibleSpace]", "24", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", "\[InvisibleSpace]", "9", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "1.203396325319866`"}], SequenceForm["H", 24, " ", "C", 9, " ", 1.203396325319866], Editable->False]], "Print", CellChangeTimes->{3.5682865497477627`*^9, 3.5682870086381683`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"H\"\>", "\[InvisibleSpace]", "25", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", "\[InvisibleSpace]", "10", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "1.2035834800504566`"}], SequenceForm["H", 25, " ", "C", 10, " ", 1.2035834800504566`], Editable->False]], "Print", CellChangeTimes->{3.5682865497477627`*^9, 3.5682870086381683`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"H\"\>", "\[InvisibleSpace]", "26", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "\<\"C\"\>", "\[InvisibleSpace]", "10", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "1.2034936068866704`"}], SequenceForm["H", 26, " ", "C", 10, " ", 1.2034936068866704`], Editable->False]], "Print", CellChangeTimes->{3.5682865497477627`*^9, 3.5682870086381683`*^9}] }, Open ]] }, Open ]], Cell[BoxData[{ RowBox[{ RowBox[{"DeleteFile", "[", " ", "\"\\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "2"}], ",", RowBox[{"i", "<", RowBox[{"n", "+", "2"}]}], ",", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"j", "=", "2"}], ",", RowBox[{"j", "\[LessEqual]", "i"}], ",", RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"Pij", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], ">", "0"}], ",", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"i", "-", "1"}], ",", RowBox[{"j", "-", "1"}], ",", RowBox[{"Extract", "[", RowBox[{"Pij", ",", RowBox[{"{", RowBox[{"i", ",", "j"}], "}"}]}], "]"}]}], "}"}], ">>>", " ", "\"\\""}]}], "]"}]}], "]"}]}], "]"}]}], "Input", CellChangeTimes->{{3.5681154444347754`*^9, 3.568115497147268*^9}, { 3.568127214924899*^9, 3.568127265796589*^9}, {3.568128215276656*^9, 3.568128238536297*^9}, {3.568130718566253*^9, 3.5681307485339055`*^9}, 3.5681308175952272`*^9, {3.5681308916497574`*^9, 3.5681309458598523`*^9}, { 3.5681309983383446`*^9, 3.5681310118947687`*^9}, {3.5681310723448744`*^9, 3.568131074357278*^9}, {3.5681315215476637`*^9, 3.568131584509374*^9}, { 3.56813162789305*^9, 3.5681316427754765`*^9}, {3.5681317113063965`*^9, 3.568131715268804*^9}, {3.5681323203938665`*^9, 3.5681323260098763`*^9}, { 3.568132387567585*^9, 3.568132440014877*^9}, {3.5681325121026034`*^9, 3.568132590945142*^9}, {3.56813274298981*^9, 3.5681327548302307`*^9}, { 3.568132877867647*^9, 3.568132888335265*^9}, {3.568133052345154*^9, 3.568133072344389*^9}, {3.5681331740533676`*^9, 3.5681331951446047`*^9}, { 3.568133267793932*^9, 3.5681332684803333`*^9}, {3.5681335079111547`*^9, 3.568133536833606*^9}, {3.5681336272201643`*^9, 3.5681336314945717`*^9}, { 3.568133737994359*^9, 3.568133744000369*^9}, {3.568133978515581*^9, 3.5681339917912045`*^9}, {3.5681340324916763`*^9, 3.5681340418204927`*^9}, {3.568134244262048*^9, 3.5681342454944506`*^9}, { 3.5681343354754086`*^9, 3.568134336099409*^9}, {3.5681345828762426`*^9, 3.568134583952645*^9}, {3.5681348406666956`*^9, 3.5681348474683075`*^9}, 3.568134889151581*^9, {3.5681352914606876`*^9, 3.568135296265496*^9}, { 3.568135518877887*^9, 3.568135519595488*^9}, 3.5681357593679094`*^9, { 3.5681787035360527`*^9, 3.568178728449296*^9}, 3.5682779128076954`*^9}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"a", "=", RowBox[{ "ReadList", "[", "\"\\"", "]"}]}]], "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.2045887149718768`"}], "}"}], ",", RowBox[{"{", RowBox[{"5", ",", "2", ",", "1.1971457193324913`"}], "}"}], ",", RowBox[{"{", RowBox[{"6", ",", "1", ",", "1.197088877556879`"}], "}"}], ",", RowBox[{"{", RowBox[{"6", ",", "3", ",", "1.2099808163272308`"}], "}"}], ",", RowBox[{"{", RowBox[{"7", ",", "1", ",", "1.2093894227973998`"}], "}"}], ",", RowBox[{"{", RowBox[{"7", ",", "4", ",", "1.1969339007869289`"}], "}"}], ",", RowBox[{"{", RowBox[{"8", ",", "2", ",", "1.2051327110449268`"}], "}"}], ",", RowBox[{"{", RowBox[{"8", ",", "3", ",", "1.190358261769134`"}], "}"}], ",", RowBox[{"{", RowBox[{"9", ",", "2", ",", "1.1957722875393675`"}], "}"}], ",", RowBox[{"{", RowBox[{"9", ",", "4", ",", "1.1979175081018865`"}], "}"}], ",", RowBox[{"{", RowBox[{"10", ",", "3", ",", "1.2037371872285063`"}], "}"}], ",", RowBox[{"{", RowBox[{"10", ",", "4", ",", "1.1829211138696971`"}], "}"}], ",", RowBox[{"{", RowBox[{"11", ",", "1", ",", "1.2034203768749958`"}], "}"}], ",", RowBox[{"{", RowBox[{"12", ",", "2", ",", "1.2034975108703974`"}], "}"}], ",", RowBox[{"{", RowBox[{"13", ",", "3", ",", "1.2034660875171552`"}], "}"}], ",", RowBox[{"{", RowBox[{"14", ",", "4", ",", "1.203333069890158`"}], "}"}], ",", RowBox[{"{", RowBox[{"15", ",", "5", ",", "1.203564180530755`"}], "}"}], ",", RowBox[{"{", RowBox[{"16", ",", "5", ",", "1.2034569571985807`"}], "}"}], ",", RowBox[{"{", RowBox[{"17", ",", "6", ",", "1.2036374034442143`"}], "}"}], ",", RowBox[{"{", RowBox[{"18", ",", "6", ",", "1.2035576768401175`"}], "}"}], ",", RowBox[{"{", RowBox[{"19", ",", "7", ",", "1.2034860485372871`"}], "}"}], ",", RowBox[{"{", RowBox[{"20", ",", "7", ",", "1.2033924679450376`"}], "}"}], ",", RowBox[{"{", RowBox[{"21", ",", "8", ",", "1.2034748389037786`"}], "}"}], ",", RowBox[{"{", RowBox[{"22", ",", "8", ",", "1.203449058905057`"}], "}"}], ",", RowBox[{"{", RowBox[{"23", ",", "9", ",", "1.203394434694922`"}], "}"}], ",", RowBox[{"{", RowBox[{"24", ",", "9", ",", "1.203396325319866`"}], "}"}], ",", RowBox[{"{", RowBox[{"25", ",", "10", ",", "1.2035834800504566`"}], "}"}], ",", RowBox[{"{", RowBox[{"26", ",", "10", ",", "1.2034936068866704`"}], "}"}]}], "}"}]], "Output", CellChangeTimes->{3.5682865498101625`*^9, 3.5682870086693687`*^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", ",", "12"}], "}"}]}], "]"}]}], ";"}], "\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", ",", "12"}], "}"}]}], "]"}]}], ";"}], "\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", ",", "12"}], "}"}]}], "]"}]}], ";"}], "\[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", ",", "12"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"wx", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "3"}], "]"}], "]"}], "+", RowBox[{"u2", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ")"}], "/", RowBox[{"(", RowBox[{ RowBox[{"u2", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "+", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "3"}], "]"}], "]"}], "+", RowBox[{"u2", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ")"}]}], "*", RowBox[{"(", RowBox[{ RowBox[{"x", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "-", RowBox[{"x", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ")"}]}], "+", RowBox[{"x", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "13", ",", "28"}], "}"}]}], "]"}]}], ";"}], "\n", " ", RowBox[{ RowBox[{"wy", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "3"}], "]"}], "]"}], "+", RowBox[{"u2", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ")"}], "/", RowBox[{"(", RowBox[{ RowBox[{"u2", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "+", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "3"}], "]"}], "]"}], "+", RowBox[{"u2", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ")"}]}], "*", RowBox[{"(", RowBox[{ RowBox[{"y", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "-", RowBox[{"y", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ")"}]}], "+", RowBox[{"y", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "13", ",", "28"}], "}"}]}], "]"}]}], ";"}], "\n", " ", RowBox[{ RowBox[{"wz", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "3"}], "]"}], "]"}], "+", RowBox[{"u2", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ")"}], "/", RowBox[{"(", RowBox[{ RowBox[{"u2", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "+", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "3"}], "]"}], "]"}], "+", RowBox[{"u2", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ")"}]}], "*", RowBox[{"(", RowBox[{ RowBox[{"z", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "-", RowBox[{"z", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ")"}]}], "+", RowBox[{"z", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "13", ",", "28"}], "}"}]}], "]"}]}], ";"}], "\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", ",", "16"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"Show", "[", RowBox[{"Plot1", ",", "Plot7", ",", "Plot8", ",", RowBox[{"{", RowBox[{ RowBox[{"AspectRatio", "\[Rule]", "Automatic"}], ",", RowBox[{"Boxed", "\[Rule]", "False"}]}], "}"}]}], "]"}]}]}]], "Input", CellChangeTimes->{{3.5681365021632137`*^9, 3.5681365171392403`*^9}, { 3.568136560351316*^9, 3.5681365863097615`*^9}, {3.5681369135895367`*^9, 3.5681369265999594`*^9}, {3.568137821916732*^9, 3.5681378378443604`*^9}, { 3.5682072492149415`*^9, 3.5682072529901485`*^9}, 3.568209802830227*^9, { 3.5682216503836985`*^9, 3.5682217583046885`*^9}, {3.5682217912831464`*^9, 3.568221843465238*^9}, {3.5682219760654707`*^9, 3.5682221133457117`*^9}, { 3.568222159615393*^9, 3.5682221693498106`*^9}, {3.5682222565227633`*^9, 3.5682222578487663`*^9}, 3.5682224736283445`*^9, 3.5682227227607822`*^9, 3.5682227558640404`*^9, {3.568222818810151*^9, 3.5682228486062036`*^9}, { 3.5682230057764797`*^9, 3.5682230079916835`*^9}, {3.568223041952943*^9, 3.5682230548229656`*^9}, 3.5682232619601297`*^9, 3.5682233000709963`*^9, { 3.568223536645412*^9, 3.5682235381118145`*^9}, {3.568223591151908*^9, 3.568223591526308*^9}, {3.5682236491528096`*^9, 3.568223657483224*^9}, { 3.568223713222122*^9, 3.5682237532361927`*^9}, {3.568223800379475*^9, 3.568223827648323*^9}, {3.5682238887380304`*^9, 3.5682238893152313`*^9}, { 3.5682240287170763`*^9, 3.5682240612431335`*^9}, 3.5682241088076167`*^9, { 3.568224198772975*^9, 3.5682242396450467`*^9}, {3.568224282170721*^9, 3.5682243185031853`*^9}, {3.568224437890195*^9, 3.568224439590598*^9}, { 3.568224572253231*^9, 3.5682245867144566`*^9}, {3.5682247081762695`*^9, 3.568224733557514*^9}, {3.56822480520844*^9, 3.5682248433505073`*^9}, { 3.568224884206979*^9, 3.5682248959069996`*^9}, {3.568224982830352*^9, 3.5682250047639904`*^9}, {3.56822503870965*^9, 3.5682250803305235`*^9}, { 3.5682251223413973`*^9, 3.568225141451431*^9}, {3.568225234177994*^9, 3.5682252929588966`*^9}, {3.5682253775890455`*^9, 3.568225482140429*^9}, { 3.56822551694409*^9, 3.5682255297205124`*^9}, {3.568225573197789*^9, 3.5682255809666023`*^9}, 3.568225660698343*^9, 3.5682257301964645`*^9, { 3.568225765078126*^9, 3.5682257828933573`*^9}, {3.568225927026843*^9, 3.568225927416844*^9}, {3.568267062420238*^9, 3.568267081358671*^9}, { 3.568267188328059*^9, 3.5682672830358257`*^9}, {3.5682676367040462`*^9, 3.568267747979042*^9}, 3.5682678404716043`*^9, 3.568267882030078*^9, { 3.5682679252577534`*^9, 3.568267996612279*^9}, {3.5682680488879704`*^9, 3.568268061133992*^9}, {3.568268100914062*^9, 3.5682681150944867`*^9}, { 3.5682681672609787`*^9, 3.5682681709113846`*^9}, {3.568268272904364*^9, 3.5682683182536435`*^9}, {3.5682685202739983`*^9, 3.568268599194537*^9}, { 3.5682686348873997`*^9, 3.568268637773405*^9}, {3.5682686944951043`*^9, 3.568268872273017*^9}, {3.5682689420207396`*^9, 3.568269016947671*^9}, { 3.5682695040741262`*^9, 3.568269601153097*^9}, {3.5682696652380095`*^9, 3.5682696728352227`*^9}, {3.5682697489633565`*^9, 3.5682698308635006`*^9}, {3.568271156491429*^9, 3.568271160063835*^9}, { 3.568271244927984*^9, 3.568271369010602*^9}, {3.5682716250070515`*^9, 3.5682717312276382`*^9}, {3.5682717843301315`*^9, 3.5682718576814604`*^9}, {3.568271940673606*^9, 3.5682720536490045`*^9}, { 3.5682722259669075`*^9, 3.5682723223126764`*^9}, {3.568272370813162*^9, 3.568272411466833*^9}, {3.5682725893539457`*^9, 3.568272622238803*^9}, { 3.5682726936713285`*^9, 3.5682726960113325`*^9}, {3.568272741079812*^9, 3.56827275143823*^9}, {3.5682729218373294`*^9, 3.568272991070251*^9}, { 3.5682734285574193`*^9, 3.5682734302890224`*^9}, {3.5682734898811274`*^9, 3.5682735111595645`*^9}, {3.5682736115925407`*^9, 3.5682736313577757`*^9}, {3.5682738601945777`*^9, 3.568273987272401*^9}, { 3.5682763944346285`*^9, 3.5682764313910933`*^9}, 3.5682765128232365`*^9, { 3.5682766527398825`*^9, 3.568276660976697*^9}, {3.5682774496048822`*^9, 3.568277461663703*^9}, {3.5682852978767633`*^9, 3.568285382460112*^9}, { 3.568285675366227*^9, 3.568285707861084*^9}, {3.568285985822372*^9, 3.568285991485182*^9}, {3.568286358600627*^9, 3.568286369505046*^9}, { 3.5682869898869357`*^9, 3.568287001134555*^9}}], Cell[BoxData[ Graphics3DBox[{ SphereBox[{1.19719811278, -2.0715176412, 1.69614246856}, 0.2555], SphereBox[{1.7163625367599997`, 0.49068625316, -2.2937494187999996`}, 0.2555], SphereBox[{-0.30091996824, 2.4556989369999997`, 1.5222687793}, 0.2555], SphereBox[{-2.6074439347999996`, -0.87883597356, -0.92035325378}, 0.2555], SphereBox[{2.9186818069999996`, -1.5428667927, -0.6026525186599999}, 0.2555], SphereBox[{0.89495533634, 0.38508836428, 3.2172585149999997`}, 0.2555], SphereBox[{-1.42978558886, -2.9498622859999997`, 0.74164186596}, 0.2555], SphereBox[{1.43243120526, 2.9477835874, -0.7394497838}, 0.2555], SphereBox[{-0.89230971994, -0.38712926836, -3.2151798163999996`}, 0.2555], SphereBox[{-2.9160361906, 1.5408258886200001`, 0.6048446008199999}, 0.2555], SphereBox[{2.0104794914, -3.5094101546, 2.8619900269999996`}, 0.1], SphereBox[{2.8937374237999998`, 0.81681516624, -3.9049298063999998`}, 0.1], SphereBox[{-0.52574067046, 4.1713811724, 2.5683266066}, 0.1], SphereBox[{-4.4259272646, -1.4875734099399998`, -1.56095147052}, 0.1], SphereBox[{4.716756095999999, -0.88716966522, 0.04983018489399999}, 0.1], SphereBox[{3.1564093378, -3.2452264597999996`, -1.6674375306199998`}, 0.1], SphereBox[{2.7110009196, 1.00671373198, 3.8533402866}, 0.1], SphereBox[{-0.30290418053999996`, 0.04616033700199999, 4.810675478199999}, 0.1], SphereBox[{-2.6350339344, -3.3039969384, 2.3260637334}, 0.1], SphereBox[{-1.21964805766, -4.6436236998, -0.34236165942}, 0.1], SphereBox[{3.2486279666, 3.546070839, -0.08228622894399999}, 0.1], SphereBox[{0.62986457306, 4.39361295, -1.902954082}, 0.1], SphereBox[{-0.6935483392599999, -2.0889031204, -4.288733156999999}, 0.1], SphereBox[{-1.7154743655399998`, 1.04730504646, -4.3783061694}, 0.1], SphereBox[{-4.1024061733999995`, 1.18260942806, 2.2026646255999998`}, 0.1], SphereBox[{-3.7592319318, 2.9657359843999997`, -0.55580621112}, 0.1], SphereBox[{2.0579399598899997`, -1.8071922169499999`, 0.5467449749500001}, 1.2045887149718768`], SphereBox[{2.3175221718799994`, -0.5260902697700001, -1.4482009687299997`}, 1.1971457193324913`], SphereBox[{1.04607672456, -0.8432146384599997, 2.45670049178}, 1.197088877556879], SphereBox[{0.29701768405000006`, 1.42039365064, 2.36976364715}, 1.2099808163272308`], SphereBox[{-0.11629373803999998`, -2.5106899636, 1.21889216726}, 1.2093894227973998`], SphereBox[{-2.01861476183, -1.9143491297799997`, -0.08935569390999998}, 1.1969339007869289`], SphereBox[{1.5743968710099998`, 1.71923492028, -1.5165996012999998`}, 1.2051327110449268`], SphereBox[{0.5657556185099999, 2.7017412621999997`, 0.39140949775}, 1.190358261769134], SphereBox[{0.41202640840999993`, 0.05177849239999999, -2.7544646175999996`}, 1.1957722875393675`], SphereBox[{-1.7498768273699996`, -0.63298262096, -2.0677665350899996`}, 1.1979175081018865`], SphereBox[{-1.6084780794200002`, 1.99826241281, 1.06355669006}, 1.2037371872285063`], SphereBox[{-2.7617400626999995`, 0.3309949575300002, -0.1577543264800001}, 1.1829211138696971`], {Opacity[0.5], SphereBox[{1.7840227914231503`, -3.1090316360871273`, 2.537361924474922}, 1.2045887149718768`]}, {Opacity[0.5], SphereBox[{2.5659120640555253`, 0.7260086360257257, -3.456316712990044}, 1.1971457193324913`]}, {Opacity[0.5], SphereBox[{-0.463141175774466, 3.6936634181342503`, 2.2770603186015355`}, 1.197088877556879]}, {Opacity[0.5], SphereBox[{-3.9195520697709725`, -1.3180642740049748`, \ -1.3825703759599959`}, 1.2099808163272308`]}, {Opacity[0.5], SphereBox[{4.216121232648983, -1.0697343366709047`, -0.1318394997964631}, 1.2093894227973998`]}, {Opacity[0.5], SphereBox[{3.0902157428454187`, -2.771216125818805, -1.3709554716922714`}, 1.1969339007869289`]}, {Opacity[0.5], SphereBox[{2.2053806478845788`, 0.833641871769609, 3.676243499407972}, 1.2051327110449268`]}, {Opacity[0.5], SphereBox[{0.030614944917054965`, 0.14052781303547504`, 4.367021590105425}, 1.190358261769134]}, {Opacity[0.5], SphereBox[{-2.2994456552277036`, -3.205392000227736, 1.8848987154642836`}, 1.1957722875393675`]}, {Opacity[0.5], SphereBox[{-1.2781612733434444`, -4.171992447847919, \ -0.04051871056663403}, 1.1979175081018865`]}, {Opacity[0.5], SphereBox[{2.7429249494502996`, 3.379483344353718, -0.26526728007543526`}, 1.2037371872285063`]}, {Opacity[0.5], SphereBox[{0.8533346016141747, 3.9910301413190505`, -1.5789830516708137`}, 1.1829211138696971`]}, {Opacity[0.5], SphereBox[{-0.7488937892511736, -1.6150412499393472`, \ -3.9898003737713843`}, 1.2034203768749958`]}, {Opacity[0.5], SphereBox[{-1.486262961538124, 0.6478847056985126, -4.05443192889689}, 1.2034975108703974`]}, {Opacity[0.5], SphereBox[{-3.772090285403262, 1.2823460953563708`, 1.757790473357436}, 1.2034660875171552`]}, {Opacity[0.5], SphereBox[{-3.5244541359262826`, 2.568986765613225, -0.2326368376379776}, 1.203333069890158]}}, AspectRatio->Automatic, Boxed->False, ImageSize->{416.68948410995415`, 458.8889010452933}, ViewPoint->{-3.0249433707382236`, 1.5162484141180566`, -0.02661485508780986}, ViewVertical->{-0.9084979034465823, -0.4593067010588772, 0.11563741355385948`}]], "Output", CellChangeTimes->{3.568286549841363*^9, 3.5682870087005687`*^9}] }, Open ]] }, WindowSize->{963, 775}, WindowMargins->{{150, Automatic}, {Automatic, 31}}, 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, 686, 15, 50, "Input"], Cell[1268, 39, 3478, 61, 226, "Output"] }, Open ]], Cell[4761, 103, 644, 13, 29, "Input"], Cell[CellGroupData[{ Cell[5430, 120, 658, 9, 29, "Input"], Cell[6091, 131, 505, 17, 28, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[6633, 153, 827, 13, 29, "Input"], Cell[7463, 168, 609, 9, 28, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[8109, 182, 213, 5, 50, "Input"], Cell[8325, 189, 95, 1, 28, "Output"] }, Open ]], Cell[8435, 193, 1100, 19, 50, "Input"], Cell[9538, 214, 827, 22, 50, "Input"], Cell[CellGroupData[{ Cell[10390, 240, 203, 4, 29, "Input"], Cell[10596, 246, 105, 1, 28, "Output"] }, Open ]], Cell[10716, 250, 720, 13, 29, "Input"], Cell[11439, 265, 478, 10, 29, "Input"], Cell[11920, 277, 3788, 79, 50, "Input"], Cell[15711, 358, 3805, 74, 70, "Input"], Cell[CellGroupData[{ Cell[19541, 436, 2500, 35, 26, "Input"], Cell[22044, 473, 2101, 31, 329, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[24182, 509, 672, 10, 26, "Input"], Cell[24857, 521, 19796, 485, 208, "Output"] }, Open ]], Cell[44668, 1009, 1254, 35, 62, "Input"], Cell[45925, 1046, 4859, 112, 98, "Input"], Cell[50787, 1160, 1985, 50, 62, "Input"], Cell[52775, 1212, 808, 14, 44, "Input"], Cell[53586, 1228, 1708, 43, 80, "Input"], Cell[55297, 1273, 511, 11, 44, "Input"], Cell[55811, 1286, 2048, 48, 80, "Input"], Cell[57862, 1336, 1771, 42, 44, "Input"], Cell[59636, 1380, 1261, 33, 44, "Input"], Cell[CellGroupData[{ Cell[60922, 1417, 1432, 36, 44, "Input"], Cell[62357, 1455, 113, 1, 25, "Output"] }, Open ]], Cell[62485, 1459, 1971, 35, 80, "Input"], Cell[64459, 1496, 1638, 24, 26, "Input"], Cell[66100, 1522, 1655, 24, 26, "Input"], Cell[CellGroupData[{ Cell[67780, 1550, 2749, 55, 62, "Input"], Cell[CellGroupData[{ Cell[70554, 1609, 431, 8, 18, "Print"], Cell[70988, 1619, 431, 8, 18, "Print"], Cell[71422, 1629, 428, 8, 18, "Print"], Cell[71853, 1639, 431, 8, 18, "Print"], Cell[72287, 1649, 429, 8, 18, "Print"], Cell[72719, 1659, 429, 8, 18, "Print"], Cell[73151, 1669, 429, 8, 18, "Print"], Cell[73583, 1679, 426, 8, 18, "Print"], Cell[74012, 1689, 429, 8, 18, "Print"], Cell[74444, 1699, 429, 8, 18, "Print"], Cell[74876, 1709, 431, 8, 18, "Print"], Cell[75310, 1719, 431, 8, 18, "Print"], Cell[75744, 1729, 431, 8, 18, "Print"], Cell[76178, 1739, 431, 8, 18, "Print"], Cell[76612, 1749, 431, 8, 18, "Print"], Cell[77046, 1759, 428, 8, 18, "Print"], Cell[77477, 1769, 428, 8, 18, "Print"], Cell[77908, 1779, 433, 8, 18, "Print"], Cell[78344, 1789, 433, 8, 18, "Print"], Cell[78780, 1799, 433, 8, 18, "Print"], Cell[79216, 1809, 433, 8, 18, "Print"], Cell[79652, 1819, 433, 8, 18, "Print"], Cell[80088, 1829, 433, 8, 18, "Print"], Cell[80524, 1839, 430, 8, 18, "Print"], Cell[80957, 1849, 430, 8, 18, "Print"], Cell[81390, 1859, 430, 8, 18, "Print"], Cell[81823, 1869, 435, 8, 18, "Print"], Cell[82261, 1879, 435, 8, 18, "Print"] }, Open ]] }, Open ]], Cell[82723, 1891, 2696, 53, 80, "Input"], Cell[CellGroupData[{ Cell[85444, 1948, 675, 11, 26, "Input"], Cell[86122, 1961, 2573, 60, 77, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[88732, 2026, 17024, 443, 404, "Input"], Cell[105759, 2471, 5567, 108, 383, "Output"] }, Open ]] } ] *) (* End of internal cache information *)