(* 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[ 181594, 4027] NotebookOptionsPosition[ 176263, 3862] NotebookOutlinePosition[ 176798, 3883] CellTagsIndexPosition[ 176755, 3880] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ RowBox[{ RowBox[{ "Kimball", " ", "model", " ", "and", " ", "energy", " ", "derived", " ", "from", " ", RowBox[{"MP2", "/", "6"}]}], "-", RowBox[{"31", "G", RowBox[{"(", "d", ")"}], " ", "data", " ", "source"}]}], ",", " ", RowBox[{"ES", " ", "19", " ", "June", " ", "2013"}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Clear", "[", RowBox[{"a", ",", "u", ",", "u1", ",", "u2", ",", "t"}], "]"}], ";"}]}]], "Input", CellChangeTimes->{{3.568278170270548*^9, 3.568278174310955*^9}, { 3.5682782109554195`*^9, 3.5682782183342323`*^9}, 3.568284752749406*^9, 3.5682847957898817`*^9, {3.5682872948362713`*^9, 3.568287369232802*^9}, 3.5682960154833097`*^9, {3.568305453232698*^9, 3.5683054626863146`*^9}, { 3.568307776310778*^9, 3.56830777732478*^9}, {3.5804443791490765`*^9, 3.58044438139548*^9}, {3.5804469366841097`*^9, 3.5804469795061855`*^9}, { 3.580459588972727*^9, 3.580459618113578*^9}, {3.5805250613979683`*^9, 3.5805250682307806`*^9}, {3.580533517520193*^9, 3.5805335275198107`*^9}, { 3.580623964920255*^9, 3.5806239657314568`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"t", "=", RowBox[{"1.889726", "*", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0.0000000000", ",", RowBox[{"-", "0.07962744727"}], ",", "0.0000000000"}], "}"}], ",", RowBox[{"{", RowBox[{"0.7632077729", ",", "0.5169595055", ",", "0.0000000000"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "0.7632077729"}], ",", "0.5169595055", ",", "0.0000000000"}], "}"}]}], "}"}]}]}], ";"}]], "Input", CellChangeTimes->{{3.567862912170579*^9, 3.567862949329844*^9}, { 3.5678630095303497`*^9, 3.5678630107471523`*^9}, {3.5679573013548217`*^9, 3.5679573132576427`*^9}, 3.568114367533684*^9, {3.5682782383334675`*^9, 3.5682782435438766`*^9}, 3.5682847420789876`*^9, 3.568284795821082*^9, { 3.568287382898426*^9, 3.5682873928356433`*^9}, 3.5682960154833097`*^9, 3.5683055034491863`*^9, {3.5804444043275204`*^9, 3.580444406558325*^9}, { 3.5804469063576565`*^9, 3.58044691935248*^9}, 3.580447101826*^9, 3.58048985172836*^9, {3.5805250473735437`*^9, 3.58052505082115*^9}, { 3.5805257978758163`*^9, 3.5805257998726196`*^9}, 3.580533417945218*^9, { 3.5805334506428757`*^9, 3.5805334825293317`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"u", "=", RowBox[{"{", 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, 3.56829601549891*^9, 3.5683055034647865`*^9, {3.5804444900184712`*^9, 3.5804444924208755`*^9}, {3.5804469928442087`*^9, 3.580447008101035*^9}, { 3.5805253295468397`*^9, 3.580525331153642*^9}}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{"n", "=", RowBox[{"Length", "[", "u", "]"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"num1", "=", RowBox[{"Table", "[", RowBox[{"i", ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"u1", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\\""}], ",", "0.255", ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\\""}], ",", "0.12", ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\\""}], ",", "0.24", ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\\""}], ",", "0.18", ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\\""}], ",", "0.5"}], "]"}]}], "]"}]}], "]"}]}], "]"}]}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], "\[IndentingNewLine]", RowBox[{"u2", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\\""}], ",", RowBox[{"(*", "0.25618", "*)"}], "0.2562", ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\\""}], ",", "0.6", ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\\""}], ",", "0.22522", ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\\""}], ",", StyleBox["0.1792", FontColor->RGBColor[1, 0, 0]], ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\\""}], ",", "0.48698"}], "]"}]}], "]"}]}], "]"}]}], "]"}]}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], "\[IndentingNewLine]", RowBox[{"cn", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\\""}], ",", "6.0", ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\\""}], ",", "1.0", ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\\""}], ",", "7.0", ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\\""}], ",", "8.0", ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\\""}], ",", "7.0"}], "]"}]}], "]"}]}], "]"}]}], "]"}]}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"bc", "=", "0"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", " ", RowBox[{"i", "<", RowBox[{"n", "+", "1"}]}], ",", RowBox[{"i", "++"}], ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"[", "i", "]"}], "]"}], "\[NotEqual]", "\"\\""}], ",", " ", RowBox[{"bc", "=", RowBox[{"bc", "+", "1"}]}]}], "]"}]}], "]"}], ";"}], "\[IndentingNewLine]", "bc", "\[IndentingNewLine]", RowBox[{"bh", "=", RowBox[{"n", "-", "bc"}]}]}], "Input", CellChangeTimes->{{3.568277663590458*^9, 3.568277712652544*^9}, 3.5682785838808746`*^9, 3.580447125101241*^9, {3.580447338150815*^9, 3.580447473683853*^9}, {3.580447529251151*^9, 3.5804475516215897`*^9}, { 3.580448220971966*^9, 3.580448222141968*^9}, {3.5804483597810097`*^9, 3.5804483678618236`*^9}, {3.5804484008402815`*^9, 3.5804484262371264`*^9}, {3.5804484577023816`*^9, 3.5804484578895817`*^9}, {3.5804485097752733`*^9, 3.5804485333157144`*^9}, 3.58044856532697*^9, {3.580449794250329*^9, 3.5804497945467296`*^9}, { 3.580449861034046*^9, 3.580449943339791*^9}, {3.580450110182084*^9, 3.5804501195421004`*^9}, {3.5804595061521816`*^9, 3.580459540721842*^9}, 3.5804625847395887`*^9, {3.5804626281388645`*^9, 3.580462631758071*^9}, { 3.5804626688549366`*^9, 3.5804627184318237`*^9}, {3.5804763932858744`*^9, 3.5804763978722825`*^9}, 3.5804765357297244`*^9, 3.5804765911254215`*^9, { 3.580476629189489*^9, 3.5804766297510896`*^9}, {3.580476668641958*^9, 3.5804766699055605`*^9}, {3.580477336588331*^9, 3.5804773389439354`*^9}, { 3.58047738137601*^9, 3.580477410969262*^9}, 3.580623917574172*^9}], Cell[BoxData["3"], "Output", CellChangeTimes->{ 3.5806206424420195`*^9, 3.5806213346464357`*^9, 3.5806222658588715`*^9, 3.580622297667327*^9, 3.580622353375025*^9, {3.580622387039884*^9, 3.580622429830759*^9}, {3.5806224688776274`*^9, 3.580622518173714*^9}, { 3.580622550777772*^9, 3.580622580573824*^9}, 3.5806227155608606`*^9, 3.5806228001130095`*^9, 3.580622873417538*^9, 3.5806229089232006`*^9, { 3.580622941901658*^9, 3.580622966674502*^9}, 3.5806230101829786`*^9, { 3.580623063441472*^9, 3.5806230931907244`*^9}, {3.5806231448268147`*^9, 3.580623172002063*^9}, {3.580623203482918*^9, 3.580623223918954*^9}, 3.5806233379707537`*^9, 3.580623395207255*^9, {3.58062342665691*^9, 3.5806234736909924`*^9}, 3.580623536434303*^9, {3.580623593608403*^9, 3.580623622593254*^9}, 3.580623682107359*^9, 3.580623718595823*^9, 3.5806237837571373`*^9, 3.580623816017994*^9, 3.5806238611020727`*^9, 3.580623922675381*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{"0.18`", ",", "0.12`", ",", "0.12`"}], "}"}]], "Output", CellChangeTimes->{ 3.5806206424420195`*^9, 3.5806213346464357`*^9, 3.5806222658588715`*^9, 3.580622297667327*^9, 3.580622353375025*^9, {3.580622387039884*^9, 3.580622429830759*^9}, {3.5806224688776274`*^9, 3.580622518173714*^9}, { 3.580622550777772*^9, 3.580622580573824*^9}, 3.5806227155608606`*^9, 3.5806228001130095`*^9, 3.580622873417538*^9, 3.5806229089232006`*^9, { 3.580622941901658*^9, 3.580622966674502*^9}, 3.5806230101829786`*^9, { 3.580623063441472*^9, 3.5806230931907244`*^9}, {3.5806231448268147`*^9, 3.580623172002063*^9}, {3.580623203482918*^9, 3.580623223918954*^9}, 3.5806233379707537`*^9, 3.580623395207255*^9, {3.58062342665691*^9, 3.5806234736909924`*^9}, 3.580623536434303*^9, {3.580623593608403*^9, 3.580623622593254*^9}, 3.580623682107359*^9, 3.580623718595823*^9, 3.5806237837571373`*^9, 3.580623816017994*^9, 3.5806238611020727`*^9, 3.580623922675381*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{"0.1792`", ",", "0.6`", ",", "0.6`"}], "}"}]], "Output", CellChangeTimes->{ 3.5806206424420195`*^9, 3.5806213346464357`*^9, 3.5806222658588715`*^9, 3.580622297667327*^9, 3.580622353375025*^9, {3.580622387039884*^9, 3.580622429830759*^9}, {3.5806224688776274`*^9, 3.580622518173714*^9}, { 3.580622550777772*^9, 3.580622580573824*^9}, 3.5806227155608606`*^9, 3.5806228001130095`*^9, 3.580622873417538*^9, 3.5806229089232006`*^9, { 3.580622941901658*^9, 3.580622966674502*^9}, 3.5806230101829786`*^9, { 3.580623063441472*^9, 3.5806230931907244`*^9}, {3.5806231448268147`*^9, 3.580623172002063*^9}, {3.580623203482918*^9, 3.580623223918954*^9}, 3.5806233379707537`*^9, 3.580623395207255*^9, {3.58062342665691*^9, 3.5806234736909924`*^9}, 3.580623536434303*^9, {3.580623593608403*^9, 3.580623622593254*^9}, 3.580623682107359*^9, 3.580623718595823*^9, 3.5806237837571373`*^9, 3.580623816017994*^9, 3.5806238611020727`*^9, 3.580623922675381*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{"8.`", ",", "1.`", ",", "1.`"}], "}"}]], "Output", CellChangeTimes->{ 3.5806206424420195`*^9, 3.5806213346464357`*^9, 3.5806222658588715`*^9, 3.580622297667327*^9, 3.580622353375025*^9, {3.580622387039884*^9, 3.580622429830759*^9}, {3.5806224688776274`*^9, 3.580622518173714*^9}, { 3.580622550777772*^9, 3.580622580573824*^9}, 3.5806227155608606`*^9, 3.5806228001130095`*^9, 3.580622873417538*^9, 3.5806229089232006`*^9, { 3.580622941901658*^9, 3.580622966674502*^9}, 3.5806230101829786`*^9, { 3.580623063441472*^9, 3.5806230931907244`*^9}, {3.5806231448268147`*^9, 3.580623172002063*^9}, {3.580623203482918*^9, 3.580623223918954*^9}, 3.5806233379707537`*^9, 3.580623395207255*^9, {3.58062342665691*^9, 3.5806234736909924`*^9}, 3.580623536434303*^9, {3.580623593608403*^9, 3.580623622593254*^9}, 3.580623682107359*^9, 3.580623718595823*^9, 3.5806237837571373`*^9, 3.580623816017994*^9, 3.5806238611020727`*^9, 3.580623922690981*^9}], Cell[BoxData["1"], "Output", CellChangeTimes->{ 3.5806206424420195`*^9, 3.5806213346464357`*^9, 3.5806222658588715`*^9, 3.580622297667327*^9, 3.580622353375025*^9, {3.580622387039884*^9, 3.580622429830759*^9}, {3.5806224688776274`*^9, 3.580622518173714*^9}, { 3.580622550777772*^9, 3.580622580573824*^9}, 3.5806227155608606`*^9, 3.5806228001130095`*^9, 3.580622873417538*^9, 3.5806229089232006`*^9, { 3.580622941901658*^9, 3.580622966674502*^9}, 3.5806230101829786`*^9, { 3.580623063441472*^9, 3.5806230931907244`*^9}, {3.5806231448268147`*^9, 3.580623172002063*^9}, {3.580623203482918*^9, 3.580623223918954*^9}, 3.5806233379707537`*^9, 3.580623395207255*^9, {3.58062342665691*^9, 3.5806234736909924`*^9}, 3.580623536434303*^9, {3.580623593608403*^9, 3.580623622593254*^9}, 3.580623682107359*^9, 3.580623718595823*^9, 3.5806237837571373`*^9, 3.580623816017994*^9, 3.5806238611020727`*^9, 3.580623922690981*^9}], Cell[BoxData["2"], "Output", CellChangeTimes->{ 3.5806206424420195`*^9, 3.5806213346464357`*^9, 3.5806222658588715`*^9, 3.580622297667327*^9, 3.580622353375025*^9, {3.580622387039884*^9, 3.580622429830759*^9}, {3.5806224688776274`*^9, 3.580622518173714*^9}, { 3.580622550777772*^9, 3.580622580573824*^9}, 3.5806227155608606`*^9, 3.5806228001130095`*^9, 3.580622873417538*^9, 3.5806229089232006`*^9, { 3.580622941901658*^9, 3.580622966674502*^9}, 3.5806230101829786`*^9, { 3.580623063441472*^9, 3.5806230931907244`*^9}, {3.5806231448268147`*^9, 3.580623172002063*^9}, {3.580623203482918*^9, 3.580623223918954*^9}, 3.5806233379707537`*^9, 3.580623395207255*^9, {3.58062342665691*^9, 3.5806234736909924`*^9}, 3.580623536434303*^9, {3.580623593608403*^9, 3.580623622593254*^9}, 3.580623682107359*^9, 3.580623718595823*^9, 3.5806237837571373`*^9, 3.580623816017994*^9, 3.5806238611020727`*^9, 3.580623922690981*^9}] }, Open ]], Cell[BoxData[{ RowBox[{ RowBox[{"Plot1", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"Graphics3D", "[", RowBox[{"{", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{"t", "[", RowBox[{"[", "i", "]"}], "]"}], "}"}], ",", RowBox[{"u1", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}], "}"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Plot1o", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"Graphics3D", "[", RowBox[{"{", RowBox[{ RowBox[{"Opacity", "[", "0.8", "]"}], ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\\""}], ",", RowBox[{"Lighter", "[", "Blue", "]"}], ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\\""}], ",", RowBox[{"Lighter", "[", "White", "]"}], ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\\""}], ",", RowBox[{"Lighter", "[", RowBox[{"Red", ",", "1"}], "]"}], ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\\""}], ",", RowBox[{"Lighter", "[", "Green", "]"}], ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\\""}], ",", RowBox[{"Lighter", "[", "Orange", "]"}], ",", RowBox[{"Darker", "[", "Brown", "]"}]}], "]"}]}], "]"}]}], "]"}]}], "]"}]}], "]"}], ",", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{"t", "[", RowBox[{"[", "i", "]"}], "]"}], "}"}], ",", RowBox[{"u1", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}], "}"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}], "\[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.5682659184078283`*^9, 3.5682659267694435`*^9}, {3.568265972477524*^9, 3.568266002725977*^9}, 3.568277802649102*^9, {3.5804449976745625`*^9, 3.580445016550596*^9}, {3.5804471386732645`*^9, 3.580447144570075*^9}, { 3.5804597731466503`*^9, 3.580459773926652*^9}, {3.580459890942457*^9, 3.5804599070572853`*^9}, {3.5804612382856236`*^9, 3.5804612398300266`*^9}, {3.580461274914488*^9, 3.580461275585289*^9}, { 3.5804613451302114`*^9, 3.580461345504612*^9}, {3.580465262094291*^9, 3.5804652627026925`*^9}, {3.580465328893608*^9, 3.5804653540408525`*^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[{"Distance", " ", "Matrix", " ", "Bohr"}], "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"dc", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"Norm", "[", RowBox[{ RowBox[{"t", "[", RowBox[{"[", "i", "]"}], "]"}], "-", RowBox[{"t", "[", RowBox[{"[", "j", "]"}], "]"}]}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}]}], "Input", CellChangeTimes->{{3.5679440895764875`*^9, 3.5679441000441055`*^9}, { 3.5679441319773617`*^9, 3.5679441523041973`*^9}, {3.567944240210352*^9, 3.5679442813476243`*^9}, {3.567944323966899*^9, 3.567944345931738*^9}, { 3.567944397567828*^9, 3.567944456941533*^9}, {3.567944491557994*^9, 3.5679445468444905`*^9}, {3.5679446505690727`*^9, 3.567944683875131*^9}, { 3.567944736244423*^9, 3.5679448177233667`*^9}, {3.567945549707852*^9, 3.5679455499730525`*^9}, {3.5679456039959474`*^9, 3.567945606195551*^9}, { 3.567945704460124*^9, 3.5679457164877453`*^9}, {3.5679457857362666`*^9, 3.5679458050023003`*^9}, {3.567946353670864*^9, 3.5679464127169676`*^9}, { 3.5679466013368993`*^9, 3.5679466311173515`*^9}, {3.5679467590999765`*^9, 3.567946915739851*^9}, {3.5679469466123056`*^9, 3.567946946815106*^9}, { 3.5679471385082426`*^9, 3.5679471497870626`*^9}, {3.5679483907224417`*^9, 3.5679484036080647`*^9}, {3.567948441937332*^9, 3.5679484421869326`*^9}, { 3.567948496740228*^9, 3.56794849772303*^9}, {3.567948847912445*^9, 3.5679488949933276`*^9}, {3.568011728877214*^9, 3.5680117308213253`*^9}, { 3.568011766325356*^9, 3.568011773436763*^9}, {3.5681173314452896`*^9, 3.5681173370612993`*^9}, {3.568117919800323*^9, 3.5681179678952074`*^9}, { 3.56811809459863*^9, 3.568118105409449*^9}, {3.5681181818963833`*^9, 3.5681182389144835`*^9}, {3.568118274888147*^9, 3.5681182883977704`*^9}, { 3.568118518513775*^9, 3.56811852739019*^9}, {3.5681188250543127`*^9, 3.568118862634779*^9}, {3.568118972396572*^9, 3.5681189862493963`*^9}, { 3.568119020694257*^9, 3.568119047136303*^9}, {3.5681190903015785`*^9, 3.568119136930061*^9}, {3.5681191671005135`*^9, 3.568119202808976*^9}, { 3.5681194788046613`*^9, 3.568119487603077*^9}, {3.5681195353235607`*^9, 3.5681195410331707`*^9}, 3.56811959201406*^9, {3.5681196426985493`*^9, 3.568119659702579*^9}, {3.568119726813897*^9, 3.56811974027672*^9}, { 3.5681199008322024`*^9, 3.568119920691037*^9}, {3.5681199536070952`*^9, 3.568119958521104*^9}, {3.568120000641178*^9, 3.568120067206495*^9}, { 3.5681211528280015`*^9, 3.5681211643408213`*^9}, 3.5681213201694956`*^9, { 3.568121503313817*^9, 3.568121508961027*^9}, {3.568266015439999*^9, 3.5682660461096525`*^9}, 3.568277806081108*^9, {3.5804471990141706`*^9, 3.5804472170010023`*^9}, {3.5804476007616763`*^9, 3.5804476229917154`*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"x", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"y", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"z", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{"i", ",", "3"}], "]"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}]}], "Input", CellChangeTimes->{{3.5678630319943895`*^9, 3.5678630956269007`*^9}, { 3.5678632280087337`*^9, 3.567863228336334*^9}, {3.5678632621727934`*^9, 3.5678633509213495`*^9}, {3.5678644293824434`*^9, 3.567864451908883*^9}, { 3.567864923404111*^9, 3.56786494544695*^9}, {3.5678650143210707`*^9, 3.5678650845367947`*^9}, {3.567869231084096*^9, 3.567869278118179*^9}, { 3.5682778102931156`*^9, 3.5682778156907253`*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"(*", " ", RowBox[{"Distance", " ", "matrix", " ", "Bohr"}], " ", "*)"}], "Text"], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Dij", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"dc", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], "<", "3.0"}], "||", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\\""}], " ", "&&", " ", RowBox[{ RowBox[{"dc", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], "<", "3.24"}]}], ")"}]}], ")"}], ",", RowBox[{"dc", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], ",", "0"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"MatrixForm", "[", RowBox[{"Prepend", "[", RowBox[{"Dij", ",", "num1"}], "]"}], "]"}]}]}]], "Input", CellChangeTimes->{{3.567863535157673*^9, 3.5678635412104836`*^9}, { 3.5678638182825704`*^9, 3.567863827283786*^9}, {3.5679503910983553`*^9, 3.5679504266352177`*^9}, {3.5679504653700857`*^9, 3.5679505048069553`*^9}, 3.567950534915008*^9, {3.567951107623214*^9, 3.567951136046464*^9}, 3.567951302545556*^9, {3.5679562986474605`*^9, 3.567956305355472*^9}, { 3.5679563590975666`*^9, 3.567956370688387*^9}, {3.567956414930065*^9, 3.567956430467692*^9}, {3.567956484740187*^9, 3.567956500667815*^9}, { 3.5679565456270943`*^9, 3.5679565868267665`*^9}, {3.5679569267045636`*^9, 3.5679569448941956`*^9}, {3.567956976562251*^9, 3.567957000711094*^9}, { 3.5679570348595533`*^9, 3.5679570632672033`*^9}, {3.5679573263928657`*^9, 3.5679573286548696`*^9}, {3.5679573611029263`*^9, 3.5679573637393312`*^9}, {3.5679592820903006`*^9, 3.567959317127962*^9}, { 3.567959397655304*^9, 3.5679594004945087`*^9}, {3.5679594704762316`*^9, 3.5679594724574347`*^9}, {3.567959655898157*^9, 3.567959673073787*^9}, 3.5681143121223865`*^9, 3.568121866388855*^9, 3.568220668221574*^9, { 3.568220718968463*^9, 3.568220723742071*^9}, {3.568220862051914*^9, 3.5682208641111183`*^9}, {3.568220906761593*^9, 3.5682209083527956`*^9}, { 3.5682209451844606`*^9, 3.5682209651680956`*^9}, {3.568221147782016*^9, 3.5682211695908546`*^9}, {3.568221237279373*^9, 3.568221253612602*^9}, { 3.568221305420293*^9, 3.568221306247094*^9}, {3.5682215161142635`*^9, 3.568221546206716*^9}, 3.56822158247678*^9, {3.568277819185131*^9, 3.5682778284359474`*^9}, {3.568286391407484*^9, 3.568286398053096*^9}, { 3.5683070136254387`*^9, 3.5683070186174474`*^9}, 3.5804476741754055`*^9}], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"1", "2", "3"}, {"0.`", "1.830599431530426`", "1.830599431530426`"}, {"1.830599431530426`", "0.`", "2.8845071437024505`"}, {"1.830599431530426`", "2.8845071437024505`", "0.`"} }, GridBoxAlignment->{ "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Output", CellChangeTimes->{ 3.58062064247322*^9, 3.5806213346776357`*^9, 3.580622265890071*^9, 3.5806222976985273`*^9, 3.5806223534062247`*^9, {3.580622387071084*^9, 3.5806224298619595`*^9}, {3.5806224689244275`*^9, 3.580622518204914*^9}, { 3.5806225508089714`*^9, 3.5806225806050243`*^9}, 3.580622715592061*^9, 3.58062280014421*^9, 3.580622873448738*^9, 3.5806229089544005`*^9, { 3.5806229419484587`*^9, 3.580622966705702*^9}, 3.580623010214178*^9, { 3.580623063472672*^9, 3.5806230932219243`*^9}, {3.580623144858015*^9, 3.5806231720332627`*^9}, {3.580623203514118*^9, 3.5806232239501543`*^9}, 3.580623338001954*^9, 3.580623395238455*^9, {3.5806234266881104`*^9, 3.580623473722193*^9}, 3.5806235364655027`*^9, {3.580623593639603*^9, 3.5806236226400537`*^9}, 3.5806236821385584`*^9, 3.5806237186270227`*^9, 3.580623783788337*^9, 3.580623816033594*^9, 3.580623861133273*^9, 3.580623922706581*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], " ", ",", RowBox[{"i", "<", RowBox[{"n", "+", "1"}]}], ",", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"j", "=", "1"}], ",", " ", RowBox[{"j", "<", RowBox[{"bc", "+", "1"}]}], ",", RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"Dij", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], "\[NotEqual]", "0"}], " ", "&&", " ", RowBox[{ RowBox[{"u", "[", RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", " ", "\"\\""}]}], ")"}], ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"[", "j", "]"}], "]"}], "\[Equal]", "\"\\""}], ",", RowBox[{ RowBox[{"u2", "[", RowBox[{"[", "i", "]"}], "]"}], "=", "0.59"}], ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"[", "j", "]"}], "]"}], "\[Equal]", "\"\\""}], ",", RowBox[{ RowBox[{"u2", "[", RowBox[{"[", "i", "]"}], "]"}], "=", "0.5"}], ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"[", "j", "]"}], "]"}], "\[Equal]", "\"\\""}], ",", RowBox[{ RowBox[{"u2", "[", RowBox[{"[", "i", "]"}], "]"}], "=", "0.5541"}], ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"[", "j", "]"}], "]"}], "\[Equal]", "\"\\""}], ",", RowBox[{ RowBox[{"u2", "[", RowBox[{"[", "i", "]"}], "]"}], "=", StyleBox["0.52", FontColor->RGBColor[1, 0, 0]]}]}], "]"}]}], "]"}]}], "]"}]}], "]"}]}], "]"}]}], "]"}]}], "]"}], "\[IndentingNewLine]", "u2"}], "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, { 3.568296148255143*^9, 3.5682961775675945`*^9}, {3.568305693114319*^9, 3.5683057123647532`*^9}, {3.5683062281640587`*^9, 3.568306232532067*^9}, { 3.568306348658671*^9, 3.5683063634318967`*^9}, {3.5683064594344654`*^9, 3.5683064674060793`*^9}, {3.5683072458942466`*^9, 3.568307254973462*^9}, { 3.56830729910594*^9, 3.56830731050956*^9}, {3.568307613368492*^9, 3.568307617174899*^9}, 3.580447736825115*^9, 3.5804627951995583`*^9, { 3.5804628372416325`*^9, 3.5804628374600325`*^9}, 3.580476792802576*^9, 3.580622794013399*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{"0.1792`", ",", "0.52`", ",", "0.52`"}], "}"}]], "Output", CellChangeTimes->{ 3.58062064248882*^9, 3.5806213346932354`*^9, 3.580622265890071*^9, 3.580622297714127*^9, 3.5806223534218245`*^9, {3.580622387086684*^9, 3.580622429877559*^9}, {3.5806224689400277`*^9, 3.5806225182205143`*^9}, { 3.5806225508245716`*^9, 3.580622580620624*^9}, 3.5806227156076612`*^9, 3.5806228001598096`*^9, 3.580622873448738*^9, 3.5806229089700003`*^9, { 3.5806229419484587`*^9, 3.580622966721302*^9}, 3.5806230102297783`*^9, { 3.5806230634882717`*^9, 3.5806230932375245`*^9}, {3.5806231448736153`*^9, 3.5806231720488625`*^9}, {3.580623203529718*^9, 3.580623223965754*^9}, 3.5806233380175543`*^9, 3.580623395254055*^9, {3.58062342670371*^9, 3.5806234737377925`*^9}, 3.580623536481103*^9, {3.580623593655203*^9, 3.580623622655654*^9}, 3.5806236821541586`*^9, 3.580623718642623*^9, 3.580623783803937*^9, 3.580623816049194*^9, 3.5806238611488733`*^9, 3.5806239227221813`*^9}] }, Open ]], Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"Fij", "=", RowBox[{"Dij", "-", "u2"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{"MatrixForm", "[", "Fij", "]"}], ";"}], "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Gij", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{"i", "\[GreaterEqual]", "j"}], ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"Fij", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], ">", "0"}], ",", RowBox[{"Fij", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], ",", "0"}], "]"}], ",", "\"\< \>\""}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"(*", RowBox[{ RowBox[{"MatrixForm", "[", "Gij", "]"}], ";"}], "*)"}]}]}], "Input", CellChangeTimes->{{3.567958596890297*^9, 3.567958603270708*^9}, { 3.5679587604565845`*^9, 3.5679587707682023`*^9}, {3.5679588272871017`*^9, 3.567958828254304*^9}, {3.567959580066424*^9, 3.5679595814704266`*^9}, { 3.567959624994503*^9, 3.5679596300177116`*^9}, {3.5679597079086485`*^9, 3.567959902877791*^9}, {3.5679599961659546`*^9, 3.5679600446196404`*^9}, { 3.5679609611992497`*^9, 3.567960974006872*^9}, 3.567961975154231*^9, { 3.5679622637079377`*^9, 3.5679622945179915`*^9}, 3.5681142413606625`*^9, 3.568121878697276*^9, {3.5681219332193723`*^9, 3.56812197189184*^9}, { 3.5682778370003624`*^9, 3.5682778397927675`*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"Hij", "=", RowBox[{"Transpose", "[", "Gij", "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"(*", RowBox[{ RowBox[{"MatrixForm", "[", "Hij", "]"}], ";"}], "*)"}]}]}], "Input", CellChangeTimes->{{3.5679601758938704`*^9, 3.5679601965327067`*^9}, { 3.567960427147912*^9, 3.567960432217921*^9}, {3.5679605764089737`*^9, 3.5679605968294096`*^9}, {3.5679607389456596`*^9, 3.5679607643581038`*^9}, 3.5681142501278777`*^9, 3.5681218877140923`*^9}], Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"Kij", "=", RowBox[{"Hij", "-", "u2"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{"MatrixForm", "[", "Kij", "]"}], ";"}], "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Lij", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{"i", "\[LessEqual]", "j"}], ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"Kij", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], ">", "0"}], ",", RowBox[{"Kij", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], ",", "0"}], "]"}], ",", "\"\< \>\""}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"(*", RowBox[{ RowBox[{"MatrixForm", "[", "Lij", "]"}], ";"}], "*)"}]}]}], "Input", CellChangeTimes->{{3.5679602339571724`*^9, 3.567960257419614*^9}, { 3.567960343406965*^9, 3.567960353500182*^9}, {3.5679604155258913`*^9, 3.567960421422702*^9}, {3.567960854323462*^9, 3.567960874197897*^9}, { 3.5679611596939983`*^9, 3.567961170270817*^9}, {3.5679614000124207`*^9, 3.567961407016833*^9}, {3.5679615066074076`*^9, 3.567961509649413*^9}, 3.5679616124379935`*^9, 3.5679616766477065`*^9, {3.5679617424486217`*^9, 3.5679617456310277`*^9}, 3.567961801276325*^9, {3.5679618582788253`*^9, 3.5679618970760937`*^9}, {3.5679620279135237`*^9, 3.567962030113127*^9}, 3.567962107894864*^9, {3.567962149468937*^9, 3.5679621572533503`*^9}, 3.5679623389156694`*^9, {3.567962645815009*^9, 3.567962648888214*^9}, { 3.567962691647889*^9, 3.567962693020692*^9}, {3.568048161486014*^9, 3.568048162827616*^9}, 3.5680482096900983`*^9, 3.5681142588638926`*^9, 3.568121899632513*^9, 3.5681219958534822`*^9, {3.568277844675576*^9, 3.5682778463603787`*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"Lij", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{"i", "\[LessEqual]", "j"}], ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"Lij", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], ">", "0"}], ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"[", "j", "]"}], "]"}], "\[NotEqual]", "\"\\""}], ",", RowBox[{"0.5", "*", RowBox[{"Lij", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}]}], ",", RowBox[{"Lij", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}]}], "]"}], ",", "0"}], "]"}], ",", "\"\< \>\""}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}], " "}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Qij", "=", "Lij"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"(*", RowBox[{ RowBox[{"MatrixForm", "[", "Lij", "]"}], ";"}], "*)"}]}]}], "Input", CellChangeTimes->{{3.5679627617388124`*^9, 3.567962837461345*^9}, { 3.567963118121438*^9, 3.567963118698639*^9}, {3.5679631520670977`*^9, 3.5679631541731014`*^9}, {3.5679631898191643`*^9, 3.5679632154656096`*^9}, {3.5679635823158536`*^9, 3.56796361468591*^9}, { 3.567963691516045*^9, 3.567963693590849*^9}, {3.567963737957327*^9, 3.567963807299449*^9}, {3.568048114904332*^9, 3.5680481161835346`*^9}, 3.5680482221389203`*^9, {3.5680485003498087`*^9, 3.5680485123618298`*^9}, 3.5681142659931054`*^9, 3.5681219084621286`*^9, {3.568277848107582*^9, 3.5682778496519847`*^9}, {3.568285867917365*^9, 3.568285868634966*^9}, { 3.568297947687104*^9, 3.568297949137906*^9}, 3.5683030673497076`*^9, { 3.5683057620820403`*^9, 3.5683057631428423`*^9}, 3.5804497010089655`*^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[{ RowBox[{"Tval", "=", RowBox[{"Sum", "[", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{"i", "\[LessEqual]", "j"}], ",", RowBox[{"Nij", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], ",", "0"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", "n"}], "}"}]}], "]"}]}], "\n", RowBox[{"Tcor", "=", RowBox[{"Sum", "[", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"[", "i", "]"}], "]"}], "\[NotEqual]", "\"\\""}], ",", RowBox[{"2.25", "/", RowBox[{ RowBox[{"u2", "[", RowBox[{"[", "i", "]"}], "]"}], "^", "2"}]}], ",", "0"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], "\n", RowBox[{"Ttot", "=", RowBox[{"Tval", "+", "Tcor"}]}], " ", RowBox[{"(*", " ", RowBox[{"plus", " ", "\[Pi]", " ", "and", " ", "lone", " ", RowBox[{"pairs", ":", " ", RowBox[{"+", " ", "8.7616"}]}]}], " ", "*)"}], "\n", RowBox[{"TvalH", "=", RowBox[{"Sum", "[", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"i", "\[LessEqual]", "j"}], " ", "&&", " ", RowBox[{ RowBox[{"u", "[", RowBox[{"[", "j", "]"}], "]"}], "\[Equal]", " ", "\"\\""}]}], ")"}], ",", RowBox[{"Nij", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], ",", "0"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", "n"}], "}"}]}], "]"}]}], "\n", RowBox[{"TvalC", "=", RowBox[{"Tval", "-", "TvalH"}]}]}]}]], "Input", CellChangeTimes->{{3.567963992456174*^9, 3.567964003890994*^9}, { 3.5679640974287586`*^9, 3.5679642038209453`*^9}, {3.56796433180357*^9, 3.567964353487608*^9}, {3.567964403641696*^9, 3.567964453046983*^9}, { 3.5679646160984693`*^9, 3.5679646168316708`*^9}, {3.567964688747797*^9, 3.567964762255126*^9}, {3.568277859480002*^9, 3.568277863598409*^9}, { 3.568286113290196*^9, 3.5682861137113967`*^9}, {3.5682870401502237`*^9, 3.5682870902263117`*^9}, {3.568297787256422*^9, 3.5682977875216227`*^9}, 3.568305788461687*^9, 3.580447849972114*^9}], Cell[BoxData["4.176782875969248`"], "Output", CellChangeTimes->{ 3.5806206425200195`*^9, 3.580621334724436*^9, 3.5806222659212713`*^9, 3.580622297745327*^9, 3.580622353453025*^9, {3.580622387102284*^9, 3.580622429908759*^9}, {3.5806224689712276`*^9, 3.5806225182517147`*^9}, { 3.5806225508557715`*^9, 3.580622580636224*^9}, 3.580622715623261*^9, 3.5806228002066097`*^9, 3.5806228734799385`*^9, 3.5806229090012007`*^9, { 3.5806229419796586`*^9, 3.580622966752502*^9}, 3.5806230102609787`*^9, { 3.580623063519472*^9, 3.5806230932687244`*^9}, {3.5806231449048147`*^9, 3.580623172080063*^9}, {3.580623203560918*^9, 3.5806232239813538`*^9}, 3.5806233380331545`*^9, 3.5806233952696548`*^9, {3.58062342673491*^9, 3.5806234737533927`*^9}, 3.580623536512303*^9, {3.580623593670803*^9, 3.580623622671254*^9}, 3.580623682185359*^9, 3.580623718658223*^9, 3.5806237838351374`*^9, 3.580623816080394*^9, 3.5806238611800733`*^9, 3.5806239227533813`*^9}], Cell[BoxData["70.06586814413265`"], "Output", CellChangeTimes->{ 3.5806206425200195`*^9, 3.580621334724436*^9, 3.5806222659212713`*^9, 3.580622297745327*^9, 3.580622353453025*^9, {3.580622387102284*^9, 3.580622429908759*^9}, {3.5806224689712276`*^9, 3.5806225182517147`*^9}, { 3.5806225508557715`*^9, 3.580622580636224*^9}, 3.580622715623261*^9, 3.5806228002066097`*^9, 3.5806228734799385`*^9, 3.5806229090012007`*^9, { 3.5806229419796586`*^9, 3.580622966752502*^9}, 3.5806230102609787`*^9, { 3.580623063519472*^9, 3.5806230932687244`*^9}, {3.5806231449048147`*^9, 3.580623172080063*^9}, {3.580623203560918*^9, 3.5806232239813538`*^9}, 3.5806233380331545`*^9, 3.5806233952696548`*^9, {3.58062342673491*^9, 3.5806234737533927`*^9}, 3.580623536512303*^9, {3.580623593670803*^9, 3.580623622671254*^9}, 3.580623682185359*^9, 3.580623718658223*^9, 3.5806237838351374`*^9, 3.580623816080394*^9, 3.5806238611800733`*^9, 3.5806239227533813`*^9}], Cell[BoxData["74.2426510201019`"], "Output", CellChangeTimes->{ 3.5806206425200195`*^9, 3.580621334724436*^9, 3.5806222659212713`*^9, 3.580622297745327*^9, 3.580622353453025*^9, {3.580622387102284*^9, 3.580622429908759*^9}, {3.5806224689712276`*^9, 3.5806225182517147`*^9}, { 3.5806225508557715`*^9, 3.580622580636224*^9}, 3.580622715623261*^9, 3.5806228002066097`*^9, 3.5806228734799385`*^9, 3.5806229090012007`*^9, { 3.5806229419796586`*^9, 3.580622966752502*^9}, 3.5806230102609787`*^9, { 3.580623063519472*^9, 3.5806230932687244`*^9}, {3.5806231449048147`*^9, 3.580623172080063*^9}, {3.580623203560918*^9, 3.5806232239813538`*^9}, 3.5806233380331545`*^9, 3.5806233952696548`*^9, {3.58062342673491*^9, 3.5806234737533927`*^9}, 3.580623536512303*^9, {3.580623593670803*^9, 3.580623622671254*^9}, 3.580623682185359*^9, 3.580623718658223*^9, 3.5806237838351374`*^9, 3.580623816080394*^9, 3.5806238611800733`*^9, 3.580623922768981*^9}], Cell[BoxData["4.176782875969248`"], "Output", CellChangeTimes->{ 3.5806206425200195`*^9, 3.580621334724436*^9, 3.5806222659212713`*^9, 3.580622297745327*^9, 3.580622353453025*^9, {3.580622387102284*^9, 3.580622429908759*^9}, {3.5806224689712276`*^9, 3.5806225182517147`*^9}, { 3.5806225508557715`*^9, 3.580622580636224*^9}, 3.580622715623261*^9, 3.5806228002066097`*^9, 3.5806228734799385`*^9, 3.5806229090012007`*^9, { 3.5806229419796586`*^9, 3.580622966752502*^9}, 3.5806230102609787`*^9, { 3.580623063519472*^9, 3.5806230932687244`*^9}, {3.5806231449048147`*^9, 3.580623172080063*^9}, {3.580623203560918*^9, 3.5806232239813538`*^9}, 3.5806233380331545`*^9, 3.5806233952696548`*^9, {3.58062342673491*^9, 3.5806234737533927`*^9}, 3.580623536512303*^9, {3.580623593670803*^9, 3.580623622671254*^9}, 3.580623682185359*^9, 3.580623718658223*^9, 3.5806237838351374`*^9, 3.580623816080394*^9, 3.5806238611800733`*^9, 3.580623922768981*^9}], Cell[BoxData["0.`"], "Output", CellChangeTimes->{ 3.5806206425200195`*^9, 3.580621334724436*^9, 3.5806222659212713`*^9, 3.580622297745327*^9, 3.580622353453025*^9, {3.580622387102284*^9, 3.580622429908759*^9}, {3.5806224689712276`*^9, 3.5806225182517147`*^9}, { 3.5806225508557715`*^9, 3.580622580636224*^9}, 3.580622715623261*^9, 3.5806228002066097`*^9, 3.5806228734799385`*^9, 3.5806229090012007`*^9, { 3.5806229419796586`*^9, 3.580622966752502*^9}, 3.5806230102609787`*^9, { 3.580623063519472*^9, 3.5806230932687244`*^9}, {3.5806231449048147`*^9, 3.580623172080063*^9}, {3.580623203560918*^9, 3.5806232239813538`*^9}, 3.5806233380331545`*^9, 3.5806233952696548`*^9, {3.58062342673491*^9, 3.5806234737533927`*^9}, 3.580623536512303*^9, {3.580623593670803*^9, 3.580623622671254*^9}, 3.580623682185359*^9, 3.580623718658223*^9, 3.5806237838351374`*^9, 3.580623816080394*^9, 3.5806238611800733`*^9, 3.580623922768981*^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[{"Eij", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{"i", "\[GreaterEqual]", "j"}], ",", RowBox[{"Dij", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], ",", "0"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", RowBox[{"bc", "-", "1"}]}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", RowBox[{"bc", "-", "1"}]}], "}"}]}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.568045828002715*^9, 3.56804588761042*^9}, 3.5680459828797874`*^9, {3.568046033938677*^9, 3.5680460518787084`*^9}, { 3.568046989190755*^9, 3.5680470284404235`*^9}, {3.568047061028881*^9, 3.568047085567724*^9}, {3.5680471562982483`*^9, 3.568047169199471*^9}, { 3.568047223393966*^9, 3.568047364059413*^9}, {3.568047417286707*^9, 3.5680474225127163`*^9}, {3.5680474583459787`*^9, 3.568047487861231*^9}, { 3.5680475345521126`*^9, 3.5680475468605347`*^9}, {3.5680479168775845`*^9, 3.5680479311516094`*^9}, {3.568048066918648*^9, 3.568048067667449*^9}, { 3.5680483816336*^9, 3.5680483930216208`*^9}, {3.568048424159275*^9, 3.568048450211321*^9}, {3.5680486326224413`*^9, 3.5680486370996494`*^9}, { 3.5680487358166227`*^9, 3.5680487967191296`*^9}, 3.568048834377596*^9, 3.568048984418659*^9, {3.5680490255559316`*^9, 3.5680490790796256`*^9}, { 3.5680491274553103`*^9, 3.5680491280949116`*^9}, {3.5680491782333994`*^9, 3.568049231944294*^9}, {3.568049361533722*^9, 3.5680493691933355`*^9}, { 3.568049422171028*^9, 3.568049427490638*^9}, {3.568049476427923*^9, 3.5680495489368505`*^9}, {3.568049620634577*^9, 3.5680496274361887`*^9}, { 3.5680496655314555`*^9, 3.568049693876705*^9}, 3.568049746573598*^9, { 3.568049846460573*^9, 3.568049932759925*^9}, {3.568050009059659*^9, 3.5680500167348723`*^9}, {3.5680501572755194`*^9, 3.568050179911159*^9}, { 3.5681087021497335`*^9, 3.5681087223049684`*^9}, 3.568114284401138*^9, 3.5681220092695055`*^9, {3.5804479260534477`*^9, 3.580447971730328*^9}, { 3.580449452781329*^9, 3.5804494665249534`*^9}, 3.5804629573930435`*^9, 3.5804651851393557`*^9}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"nc", "=", "0"}], ";"}], " "}], "\n", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "<", RowBox[{"bc", "+", "1"}]}], ",", RowBox[{"i", "++"}], ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\\""}], ",", RowBox[{"nc", "=", RowBox[{"nc", "+", "1"}]}], ",", "nc"}], "]"}]}], "]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"cnt", "=", RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{"{", "nc", "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", "nc", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"no", "=", "0"}], ";"}], " "}], "\n", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "<", RowBox[{"bc", "+", "1"}]}], ",", RowBox[{"i", "++"}], ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\\""}], ",", RowBox[{"no", "=", RowBox[{"no", "+", "1"}]}], ",", "no"}], "]"}]}], "]"}], "\[IndentingNewLine]", "no"}], "Input", CellChangeTimes->{{3.5804503559293156`*^9, 3.5804503572241178`*^9}, { 3.580450419593027*^9, 3.5804504349902544`*^9}, 3.5804645273330007`*^9, 3.580465176465741*^9}], Cell[BoxData["0"], "Output", CellChangeTimes->{ 3.5806206425356197`*^9, 3.5806213347556353`*^9, 3.5806222659524713`*^9, 3.580622297760927*^9, 3.580622353468625*^9, {3.580622387133484*^9, 3.5806224299243593`*^9}, {3.580622468986828*^9, 3.5806225182673144`*^9}, { 3.5806225508713713`*^9, 3.580622580667424*^9}, 3.580622715654461*^9, 3.58062280022221*^9, 3.5806228735111384`*^9, 3.580622909016801*^9, { 3.5806229420108585`*^9, 3.580622966783702*^9}, 3.580623010276579*^9, { 3.5806230635350723`*^9, 3.5806230932843246`*^9}, {3.580623144936015*^9, 3.580623172111263*^9}, {3.580623203576518*^9, 3.580623224012554*^9}, 3.5806233380643544`*^9, 3.5806233953008547`*^9, {3.58062342675051*^9, 3.5806234737845926`*^9}, 3.580623536543503*^9, {3.5806235937020035`*^9, 3.580623622702454*^9}, 3.5806236822165585`*^9, 3.5806237186894226`*^9, 3.5806237838663373`*^9, 3.580623816111594*^9, 3.580623861211273*^9, 3.580623922784581*^9}], Cell[BoxData["1"], "Output", CellChangeTimes->{ 3.5806206425356197`*^9, 3.5806213347556353`*^9, 3.5806222659524713`*^9, 3.580622297760927*^9, 3.580622353468625*^9, {3.580622387133484*^9, 3.5806224299243593`*^9}, {3.580622468986828*^9, 3.5806225182673144`*^9}, { 3.5806225508713713`*^9, 3.580622580667424*^9}, 3.580622715654461*^9, 3.58062280022221*^9, 3.5806228735111384`*^9, 3.580622909016801*^9, { 3.5806229420108585`*^9, 3.580622966783702*^9}, 3.580623010276579*^9, { 3.5806230635350723`*^9, 3.5806230932843246`*^9}, {3.580623144936015*^9, 3.580623172111263*^9}, {3.580623203576518*^9, 3.580623224012554*^9}, 3.5806233380643544`*^9, 3.5806233953008547`*^9, {3.58062342675051*^9, 3.5806234737845926`*^9}, 3.580623536543503*^9, {3.5806235937020035`*^9, 3.580623622702454*^9}, 3.5806236822165585`*^9, 3.5806237186894226`*^9, 3.5806237838663373`*^9, 3.580623816111594*^9, 3.580623861211273*^9, 3.580623922784581*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"j", "=", "1"}], ",", RowBox[{"j", "\[LessEqual]", "nc"}], ",", RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", RowBox[{"bc", "+", "1"}]}], ",", RowBox[{"i", "\[LessEqual]", "n"}], ",", RowBox[{"i", "++"}], ",", RowBox[{"If", "[", " ", RowBox[{ RowBox[{ RowBox[{"Dij", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], "\[NotEqual]", "0"}], " ", ",", RowBox[{ RowBox[{"cnt", "[", RowBox[{"[", "j", "]"}], "]"}], "=", RowBox[{ RowBox[{"cnt", "[", RowBox[{"[", "j", "]"}], "]"}], "+", "1"}]}], ",", "0"}], "]"}]}], "]"}]}], "]"}], ";"}], "\[IndentingNewLine]", "cnt"}], "Input"], Cell[BoxData[ RowBox[{"{", "}"}]], "Output", CellChangeTimes->{ 3.58062064255122*^9, 3.5806213347712355`*^9, 3.580622265968071*^9, 3.5806222977765274`*^9, 3.5806223534842253`*^9, {3.580622387149084*^9, 3.5806224299399595`*^9}, {3.580622469002428*^9, 3.580622518282914*^9}, { 3.5806225508869715`*^9, 3.580622580683024*^9}, 3.580622715654461*^9, 3.5806228002378097`*^9, 3.580622873526738*^9, 3.580622909032401*^9, { 3.5806229420264587`*^9, 3.580622966799302*^9}, 3.5806230102921786`*^9, { 3.580623063550672*^9, 3.5806230932999244`*^9}, {3.5806231449516153`*^9, 3.580623172111263*^9}, {3.5806232035921183`*^9, 3.580623224028154*^9}, 3.580623338079954*^9, 3.580623395316455*^9, {3.58062342676611*^9, 3.580623473800193*^9}, 3.580623536559103*^9, {3.580623593717603*^9, 3.5806236227180543`*^9}, 3.5806236822165585`*^9, 3.580623718705023*^9, 3.580623783881937*^9, 3.580623816127194*^9, 3.5806238612268734`*^9, 3.5806239228001814`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"con", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"{", "i", "}"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";", " ", RowBox[{"conj", "=", RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{"{", "nc", "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "\[LessEqual]", "n"}], ",", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"j", "=", "1"}], ",", RowBox[{"j", "\[LessEqual]", "n"}], ",", RowBox[{"j", "++"}], ",", RowBox[{"If", "[", " ", RowBox[{ RowBox[{ RowBox[{"Dij", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], "\[NotEqual]", "0"}], " ", ",", RowBox[{ RowBox[{"con", "[", RowBox[{"[", "i", "]"}], "]"}], "=", RowBox[{"{", RowBox[{ RowBox[{"con", "[", RowBox[{"[", "i", "]"}], "]"}], ",", "j"}], "}"}]}]}], "]"}]}], "]"}]}], "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"conn", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"Flatten", "[", RowBox[{"con", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ "identify", " ", "C", " ", "atoms", " ", "at", " ", "double", " ", "bonds"}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"j", "=", "0"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "\[LessEqual]", "nc"}], ",", RowBox[{"i", "++"}], ",", RowBox[{"If", "[", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"Length", "[", RowBox[{"conn", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}], "\[Equal]", "4"}], " ", "&&", " ", RowBox[{ RowBox[{"u", "[", RowBox[{"[", RowBox[{"conn", "[", RowBox[{"[", RowBox[{"i", ",", "4"}], "]"}], "]"}], "]"}], "]"}], "\[NotEqual]", "\"\\""}]}], ")"}], ",", RowBox[{ RowBox[{"j", "=", RowBox[{"j", "+", "1"}]}], ";", RowBox[{ RowBox[{"conj", "[", RowBox[{"[", "j", "]"}], "]"}], "=", "i"}]}]}], "]"}]}], "]"}], ";"}], "\[IndentingNewLine]", RowBox[{"conj", "=", RowBox[{"Take", "[", RowBox[{"conj", ",", RowBox[{"{", RowBox[{"1", ",", "j"}], "}"}]}], "]"}]}], "\[IndentingNewLine]", RowBox[{"nj", "=", RowBox[{"Length", "[", "conj", "]"}]}], "\[IndentingNewLine]", RowBox[{"MatrixForm", "[", RowBox[{"conn", ",", " ", RowBox[{"TableAlignments", "\[Rule]", "Left"}]}], "]"}]}], "Input", CellChangeTimes->{{3.5804504793099327`*^9, 3.580450480916735*^9}}], Cell[BoxData[ RowBox[{"{", "}"}]], "Output", CellChangeTimes->{ 3.58062064256682*^9, 3.5806213347712355`*^9, 3.580622265983671*^9, 3.580622297792127*^9, 3.580622353499825*^9, {3.580622387149084*^9, 3.5806224299555593`*^9}, {3.5806224690180283`*^9, 3.5806225182985144`*^9}, {3.5806225509025717`*^9, 3.5806225806986237`*^9}, 3.580622715670061*^9, 3.5806228002534094`*^9, 3.5806228735423384`*^9, 3.580622909048001*^9, {3.580622942042059*^9, 3.5806229668149023`*^9}, 3.5806230103077784`*^9, {3.5806230635662723`*^9, 3.580623093315524*^9}, { 3.5806231449672155`*^9, 3.5806231721268625`*^9}, {3.580623203607718*^9, 3.580623224043754*^9}, 3.5806233380955544`*^9, 3.580623395332055*^9, { 3.5806234267817097`*^9, 3.580623473815793*^9}, 3.580623536574703*^9, { 3.5806235937332034`*^9, 3.5806236227336545`*^9}, 3.5806236822321587`*^9, 3.580623718720623*^9, 3.580623783897537*^9, 3.5806238161427937`*^9, 3.580623861242473*^9, 3.5806239228157816`*^9}], Cell[BoxData["0"], "Output", CellChangeTimes->{ 3.58062064256682*^9, 3.5806213347712355`*^9, 3.580622265983671*^9, 3.580622297792127*^9, 3.580622353499825*^9, {3.580622387149084*^9, 3.5806224299555593`*^9}, {3.5806224690180283`*^9, 3.5806225182985144`*^9}, {3.5806225509025717`*^9, 3.5806225806986237`*^9}, 3.580622715670061*^9, 3.5806228002534094`*^9, 3.5806228735423384`*^9, 3.580622909048001*^9, {3.580622942042059*^9, 3.5806229668149023`*^9}, 3.5806230103077784`*^9, {3.5806230635662723`*^9, 3.580623093315524*^9}, { 3.5806231449672155`*^9, 3.5806231721268625`*^9}, {3.580623203607718*^9, 3.580623224043754*^9}, 3.5806233380955544`*^9, 3.580623395332055*^9, { 3.5806234267817097`*^9, 3.580623473815793*^9}, 3.580623536574703*^9, { 3.5806235937332034`*^9, 3.5806236227336545`*^9}, 3.5806236822321587`*^9, 3.580623718720623*^9, 3.580623783897537*^9, 3.5806238161427937`*^9, 3.580623861242473*^9, 3.5806239228157816`*^9}], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"1", "2", "3"}, {"2", "1", "3"}, {"3", "1", "2"} }, GridBoxAlignment->{ "Columns" -> {{Left}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$, TableAlignments -> Left]]]], "Output", CellChangeTimes->{ 3.58062064256682*^9, 3.5806213347712355`*^9, 3.580622265983671*^9, 3.580622297792127*^9, 3.580622353499825*^9, {3.580622387149084*^9, 3.5806224299555593`*^9}, {3.5806224690180283`*^9, 3.5806225182985144`*^9}, {3.5806225509025717`*^9, 3.5806225806986237`*^9}, 3.580622715670061*^9, 3.5806228002534094`*^9, 3.5806228735423384`*^9, 3.580622909048001*^9, {3.580622942042059*^9, 3.5806229668149023`*^9}, 3.5806230103077784`*^9, {3.5806230635662723`*^9, 3.580623093315524*^9}, { 3.5806231449672155`*^9, 3.5806231721268625`*^9}, {3.580623203607718*^9, 3.580623224043754*^9}, 3.5806233380955544`*^9, 3.580623395332055*^9, { 3.5806234267817097`*^9, 3.580623473815793*^9}, 3.580623536574703*^9, { 3.5806235937332034`*^9, 3.5806236227336545`*^9}, 3.5806236822321587`*^9, 3.580623718720623*^9, 3.580623783897537*^9, 3.5806238161427937`*^9, 3.580623861242473*^9, 3.5806239228157816`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"j", "=", "0"}], ";", " ", RowBox[{"h", "=", "0"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"coco", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0", ",", "0", ",", "0"}], "}"}], ",", RowBox[{"{", "no", "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"coch", "=", "coco"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", RowBox[{"nc", "+", "1"}]}], ",", " ", RowBox[{"i", "\[LessEqual]", RowBox[{"nc", "+", "no"}]}], ",", " ", RowBox[{"i", "++"}], ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"Length", "[", RowBox[{"conn", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}], "\[Equal]", "2"}], ",", RowBox[{ RowBox[{"j", "=", RowBox[{"j", "+", "1"}]}], ";", RowBox[{ RowBox[{"coco", "[", RowBox[{"[", "j", "]"}], "]"}], "=", RowBox[{"Join", "[", RowBox[{ RowBox[{"{", RowBox[{"conn", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "}"}], ",", RowBox[{"Take", "[", RowBox[{ RowBox[{"conn", "[", RowBox[{"[", RowBox[{"conn", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}], ",", "3"}], "]"}]}], "]"}]}]}], ",", RowBox[{ RowBox[{"h", "=", RowBox[{"h", "+", "1"}]}], ";", RowBox[{ RowBox[{"coch", "[", RowBox[{"[", "h", "]"}], "]"}], "=", RowBox[{"conn", "[", RowBox[{"[", "i", "]"}], "]"}]}]}]}], "]"}]}], "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"coco", "=", RowBox[{"Take", "[", RowBox[{"coco", ",", RowBox[{"{", RowBox[{"1", ",", "j"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"coch", "=", RowBox[{"Take", "[", RowBox[{"coch", ",", RowBox[{"{", RowBox[{"1", ",", "h"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"oj", "=", RowBox[{"Length", "[", "coco", "]"}]}], "\[IndentingNewLine]", RowBox[{"hj", "=", RowBox[{"Length", "[", "coch", "]"}]}]}], "Input", CellChangeTimes->{{3.5804505432700443`*^9, 3.5804506082285585`*^9}, { 3.5804506582734466`*^9, 3.5804506742478747`*^9}, 3.5804507096443367`*^9, { 3.5804508095313125`*^9, 3.5804508306849494`*^9}, 3.5804510247492905`*^9, 3.5804511224366617`*^9, {3.580460502510331*^9, 3.580460514147952*^9}, { 3.58046055297642*^9, 3.5804605556908245`*^9}, 3.580460689991461*^9, { 3.5804610478404894`*^9, 3.5804610672625237`*^9}, {3.5804611544354763`*^9, 3.5804611697547035`*^9}, {3.580464611183148*^9, 3.5804646160191565`*^9}}], Cell[BoxData["0"], "Output", CellChangeTimes->{ 3.58062064259802*^9, 3.5806213347868357`*^9, 3.5806222659992714`*^9, 3.580622297807727*^9, 3.5806223535154247`*^9, {3.5806223871646843`*^9, 3.580622429971159*^9}, {3.580622469033628*^9, 3.5806225183141146`*^9}, { 3.580622550918172*^9, 3.580622580714224*^9}, 3.5806227156856613`*^9, 3.5806228002690096`*^9, 3.5806228735579386`*^9, 3.5806229090636005`*^9, { 3.580622942057659*^9, 3.5806229668305025`*^9}, 3.5806230103233786`*^9, { 3.580623063581872*^9, 3.5806230933311243`*^9}, {3.5806231449828153`*^9, 3.5806231721424627`*^9}, {3.580623203638918*^9, 3.580623224059354*^9}, 3.5806233381111546`*^9, 3.580623395347655*^9, {3.58062342679731*^9, 3.580623473831393*^9}, 3.580623536590303*^9, {3.580623593748803*^9, 3.580623622749254*^9}, 3.580623682263359*^9, 3.5806237187362227`*^9, 3.5806237839131374`*^9, 3.580623816158394*^9, 3.580623861258073*^9, 3.5806239228313813`*^9}], Cell[BoxData["1"], "Output", CellChangeTimes->{ 3.58062064259802*^9, 3.5806213347868357`*^9, 3.5806222659992714`*^9, 3.580622297807727*^9, 3.5806223535154247`*^9, {3.5806223871646843`*^9, 3.580622429971159*^9}, {3.580622469033628*^9, 3.5806225183141146`*^9}, { 3.580622550918172*^9, 3.580622580714224*^9}, 3.5806227156856613`*^9, 3.5806228002690096`*^9, 3.5806228735579386`*^9, 3.5806229090636005`*^9, { 3.580622942057659*^9, 3.5806229668305025`*^9}, 3.5806230103233786`*^9, { 3.580623063581872*^9, 3.5806230933311243`*^9}, {3.5806231449828153`*^9, 3.5806231721424627`*^9}, {3.580623203638918*^9, 3.580623224059354*^9}, 3.5806233381111546`*^9, 3.580623395347655*^9, {3.58062342679731*^9, 3.580623473831393*^9}, 3.580623536590303*^9, {3.580623593748803*^9, 3.580623622749254*^9}, 3.580623682263359*^9, 3.5806237187362227`*^9, 3.5806237839131374`*^9, 3.580623816158394*^9, 3.580623861258073*^9, 3.5806239228313813`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"(*", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{"--", RowBox[{ "--", "-"}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}\ ]}]}]}]}]}]}]}]}]}]}]}], "*)"}], "Text"], "\[IndentingNewLine]", StyleBox[ RowBox[{"(*", RowBox[{ RowBox[{"next", ":", RowBox[{ RowBox[{"O", " ", "atoms", " ", "in", " ", "C"}], "-", "O", "-", "H"}]}], ",", RowBox[{"C", "-", "O", "-", "C"}], ",", RowBox[{"O", "-", "O", "-", RowBox[{"H", ":", RowBox[{"two", " ", "doubly", " ", "occ", " ", "LP"}]}]}]}], "*)"}], "Text"], "\n", RowBox[{ RowBox[{ RowBox[{"tpso", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", RowBox[{"{", "hj", "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"tpsu", "=", "tpso"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"P", "=", "3.0"}], ";"}], "\[IndentingNewLine]", StyleBox[ RowBox[{ RowBox[{"rlos", "=", RowBox[{"{", "1.125", "}"}]}], ";"}], "Input"], "\[IndentingNewLine]", "coch", "\[IndentingNewLine]", RowBox[{ RowBox[{"Ok", "=", RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{"-", "1"}], ",", RowBox[{"-", "1"}], ",", "1"}], "}"}], "}"}]}], ";"}]}]}]], "Input", CellChangeTimes->{{3.5804505432700443`*^9, 3.5804506082285585`*^9}, { 3.5804506582734466`*^9, 3.5804506742478747`*^9}, 3.5804507096443367`*^9, { 3.5804508095313125`*^9, 3.5804508306849494`*^9}, 3.5804510247492905`*^9, 3.5804511224366617`*^9, {3.580460502510331*^9, 3.580460514147952*^9}, { 3.58046055297642*^9, 3.5804605556908245`*^9}, 3.580460689991461*^9, { 3.5804610478404894`*^9, 3.5804610672625237`*^9}, {3.5804611544354763`*^9, 3.5804611697547035`*^9}, 3.580464611183148*^9, {3.5804885900917444`*^9, 3.5804886187021947`*^9}, {3.5804887019907413`*^9, 3.580488708293152*^9}, { 3.580489378080329*^9, 3.5804893866759434`*^9}, {3.580525866734337*^9, 3.5805259166232247`*^9}, {3.5806227077140474`*^9, 3.5806227086032486`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{"1", ",", "2", ",", "3"}], "}"}], "}"}]], "Output", CellChangeTimes->{ 3.58062064261362*^9, 3.580621334802436*^9, 3.580622266014871*^9, 3.580622297823327*^9, 3.580622353531025*^9, {3.5806223871802845`*^9, 3.580622429986759*^9}, {3.5806224690492277`*^9, 3.5806225183297143`*^9}, { 3.580622550933772*^9, 3.580622580729824*^9}, 3.5806227157012615`*^9, 3.58062280028461*^9, 3.5806228735735383`*^9, 3.5806229090792007`*^9, { 3.580622942073259*^9, 3.5806229668461027`*^9}, 3.580623010338979*^9, { 3.580623063597472*^9, 3.5806230933467245`*^9}, {3.580623144998415*^9, 3.580623172158063*^9}, {3.580623203638918*^9, 3.580623224074954*^9}, 3.5806233381267543`*^9, 3.5806233953632545`*^9, {3.58062342681291*^9, 3.5806234738469925`*^9}, 3.5806235366059027`*^9, {3.5806235937644033`*^9, 3.5806236227648544`*^9}, 3.580623682278959*^9, 3.5806237187518225`*^9, 3.5806237839287376`*^9, 3.580623816173994*^9, 3.580623861273673*^9, 3.580623922846981*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"<<", " ", "\"\\""}]], "Input", CellChangeTimes->{{3.580460580931669*^9, 3.5804605820704713`*^9}, 3.580460952617922*^9, {3.5805326737704744`*^9, 3.5805326831460915`*^9}, { 3.5805328566027956`*^9, 3.5805328613452044`*^9}}], Cell[BoxData[ RowBox[{"{", Graphics3DBox[{ {RGBColor[1, NCache[ Rational[1, 3], 0.3333333333333333], NCache[ Rational[1, 3], 0.3333333333333333]], Opacity[0.5], SphereBox[{-0.425494998113922, -0.5443315318627366, 1.153125}, 1.125]}, {RGBColor[1, NCache[ Rational[1, 3], 0.3333333333333333], NCache[ Rational[1, 3], 0.3333333333333333]], Opacity[0.5], SphereBox[{-0.425494998113922, -0.5443315318627366, -1.153125}, 1.125]}, {Opacity[0.5], SphereBox[{1.830599431530426, 2.220446049250313*^-16, 0.}, 0.52]}, {Opacity[0.5], SphereBox[{0., 0., 0.}, 0.1792]}, {Opacity[0.5], SphereBox[{-0.44198443328439774`, 1.7764413976975961`, 0.}, 0.52]}, Text3DBox["2", {1.830599431530426, 2.220446049250313*^-16, 0.}, {1, 0}], Text3DBox["1", {0., 0., 0.}, {1, 0}], Text3DBox["3", {-0.44198443328439774`, 1.7764413976975961`, 0.}, {1, 0}], TubeBox[{{1.830599431530426, 2.220446049250313*^-16, 0.}, {0., 0., 0.}}, 0.02], TubeBox[{{1.830599431530426, 2.220446049250313*^-16, 0.}, {-0.44198443328439774`, 1.7764413976975961`, 0.}}, 0.02], TubeBox[{{0., 0., 0.}, {-0.44198443328439774`, 1.7764413976975961`, 0.}}, 0.02]}, AspectRatio->Automatic, Axes->Automatic, AxesLabel->{ FormBox["xl", TraditionalForm], FormBox["yl", TraditionalForm], FormBox["zl", TraditionalForm]}, Boxed->False, PlotRange->Automatic, SphericalRegion->True, ImageCache->GraphicsData["CompressedBitmap", "\<\ eJzVWnlQm2d6/y1gO4k3sTf3JtnUSZo23e1mpu3sttkk22M7zXR6bo+0/3Uy 7nay06mTtZ1pprvdXTfJxuBNsvbGHQE2NtgcBmNAF5IAgYyEEEhISEKWxGEO cwiJ04BtLvV53+99xSdZ0Oy0/aOaefn0HXqO33O/H3/2xtE3v/Mvbxz97sE3 Drx2+I233/zuwSMH/uitw3Qp93MAymjlHwD7nqSv4k8Zkskku/m/fzjIDnmM y8vKhX/kh1xwxknOX976TsYttr6h3PqnLLde+m9vcebIS7uVIUkOu/CKcuFt 8fjt27dhs9nw+uuvo7y8/EXl7iF+2INYLIb8/Hy89tprMBqN7AnOgH3/qvLo n7IDl+gpfmFvGkn2Q6/Xi6WlJX53YWEB4+Pj/B679usKka+lxCPlJFLKGX32 Z9H6cSFhJjMpJbuWJ3SU9w4ePMiYfkXh8PUsxsnJwmlflmuP8L+7OXeJEfsw ziMjI5IM01ZqniuEkefsURJWCPPb/LArjYkU6gtZBHiI/70nK9YNDQ1S1ftU grBH2O0v/8/V369SX3JnHwYDcc0VP2EqBoNBfnavUJ9JwB4XAu5Nec+Td3nC rjRP2AmN+1VoZLosw5upf/Lkyf0qGchK0gSfz0LxCYXfy1luPZCS8P4sMDC2 7ENsH80Cfl4qReRuz/UbKQ5PCeyyebl0JGZu9v1RVUTX1NRkukAWPl9P8fmS KpqYmGoAZag+KegzOd59912Onkaj2YH+17KnpCdUSjEijAEjyEyizhPs+9OC usRPsfMvyvGxHZIFg49xf05QZf4joKNfpBtiO/oPZziCTAeMrioepHJs/ZII eBkPTBxyUc76+R2ZPahy9kz4ZGolqs8LmZg87JakIG7/8o489m3j2ewYDoel QsxKDEN29oLq2rFjx7iyz2X62f0Z/pwZqGrTv6gyvczeguC3VPZUQ6DKfPjN LPmWWD6rUHgmlXee2iaXsefptzIlqiH8rYzAZiyeUeg+m8U3H98mV0ulmdTN zc0vCYEzqtUBhfCTOyaOR1VKqAERNF7NggVxP7C9zA99hsQj6/rvZoFjV8rk WzG0b3t2+3eASGXUbwk9mGoqiGRpzXDiPdtWiqwyPCAikiUABtChQ4f4fVbB JOZMFnaPnf2B4Co9Q4C6jTDZkhY/u1flkKoY2d5d2DVyF6WPuw8z4wkMuyO4 4Y4iqnUioNEjUqjHULEe00YnlnqjWPZH8jKMxIyXs5Oo2XsRJduL3i97QmXB z/KQ4nSPIUYCXiPhdCRYs8YAu6YR3RoT/IVmhIstGCppwliZBVMVFiQumTB3 uRGLV/S41e7EajjCqWQYPCVOtoZlV8a17Z5T1MRXVGZn6UjGPzOruo+TP3tZ xLRV24FLpFK9xgiTxoI2TQs6NW3wFl1F6Ew7BkrtGK1wYLLajsQVO+YbrmJJ 34YVYwvumMxYNRmw1tWBjZlEyn2onEqRs0mcl91eOyj3wg52YjZiPq6E7xPc TlUaHS0DqjVNZK82sped7OWEu8iFwLluRC56cL2qB+O1PZhu6MGcwYObZjdu NXdh1erERpsdG7ZWbNrM2GjVY3M2wYFlH1G+uY/vTn1TvkL9EefC1Z7fRgVZ LJk7/D5/5kmVCiayTStqSXyjxoVWjQfOYi96zvnRVx7AQHUfRutCmNKHMGMK YbG5DyutQay2+7Hu8CHp9CDpctGy4/bVRuS/+Q94/S//Qp1Ss/bFAvXnsriU mAyQrZ78nki2iuhmWjayQCe5lgdmTS98uuu44Y3j5uSyYEHYjMeRjA4j6fPT 8mDB3gRvZSEW3Vb+BIOKRN1qVgWcB7bJLsw6zCNE9P6hQD1ImYWhWaW5SqsL NRovIu7YliPeJxWPTyN5jSwccpI8Fn6V8RCW/6LC/UtZii3jzrimcqcSb93d 3fxM5nzmllXkkkwKq3YgLwN+9tkrZVmhDDEeRTJsx1RPN/dymZnIdI+rKio9 nWV+43eYEPT9jwWjKorxKjIK8RY+m277zwvVkhNhxPy9qevCEILpPSl7ZJZx JZJlOG99y5jRCa0cbqxcHgGSy5+kpLSRhFFhvbul3CvuTA9NIBYKpkkpg1Tp llWFWzWJqY8K5Xu2i19+LCkpkWM9YyHZ/bkU96QldU0trv3CECr+1ZKyqFUb wnV/NPWszF+i3VDaHwWve7GFrvzkptDMy9BHfliRPyAQinsjGDM5MHBei+ul WoyWaXHjghaTtGIXtUiUayljOHA7EBHZKRcBQ1NGmt7N6XZUjCD/nxVs3zs5 Btv1dTkwMTNVUYIKugcgvVk4yxYZLvkDQtx0pbbdLvnM1+6GYWvk2oeJnjB8 JQ0IntPiGkERJSgGaQ0THGMExThBMUVrmuBIVGgxS2u+sgF3AmGRASm2HATL 5JbdrnWH4awax08PtcOkvy1VqCoNocaVFGlQUZ7FGYu3Kk0zQRSVbRWLVlmi yQsyiz1bosfjQbp3B+VlYWIfmdWVBPUg+gwOdJ1tQE+JFr0EQB8BECbl+0n5 IVrDBMAYKT/BQCAApisECJUKCItVDViLxZ8RrjBp1SNcV4fOsgF8fNgBU+Mq AuNJmQMC3bM4cX5WVUN2Uy1bEhh0ibxnZPUtV9QWOSWpm+7MnJgJhgL3nh0D dnuvymOTH8+lYjRSsHoYPr0D9jMN6CS83ISXl/DyE14hwitSxjDTYeiCDiMX dRgjrCZoTVUIzCoVzOYIrwUFs6dFNoh098NVG0fV9y/DVBeH3zdK1YUqzHQM W8CFcfQnNzC3KvPHHmqIl2DXDlEF9RJoLp4SlaqqE4lbyajUQ6cPENsDp+q5 UsGXwjMvDSuZY54RRWGG6nRLcT1sBJCDAHIRQG4CyJp/Bq35Z2ErOIurbJ04 i/aCszlopy97ME7YTNKKETZxWjOEzSytuUv1onDtQ7R7EF31szj1rhumlnX0 TiaxPp/Aposm8h5qA7wdtNz0napR5Do2b8Tx0adLmFqSOWovlmPLmPLHcaNn GrbCALUbPtRp3NQ9dQivU4BjnZWEgtWILEV7O/B2fcaclLMF2270e8KwEGxW gq3lVAUM/34KdUc+zEH9keOUD+uO5OMKrdrD+bhMq5pW1ffyUf72cejpUTuB OEWwTRNkcVqJqnqsxuKPieCpp1p+7kdd+PQHPhhbN+Cbkog8ho3hCDYdBiQ7 m2m1Uy9IDY3DjfcPueDpmMJSsx9zpgBi+iBu1AUxWB1A6KIf3hIfnIU91Gt2 U8/ppN6znZBrodWYQvBeVX+TPoT9YhksDaoEeZixqJ70Pgnd0Q+hP0oYHM3P geFoPkGlP1pAlwugPVJA4BWg7nABh41BVikgK33rOPqKKjBFUMUIquVg+GEB VQOl3/P/4cbpHwehI6g8UzJFPYKN6xFsUGbbbDVho9WKimMm1DbNY9nUiQWD C4mGLkzUdmGkqgvRC10InnXBXeikicJOk4UNWoKnmjuYPgOiXWk9j0zO6l2E 3alIy7urX8qIVTns6wggAwHU+M5xetT0Tj77U0AQNdLfRoLISEuvguoKQXX5 ewW4JGAqO3Sc4vQMxgmiO7H4F4SoOo0VpR948Z/vhRHwJhG7KbvBB/l0t+py YvRMLTmlHku6JizUt2CmthVTl9ooH9oweN6GUHEbfIWtNEG20CRpoYnSRHFo oDjUsR58r2AlO1nV7rXYEFe1a1t7L7mpb88KAvPBawgVlcNETmJ+Jx8WWjn0 t4D9OUGOZ6a/pqMnCI8TBNYJ6I6cQMPhE9xtagiLqrcJC3KXUsKit7Bc5bx5 GKJhYYpWgtbSeGKP4KknU1847oPmJwMI9iZ5/tkvcligO4IP3gtgsS+KuWo9 BaqBmgojbpQ2YuhsIyJFjfBT5eumZSc8mslVou4oQsRCvhZSD5jKtdTrlIS1 nVsrWlxB+h1HE1e3mVu9mTRtEstyVKX1EUXrehYopHE1aVz5Vj4ukMa1n5Tn iDLDPkaSqKewEcPnTFTJjKSBDmvTCcgibSxsR/kJPzT5w9RhKpo/JKT+pHAF seWkKLpKKN+eSGDeE8EMaTdNa5LWmDvysPBrNqWoRtB9al13Y85K0w7pOklB PFBciY5jPycd87m6OVzbL3I9zWo9SccGWldIzxrSs/bfTsHwM0q3lE9YXul1 Bbii9qJu9Ff1Ysbox522Hmx22DF+SQtTkQOVHwdQ+NEYgqEkJpflu7FcfPBD Cocp9W53LvdXucnxgFr8A1gJhnmajvMc1MDT90BRFXo/LkEzZbJmCswmWhZa ZgrORloGWooK+TD84OdoPVUFGxXZluIGmIUK5R/VcBW6K0cxapnCTVcMA6Uu BMucaD3TA5cthonl5GNCwksXo9D2JoVQyjX51pM5mJD6d/jhaV6PWV1m9ZnV aVavWd1m9XusXAfPT0vgptV1ogSdBSVw0PJTX+Q9r7QBndQO2EniNqpzzSQx o1qZb+HsPDVzGDGtwHZhHsNty4jaYwg29iJ2LSylDXTPo6T+Dj97RlWe5SRP XrL1xmlXuujUavGWi7VerAVjrRhryVhrxlo01t6ylo21bqztZe1vX6nS2vWQ 6N0kupNEbyfR5ybinGpLmU/JvqFBrHluQN92CyOtw+isHELAFsRE/7VHU6L3 4/sfT2JhdSP17oUlfZoYuYNn28YTGrwkkivTYNnq4A0310A0lawhZ40504A1 6oOsCS1VppjgeaWh99Bk4yLp5yfiql5uN9ZtdiyZ3Wg7F0Y0NEEhfhXmyy0S 8J+dvYPpW4p7PCcAZ20xC0tR0bONJKJHfFlkPTYiLbU4MFehTJNshGKjFBup Ri8oIxYbtSKsiyZhA7QGjA4sEsx5ss21dNBvdRg/b6CMbsD59/Wo/rSGF1GJ 8Yc/7kff1Co/+xVV7yHfz2QTVXQa3+SHe3GLhtwVfwRzTQ4+ALNBeITEmzQ7 MOOjXE9LbhjfJNYs/4/QsrVEEJ+dTYlyqWIAxqDC4lfFNWZquau7gyh/LU7a B+jysuJoMosGuhdxXq9o+IJIhDMzM1w7+W8AwnuyNVoZHELjSczNxIR7SQ4D +OEp5dqXVc0Je8XKPnI7emtfE/gb8duuEWKzMMHP5HvrU6VrSNxW2P+akFh2 1BmvO+5JpxaNJbE4PcrPZO07fmwIY7MKta+qokg2CoLE34pbvvENbCaGVPk3 F9WVgzCHFRIvCm9lsLHtJTFZ7kmnMpQgQcYiqlaVwXQTF8xr/Ow3VK7GAkNR kpN4XSIdW8edsTA/25siMYhjp6f5mXzFpn5luiudxAipfWt1NY3E6cp1zNxR b9Tj78StevNMipLsw/PfG8H4vHLtlfTHze13cGf1Fj+T725+9I4ZLVHl8VfV wtyPaXL9fupLDLYFOKlH8VJnECo0YPCMAaPnDJgs0yNeocdt6nN8pjGUW9fF xrsw736+1cravXr6qYXav7Oafow6/AgWWzFwzkqJwUr13YqZyy3UQjah+2Q1 PtDEcNugx/pA5JtqV36Iv55gfXUN9dcGalBbKY+5Ch2oO91DaaUTwxUuTNRQ j15PvbreiRVTB86+a0Wiy4eN1iZaRqxbddiYTUho/oofHhGkG0lYK65QK2/W dKK90A33GS/0H/nQf9GFsSt+xHRsUvJjubkX1vevYr7FQ42Di6YqB5KuNjqa se4g4ecS0ljfFoZWtp0b+exeQ+S1mp7UzBgPJmB1J0TQKeY8ddiIvh4KtBvk Pp4Ajb1eJHu7aNEA523CWsh5N4v+1H5ym3aAusppkQ6Unbt2MvTmTcUbA54V VF3dkG1dkjnl6CCSYQ+S1xzY6DWlU6eW4ZNqTjnonhJ+u/U+ze6dIF8fwenS 2dROcnLlJpJDRC5qyyRVcfwiJfXFrXyVo85mncNJOJwjmFvb8tNYRxPGfS5+ 9lI6saaiy0jGIqqidPd7vsmbKlLEWFfShvGAW/572rfFj73VtUj2kyknomzH Pdvmk6TCUFD2sVoyzbDBdgcDNGqHOgjLHiSHBxUich9sdnyJbBMjN/CptnOM aWQe5fse6100tfcQKa+dFhnfQ+S6/dgMXcf6yDRWr8dxayCOpWgc89fiSATj mPJNw1kcQAtR12rcNI46iXordz7y85fVnv8wVn0d2CSPTTotSHaQBzvs2Gzv xFpbN1aaPVgweRHXUftV68VQpReh0h54it1oJ7nNmg6Klqt8P5Pt5RH1V9Qh +yB/WbbWosN6ixFrzRasWqxYabRhUd+OmTo7JmscGCl3IHrOAX+RHS6i1kqy GjRNhI6yy0BUX1Un631Y76cCrqXJs8GAhSuNmKk28Veto2UWDJyxoK/QAg8l CfZa1kKpp56SRitFOOsjRE75exFlK9QLsNfJi14a1jxRjGidVO+j/NUzewXN XkUrr+L+T/698v/1AZ/7L/RjegM=\ \>"]], "}"}]], "Output", CellChangeTimes->{ 3.58062064262922*^9, 3.580621334818036*^9, 3.5806222660304713`*^9, 3.5806222978545275`*^9, 3.580622353546625*^9, {3.580622387195884*^9, 3.5806224300023594`*^9}, {3.580622469064828*^9, 3.5806225183609147`*^9}, { 3.580622550949372*^9, 3.5806225807454243`*^9}, 3.5806227157168612`*^9, 3.5806228003158097`*^9, 3.5806228735891385`*^9, 3.580622909094801*^9, { 3.5806229420888586`*^9, 3.5806229668617024`*^9}, 3.5806230103701787`*^9, { 3.5806230636130724`*^9, 3.5806230933623247`*^9}, {3.580623145014015*^9, 3.5806231721892633`*^9}, {3.5806232036701183`*^9, 3.580623224090554*^9}, 3.580623338157954*^9, 3.5806233953788548`*^9, {3.5806234268441105`*^9, 3.580623473878193*^9}, 3.580623536621503*^9, {3.5806235937956038`*^9, 3.5806236227960544`*^9}, 3.5806236822945585`*^9, 3.5806237187674227`*^9, 3.5806237839443374`*^9, 3.5806238161895943`*^9, 3.5806238612892733`*^9, 3.5806239228625813`*^9}] }, Open ]], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{"sigma", " ", "bond", " ", "cloud", " ", "radii"}], " ", "*)"}]], "Input", CellChangeTimes->{{3.568045828002715*^9, 3.56804588761042*^9}, 3.5680459828797874`*^9, {3.568046033938677*^9, 3.5680460518787084`*^9}, { 3.568046989190755*^9, 3.5680470284404235`*^9}, {3.568047061028881*^9, 3.568047085567724*^9}, {3.5680471562982483`*^9, 3.568047169199471*^9}, { 3.568047223393966*^9, 3.568047364059413*^9}, {3.568047417286707*^9, 3.5680474225127163`*^9}, {3.5680474583459787`*^9, 3.568047487861231*^9}, { 3.5680475345521126`*^9, 3.5680475468605347`*^9}, {3.5680479168775845`*^9, 3.5680479311516094`*^9}, {3.568048066918648*^9, 3.568048067667449*^9}, { 3.5680483816336*^9, 3.5680483930216208`*^9}, {3.568048424159275*^9, 3.568048450211321*^9}, {3.5680486326224413`*^9, 3.5680486370996494`*^9}, { 3.5680487358166227`*^9, 3.5680487967191296`*^9}, 3.568048834377596*^9, 3.568048984418659*^9, {3.5680490255559316`*^9, 3.5680490790796256`*^9}, { 3.5680491274553103`*^9, 3.5680491280949116`*^9}, {3.5680491782333994`*^9, 3.568049231944294*^9}, {3.568049361533722*^9, 3.5680493691933355`*^9}, { 3.568049422171028*^9, 3.568049427490638*^9}, {3.568049476427923*^9, 3.5680495489368505`*^9}, {3.568049620634577*^9, 3.5680496274361887`*^9}, { 3.5680496655314555`*^9, 3.568049693876705*^9}, 3.568049746573598*^9, { 3.568049846460573*^9, 3.568049932759925*^9}, {3.568050009059659*^9, 3.5680500167348723`*^9}, {3.5680501572755194`*^9, 3.568050179911159*^9}, { 3.5681087021497335`*^9, 3.5681087223049684`*^9}, {3.5681133837179556`*^9, 3.5681133888191648`*^9}}], Cell[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, 3.5804489138003826`*^9}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{"a", "=", RowBox[{ "ReadList", "[", "\"\\"", "]"}]}], "\[IndentingNewLine]", RowBox[{"nb", "=", RowBox[{ RowBox[{"Length", "[", "a", "]"}], "-", "1"}]}], "\[IndentingNewLine]", RowBox[{"cw", "=", RowBox[{"nb", "-", "bh"}]}]}], "Input", CellChangeTimes->{{3.5681343270045934`*^9, 3.568134331481801*^9}, 3.5681343823222904`*^9, {3.5681345934218616`*^9, 3.5681345986322703`*^9}, { 3.568134635479535*^9, 3.5681346521559644`*^9}, {3.5681351451792307`*^9, 3.568135145428831*^9}, {3.5681351781888885`*^9, 3.5681352349417877`*^9}, { 3.568135340163973*^9, 3.568135340335573*^9}, {3.568135523495495*^9, 3.5681355244314966`*^9}, {3.5681357674799232`*^9, 3.5681357771675406`*^9}, {3.568136284433632*^9, 3.5681362987700567`*^9}, { 3.568136395053426*^9, 3.5681363965510283`*^9}, {3.568303192134327*^9, 3.568303196237134*^9}, {3.568307596598462*^9, 3.5683075983924656`*^9}, 3.5804489330664167`*^9, {3.580526129501199*^9, 3.580526130624401*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1", ",", "1.1313994315304259`"}], "}"}], ",", RowBox[{"{", RowBox[{"3", ",", "1", ",", "1.1313994315304259`"}], "}"}], ",", RowBox[{"{", RowBox[{"3", ",", "2", ",", "1.8445071437024505`"}], "}"}]}], "}"}]], "Output", CellChangeTimes->{ 3.5806206427228203`*^9, 3.5806213348492355`*^9, 3.5806222660616713`*^9, 3.580622297870127*^9, 3.5806223535778255`*^9, {3.580622387227084*^9, 3.5806224300335593`*^9}, {3.5806224690960283`*^9, 3.5806225183765144`*^9}, {3.5806225509805717`*^9, 3.580622580776624*^9}, 3.580622715810461*^9, 3.5806228003314095`*^9, 3.5806228736203384`*^9, 3.580622909126001*^9, {3.580622942120059*^9, 3.5806229668929024`*^9}, 3.580623010385779*^9, {3.5806230636442723`*^9, 3.5806230933935246`*^9}, { 3.580623145045215*^9, 3.580623172204863*^9}, {3.5806232036857185`*^9, 3.580623224121754*^9}, 3.5806233381735544`*^9, 3.580623395410055*^9, { 3.5806234268597107`*^9, 3.580623473893793*^9}, 3.5806235366527033`*^9, { 3.580623593811204*^9, 3.580623622811654*^9}, 3.5806236823101587`*^9, 3.580623718798623*^9, 3.5806237839755373`*^9, 3.5806238162207937`*^9, 3.580623861320473*^9, 3.5806239228937817`*^9}], Cell[BoxData["2"], "Output", CellChangeTimes->{ 3.5806206427228203`*^9, 3.5806213348492355`*^9, 3.5806222660616713`*^9, 3.580622297870127*^9, 3.5806223535778255`*^9, {3.580622387227084*^9, 3.5806224300335593`*^9}, {3.5806224690960283`*^9, 3.5806225183765144`*^9}, {3.5806225509805717`*^9, 3.580622580776624*^9}, 3.580622715810461*^9, 3.5806228003314095`*^9, 3.5806228736203384`*^9, 3.580622909126001*^9, {3.580622942120059*^9, 3.5806229668929024`*^9}, 3.580623010385779*^9, {3.5806230636442723`*^9, 3.5806230933935246`*^9}, { 3.580623145045215*^9, 3.580623172204863*^9}, {3.5806232036857185`*^9, 3.580623224121754*^9}, 3.5806233381735544`*^9, 3.580623395410055*^9, { 3.5806234268597107`*^9, 3.580623473893793*^9}, 3.5806235366527033`*^9, { 3.580623593811204*^9, 3.580623622811654*^9}, 3.5806236823101587`*^9, 3.580623718798623*^9, 3.5806237839755373`*^9, 3.5806238162207937`*^9, 3.580623861320473*^9, 3.5806239228937817`*^9}], Cell[BoxData["0"], "Output", CellChangeTimes->{ 3.5806206427228203`*^9, 3.5806213348492355`*^9, 3.5806222660616713`*^9, 3.580622297870127*^9, 3.5806223535778255`*^9, {3.580622387227084*^9, 3.5806224300335593`*^9}, {3.5806224690960283`*^9, 3.5806225183765144`*^9}, {3.5806225509805717`*^9, 3.580622580776624*^9}, 3.580622715810461*^9, 3.5806228003314095`*^9, 3.5806228736203384`*^9, 3.580622909126001*^9, {3.580622942120059*^9, 3.5806229668929024`*^9}, 3.580623010385779*^9, {3.5806230636442723`*^9, 3.5806230933935246`*^9}, { 3.580623145045215*^9, 3.580623172204863*^9}, {3.5806232036857185`*^9, 3.580623224121754*^9}, 3.5806233381735544`*^9, 3.580623395410055*^9, { 3.5806234268597107`*^9, 3.580623473893793*^9}, 3.5806235366527033`*^9, { 3.580623593811204*^9, 3.580623622811654*^9}, 3.5806236823101587`*^9, 3.580623718798623*^9, 3.5806237839755373`*^9, 3.5806238162207937`*^9, 3.580623861320473*^9, 3.5806239228937817`*^9}] }, Open ]], 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", ",", "cw"}], "}"}]}], "]"}]}], ";", " ", 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", ",", "cw"}], "}"}]}], "]"}]}], ";", " ", 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", ",", "cw"}], "}"}]}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"Plot7", "=", 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", ",", "cw"}], "}"}]}], "]"}]}], ";"}]}]}]], "Input", CellChangeTimes->{{3.5681365021632137`*^9, 3.5681365171392403`*^9}, { 3.568136560351316*^9, 3.5681365863097615`*^9}, {3.5681369135895367`*^9, 3.5681369265999594`*^9}, {3.568137821916732*^9, 3.5681378378443604`*^9}, { 3.5682072492149415`*^9, 3.5682072529901485`*^9}, 3.568209802830227*^9, { 3.5682216503836985`*^9, 3.5682217583046885`*^9}, {3.5682217912831464`*^9, 3.568221843465238*^9}, {3.5682219760654707`*^9, 3.5682221133457117`*^9}, { 3.568222159615393*^9, 3.5682221693498106`*^9}, {3.5682222565227633`*^9, 3.5682222578487663`*^9}, 3.5682224736283445`*^9, 3.5682227227607822`*^9, 3.5682227558640404`*^9, {3.568222818810151*^9, 3.5682228486062036`*^9}, { 3.5682230057764797`*^9, 3.5682230079916835`*^9}, {3.568223041952943*^9, 3.5682230548229656`*^9}, 3.5682232619601297`*^9, 3.5682233000709963`*^9, { 3.568223536645412*^9, 3.5682235381118145`*^9}, {3.568223591151908*^9, 3.568223591526308*^9}, {3.5682236491528096`*^9, 3.568223657483224*^9}, { 3.568223713222122*^9, 3.5682237532361927`*^9}, {3.568223800379475*^9, 3.568223827648323*^9}, {3.5682238887380304`*^9, 3.5682238893152313`*^9}, { 3.5682240287170763`*^9, 3.5682240612431335`*^9}, 3.5682241088076167`*^9, { 3.568224198772975*^9, 3.5682242396450467`*^9}, {3.568224282170721*^9, 3.5682243185031853`*^9}, {3.568224437890195*^9, 3.568224439590598*^9}, { 3.568224572253231*^9, 3.5682245867144566`*^9}, {3.5682247081762695`*^9, 3.568224733557514*^9}, {3.56822480520844*^9, 3.5682248433505073`*^9}, { 3.568224884206979*^9, 3.5682248959069996`*^9}, {3.568224982830352*^9, 3.5682250047639904`*^9}, {3.56822503870965*^9, 3.5682250803305235`*^9}, { 3.5682251223413973`*^9, 3.568225141451431*^9}, {3.568225234177994*^9, 3.5682252929588966`*^9}, {3.5682253775890455`*^9, 3.568225482140429*^9}, { 3.56822551694409*^9, 3.5682255297205124`*^9}, {3.568225573197789*^9, 3.5682255809666023`*^9}, 3.568225660698343*^9, 3.5682257301964645`*^9, { 3.568225765078126*^9, 3.5682257828933573`*^9}, {3.568225927026843*^9, 3.568225927416844*^9}, {3.568267062420238*^9, 3.568267081358671*^9}, { 3.568267188328059*^9, 3.5682672830358257`*^9}, {3.5682676367040462`*^9, 3.568267747979042*^9}, 3.5682678404716043`*^9, 3.568267882030078*^9, { 3.5682679252577534`*^9, 3.568267996612279*^9}, {3.5682680488879704`*^9, 3.568268061133992*^9}, {3.568268100914062*^9, 3.5682681150944867`*^9}, { 3.5682681672609787`*^9, 3.5682681709113846`*^9}, {3.568268272904364*^9, 3.5682683182536435`*^9}, {3.5682685202739983`*^9, 3.568268599194537*^9}, { 3.5682686348873997`*^9, 3.568268637773405*^9}, {3.5682686944951043`*^9, 3.568268872273017*^9}, {3.5682689420207396`*^9, 3.568269016947671*^9}, { 3.5682695040741262`*^9, 3.568269601153097*^9}, {3.5682696652380095`*^9, 3.5682696728352227`*^9}, {3.5682697489633565`*^9, 3.5682698308635006`*^9}, {3.568271156491429*^9, 3.568271160063835*^9}, { 3.568271244927984*^9, 3.568271369010602*^9}, {3.5682716250070515`*^9, 3.5682717312276382`*^9}, {3.5682717843301315`*^9, 3.5682718576814604`*^9}, {3.568271940673606*^9, 3.5682720536490045`*^9}, { 3.5682722259669075`*^9, 3.5682723223126764`*^9}, {3.568272370813162*^9, 3.568272411466833*^9}, {3.5682725893539457`*^9, 3.568272622238803*^9}, { 3.5682726936713285`*^9, 3.5682726960113325`*^9}, {3.568272741079812*^9, 3.56827275143823*^9}, {3.5682729218373294`*^9, 3.568272991070251*^9}, { 3.5682734285574193`*^9, 3.5682734302890224`*^9}, {3.5682734898811274`*^9, 3.5682735111595645`*^9}, {3.5682736115925407`*^9, 3.5682736313577757`*^9}, {3.5682738601945777`*^9, 3.568273987272401*^9}, { 3.5682763944346285`*^9, 3.5682764313910933`*^9}, 3.5682765128232365`*^9, { 3.5682766527398825`*^9, 3.568276660976697*^9}, {3.5682774496048822`*^9, 3.568277461663703*^9}, {3.5682852978767633`*^9, 3.568285382460112*^9}, { 3.568285675366227*^9, 3.568285707861084*^9}, {3.568285985822372*^9, 3.568285991485182*^9}, {3.568286358600627*^9, 3.568286369505046*^9}, { 3.5682869898869357`*^9, 3.568287001134555*^9}, {3.5682977174462996`*^9, 3.568297735994732*^9}, {3.568305934649544*^9, 3.5683059800144234`*^9}, { 3.5804489722380853`*^9, 3.580449074449465*^9}, {3.580449125945155*^9, 3.580449181777653*^9}, 3.580449245191765*^9, 3.580449330789115*^9, { 3.5804589977940884`*^9, 3.580459000102893*^9}, {3.5804590597729974`*^9, 3.5804590762778263`*^9}, 3.5804591100518856`*^9, 3.580459208191658*^9, { 3.5804592899358015`*^9, 3.580459292322606*^9}, {3.580459332664277*^9, 3.5804593337250786`*^9}, {3.5804599640909853`*^9, 3.5804599656041884`*^9}, {3.580460783217224*^9, 3.580460832840912*^9}, { 3.580460867847373*^9, 3.5804608846018023`*^9}, {3.580465026299877*^9, 3.5804650871087837`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"wxh", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "3"}], "]"}], "]"}], "+", RowBox[{"u2", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ")"}], "/", RowBox[{"(", RowBox[{ RowBox[{"u2", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "+", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "3"}], "]"}], "]"}], "+", RowBox[{"u2", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ")"}]}], "*", RowBox[{"(", RowBox[{ RowBox[{"x", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "-", RowBox[{"x", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ")"}]}], "+", RowBox[{"x", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ",", RowBox[{"{", RowBox[{"i", ",", RowBox[{"cw", "+", "1"}], ",", "nb"}], "}"}]}], "]"}]}], ";", " ", RowBox[{"wyh", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "3"}], "]"}], "]"}], "+", RowBox[{"u2", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ")"}], "/", RowBox[{"(", RowBox[{ RowBox[{"u2", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "+", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "3"}], "]"}], "]"}], "+", RowBox[{"u2", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ")"}]}], "*", RowBox[{"(", RowBox[{ RowBox[{"y", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "-", RowBox[{"y", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ")"}]}], "+", RowBox[{"y", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ",", RowBox[{"{", RowBox[{"i", ",", RowBox[{"cw", "+", "1"}], ",", "nb"}], "}"}]}], "]"}]}], ";", " ", RowBox[{"wzh", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "3"}], "]"}], "]"}], "+", RowBox[{"u2", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ")"}], "/", RowBox[{"(", RowBox[{ RowBox[{"u2", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "+", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "3"}], "]"}], "]"}], "+", RowBox[{"u2", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ")"}]}], "*", RowBox[{"(", RowBox[{ RowBox[{"z", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "-", RowBox[{"z", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ")"}]}], "+", RowBox[{"z", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ",", RowBox[{"{", RowBox[{"i", ",", RowBox[{"cw", "+", "1"}], ",", "nb"}], "}"}]}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.5681365021632137`*^9, 3.5681365171392403`*^9}, { 3.568136560351316*^9, 3.5681365863097615`*^9}, {3.5681369135895367`*^9, 3.5681369265999594`*^9}, {3.568137821916732*^9, 3.5681378378443604`*^9}, { 3.5682072492149415`*^9, 3.5682072529901485`*^9}, 3.568209802830227*^9, { 3.5682216503836985`*^9, 3.5682217583046885`*^9}, {3.5682217912831464`*^9, 3.568221843465238*^9}, {3.5682219760654707`*^9, 3.5682221133457117`*^9}, { 3.568222159615393*^9, 3.5682221693498106`*^9}, {3.5682222565227633`*^9, 3.5682222578487663`*^9}, 3.5682224736283445`*^9, 3.5682227227607822`*^9, 3.5682227558640404`*^9, {3.568222818810151*^9, 3.5682228486062036`*^9}, { 3.5682230057764797`*^9, 3.5682230079916835`*^9}, {3.568223041952943*^9, 3.5682230548229656`*^9}, 3.5682232619601297`*^9, 3.5682233000709963`*^9, { 3.568223536645412*^9, 3.5682235381118145`*^9}, {3.568223591151908*^9, 3.568223591526308*^9}, {3.5682236491528096`*^9, 3.568223657483224*^9}, { 3.568223713222122*^9, 3.5682237532361927`*^9}, {3.568223800379475*^9, 3.568223827648323*^9}, {3.5682238887380304`*^9, 3.5682238893152313`*^9}, { 3.5682240287170763`*^9, 3.5682240612431335`*^9}, 3.5682241088076167`*^9, { 3.568224198772975*^9, 3.5682242396450467`*^9}, {3.568224282170721*^9, 3.5682243185031853`*^9}, {3.568224437890195*^9, 3.568224439590598*^9}, { 3.568224572253231*^9, 3.5682245867144566`*^9}, {3.5682247081762695`*^9, 3.568224733557514*^9}, {3.56822480520844*^9, 3.5682248433505073`*^9}, { 3.568224884206979*^9, 3.5682248959069996`*^9}, {3.568224982830352*^9, 3.5682250047639904`*^9}, {3.56822503870965*^9, 3.5682250803305235`*^9}, { 3.5682251223413973`*^9, 3.568225141451431*^9}, {3.568225234177994*^9, 3.5682252929588966`*^9}, {3.5682253775890455`*^9, 3.568225482140429*^9}, { 3.56822551694409*^9, 3.5682255297205124`*^9}, {3.568225573197789*^9, 3.5682255809666023`*^9}, 3.568225660698343*^9, 3.5682257301964645`*^9, { 3.568225765078126*^9, 3.5682257828933573`*^9}, {3.568225927026843*^9, 3.568225927416844*^9}, {3.568267062420238*^9, 3.568267081358671*^9}, { 3.568267188328059*^9, 3.5682672830358257`*^9}, {3.5682676367040462`*^9, 3.568267747979042*^9}, 3.5682678404716043`*^9, 3.568267882030078*^9, { 3.5682679252577534`*^9, 3.568267996612279*^9}, {3.5682680488879704`*^9, 3.568268061133992*^9}, {3.568268100914062*^9, 3.5682681150944867`*^9}, { 3.5682681672609787`*^9, 3.5682681709113846`*^9}, {3.568268272904364*^9, 3.5682683182536435`*^9}, {3.5682685202739983`*^9, 3.568268599194537*^9}, { 3.5682686348873997`*^9, 3.568268637773405*^9}, {3.5682686944951043`*^9, 3.568268872273017*^9}, {3.5682689420207396`*^9, 3.568269016947671*^9}, { 3.5682695040741262`*^9, 3.568269601153097*^9}, {3.5682696652380095`*^9, 3.5682696728352227`*^9}, {3.5682697489633565`*^9, 3.5682698308635006`*^9}, {3.568271156491429*^9, 3.568271160063835*^9}, { 3.568271244927984*^9, 3.568271369010602*^9}, {3.5682716250070515`*^9, 3.5682717312276382`*^9}, {3.5682717843301315`*^9, 3.5682718576814604`*^9}, {3.568271940673606*^9, 3.5682720536490045`*^9}, { 3.5682722259669075`*^9, 3.5682723223126764`*^9}, {3.568272370813162*^9, 3.568272411466833*^9}, {3.5682725893539457`*^9, 3.568272622238803*^9}, { 3.5682726936713285`*^9, 3.5682726960113325`*^9}, {3.568272741079812*^9, 3.56827275143823*^9}, {3.5682729218373294`*^9, 3.568272991070251*^9}, { 3.5682734285574193`*^9, 3.5682734302890224`*^9}, {3.5682734898811274`*^9, 3.5682735111595645`*^9}, {3.5682736115925407`*^9, 3.5682736313577757`*^9}, {3.5682738601945777`*^9, 3.568273987272401*^9}, { 3.5682763944346285`*^9, 3.5682764313910933`*^9}, 3.5682765128232365`*^9, { 3.5682766527398825`*^9, 3.568276660976697*^9}, {3.5682774496048822`*^9, 3.568277461663703*^9}, {3.5682852978767633`*^9, 3.568285382460112*^9}, { 3.568285675366227*^9, 3.568285707861084*^9}, {3.568285985822372*^9, 3.568285991485182*^9}, {3.568286358600627*^9, 3.568286369505046*^9}, { 3.5682869898869357`*^9, 3.568287001134555*^9}, {3.5682977174462996`*^9, 3.568297735994732*^9}, {3.568305934649544*^9, 3.5683059800144234`*^9}, { 3.5804489722380853`*^9, 3.580449074449465*^9}, {3.580449125945155*^9, 3.580449181777653*^9}, 3.580449245191765*^9, 3.580449330789115*^9, { 3.5804589977940884`*^9, 3.580459000102893*^9}, {3.5804590597729974`*^9, 3.5804590762778263`*^9}, 3.5804591100518856`*^9, 3.580459208191658*^9, { 3.5804592899358015`*^9, 3.580459292322606*^9}, {3.580459332664277*^9, 3.5804593337250786`*^9}, {3.5804599640909853`*^9, 3.5804599656041884`*^9}, {3.580460783217224*^9, 3.580460832840912*^9}, { 3.580460867847373*^9, 3.5804608846018023`*^9}, {3.580465026299877*^9, 3.5804650896203885`*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"Plot8", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"Graphics3D", "[", RowBox[{"{", RowBox[{ RowBox[{"Opacity", "[", "0.5", "]"}], ",", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"wxh", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"wyh", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"wzh", "[", RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", RowBox[{"a", "[", RowBox[{"[", RowBox[{ RowBox[{"i", "+", "cw"}], ",", "3"}], "]"}], "]"}]}], "]"}]}], "}"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "bh"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Plot8ib", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"Graphics3D", "[", RowBox[{"{", RowBox[{ RowBox[{"Opacity", "[", "0.4", "]"}], ",", RowBox[{"Lighter", "[", "Red", "]"}], ",", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"tpso", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"rlos", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}], "}"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "hj"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Plot8jb", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"Graphics3D", "[", RowBox[{"{", RowBox[{ RowBox[{"Opacity", "[", "0.4", "]"}], ",", RowBox[{"Lighter", "[", "Red", "]"}], ",", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"tpsu", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"rlos", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}], "}"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "hj"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Plot9", "=", " ", RowBox[{"Table", "[", RowBox[{ RowBox[{"Graphics3D", "[", RowBox[{"{", RowBox[{ RowBox[{"Lighter", "[", RowBox[{"Gray", ",", "0.8"}], "]"}], ",", RowBox[{"Tube", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"x", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], ",", RowBox[{"y", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], ",", RowBox[{"z", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"x", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}], ",", RowBox[{"y", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}], ",", RowBox[{"z", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], "}"}]}], "}"}], ",", "0.07"}], "]"}]}], "}"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "nb"}], "}"}]}], "]"}]}], ";"}]}], "Input", CellChangeTimes->{{3.5681365021632137`*^9, 3.5681365171392403`*^9}, { 3.568136560351316*^9, 3.5681365863097615`*^9}, {3.5681369135895367`*^9, 3.5681369265999594`*^9}, {3.568137821916732*^9, 3.5681378378443604`*^9}, { 3.5682072492149415`*^9, 3.5682072529901485`*^9}, 3.568209802830227*^9, { 3.5682216503836985`*^9, 3.5682217583046885`*^9}, {3.5682217912831464`*^9, 3.568221843465238*^9}, {3.5682219760654707`*^9, 3.5682221133457117`*^9}, { 3.568222159615393*^9, 3.5682221693498106`*^9}, {3.5682222565227633`*^9, 3.5682222578487663`*^9}, 3.5682224736283445`*^9, 3.5682227227607822`*^9, 3.5682227558640404`*^9, {3.568222818810151*^9, 3.5682228486062036`*^9}, { 3.5682230057764797`*^9, 3.5682230079916835`*^9}, {3.568223041952943*^9, 3.5682230548229656`*^9}, 3.5682232619601297`*^9, 3.5682233000709963`*^9, { 3.568223536645412*^9, 3.5682235381118145`*^9}, {3.568223591151908*^9, 3.568223591526308*^9}, {3.5682236491528096`*^9, 3.568223657483224*^9}, { 3.568223713222122*^9, 3.5682237532361927`*^9}, {3.568223800379475*^9, 3.568223827648323*^9}, {3.5682238887380304`*^9, 3.5682238893152313`*^9}, { 3.5682240287170763`*^9, 3.5682240612431335`*^9}, 3.5682241088076167`*^9, { 3.568224198772975*^9, 3.5682242396450467`*^9}, {3.568224282170721*^9, 3.5682243185031853`*^9}, {3.568224437890195*^9, 3.568224439590598*^9}, { 3.568224572253231*^9, 3.5682245867144566`*^9}, {3.5682247081762695`*^9, 3.568224733557514*^9}, {3.56822480520844*^9, 3.5682248433505073`*^9}, { 3.568224884206979*^9, 3.5682248959069996`*^9}, {3.568224982830352*^9, 3.5682250047639904`*^9}, {3.56822503870965*^9, 3.5682250803305235`*^9}, { 3.5682251223413973`*^9, 3.568225141451431*^9}, {3.568225234177994*^9, 3.5682252929588966`*^9}, {3.5682253775890455`*^9, 3.568225482140429*^9}, { 3.56822551694409*^9, 3.5682255297205124`*^9}, {3.568225573197789*^9, 3.5682255809666023`*^9}, 3.568225660698343*^9, 3.5682257301964645`*^9, { 3.568225765078126*^9, 3.5682257828933573`*^9}, {3.568225927026843*^9, 3.568225927416844*^9}, {3.568267062420238*^9, 3.568267081358671*^9}, { 3.568267188328059*^9, 3.5682672830358257`*^9}, {3.5682676367040462`*^9, 3.568267747979042*^9}, 3.5682678404716043`*^9, 3.568267882030078*^9, { 3.5682679252577534`*^9, 3.568267996612279*^9}, {3.5682680488879704`*^9, 3.568268061133992*^9}, {3.568268100914062*^9, 3.5682681150944867`*^9}, { 3.5682681672609787`*^9, 3.5682681709113846`*^9}, {3.568268272904364*^9, 3.5682683182536435`*^9}, {3.5682685202739983`*^9, 3.568268599194537*^9}, { 3.5682686348873997`*^9, 3.568268637773405*^9}, {3.5682686944951043`*^9, 3.568268872273017*^9}, {3.5682689420207396`*^9, 3.568269016947671*^9}, { 3.5682695040741262`*^9, 3.568269601153097*^9}, {3.5682696652380095`*^9, 3.5682696728352227`*^9}, {3.5682697489633565`*^9, 3.5682698308635006`*^9}, {3.568271156491429*^9, 3.568271160063835*^9}, { 3.568271244927984*^9, 3.568271369010602*^9}, {3.5682716250070515`*^9, 3.5682717312276382`*^9}, {3.5682717843301315`*^9, 3.5682718576814604`*^9}, {3.568271940673606*^9, 3.5682720536490045`*^9}, { 3.5682722259669075`*^9, 3.5682723223126764`*^9}, {3.568272370813162*^9, 3.568272411466833*^9}, {3.5682725893539457`*^9, 3.568272622238803*^9}, { 3.5682726936713285`*^9, 3.5682726960113325`*^9}, {3.568272741079812*^9, 3.56827275143823*^9}, {3.5682729218373294`*^9, 3.568272991070251*^9}, { 3.5682734285574193`*^9, 3.5682734302890224`*^9}, {3.5682734898811274`*^9, 3.5682735111595645`*^9}, {3.5682736115925407`*^9, 3.5682736313577757`*^9}, {3.5682738601945777`*^9, 3.568273987272401*^9}, { 3.5682763944346285`*^9, 3.5682764313910933`*^9}, 3.5682765128232365`*^9, { 3.5682766527398825`*^9, 3.568276660976697*^9}, {3.5682774496048822`*^9, 3.568277461663703*^9}, {3.5682852978767633`*^9, 3.568285382460112*^9}, { 3.568285675366227*^9, 3.568285707861084*^9}, {3.568285985822372*^9, 3.568285991485182*^9}, {3.568286358600627*^9, 3.568286369505046*^9}, { 3.5682869898869357`*^9, 3.568287001134555*^9}, {3.5682977174462996`*^9, 3.568297735994732*^9}, {3.568305934649544*^9, 3.5683059800144234`*^9}, { 3.5804489722380853`*^9, 3.580449074449465*^9}, {3.580449125945155*^9, 3.580449181777653*^9}, 3.580449245191765*^9, 3.580449330789115*^9, { 3.5804589977940884`*^9, 3.580459000102893*^9}, {3.5804590597729974`*^9, 3.5804590762778263`*^9}, 3.5804591100518856`*^9, 3.580459208191658*^9, { 3.5804592899358015`*^9, 3.580459292322606*^9}, {3.580459332664277*^9, 3.5804593337250786`*^9}, {3.5804599640909853`*^9, 3.5804599656041884`*^9}, {3.580460783217224*^9, 3.580460832840912*^9}, { 3.580460867847373*^9, 3.5804608846018023`*^9}, {3.580465026299877*^9, 3.5804651023032103`*^9}, {3.5804653879553123`*^9, 3.580465390466917*^9}, { 3.5804654257853785`*^9, 3.580465427017781*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Show", "[", RowBox[{ "Plot1o", ",", "Plot2", ",", "Plot7", ",", "Plot8", ",", "Plot8ib", ",", "Plot8jb", ",", "Plot9", ",", RowBox[{"{", RowBox[{ RowBox[{"AspectRatio", "\[Rule]", "Automatic"}], ",", RowBox[{"Boxed", "\[Rule]", "False"}], ",", RowBox[{"SphericalRegion", "->", "True"}]}], "}"}]}], "]"}]], "Input", CellChangeTimes->{{3.5681365021632137`*^9, 3.5681365171392403`*^9}, { 3.568136560351316*^9, 3.5681365863097615`*^9}, {3.5681369135895367`*^9, 3.5681369265999594`*^9}, {3.568137821916732*^9, 3.5681378378443604`*^9}, { 3.5682072492149415`*^9, 3.5682072529901485`*^9}, 3.568209802830227*^9, { 3.5682216503836985`*^9, 3.5682217583046885`*^9}, {3.5682217912831464`*^9, 3.568221843465238*^9}, {3.5682219760654707`*^9, 3.5682221133457117`*^9}, { 3.568222159615393*^9, 3.5682221693498106`*^9}, {3.5682222565227633`*^9, 3.5682222578487663`*^9}, 3.5682224736283445`*^9, 3.5682227227607822`*^9, 3.5682227558640404`*^9, {3.568222818810151*^9, 3.5682228486062036`*^9}, { 3.5682230057764797`*^9, 3.5682230079916835`*^9}, {3.568223041952943*^9, 3.5682230548229656`*^9}, 3.5682232619601297`*^9, 3.5682233000709963`*^9, { 3.568223536645412*^9, 3.5682235381118145`*^9}, {3.568223591151908*^9, 3.568223591526308*^9}, {3.5682236491528096`*^9, 3.568223657483224*^9}, { 3.568223713222122*^9, 3.5682237532361927`*^9}, {3.568223800379475*^9, 3.568223827648323*^9}, {3.5682238887380304`*^9, 3.5682238893152313`*^9}, { 3.5682240287170763`*^9, 3.5682240612431335`*^9}, 3.5682241088076167`*^9, { 3.568224198772975*^9, 3.5682242396450467`*^9}, {3.568224282170721*^9, 3.5682243185031853`*^9}, {3.568224437890195*^9, 3.568224439590598*^9}, { 3.568224572253231*^9, 3.5682245867144566`*^9}, {3.5682247081762695`*^9, 3.568224733557514*^9}, {3.56822480520844*^9, 3.5682248433505073`*^9}, { 3.568224884206979*^9, 3.5682248959069996`*^9}, {3.568224982830352*^9, 3.5682250047639904`*^9}, {3.56822503870965*^9, 3.5682250803305235`*^9}, { 3.5682251223413973`*^9, 3.568225141451431*^9}, {3.568225234177994*^9, 3.5682252929588966`*^9}, {3.5682253775890455`*^9, 3.568225482140429*^9}, { 3.56822551694409*^9, 3.5682255297205124`*^9}, {3.568225573197789*^9, 3.5682255809666023`*^9}, 3.568225660698343*^9, 3.5682257301964645`*^9, { 3.568225765078126*^9, 3.5682257828933573`*^9}, {3.568225927026843*^9, 3.568225927416844*^9}, {3.568267062420238*^9, 3.568267081358671*^9}, { 3.568267188328059*^9, 3.5682672830358257`*^9}, {3.5682676367040462`*^9, 3.568267747979042*^9}, 3.5682678404716043`*^9, 3.568267882030078*^9, { 3.5682679252577534`*^9, 3.568267996612279*^9}, {3.5682680488879704`*^9, 3.568268061133992*^9}, {3.568268100914062*^9, 3.5682681150944867`*^9}, { 3.5682681672609787`*^9, 3.5682681709113846`*^9}, {3.568268272904364*^9, 3.5682683182536435`*^9}, {3.5682685202739983`*^9, 3.568268599194537*^9}, { 3.5682686348873997`*^9, 3.568268637773405*^9}, {3.5682686944951043`*^9, 3.568268872273017*^9}, {3.5682689420207396`*^9, 3.568269016947671*^9}, { 3.5682695040741262`*^9, 3.568269601153097*^9}, {3.5682696652380095`*^9, 3.5682696728352227`*^9}, {3.5682697489633565`*^9, 3.5682698308635006`*^9}, {3.568271156491429*^9, 3.568271160063835*^9}, { 3.568271244927984*^9, 3.568271369010602*^9}, {3.5682716250070515`*^9, 3.5682717312276382`*^9}, {3.5682717843301315`*^9, 3.5682718576814604`*^9}, {3.568271940673606*^9, 3.5682720536490045`*^9}, { 3.5682722259669075`*^9, 3.5682723223126764`*^9}, {3.568272370813162*^9, 3.568272411466833*^9}, {3.5682725893539457`*^9, 3.568272622238803*^9}, { 3.5682726936713285`*^9, 3.5682726960113325`*^9}, {3.568272741079812*^9, 3.56827275143823*^9}, {3.5682729218373294`*^9, 3.568272991070251*^9}, { 3.5682734285574193`*^9, 3.5682734302890224`*^9}, {3.5682734898811274`*^9, 3.5682735111595645`*^9}, {3.5682736115925407`*^9, 3.5682736313577757`*^9}, {3.5682738601945777`*^9, 3.568273987272401*^9}, { 3.5682763944346285`*^9, 3.5682764313910933`*^9}, 3.5682765128232365`*^9, { 3.5682766527398825`*^9, 3.568276660976697*^9}, {3.5682774496048822`*^9, 3.568277461663703*^9}, {3.5682852978767633`*^9, 3.568285382460112*^9}, { 3.568285675366227*^9, 3.568285707861084*^9}, {3.568285985822372*^9, 3.568285991485182*^9}, {3.568286358600627*^9, 3.568286369505046*^9}, { 3.5682869898869357`*^9, 3.568287001134555*^9}, {3.5682977174462996`*^9, 3.568297735994732*^9}, {3.568305934649544*^9, 3.5683059800144234`*^9}, { 3.5804489722380853`*^9, 3.580449074449465*^9}, {3.580449125945155*^9, 3.580449181777653*^9}, 3.580449245191765*^9, 3.580449330789115*^9, { 3.5804589977940884`*^9, 3.580459000102893*^9}, {3.5804590597729974`*^9, 3.5804590762778263`*^9}, 3.5804591100518856`*^9, 3.580459208191658*^9, { 3.5804592899358015`*^9, 3.580459292322606*^9}, {3.580459332664277*^9, 3.5804593337250786`*^9}, {3.5804599640909853`*^9, 3.5804599656041884`*^9}, {3.580460783217224*^9, 3.580460832840912*^9}, { 3.580460867847373*^9, 3.5804608846018023`*^9}, {3.580465026299877*^9, 3.580465098200403*^9}, 3.580465296975952*^9, {3.580488471874737*^9, 3.580488483902358*^9}}], Cell[BoxData[ Graphics3DBox[{ {RGBColor[1, 1, 1], Opacity[0.8], SphereBox[{{0., -0.150474057419748, 0.}}, 0.18]}, {RGBColor[1, 1, 1], Opacity[0.8], SphereBox[{{1.4422535718512253`, 0.976911818490493, 0.}}, 0.12]}, {RGBColor[1, 1, 1], Opacity[0.8], SphereBox[{{-1.4422535718512253`, 0.976911818490493, 0.}}, 0.12]}, Text3DBox["1", {-0.5, -0.150474057419748, 0.}, {1, 0}], Text3DBox["2", {0.9422535718512253, 0.976911818490493, 0.}, {1, 0}], Text3DBox["3", {-1.9422535718512253`, 0.976911818490493, 0.}, {1, 0}], {Opacity[0.5], SphereBox[{1.0325670809428116`, 0.6566666324733568, 0.}, 1.1313994315304259`]}, {Opacity[0.5], SphereBox[{-1.0325670809428116`, 0.6566666324733568, 0.}, 1.1313994315304259`]}, {RGBColor[1, NCache[ Rational[1, 3], 0.3333333333333333], NCache[ Rational[1, 3], 0.3333333333333333]], Opacity[0.4], SphereBox[{1.1102230246251565`*^-16, -0.8413740574197479, 1.153125}, 1.125]}, {RGBColor[1, NCache[ Rational[1, 3], 0.3333333333333333], NCache[ Rational[1, 3], 0.3333333333333333]], Opacity[0.4], SphereBox[{1.1102230246251565`*^-16, -0.8413740574197479, -1.153125}, 1.125]}, {RGBColor[0.9, 0.9, 0.9], TubeBox[{{1.4422535718512253`, 0.976911818490493, 0.}, { 0., -0.150474057419748, 0.}}, 0.07]}, {RGBColor[0.9, 0.9, 0.9], TubeBox[{{-1.4422535718512253`, 0.976911818490493, 0.}, { 0., -0.150474057419748, 0.}}, 0.07]}}, AspectRatio->Automatic, AutomaticImageSize->True, Boxed->False, ImageSize->{360., 360.}, SphericalRegion->True, ViewPoint->{-3.3785508031509934`, 0.15364580406235057`, -0.10856996556053142`}, ViewVertical->{-0.6470358401644789, -0.9347881773468355, 0.026179365228894982`}]], "Output", CellChangeTimes->{ 3.58062064275402*^9, 3.580621334896036*^9, 3.5806222660928717`*^9, 3.5806222979169273`*^9, 3.5806223536402254`*^9, {3.5806223872738843`*^9, 3.5806224300647593`*^9}, {3.580622469127228*^9, 3.580622518423315*^9}, { 3.580622551011772*^9, 3.580622580807824*^9}, 3.5806227158572617`*^9, 3.58062280037821*^9, 3.580622873651539*^9, 3.580622909157201*^9, { 3.5806229421668587`*^9, 3.5806229669241023`*^9}, 3.580623010432579*^9, { 3.5806230636754723`*^9, 3.580623093440325*^9}, {3.580623145076415*^9, 3.580623172251663*^9}, {3.580623203732518*^9, 3.580623224152954*^9}, 3.5806233382203546`*^9, 3.580623395441255*^9, {3.5806234269065104`*^9, 3.580623473940593*^9}, 3.580623536683903*^9, {3.5806235938580036`*^9, 3.5806236228584547`*^9}, 3.5806236823569593`*^9, 3.580623718845423*^9, 3.5806237840067377`*^9, 3.5806238162675943`*^9, 3.5806238613672733`*^9, 3.5806239229405813`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"tc", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"x", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"y", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"z", "[", RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "bc"}], "}"}]}], "]"}]}], ";"}], " ", RowBox[{"(*", " ", RowBox[{"core", " ", "clouds"}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"twc", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"wx", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"wy", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"wz", "[", RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "cw"}], "}"}]}], "]"}]}], ";"}], " ", RowBox[{"(*", " ", RowBox[{"Bond", " ", "c\[ODoubleDot]ouds", " ", "between", " ", "cores"}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"twh", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"wxh", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"wyh", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"wzh", "[", RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "bh"}], "}"}]}], "]"}]}], ";"}], " ", RowBox[{"(*", " ", RowBox[{"H", "-", RowBox[{"bonding", " ", "clouds"}]}], " ", "*)"}], " "}], "\[IndentingNewLine]", RowBox[{ RowBox[{"tw", "=", RowBox[{"Join", "[", RowBox[{"tc", ",", "twc", ",", "twh", ",", "tpso", ",", "tpsu"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"bt", "=", RowBox[{"Length", "[", "tw", "]"}]}]}], "Input", CellChangeTimes->{{3.5804614754216404`*^9, 3.5804614810532503`*^9}, { 3.580461527603732*^9, 3.5804615546541796`*^9}, 3.5804647573086047`*^9, 3.5804649196424894`*^9, {3.58048949832534*^9, 3.5804895043313503`*^9}, 3.580489609803136*^9}], Cell[BoxData["5"], "Output", CellChangeTimes->{ 3.5806206427696204`*^9, 3.5806213349116364`*^9, 3.5806222661552715`*^9, 3.5806222979325275`*^9, 3.580622353687025*^9, {3.580622387305084*^9, 3.5806224300803595`*^9}, {3.5806224691740284`*^9, 3.580622518438915*^9}, { 3.580622551074172*^9, 3.5806225808234243`*^9}, 3.5806227158572617`*^9, 3.5806228003938103`*^9, 3.5806228736983385`*^9, 3.580622909204001*^9, { 3.5806229422136593`*^9, 3.5806229669709024`*^9}, 3.5806230104793787`*^9, { 3.5806230636910725`*^9, 3.580623093455925*^9}, {3.5806231450920153`*^9, 3.5806231722672634`*^9}, {3.5806232037481184`*^9, 3.5806232241685543`*^9}, 3.5806233382359543`*^9, 3.5806233955036554`*^9, {3.5806234269533105`*^9, 3.580623473987393*^9}, 3.5806235367307034`*^9, {3.5806235939048033`*^9, 3.5806236229052544`*^9}, 3.580623682372559*^9, 3.580623718861023*^9, 3.5806237840535374`*^9, 3.5806238163143945`*^9, 3.5806238613828735`*^9, 3.5806239229405813`*^9}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{"cloud", "-", RowBox[{"cloud", " ", "distances"}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"dw", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"Norm", "[", RowBox[{ RowBox[{"tw", "[", RowBox[{"[", "i", "]"}], "]"}], "-", RowBox[{"tw", "[", RowBox[{"[", "j", "]"}], "]"}]}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "bt"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", "bt"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"MatrixForm", "[", "dw", "]"}], "*)"}]}]}]], "Input", CellChangeTimes->{{3.5804614754216404`*^9, 3.5804614810532503`*^9}, { 3.580461527603732*^9, 3.5804615546541796`*^9}, {3.5804647573086047`*^9, 3.58046477754184*^9}, {3.580464923870097*^9, 3.5804649422313294`*^9}, 3.5804897528397865`*^9, 3.580489783634241*^9}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{"nucleus", "-", RowBox[{"nucleus", " ", "distances"}]}], " ", "*)"}], "\n", RowBox[{ RowBox[{ RowBox[{"tv1", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"t", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "bc"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"tv2", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "cw"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"tv3", "=", " ", RowBox[{"Table", "[", RowBox[{ RowBox[{"t", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", RowBox[{"bc", "+", "1"}], ",", "n"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"tv6", " ", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", RowBox[{"2", "*", "hj"}]}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"tv", " ", "=", " ", RowBox[{"Join", "[", RowBox[{"tv1", ",", "tv2", ",", "tv3", ",", "tv6"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"MatrixForm", "[", "tv", "]"}], "*)"}], "\[IndentingNewLine]", RowBox[{"Length", "[", "tv", "]"}]}]}]], "Input", CellChangeTimes->{{3.5804614754216404`*^9, 3.5804614810532503`*^9}, { 3.580461527603732*^9, 3.5804615546541796`*^9}, {3.5804647573086047`*^9, 3.58046477754184*^9}, {3.580464923870097*^9, 3.5804649550077515`*^9}, { 3.5804899034268513`*^9, 3.5804899043784533`*^9}}], Cell[BoxData["5"], "Output", CellChangeTimes->{ 3.58062064278522*^9, 3.580621334927236*^9, 3.5806222661708717`*^9, 3.5806222979481277`*^9, 3.5806223537026253`*^9, {3.5806223873362846`*^9, 3.5806224300959597`*^9}, {3.580622469189628*^9, 3.5806225184545145`*^9}, { 3.580622551089772*^9, 3.5806225808390245`*^9}, 3.580622715872862*^9, 3.58062280040941*^9, 3.580622873729539*^9, 3.580622909235201*^9, { 3.580622942229259*^9, 3.5806229670177026`*^9}, 3.580623010494979*^9, { 3.5806230637066727`*^9, 3.5806230934715247`*^9}, {3.5806231451076155`*^9, 3.580623172282863*^9}, {3.5806232037637186`*^9, 3.5806232241841545`*^9}, 3.5806233382515545`*^9, 3.580623395519255*^9, {3.58062342696891*^9, 3.580623474002993*^9}, 3.5806235367619033`*^9, {3.5806235939204035`*^9, 3.5806236229208546`*^9}, 3.5806236823881593`*^9, 3.580623718876623*^9, 3.5806237840691376`*^9, 3.580623816329994*^9, 3.5806238613984737`*^9, 3.5806239229561815`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"ow", "=", RowBox[{"Join", "[", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"-", "2.0"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", RowBox[{"bc", "+", "nb"}]}], "}"}]}], "]"}], ",", RowBox[{"Table", "[", RowBox[{ RowBox[{"-", "2.0"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", RowBox[{"2", "*", "hj"}]}], "}"}]}], "]"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{" ", "ow"}], "\[IndentingNewLine]", RowBox[{"cnv", "=", RowBox[{"Join", "[", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"cn", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "bc"}], "}"}]}], "]"}], ",", RowBox[{"Table", "[", RowBox[{"0.0", ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "cw"}], "}"}]}], "]"}], ",", RowBox[{"Table", "[", RowBox[{ RowBox[{"cn", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", RowBox[{"bc", "+", "1"}], ",", "n"}], "}"}]}], "]"}], ",", RowBox[{"Table", "[", RowBox[{"0.0", ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", RowBox[{"2", "*", "hj"}]}], "}"}]}], "]"}]}], "]"}]}]}], "Input", CellChangeTimes->{{3.580461670421983*^9, 3.5804617339140944`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"-", "2.`"}], ",", RowBox[{"-", "2.`"}], ",", RowBox[{"-", "2.`"}], ",", RowBox[{"-", "2.`"}], ",", RowBox[{"-", "2.`"}]}], "}"}]], "Output", CellChangeTimes->{ 3.5806206428008204`*^9, 3.580621334942836*^9, 3.580622266186472*^9, 3.580622297963728*^9, 3.5806223537182255`*^9, {3.5806223873518844`*^9, 3.58062243011156*^9}, {3.580622469205228*^9, 3.5806225184701147`*^9}, { 3.580622551105372*^9, 3.5806225808546247`*^9}, 3.5806227158884616`*^9, 3.5806228004250097`*^9, 3.580622873745139*^9, 3.580622909250801*^9, { 3.580622942244859*^9, 3.5806229670333023`*^9}, 3.580623010510579*^9, { 3.5806230637222724`*^9, 3.5806230934871244`*^9}, {3.5806231451232157`*^9, 3.580623172298463*^9}, {3.5806232037793183`*^9, 3.5806232241997547`*^9}, 3.5806233382671547`*^9, 3.580623395534855*^9, {3.5806234269845104`*^9, 3.5806234740185933`*^9}, 3.580623536777503*^9, {3.5806235939360037`*^9, 3.580623622936455*^9}, 3.580623682403759*^9, 3.5806237188922234`*^9, 3.5806237841003375`*^9, 3.5806238163455944`*^9, 3.580623861414074*^9, 3.5806239229717817`*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{"8.`", ",", "1.`", ",", "1.`", ",", "0.`", ",", "0.`"}], "}"}]], "Output", CellChangeTimes->{ 3.5806206428008204`*^9, 3.580621334942836*^9, 3.580622266186472*^9, 3.580622297963728*^9, 3.5806223537182255`*^9, {3.5806223873518844`*^9, 3.58062243011156*^9}, {3.580622469205228*^9, 3.5806225184701147`*^9}, { 3.580622551105372*^9, 3.5806225808546247`*^9}, 3.5806227158884616`*^9, 3.5806228004250097`*^9, 3.580622873745139*^9, 3.580622909250801*^9, { 3.580622942244859*^9, 3.5806229670333023`*^9}, 3.580623010510579*^9, { 3.5806230637222724`*^9, 3.5806230934871244`*^9}, {3.5806231451232157`*^9, 3.580623172298463*^9}, {3.5806232037793183`*^9, 3.5806232241997547`*^9}, 3.5806233382671547`*^9, 3.580623395534855*^9, {3.5806234269845104`*^9, 3.5806234740185933`*^9}, 3.580623536777503*^9, {3.5806235939360037`*^9, 3.580623622936455*^9}, 3.580623682403759*^9, 3.5806237188922234`*^9, 3.5806237841003375`*^9, 3.5806238163455944`*^9, 3.580623861414074*^9, 3.5806239229717817`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{"nucleus", "-", RowBox[{"cloud", " ", "distances"}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"dcw", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{ RowBox[{"sq", "=", RowBox[{"Norm", "[", RowBox[{ RowBox[{"tw", "[", RowBox[{"[", "i", "]"}], "]"}], "-", RowBox[{"tv", "[", RowBox[{"[", "j", "]"}], "]"}]}], "]"}]}], ";", RowBox[{"If", "[", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"sq", ">", " ", "0.1"}], " ", "&&", " ", RowBox[{"i", "\[NotEqual]", "j"}]}], ")"}], ",", RowBox[{"1", "/", "sq"}], ",", "0"}], "]"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "bt"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", "bt"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"Length", "[", "dcw", "]"}]}]}]], "Input", CellChangeTimes->{3.5804649749445868`*^9}], Cell[BoxData["5"], "Output", CellChangeTimes->{ 3.5806206428164206`*^9, 3.580621334958436*^9, 3.5806222662020717`*^9, 3.5806222979793277`*^9, 3.5806223537338257`*^9, {3.580622387367484*^9, 3.5806224301271596`*^9}, {3.5806224692364283`*^9, 3.580622518485715*^9}, { 3.5806225511209717`*^9, 3.5806225808702245`*^9}, 3.5806227159040613`*^9, 3.58062280044061*^9, 3.580622873760739*^9, 3.5806229092664013`*^9, { 3.580622942260459*^9, 3.5806229670489025`*^9}, 3.5806230105261793`*^9, { 3.580623063737872*^9, 3.5806230935027246`*^9}, {3.5806231451388154`*^9, 3.580623172314063*^9}, {3.580623203794918*^9, 3.5806232242153544`*^9}, 3.580623338282755*^9, 3.580623395550455*^9, {3.5806234270001106`*^9, 3.580623474034193*^9}, 3.580623536793103*^9, {3.580623593951604*^9, 3.580623622952055*^9}, 3.5806236824349594`*^9, 3.5806237189234233`*^9, 3.5806237841159377`*^9, 3.5806238163611946`*^9, 3.5806238614296737`*^9, 3.580623922987382*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{"Elcharge", "=", RowBox[{"Sum", "[", RowBox[{ RowBox[{"ow", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", RowBox[{"Length", "[", "ow", "]"}]}], "}"}]}], "]"}]}], "\[IndentingNewLine]", RowBox[{"Length", "[", "ow", "]"}], "\[IndentingNewLine]", RowBox[{"Nuccharge", "=", RowBox[{"Sum", "[", RowBox[{ RowBox[{"cnv", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", RowBox[{"Length", "[", "cnv", "]"}]}], "}"}]}], "]"}]}], "\[IndentingNewLine]", RowBox[{"Length", "[", "cnv", "]"}]}], "Input", CellChangeTimes->{3.580489056610564*^9}], Cell[BoxData[ RowBox[{"-", "10.`"}]], "Output", CellChangeTimes->{ 3.5806206428164206`*^9, 3.580621334958436*^9, 3.5806222662176714`*^9, 3.5806222979949274`*^9, 3.5806223537338257`*^9, {3.5806223873830843`*^9, 3.5806224301427593`*^9}, {3.5806224692364283`*^9, 3.580622518501315*^9}, { 3.580622551136572*^9, 3.580622580885824*^9}, 3.5806227159196615`*^9, 3.58062280045621*^9, 3.5806228737763386`*^9, 3.5806229092820015`*^9, { 3.580622942276059*^9, 3.5806229670645027`*^9}, 3.580623010541779*^9, { 3.5806230637534723`*^9, 3.580623093518325*^9}, {3.5806231451544156`*^9, 3.5806231723296633`*^9}, {3.5806232038105183`*^9, 3.580623224230954*^9}, 3.580623338298355*^9, 3.5806233955660553`*^9, {3.5806234270157104`*^9, 3.5806234740497932`*^9}, 3.5806235368087034`*^9, {3.580623593967204*^9, 3.5806236229676547`*^9}, 3.580623682450559*^9, 3.580623718939023*^9, 3.580623784131538*^9, 3.5806238163767943`*^9, 3.5806238614452734`*^9, 3.5806239230029817`*^9}], Cell[BoxData["5"], "Output", CellChangeTimes->{ 3.5806206428164206`*^9, 3.580621334958436*^9, 3.5806222662176714`*^9, 3.5806222979949274`*^9, 3.5806223537338257`*^9, {3.5806223873830843`*^9, 3.5806224301427593`*^9}, {3.5806224692364283`*^9, 3.580622518501315*^9}, { 3.580622551136572*^9, 3.580622580885824*^9}, 3.5806227159196615`*^9, 3.58062280045621*^9, 3.5806228737763386`*^9, 3.5806229092820015`*^9, { 3.580622942276059*^9, 3.5806229670645027`*^9}, 3.580623010541779*^9, { 3.5806230637534723`*^9, 3.580623093518325*^9}, {3.5806231451544156`*^9, 3.5806231723296633`*^9}, {3.5806232038105183`*^9, 3.580623224230954*^9}, 3.580623338298355*^9, 3.5806233955660553`*^9, {3.5806234270157104`*^9, 3.5806234740497932`*^9}, 3.5806235368087034`*^9, {3.580623593967204*^9, 3.5806236229676547`*^9}, 3.580623682450559*^9, 3.580623718939023*^9, 3.580623784131538*^9, 3.5806238163767943`*^9, 3.5806238614452734`*^9, 3.5806239230029817`*^9}], Cell[BoxData["10.`"], "Output", CellChangeTimes->{ 3.5806206428164206`*^9, 3.580621334958436*^9, 3.5806222662176714`*^9, 3.5806222979949274`*^9, 3.5806223537338257`*^9, {3.5806223873830843`*^9, 3.5806224301427593`*^9}, {3.5806224692364283`*^9, 3.580622518501315*^9}, { 3.580622551136572*^9, 3.580622580885824*^9}, 3.5806227159196615`*^9, 3.58062280045621*^9, 3.5806228737763386`*^9, 3.5806229092820015`*^9, { 3.580622942276059*^9, 3.5806229670645027`*^9}, 3.580623010541779*^9, { 3.5806230637534723`*^9, 3.580623093518325*^9}, {3.5806231451544156`*^9, 3.5806231723296633`*^9}, {3.5806232038105183`*^9, 3.580623224230954*^9}, 3.580623338298355*^9, 3.5806233955660553`*^9, {3.5806234270157104`*^9, 3.5806234740497932`*^9}, 3.5806235368087034`*^9, {3.580623593967204*^9, 3.5806236229676547`*^9}, 3.580623682450559*^9, 3.580623718939023*^9, 3.580623784131538*^9, 3.5806238163767943`*^9, 3.5806238614452734`*^9, 3.5806239230185814`*^9}], Cell[BoxData["5"], "Output", CellChangeTimes->{ 3.5806206428164206`*^9, 3.580621334958436*^9, 3.5806222662176714`*^9, 3.5806222979949274`*^9, 3.5806223537338257`*^9, {3.5806223873830843`*^9, 3.5806224301427593`*^9}, {3.5806224692364283`*^9, 3.580622518501315*^9}, { 3.580622551136572*^9, 3.580622580885824*^9}, 3.5806227159196615`*^9, 3.58062280045621*^9, 3.5806228737763386`*^9, 3.5806229092820015`*^9, { 3.580622942276059*^9, 3.5806229670645027`*^9}, 3.580623010541779*^9, { 3.5806230637534723`*^9, 3.580623093518325*^9}, {3.5806231451544156`*^9, 3.5806231723296633`*^9}, {3.5806232038105183`*^9, 3.580623224230954*^9}, 3.580623338298355*^9, 3.5806233955660553`*^9, {3.5806234270157104`*^9, 3.5806234740497932`*^9}, 3.5806235368087034`*^9, {3.580623593967204*^9, 3.5806236229676547`*^9}, 3.580623682450559*^9, 3.580623718939023*^9, 3.580623784131538*^9, 3.5806238163767943`*^9, 3.5806238614452734`*^9, 3.5806239230185814`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"Vee", "=", "0.0"}], ";"}], "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "<", "bt"}], ",", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"j", "=", RowBox[{"i", "+", "1"}]}], ",", RowBox[{"j", "\[LessEqual]", "bt"}], ",", RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", RowBox[{"Vee", "=", RowBox[{"Vee", "+", RowBox[{ RowBox[{"ow", "[", RowBox[{"[", "i", "]"}], "]"}], "*", RowBox[{ RowBox[{"ow", "[", RowBox[{"[", "j", "]"}], "]"}], "/", RowBox[{"dw", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}]}]}]}]}]}], "]"}]}], "]"}], "\[IndentingNewLine]", "Vee", "\[IndentingNewLine]", RowBox[{ RowBox[{"sc", "=", StyleBox["0.52", FontColor->RGBColor[1, 0, 0]]}], ";", " ", RowBox[{"(*", " ", RowBox[{"core", " ", "clouds"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], " ", ",", RowBox[{"i", "\[LessEqual]", "bc"}], ",", RowBox[{"i", "++"}], ",", RowBox[{"Vee", "=", RowBox[{"Vee", "+", RowBox[{"3", "*", RowBox[{"sc", "/", RowBox[{"u2", "[", RowBox[{"[", "i", "]"}], "]"}]}]}]}]}]}], "]"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"sw", "=", StyleBox["0.52", FontColor->RGBColor[1, 0, 0]]}], ";", " ", RowBox[{"(*", " ", RowBox[{"bonding", " ", "clouds"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], " ", ",", RowBox[{"i", "\[LessEqual]", "nb"}], ",", RowBox[{"i", "++"}], ",", RowBox[{"Vee", "=", RowBox[{"Vee", "+", RowBox[{"3", "*", RowBox[{"sw", "/", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "3"}], "]"}], "]"}]}]}]}]}]}], "]"}]}], "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], " ", ",", RowBox[{"i", "\[LessEqual]", "hj"}], ",", RowBox[{"i", "++"}], ",", RowBox[{"Vee", "=", RowBox[{"Vee", "+", RowBox[{"2", "*", "3", "*", RowBox[{"sw", "/", RowBox[{"rlos", "[", RowBox[{"[", "i", "]"}], "]"}]}]}]}]}]}], "]"}], "\[IndentingNewLine]", "Vee"}], "Input", CellChangeTimes->{{3.5804618141918354`*^9, 3.580461913283209*^9}, 3.5804623255451336`*^9, {3.5804628735428963`*^9, 3.5804628958041353`*^9}, 3.580476850943878*^9, {3.580477105068324*^9, 3.5804771083443303`*^9}, { 3.5804771721796427`*^9, 3.5804771783728533`*^9}, {3.5804772740010214`*^9, 3.5804772770274267`*^9}, {3.580482892684083*^9, 3.5804828962564898`*^9}, { 3.580482954194991*^9, 3.5804829592026*^9}, {3.5804831929846106`*^9, 3.580483196182616*^9}, {3.5804872622642126`*^9, 3.5804872631222143`*^9}, { 3.5804873217939167`*^9, 3.580487322527118*^9}, {3.5804873526819715`*^9, 3.5804873772676144`*^9}, {3.580487884424505*^9, 3.5804878847053056`*^9}, 3.5804890446921434`*^9, {3.580622257965257*^9, 3.580622261210063*^9}, { 3.5806223411134033`*^9, 3.580622378725069*^9}, {3.5806224544944024`*^9, 3.5806224929172697`*^9}, {3.5806225414489555`*^9, 3.5806225448029613`*^9}, {3.5806230023205643`*^9, 3.580623005799371*^9}, { 3.580623383101633*^9, 3.580623389170044*^9}, {3.580623584887988*^9, 3.580623618568447*^9}}], Cell[BoxData["23.154435013426014`"], "Output", CellChangeTimes->{ 3.5806206428476205`*^9, 3.580621334989636*^9, 3.5806222662332716`*^9, 3.5806222980105276`*^9, 3.5806223537650256`*^9, {3.5806223873986845`*^9, 3.5806224301583595`*^9}, {3.5806224692676287`*^9, 3.580622518516915*^9}, { 3.580622551152172*^9, 3.5806225809014244`*^9}, 3.5806227159352617`*^9, 3.5806228004718103`*^9, 3.580622873791939*^9, 3.580622909297601*^9, { 3.580622942307259*^9, 3.580622967080103*^9}, 3.5806230105573792`*^9, { 3.5806230637846727`*^9, 3.580623093533925*^9}, {3.5806231451856155`*^9, 3.5806231723452635`*^9}, {3.5806232038417187`*^9, 3.5806232242621546`*^9}, 3.580623338313955*^9, 3.5806233955816555`*^9, {3.5806234270313106`*^9, 3.580623474080993*^9}, 3.5806235368243036`*^9, {3.580623593982804*^9, 3.5806236229988546`*^9}, 3.580623682466159*^9, 3.580623718954623*^9, 3.5806237841471376`*^9, 3.5806238164079943`*^9, 3.580623861476474*^9, 3.5806239230341816`*^9}], Cell[BoxData["37.3907722702761`"], "Output", CellChangeTimes->{ 3.5806206428476205`*^9, 3.580621334989636*^9, 3.5806222662332716`*^9, 3.5806222980105276`*^9, 3.5806223537650256`*^9, {3.5806223873986845`*^9, 3.5806224301583595`*^9}, {3.5806224692676287`*^9, 3.580622518516915*^9}, { 3.580622551152172*^9, 3.5806225809014244`*^9}, 3.5806227159352617`*^9, 3.5806228004718103`*^9, 3.580622873791939*^9, 3.580622909297601*^9, { 3.580622942307259*^9, 3.580622967080103*^9}, 3.5806230105573792`*^9, { 3.5806230637846727`*^9, 3.580623093533925*^9}, {3.5806231451856155`*^9, 3.5806231723452635`*^9}, {3.5806232038417187`*^9, 3.5806232242621546`*^9}, 3.580623338313955*^9, 3.5806233955816555`*^9, {3.5806234270313106`*^9, 3.580623474080993*^9}, 3.5806235368243036`*^9, {3.580623593982804*^9, 3.5806236229988546`*^9}, 3.580623682466159*^9, 3.580623718954623*^9, 3.5806237841471376`*^9, 3.5806238164079943`*^9, 3.580623861476474*^9, 3.5806239230341816`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ RowBox[{"all", " ", "nuclei", " ", "with", " ", "all", " ", "clouds"}], ",", " ", RowBox[{ RowBox[{"except", " ", "i"}], " ", "=", " ", "j"}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vne", "=", RowBox[{"cnv", ".", RowBox[{"(", RowBox[{"ow", ".", "dcw"}], ")"}]}]}], "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "\[LessEqual]", "bc"}], ",", RowBox[{"i", "++"}], ",", RowBox[{"Vne", "=", RowBox[{"Vne", "-", RowBox[{ StyleBox["3.11", FontColor->RGBColor[1, 0, 0]], "*", RowBox[{ RowBox[{"cn", "[", RowBox[{"[", "i", "]"}], "]"}], "/", RowBox[{"u2", "[", RowBox[{"[", "i", "]"}], "]"}]}]}]}]}]}], "]"}], "\[IndentingNewLine]", "Vne"}]}]], "Input", CellChangeTimes->{{3.5804618141918354`*^9, 3.580461913283209*^9}, { 3.5804623255451336`*^9, 3.5804623320503445`*^9}, 3.5804623790064273`*^9, 3.5804648933564434`*^9, {3.580486684704398*^9, 3.5804866854531994`*^9}, 3.5804867335324836`*^9, {3.580486767821344*^9, 3.5804868108774195`*^9}, { 3.5804868414066734`*^9, 3.5804868626539106`*^9}, {3.580487002024555*^9, 3.5804870027109566`*^9}, 3.580487037218217*^9, {3.580487070789476*^9, 3.5804870718814783`*^9}, {3.5804874240520964`*^9, 3.580487452054146*^9}, { 3.580487501646633*^9, 3.5804875425967045`*^9}, {3.580487596245199*^9, 3.5804878336464157`*^9}, 3.5804879072785454`*^9, 3.5806228673959274`*^9, { 3.580622903821992*^9, 3.580622904133992*^9}, {3.5806229371124496`*^9, 3.580622962181694*^9}, {3.580623418685296*^9, 3.580623468745784*^9}, { 3.5806235308650928`*^9, 3.5806235316294947`*^9}, {3.5806236760545483`*^9, 3.580623676772149*^9}, 3.5806237132606134`*^9, {3.5806237775483265`*^9, 3.580623778655928*^9}, 3.5806238112911854`*^9, 3.5806238549244623`*^9}], Cell[BoxData[ RowBox[{"-", "55.0961802314351`"}]], "Output", CellChangeTimes->{ 3.5806206428476205`*^9, 3.5806213350052357`*^9, 3.580622266248872*^9, 3.580622298026128*^9, 3.5806223537806253`*^9, {3.5806223874142847`*^9, 3.58062243018956*^9}, {3.5806224692832284`*^9, 3.5806225185325146`*^9}, { 3.5806225511677723`*^9, 3.5806225809326243`*^9}, 3.580622715950862*^9, 3.58062280048741*^9, 3.580622873807539*^9, 3.580622909313201*^9, { 3.5806229423228593`*^9, 3.580622967095703*^9}, 3.580623010572979*^9, { 3.580623063800273*^9, 3.580623093549525*^9}, {3.5806231452012157`*^9, 3.580623172360863*^9}, {3.580623203857319*^9, 3.5806232242777543`*^9}, 3.580623338345155*^9, 3.5806233956128554`*^9, {3.5806234270625105`*^9, 3.580623474096593*^9}, 3.5806235368399034`*^9, {3.5806235939984035`*^9, 3.580623623014455*^9}, 3.580623682481759*^9, 3.5806237189702234`*^9, 3.580623784162738*^9, 3.5806238164235945`*^9, 3.5806238614920735`*^9, 3.580623923049782*^9}], Cell[BoxData[ RowBox[{"-", "193.9354659457208`"}]], "Output", CellChangeTimes->{ 3.5806206428476205`*^9, 3.5806213350052357`*^9, 3.580622266248872*^9, 3.580622298026128*^9, 3.5806223537806253`*^9, {3.5806223874142847`*^9, 3.58062243018956*^9}, {3.5806224692832284`*^9, 3.5806225185325146`*^9}, { 3.5806225511677723`*^9, 3.5806225809326243`*^9}, 3.580622715950862*^9, 3.58062280048741*^9, 3.580622873807539*^9, 3.580622909313201*^9, { 3.5806229423228593`*^9, 3.580622967095703*^9}, 3.580623010572979*^9, { 3.580623063800273*^9, 3.580623093549525*^9}, {3.5806231452012157`*^9, 3.580623172360863*^9}, {3.580623203857319*^9, 3.5806232242777543`*^9}, 3.580623338345155*^9, 3.5806233956128554`*^9, {3.5806234270625105`*^9, 3.580623474096593*^9}, 3.5806235368399034`*^9, {3.5806235939984035`*^9, 3.580623623014455*^9}, 3.580623682481759*^9, 3.5806237189702234`*^9, 3.580623784162738*^9, 3.5806238164235945`*^9, 3.5806238614920735`*^9, 3.580623923049782*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", RowBox[{"cw", "+", "1"}]}], ",", RowBox[{"i", "\[LessEqual]", "nb"}], ",", RowBox[{"i", "++"}], ",", RowBox[{"Vne", "=", RowBox[{"Vne", "-", RowBox[{ RowBox[{"(", RowBox[{"3", "-", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"u2", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "/", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "3"}], "]"}], "]"}]}], ")"}], "^", "2"}]}], ")"}], "/", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "3"}], "]"}], "]"}]}]}]}]}], "]"}], "\[IndentingNewLine]", "Vne"}], "Input", CellChangeTimes->{{3.580461986228937*^9, 3.5804620091453776`*^9}, { 3.5804624463517456`*^9, 3.580462451655755*^9}, 3.5804648450587587`*^9, 3.580464883887227*^9}], Cell[BoxData[ RowBox[{"-", "198.86522036387504`"}]], "Output", CellChangeTimes->{ 3.58062064286322*^9, 3.580621335020836*^9, 3.580622266264472*^9, 3.5806222980417275`*^9, 3.5806223537962255`*^9, {3.580622387429885*^9, 3.5806224302051597`*^9}, {3.580622469298828*^9, 3.580622518563715*^9}, { 3.5806225511833725`*^9, 3.580622580948224*^9}, 3.580622715950862*^9, 3.58062280050301*^9, 3.5806228738231387`*^9, 3.580622909328801*^9, { 3.580622942338459*^9, 3.580622967111303*^9}, 3.5806230106041794`*^9, { 3.5806230638158727`*^9, 3.5806230935651245`*^9}, {3.580623145216816*^9, 3.580623172376463*^9}, {3.5806232038729186`*^9, 3.5806232242933545`*^9}, 3.580623338360755*^9, 3.580623395628455*^9, {3.5806234270781107`*^9, 3.580623474112193*^9}, 3.5806235368711033`*^9, {3.580623594029604*^9, 3.580623623030055*^9}, 3.5806236825129595`*^9, 3.5806237190014234`*^9, 3.5806237841939383`*^9, 3.5806238164391947`*^9, 3.5806238615076733`*^9, 3.580623923065382*^9}] }, Open ]], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{"nucleus", "-", "nucleus"}], " ", "*)"}]], "Input", CellChangeTimes->{{3.580461986228937*^9, 3.5804620091453776`*^9}, { 3.5804624463517456`*^9, 3.580462451655755*^9}, {3.5804648450587587`*^9, 3.5804648790824184`*^9}}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"Vnn", "=", "0.0"}], ";"}], "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "<", "n"}], ",", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"j", "=", RowBox[{"i", "+", "1"}]}], ",", RowBox[{"j", "<", RowBox[{"n", "+", "1"}]}], ",", RowBox[{"j", "++"}], ",", RowBox[{"Vnn", "=", RowBox[{"Vnn", "+", RowBox[{ RowBox[{"cn", "[", RowBox[{"[", "i", "]"}], "]"}], "*", RowBox[{ RowBox[{"cn", "[", RowBox[{"[", "j", "]"}], "]"}], "/", RowBox[{"dc", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}]}]}]}]}]}], "]"}]}], "]"}], "\[IndentingNewLine]", "Vnn"}], "Input", CellChangeTimes->{{3.580461986228937*^9, 3.5804620091453776`*^9}, { 3.5804624463517456`*^9, 3.580462451655755*^9}, {3.5804648450587587`*^9, 3.580464866945597*^9}}], Cell[BoxData["9.086986115052769`"], "Output", CellChangeTimes->{ 3.5806206428788204`*^9, 3.580621335036436*^9, 3.5806222662800717`*^9, 3.5806222980573273`*^9, 3.5806223537962255`*^9, {3.5806223874454846`*^9, 3.58062243022076*^9}, {3.5806224693144283`*^9, 3.580622518563715*^9}, { 3.580622551198972*^9, 3.580622580948224*^9}, 3.5806227159664617`*^9, 3.58062280051861*^9, 3.580622873838739*^9, 3.5806229093444014`*^9, { 3.580622942354059*^9, 3.5806229671269026`*^9}, 3.5806230106041794`*^9, { 3.5806230638314724`*^9, 3.5806230935807247`*^9}, {3.5806231452324157`*^9, 3.580623172392063*^9}, {3.5806232038885183`*^9, 3.5806232243089542`*^9}, 3.5806233383763547`*^9, 3.5806233956440554`*^9, {3.580623427093711*^9, 3.5806234741277933`*^9}, 3.5806235368711033`*^9, {3.580623594045204*^9, 3.580623623045655*^9}, 3.580623682528559*^9, 3.580623719017023*^9, 3.580623784209538*^9, 3.5806238164547944`*^9, 3.5806238615232735`*^9, 3.5806239230809817`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"ratio", "=", RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"Vne", "+", RowBox[{"2", "*", "Vnn"}]}], ")"}]}], "/", RowBox[{"(", RowBox[{"Vee", "-", "Vnn"}], ")"}]}]}]], "Input", CellChangeTimes->{{3.580461986228937*^9, 3.5804620091453776`*^9}, { 3.5804624463517456`*^9, 3.580462451655755*^9}, {3.5804648450587587`*^9, 3.580464857866381*^9}}], Cell[BoxData["6.383995665555993`"], "Output", CellChangeTimes->{ 3.5806206428944206`*^9, 3.580621335036436*^9, 3.5806222662956715`*^9, 3.5806222980729275`*^9, 3.5806223538118258`*^9, {3.5806223874610844`*^9, 3.5806224302363596`*^9}, {3.5806224693300285`*^9, 3.5806225185793147`*^9}, {3.5806225512145724`*^9, 3.5806225809638243`*^9}, 3.5806227159820614`*^9, 3.58062280053421*^9, 3.5806228738543386`*^9, 3.580622909360001*^9, {3.5806229423696594`*^9, 3.580622967142503*^9}, 3.580623010619779*^9, {3.5806230638470726`*^9, 3.580623093596325*^9}, { 3.5806231452480154`*^9, 3.5806231724076633`*^9}, {3.5806232039041185`*^9, 3.5806232243245544`*^9}, 3.580623338391955*^9, 3.5806233956596556`*^9, { 3.580623427109311*^9, 3.5806234741433935`*^9}, 3.5806235368867035`*^9, { 3.580623594060804*^9, 3.5806236230612545`*^9}, 3.580623682528559*^9, 3.5806237190326233`*^9, 3.5806237842251377`*^9, 3.5806238164703946`*^9, 3.5806238615388737`*^9, 3.580623923096582*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"ratiomp2", "=", RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{ RowBox[{"-", "198.599284"}], "+", RowBox[{"2", "*", "9.0869862"}]}], ")"}]}], "/", RowBox[{"(", RowBox[{"37.3841105", "-", "9.0869862"}], ")"}]}]}]], "Input", CellChangeTimes->{{3.580461986228937*^9, 3.5804620091453776`*^9}, { 3.5804624463517456`*^9, 3.580462451655755*^9}, {3.5804648450587587`*^9, 3.580464857866381*^9}, {3.5804657856624107`*^9, 3.580465868716957*^9}, { 3.58052513033449*^9, 3.580525195386604*^9}, {3.580526224146565*^9, 3.5805263240959406`*^9}, {3.580533563696274*^9, 3.580533628046387*^9}, { 3.5805336590436416`*^9, 3.580533659995243*^9}, {3.5806213276576233`*^9, 3.5806213282504244`*^9}}], Cell[BoxData["6.376100613163721`"], "Output", CellChangeTimes->{ 3.5806206429100204`*^9, 3.5806213350520363`*^9, 3.5806222663112717`*^9, 3.5806222980729275`*^9, 3.580622353827426*^9, {3.5806223874766846`*^9, 3.5806224302363596`*^9}, {3.5806224693300285`*^9, 3.580622518594915*^9}, { 3.580622551230172*^9, 3.5806225809794245`*^9}, 3.5806227159976616`*^9, 3.5806228005498104`*^9, 3.580622873869939*^9, 3.5806229093756013`*^9, { 3.580622942385259*^9, 3.580622967158103*^9}, 3.5806230106353793`*^9, { 3.580623063862673*^9, 3.580623093611925*^9}, {3.5806231452636156`*^9, 3.5806231724232635`*^9}, {3.5806232039197187`*^9, 3.5806232243401546`*^9}, 3.5806233384075546`*^9, 3.580623395675256*^9, {3.580623427124911*^9, 3.5806234741589932`*^9}, 3.5806235369023037`*^9, {3.580623594076404*^9, 3.5806236230768547`*^9}, 3.580623682544159*^9, 3.5806237190482235`*^9, 3.580623784240738*^9, 3.5806238164859943`*^9, 3.580623861554474*^9, 3.5806239231121817`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"\[IndentingNewLine]", RowBox[{ "Vee", "\[IndentingNewLine]", "Vnn", "\[IndentingNewLine]", "Vne", "\[IndentingNewLine]", RowBox[{"Tlp", " ", "=", RowBox[{ StyleBox["0.5482", FontColor->RGBColor[1, 0, 0]], "*", "2.25", "*", RowBox[{"(", RowBox[{"Sum", "[", RowBox[{ RowBox[{"2", "/", RowBox[{ RowBox[{"rlos", "[", RowBox[{"[", "i", "]"}], "]"}], "^", "2"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "hj"}], "}"}]}], "]"}], ")"}]}]}]}]}]], "Input", CellChangeTimes->{{3.5804620346046224`*^9, 3.5804620359462247`*^9}, { 3.580462097113932*^9, 3.5804620981435337`*^9}, 3.580462151105627*^9, { 3.5804763243493533`*^9, 3.580476329840563*^9}, {3.580476444937565*^9, 3.580476458634389*^9}, {3.58047649332885*^9, 3.580476496043255*^9}, 3.5804767076888266`*^9, {3.58047745555414*^9, 3.5804774557413406`*^9}, 3.580483254167918*^9, {3.5804834637194357`*^9, 3.580483464452637*^9}, { 3.580486296076515*^9, 3.5804864867556505`*^9}, {3.5804868827935457`*^9, 3.5804869648652897`*^9}, {3.580487096810322*^9, 3.5804871692724495`*^9}, { 3.5804886599642673`*^9, 3.5804886611810694`*^9}, {3.580488733908397*^9, 3.5804887603348436`*^9}, {3.5804887932197013`*^9, 3.5804888678658323`*^9}, {3.580488905711499*^9, 3.5804889063667*^9}, 3.5806222929561186`*^9, {3.580622404231114*^9, 3.580622425166351*^9}, 3.580622513587306*^9, {3.5806225751606145`*^9, 3.5806225758626156`*^9}, { 3.580623056265459*^9, 3.5806230569674606`*^9}, 3.580623087886715*^9, { 3.580623139132805*^9, 3.5806231663236523`*^9}, {3.580623197944908*^9, 3.580623217819343*^9}, {3.5806233326979446`*^9, 3.5806233335871463`*^9}}], Cell[BoxData["37.3907722702761`"], "Output", CellChangeTimes->{ 3.58062064292562*^9, 3.5806213350676365`*^9, 3.5806222663112717`*^9, 3.5806222980885277`*^9, 3.5806223538430257`*^9, {3.580622387492285*^9, 3.58062243025196*^9}, {3.5806224693456287`*^9, 3.580622518610515*^9}, { 3.5806225512457724`*^9, 3.5806225809950247`*^9}, 3.580622716013262*^9, 3.58062280056541*^9, 3.580622873885539*^9, 3.5806229093912015`*^9, { 3.580622942385259*^9, 3.5806229671737027`*^9}, 3.580623010650979*^9, { 3.580623063878273*^9, 3.5806230936275253`*^9}, {3.580623145279216*^9, 3.5806231724388638`*^9}, {3.580623203935319*^9, 3.580623224355755*^9}, 3.580623338423155*^9, 3.5806233956908555`*^9, {3.580623427124911*^9, 3.5806234741745934`*^9}, 3.580623536917904*^9, {3.580623594092004*^9, 3.580623623092455*^9}, 3.580623682559759*^9, 3.5806237190482235`*^9, 3.580623784256338*^9, 3.5806238165015945`*^9, 3.580623861570074*^9, 3.580623923127782*^9}], Cell[BoxData["9.086986115052769`"], "Output", CellChangeTimes->{ 3.58062064292562*^9, 3.5806213350676365`*^9, 3.5806222663112717`*^9, 3.5806222980885277`*^9, 3.5806223538430257`*^9, {3.580622387492285*^9, 3.58062243025196*^9}, {3.5806224693456287`*^9, 3.580622518610515*^9}, { 3.5806225512457724`*^9, 3.5806225809950247`*^9}, 3.580622716013262*^9, 3.58062280056541*^9, 3.580622873885539*^9, 3.5806229093912015`*^9, { 3.580622942385259*^9, 3.5806229671737027`*^9}, 3.580623010650979*^9, { 3.580623063878273*^9, 3.5806230936275253`*^9}, {3.580623145279216*^9, 3.5806231724388638`*^9}, {3.580623203935319*^9, 3.580623224355755*^9}, 3.580623338423155*^9, 3.5806233956908555`*^9, {3.580623427124911*^9, 3.5806234741745934`*^9}, 3.580623536917904*^9, {3.580623594092004*^9, 3.580623623092455*^9}, 3.580623682559759*^9, 3.5806237190482235`*^9, 3.580623784256338*^9, 3.5806238165015945`*^9, 3.580623861570074*^9, 3.580623923127782*^9}], Cell[BoxData[ RowBox[{"-", "198.86522036387504`"}]], "Output", CellChangeTimes->{ 3.58062064292562*^9, 3.5806213350676365`*^9, 3.5806222663112717`*^9, 3.5806222980885277`*^9, 3.5806223538430257`*^9, {3.580622387492285*^9, 3.58062243025196*^9}, {3.5806224693456287`*^9, 3.580622518610515*^9}, { 3.5806225512457724`*^9, 3.5806225809950247`*^9}, 3.580622716013262*^9, 3.58062280056541*^9, 3.580622873885539*^9, 3.5806229093912015`*^9, { 3.580622942385259*^9, 3.5806229671737027`*^9}, 3.580623010650979*^9, { 3.580623063878273*^9, 3.5806230936275253`*^9}, {3.580623145279216*^9, 3.5806231724388638`*^9}, {3.580623203935319*^9, 3.580623224355755*^9}, 3.580623338423155*^9, 3.5806233956908555`*^9, {3.580623427124911*^9, 3.5806234741745934`*^9}, 3.580623536917904*^9, {3.580623594092004*^9, 3.580623623092455*^9}, 3.580623682559759*^9, 3.5806237190482235`*^9, 3.580623784256338*^9, 3.5806238165015945`*^9, 3.580623861570074*^9, 3.580623923127782*^9}], Cell[BoxData["1.9491555555555553`"], "Output", CellChangeTimes->{ 3.58062064292562*^9, 3.5806213350676365`*^9, 3.5806222663112717`*^9, 3.5806222980885277`*^9, 3.5806223538430257`*^9, {3.580622387492285*^9, 3.58062243025196*^9}, {3.5806224693456287`*^9, 3.580622518610515*^9}, { 3.5806225512457724`*^9, 3.5806225809950247`*^9}, 3.580622716013262*^9, 3.58062280056541*^9, 3.580622873885539*^9, 3.5806229093912015`*^9, { 3.580622942385259*^9, 3.5806229671737027`*^9}, 3.580623010650979*^9, { 3.580623063878273*^9, 3.5806230936275253`*^9}, {3.580623145279216*^9, 3.5806231724388638`*^9}, {3.580623203935319*^9, 3.580623224355755*^9}, 3.580623338423155*^9, 3.5806233956908555`*^9, {3.580623427124911*^9, 3.5806234741745934`*^9}, 3.580623536917904*^9, {3.580623594092004*^9, 3.580623623092455*^9}, 3.580623682559759*^9, 3.5806237190482235`*^9, 3.580623784256338*^9, 3.5806238165015945`*^9, 3.580623861570074*^9, 3.580623923143382*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{"Vtot", "=", RowBox[{"Vne", "+", "Vee", "+", "Vnn"}]}], "\[IndentingNewLine]", RowBox[{"Ttot", " ", "=", RowBox[{"Ttot", "+", "Tlp"}]}], "\[IndentingNewLine]", RowBox[{"Vir", " ", "=", RowBox[{ RowBox[{"-", "Vtot"}], "/", "Ttot"}]}], "\[IndentingNewLine]", RowBox[{"Etot", "=", RowBox[{"Ttot", "+", "Vtot"}]}]}], "Input", CellChangeTimes->{ 3.580462058613065*^9, {3.5804621657228527`*^9, 3.580462173351266*^9}}], Cell[BoxData[ RowBox[{"-", "152.38746197854616`"}]], "Output", CellChangeTimes->{ 3.5806206429412203`*^9, 3.580621335083236*^9, 3.580622266326872*^9, 3.580622298104128*^9, 3.5806223538586254`*^9, {3.580622387507885*^9, 3.58062243028316*^9}, {3.5806224693612285`*^9, 3.580622518641715*^9}, { 3.5806225512613726`*^9, 3.580622581010625*^9}, 3.5806227160288615`*^9, 3.58062280059661*^9, 3.5806228739011393`*^9, 3.580622909406801*^9, { 3.580622942416459*^9, 3.5806229671893024`*^9}, 3.5806230106821795`*^9, { 3.5806230638938727`*^9, 3.580623093643125*^9}, {3.5806231452948155`*^9, 3.580623172470063*^9}, {3.580623203950919*^9, 3.5806232243713546`*^9}, 3.580623338438755*^9, 3.5806233957064557`*^9, {3.580623427156111*^9, 3.580623474190193*^9}, 3.5806235369491034`*^9, {3.580623594107604*^9, 3.580623623108055*^9}, 3.5806236825909595`*^9, 3.580623719079424*^9, 3.5806237842719383`*^9, 3.580623816532795*^9, 3.580623861585674*^9, 3.5806239231589823`*^9}], Cell[BoxData["76.19180657565745`"], "Output", CellChangeTimes->{ 3.5806206429412203`*^9, 3.580621335083236*^9, 3.580622266326872*^9, 3.580622298104128*^9, 3.5806223538586254`*^9, {3.580622387507885*^9, 3.58062243028316*^9}, {3.5806224693612285`*^9, 3.580622518641715*^9}, { 3.5806225512613726`*^9, 3.580622581010625*^9}, 3.5806227160288615`*^9, 3.58062280059661*^9, 3.5806228739011393`*^9, 3.580622909406801*^9, { 3.580622942416459*^9, 3.5806229671893024`*^9}, 3.5806230106821795`*^9, { 3.5806230638938727`*^9, 3.580623093643125*^9}, {3.5806231452948155`*^9, 3.580623172470063*^9}, {3.580623203950919*^9, 3.5806232243713546`*^9}, 3.580623338438755*^9, 3.5806233957064557`*^9, {3.580623427156111*^9, 3.580623474190193*^9}, 3.5806235369491034`*^9, {3.580623594107604*^9, 3.580623623108055*^9}, 3.5806236825909595`*^9, 3.580623719079424*^9, 3.5806237842719383`*^9, 3.580623816532795*^9, 3.580623861585674*^9, 3.5806239231589823`*^9}], Cell[BoxData["2.0000505149753525`"], "Output", CellChangeTimes->{ 3.5806206429412203`*^9, 3.580621335083236*^9, 3.580622266326872*^9, 3.580622298104128*^9, 3.5806223538586254`*^9, {3.580622387507885*^9, 3.58062243028316*^9}, {3.5806224693612285`*^9, 3.580622518641715*^9}, { 3.5806225512613726`*^9, 3.580622581010625*^9}, 3.5806227160288615`*^9, 3.58062280059661*^9, 3.5806228739011393`*^9, 3.580622909406801*^9, { 3.580622942416459*^9, 3.5806229671893024`*^9}, 3.5806230106821795`*^9, { 3.5806230638938727`*^9, 3.580623093643125*^9}, {3.5806231452948155`*^9, 3.580623172470063*^9}, {3.580623203950919*^9, 3.5806232243713546`*^9}, 3.580623338438755*^9, 3.5806233957064557`*^9, {3.580623427156111*^9, 3.580623474190193*^9}, 3.5806235369491034`*^9, {3.580623594107604*^9, 3.580623623108055*^9}, 3.5806236825909595`*^9, 3.580623719079424*^9, 3.5806237842719383`*^9, 3.580623816532795*^9, 3.580623861585674*^9, 3.5806239231589823`*^9}], Cell[BoxData[ RowBox[{"-", "76.1956554028887`"}]], "Output", CellChangeTimes->{ 3.5806206429412203`*^9, 3.580621335083236*^9, 3.580622266326872*^9, 3.580622298104128*^9, 3.5806223538586254`*^9, {3.580622387507885*^9, 3.58062243028316*^9}, {3.5806224693612285`*^9, 3.580622518641715*^9}, { 3.5806225512613726`*^9, 3.580622581010625*^9}, 3.5806227160288615`*^9, 3.58062280059661*^9, 3.5806228739011393`*^9, 3.580622909406801*^9, { 3.580622942416459*^9, 3.5806229671893024`*^9}, 3.5806230106821795`*^9, { 3.5806230638938727`*^9, 3.580623093643125*^9}, {3.5806231452948155`*^9, 3.580623172470063*^9}, {3.580623203950919*^9, 3.5806232243713546`*^9}, 3.580623338438755*^9, 3.5806233957064557`*^9, {3.580623427156111*^9, 3.580623474190193*^9}, 3.5806235369491034`*^9, {3.580623594107604*^9, 3.580623623108055*^9}, 3.5806236825909595`*^9, 3.580623719079424*^9, 3.5806237842719383`*^9, 3.580623816532795*^9, 3.580623861585674*^9, 3.5806239231589823`*^9}] }, Open ]] }, WindowSize->{850, 881}, WindowMargins->{{Automatic, 153}, {Automatic, 22}}, PrintingCopies->1, PrintingPageRange->{32000, 32000}, PrintingOptions->{"Magnification"->1., "PaperOrientation"->"Portrait", "PaperSize"->{595.2, 841.8}}, Magnification->0.8999999761581421, FrontEndVersion->"9.0 for Microsoft Windows (64-bit) (January 25, 2013)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[557, 20, 1198, 25, 50, "Input"], Cell[1758, 47, 1239, 25, 50, "Input"], Cell[3000, 74, 1018, 17, 29, "Input"], Cell[CellGroupData[{ Cell[4043, 95, 5722, 155, 350, "Input"], Cell[9768, 252, 958, 14, 28, "Output"], Cell[10729, 268, 1025, 16, 28, "Output"], Cell[11757, 286, 1025, 16, 28, "Output"], Cell[12785, 304, 1019, 16, 28, "Output"], Cell[13807, 322, 958, 14, 28, "Output"], Cell[14768, 338, 958, 14, 28, "Output"] }, Open ]], Cell[15741, 355, 5105, 104, 190, "Input"], Cell[20849, 461, 3790, 74, 70, "Input"], Cell[24642, 537, 1254, 35, 70, "Input"], Cell[CellGroupData[{ Cell[25921, 576, 3042, 63, 70, "Input"], Cell[28966, 641, 1684, 33, 82, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[30687, 679, 3060, 74, 130, "Input"], Cell[33750, 755, 1036, 16, 28, "Output"] }, Open ]], Cell[34801, 774, 1708, 43, 90, "Input"], Cell[36512, 819, 511, 11, 50, "Input"], Cell[37026, 832, 2048, 48, 90, "Input"], Cell[39077, 882, 2106, 50, 90, "Input"], Cell[41186, 934, 1261, 33, 50, "Input"], Cell[CellGroupData[{ Cell[42472, 971, 2727, 71, 130, "Input"], Cell[45202, 1044, 985, 14, 28, "Output"], Cell[46190, 1060, 985, 14, 28, "Output"], Cell[47178, 1076, 982, 14, 28, "Output"], Cell[48163, 1092, 983, 14, 28, "Output"], Cell[49149, 1108, 968, 14, 28, "Output"] }, Open ]], Cell[50132, 1125, 1971, 35, 90, "Input"], Cell[52106, 1162, 2207, 39, 29, "Input"], Cell[CellGroupData[{ Cell[54338, 1205, 1371, 43, 150, "Input"], Cell[55712, 1250, 963, 14, 28, "Output"], Cell[56678, 1266, 963, 14, 28, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[57678, 1285, 883, 26, 70, "Input"], Cell[58564, 1313, 976, 15, 28, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[59577, 1333, 3081, 94, 210, "Input"], Cell[62661, 1429, 987, 15, 28, "Output"], Cell[63651, 1446, 970, 14, 28, "Output"], Cell[64624, 1462, 1582, 32, 67, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[66243, 1499, 2861, 78, 210, "Input"], Cell[69107, 1579, 965, 14, 28, "Output"], Cell[70075, 1595, 965, 14, 28, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[71077, 1614, 3221, 87, 170, "Input"], Cell[74301, 1703, 1039, 17, 28, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[75377, 1725, 286, 4, 29, "Input"], Cell[75666, 1731, 9022, 159, 152, "Output"] }, Open ]], Cell[84703, 1893, 1655, 24, 29, "Input"], Cell[86361, 1919, 2724, 54, 90, "Input"], Cell[CellGroupData[{ Cell[89110, 1977, 1037, 19, 70, "Input"], Cell[90150, 1998, 1262, 23, 28, "Output"], Cell[91415, 2023, 973, 14, 28, "Output"], Cell[92391, 2039, 973, 14, 28, "Output"] }, Open ]], Cell[93379, 2056, 11150, 255, 210, "Input"], Cell[104532, 2313, 9944, 229, 190, "Input"], Cell[114479, 2544, 8783, 179, 170, "Input"], Cell[CellGroupData[{ Cell[123287, 2727, 5231, 75, 50, "Input"], Cell[128521, 2804, 2757, 59, 337, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[131315, 2868, 2228, 65, 110, "Input"], Cell[133546, 2935, 987, 14, 28, "Output"] }, Open ]], Cell[134548, 2952, 1001, 26, 70, "Input"], Cell[CellGroupData[{ Cell[135574, 2982, 1905, 55, 170, "Input"], Cell[137482, 3039, 970, 14, 28, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[138489, 3058, 1430, 44, 90, "Input"], Cell[139922, 3104, 1152, 21, 28, "Output"], Cell[141077, 3127, 1060, 17, 28, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[142174, 3149, 1082, 31, 70, "Input"], Cell[143259, 3182, 972, 14, 28, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[144268, 3201, 711, 21, 90, "Input"], Cell[144982, 3224, 996, 15, 28, "Output"], Cell[145981, 3241, 976, 14, 28, "Output"], Cell[146960, 3257, 979, 14, 28, "Output"], Cell[147942, 3273, 976, 14, 28, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[148955, 3292, 3487, 89, 230, "Input"], Cell[152445, 3383, 995, 14, 28, "Output"], Cell[153443, 3399, 993, 14, 28, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[154473, 3418, 1988, 43, 90, "Input"], Cell[156464, 3463, 998, 15, 28, "Output"], Cell[157465, 3480, 999, 15, 28, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[158501, 3500, 1016, 30, 50, "Input"], Cell[159520, 3532, 998, 15, 28, "Output"] }, Open ]], Cell[160533, 3550, 268, 5, 29, "Input"], Cell[CellGroupData[{ Cell[160826, 3559, 1016, 29, 90, "Input"], Cell[161845, 3590, 990, 14, 28, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[162872, 3609, 400, 11, 29, "Input"], Cell[163275, 3622, 999, 14, 28, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[164311, 3641, 748, 16, 29, "Input"], Cell[165062, 3659, 994, 14, 28, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[166093, 3678, 1767, 34, 110, "Input"], Cell[167863, 3714, 974, 14, 28, "Output"], Cell[168840, 3730, 975, 14, 28, "Output"], Cell[169818, 3746, 993, 15, 28, "Output"], Cell[170814, 3763, 976, 14, 28, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[171827, 3782, 461, 11, 90, "Input"], Cell[172291, 3795, 996, 15, 28, "Output"], Cell[173290, 3812, 978, 14, 28, "Output"], Cell[174271, 3828, 979, 14, 28, "Output"], Cell[175253, 3844, 994, 15, 60, "Output"] }, Open ]] } ] *) (* End of internal cache information *)