(* 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[ 206509, 4333] NotebookOptionsPosition[ 200501, 4141] NotebookOutlinePosition[ 201036, 4162] CellTagsIndexPosition[ 200993, 4159] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ RowBox[{ "Kimball", " ", "model", " ", "and", " ", "energy", " ", "derived", " ", "from", " ", "experimental", " ", "molecule", " ", "data", " ", "source"}], ",", " ", RowBox[{"ES", " ", "Febuary", " ", "2013"}]}], " ", "*)"}], " ", "\[IndentingNewLine]", RowBox[{ RowBox[{"Clear", "[", RowBox[{ "tw", ",", "n2", ",", "n", ",", "nb", ",", "bc", ",", "bh", ",", "u", ",", "u1", ",", "u2"}], "]"}], ";"}]}]], "Input", CellChangeTimes->{{3.568278170270548*^9, 3.568278174310955*^9}, { 3.5682782109554195`*^9, 3.5682782183342323`*^9}, {3.5688220426984916`*^9, 3.5688220506389055`*^9}, {3.568822092353379*^9, 3.568822134660653*^9}, { 3.5688728241768503`*^9, 3.568872825034852*^9}, {3.5688958900531473`*^9, 3.5688958965115585`*^9}, {3.5689027504159966`*^9, 3.568902751695199*^9}, { 3.5689109808408637`*^9, 3.568910992618884*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"ChemicalData", "[", RowBox[{"\"\\"", ",", "\"\\""}], "]"}]], "Input", CellChangeTimes->{{3.568278170270548*^9, 3.568278174310955*^9}, { 3.5682782109554195`*^9, 3.5682782183342323`*^9}, {3.5688220426984916`*^9, 3.5688220506389055`*^9}, {3.568822092353379*^9, 3.568822123693834*^9}, 3.568873215924739*^9, 3.5689108676002645`*^9, 3.5689997098237543`*^9}], Cell[BoxData[ Graphics3DBox[ {AbsoluteThickness[3], EdgeForm[None], Specularity[ GrayLevel[1], 100], GraphicsComplex3DBox[{{0.0115361, -0.00644059, -0.0023318}, {-6.98724, 109.314, 19.6968}, {-53.4044, -24.1997, -94.6041}, {-46.6211, -55.8814, 84.2109}, {107.001, -29.2263, -9.30131}, {-3.48785195, 54.653779705, 9.8472341}, {-26.69643195, -12.103070295, -47.303215900000005`}, \ {-23.30478195, -27.943920295, 42.104284099999994`}, { 53.50626805, -14.616370295, -4.651820900000001}}, { {RGBColor[0.65, 0.7, 0.7], SphereBox[2, 24.], SphereBox[3, 24.], SphereBox[4, 24.], SphereBox[5, 24.]}, {RGBColor[0.4, 0.4, 0.4], SphereBox[1, 34.]}, {RGBColor[0.65, 0.7, 0.7], CylinderBox[{6, 2}, 15.], CylinderBox[{7, 3}, 15.], CylinderBox[{8, 4}, 15.], CylinderBox[{9, 5}, 15.]}, {RGBColor[0.4, 0.4, 0.4], CylinderBox[{1, 6}, 15.], CylinderBox[{1, 7}, 15.], CylinderBox[{1, 8}, 15.], CylinderBox[{1, 9}, 15.]}}]}, Boxed->False, ImageSize->132.91425168924724`, Lighting->"Neutral", ViewPoint->{1.3, -2.4, -2.}]], "Output", CellChangeTimes->{ 3.5689997293861885`*^9, 3.569060092367221*^9, 3.569065230094085*^9, 3.5690652766601667`*^9, {3.5690653768591433`*^9, 3.5690654060935946`*^9}, 3.569065561766268*^9, 3.5690656435572114`*^9, 3.5690658373407516`*^9, 3.569066361954073*^9, 3.569066404417348*^9, {3.5690665341315756`*^9, 3.569066555410013*^9}, {3.5690666161721196`*^9, 3.5690666340497513`*^9}, { 3.5690666757954245`*^9, 3.5690667053106766`*^9}, {3.5690667504571557`*^9, 3.5690667705187907`*^9}, {3.5690668202516785`*^9, 3.5690668470681252`*^9}, 3.569066880467784*^9, 3.5690694994659567`*^9, 3.5690695660780735`*^9, 3.569069721501147*^9, 3.569069809625702*^9, 3.569069854959381*^9, 3.56906994540834*^9, 3.569070175571144*^9, 3.5690702503420753`*^9, 3.569070440116409*^9, 3.56907076525218*^9, 3.569071524966315*^9, 3.5690715823744154`*^9}] }, Open ]], Cell[BoxData[{ RowBox[{ RowBox[{"t1", "=", RowBox[{"0.01889726", "*", "0.97956684", "*", RowBox[{"ChemicalData", "[", RowBox[{"\"\\"", ",", "\"\\""}], "]"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"(*", RowBox[{"MatrixForm", "[", "t1", "]"}], "*)"}]}]}], "Input", CellChangeTimes->{{3.567862912170579*^9, 3.567862949329844*^9}, { 3.5678630095303497`*^9, 3.5678630107471523`*^9}, {3.5679573013548217`*^9, 3.5679573132576427`*^9}, 3.568114367533684*^9, {3.5682782383334675`*^9, 3.5682782435438766`*^9}, {3.568381606879155*^9, 3.5683816154123707`*^9}, 3.5688204675481253`*^9, {3.5688215945097046`*^9, 3.568821600578115*^9}, 3.568873215924739*^9, 3.5689108676002645`*^9, 3.5689997098393545`*^9, { 3.569065470786908*^9, 3.5690654833917303`*^9}, 3.569065515652587*^9, { 3.5690656236047764`*^9, 3.5690656379100018`*^9}}], Cell[BoxData[ RowBox[{"(*", RowBox[{"ChemicalData", "\[NoBreak]", "[", RowBox[{"\"\\"", ",", "\"\\""}], "]"}], "*)"}]], "Input", CellChangeTimes->{{3.5678630319943895`*^9, 3.5678630956269007`*^9}, { 3.5678632280087337`*^9, 3.567863228336334*^9}, {3.5678632621727934`*^9, 3.5678633509213495`*^9}, {3.5678644293824434`*^9, 3.567864451908883*^9}, { 3.567864923404111*^9, 3.56786494544695*^9}, {3.5678650143210707`*^9, 3.5678650845367947`*^9}, 3.567869231084096*^9, 3.5678693210962543`*^9, { 3.567869354027912*^9, 3.5678693684423375`*^9}, {3.5682782542610955`*^9, 3.5682782564138994`*^9}, {3.568821613806938*^9, 3.568821623151355*^9}, 3.5688732159403386`*^9, 3.5689108676158648`*^9, 3.5689997098393545`*^9}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"u", "=", RowBox[{"ChemicalData", "\[NoBreak]", "[", RowBox[{"\"\\"", ",", "\"\\""}], "]"}]}]], "Input", CellChangeTimes->{{3.567869377256353*^9, 3.5678694101568108`*^9}, { 3.5678695244894114`*^9, 3.56786953511303*^9}, 3.56786988243204*^9, 3.5678699252853155`*^9, {3.567870069179968*^9, 3.5678700766835814`*^9}, { 3.5678703216664114`*^9, 3.567870325738019*^9}, 3.5679271251072817`*^9, { 3.567942840466694*^9, 3.567942854381918*^9}, {3.567942932272855*^9, 3.567942950509287*^9}, {3.567943087867528*^9, 3.567943093733138*^9}, 3.567943185211699*^9, 3.5679432253505697`*^9, 3.5679433277647495`*^9, { 3.5679434430645523`*^9, 3.5679434440941534`*^9}, {3.568278260766307*^9, 3.5682782624355097`*^9}, 3.5688732159403386`*^9, 3.5689108676158648`*^9, 3.5689997098393545`*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{"\<\"C\"\>", ",", "\<\"H\"\>", ",", "\<\"H\"\>", ",", "\<\"H\"\>", ",", "\<\"H\"\>"}], "}"}]], "Output", CellChangeTimes->{ 3.568999729495389*^9, 3.569060092382821*^9, 3.569065230109685*^9, 3.569065276675767*^9, {3.569065376874743*^9, 3.5690654061091948`*^9}, 3.569065561813068*^9, 3.5690656435884113`*^9, 3.569065837371952*^9, 3.569066361969673*^9, 3.569066404432948*^9, {3.569066534147176*^9, 3.569066555425613*^9}, {3.56906661618772*^9, 3.5690666340653515`*^9}, { 3.5690666758266244`*^9, 3.5690667053418765`*^9}, {3.5690667504727554`*^9, 3.5690667705343904`*^9}, {3.569066820267278*^9, 3.569066847099325*^9}, 3.5690668804989843`*^9, 3.569069499481557*^9, 3.569069566109274*^9, 3.5690697215167465`*^9, 3.5690698096413016`*^9, 3.5690698549749813`*^9, 3.5690699454239397`*^9, 3.5690701755867443`*^9, 3.5690702503576756`*^9, 3.5690704401320086`*^9, 3.5690707652677803`*^9, 3.5690715249819145`*^9, 3.5690715823900156`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"n", "=", RowBox[{"Length", "[", "u", "]"}]}]], "Input", CellChangeTimes->{{3.568277663590458*^9, 3.568277712652544*^9}, 3.5682785838808746`*^9, {3.56880682872497*^9, 3.5688068517038097`*^9}, 3.5688069726508226`*^9}], Cell[BoxData["5"], "Output", CellChangeTimes->{ 3.5689997295109887`*^9, 3.569060092398421*^9, 3.569065230125285*^9, 3.569065276691367*^9, {3.569065376890343*^9, 3.5690654061247945`*^9}, 3.569065561828668*^9, 3.5690656436040115`*^9, 3.5690658373875523`*^9, 3.5690663619852734`*^9, 3.5690664044485483`*^9, {3.5690665341627755`*^9, 3.569066555441213*^9}, {3.5690666162033195`*^9, 3.569066634080951*^9}, { 3.569066675842224*^9, 3.5690667053418765`*^9}, {3.5690667504883556`*^9, 3.5690667705499907`*^9}, {3.569066820282878*^9, 3.5690668471149254`*^9}, 3.569066880514584*^9, 3.569069499497157*^9, 3.569069566124874*^9, 3.5690697215323467`*^9, 3.5690698096569014`*^9, 3.5690698549905815`*^9, 3.56906994543954*^9, 3.5690701756023445`*^9, 3.5690702503732758`*^9, 3.5690704401476088`*^9, 3.56907076528338*^9, 3.5690715249975147`*^9, 3.569071582405616*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"un", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\\""}], ",", "\"\\"", ",", RowBox[{"u", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.568807971988578*^9, 3.5688080897687845`*^9}, { 3.5688081769573374`*^9, 3.5688081865669546`*^9}, {3.568808274317109*^9, 3.5688082853151283`*^9}, {3.5688085365847692`*^9, 3.568808578330443*^9}, { 3.568808822704872*^9, 3.5688088761037655`*^9}, {3.56880890727262*^9, 3.568808937365073*^9}, {3.568809097608555*^9, 3.5688091594002633`*^9}, 3.5688092418308077`*^9, {3.5688092753552666`*^9, 3.568809302577315*^9}, { 3.5688096257786827`*^9, 3.568809703201618*^9}, {3.5688097456336927`*^9, 3.568809788268568*^9}, 3.568809850575077*^9, {3.568809946562046*^9, 3.5688099535664577`*^9}, {3.5688100003197403`*^9, 3.568810008556555*^9}, { 3.568810075371472*^9, 3.56881007981748*^9}, {3.568813431920168*^9, 3.568813533086345*^9}, {3.5688135960324554`*^9, 3.5688136045344706`*^9}, { 3.5688136407889347`*^9, 3.568813707494652*^9}, {3.568813812529636*^9, 3.5688138275524626`*^9}, {3.5688140611940727`*^9, 3.568814062738476*^9}, { 3.568814121113778*^9, 3.5688141965555105`*^9}, {3.568814326675339*^9, 3.568814331074547*^9}, 3.568814589473401*^9, 3.5688146762563534`*^9, 3.56881472549004*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{"\<\"C\"\>", ",", "\<\"Z\"\>", ",", "\<\"Z\"\>", ",", "\<\"Z\"\>", ",", "\<\"Z\"\>"}], "}"}]], "Output", CellChangeTimes->{ 3.568999729526589*^9, 3.569060092398421*^9, 3.5690652301408854`*^9, 3.5690652767069674`*^9, {3.569065376905943*^9, 3.569065406140394*^9}, 3.569065561828668*^9, 3.5690656436196117`*^9, 3.569065837403152*^9, 3.569066362000873*^9, 3.5690664044485483`*^9, {3.5690665341783752`*^9, 3.569066555456813*^9}, {3.5690666162189198`*^9, 3.569066634096551*^9}, { 3.569066675842224*^9, 3.569066705357476*^9}, {3.569066750503956*^9, 3.569066770565591*^9}, {3.569066820298478*^9, 3.5690668471305256`*^9}, 3.569066880530184*^9, 3.569069499512757*^9, 3.569069566124874*^9, 3.569069721547947*^9, 3.5690698096725016`*^9, 3.5690698550061812`*^9, 3.56906994545514*^9, 3.5690701756179447`*^9, 3.5690702503888755`*^9, 3.569070440163209*^9, 3.5690707652989798`*^9, 3.569071525013115*^9, 3.5690715824212155`*^9}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"num", "=", RowBox[{"Table", "[", RowBox[{"i", ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.568807971988578*^9, 3.5688080897687845`*^9}, { 3.5688081769573374`*^9, 3.5688081865669546`*^9}, {3.568808274317109*^9, 3.5688082853151283`*^9}, {3.5688085365847692`*^9, 3.568808578330443*^9}, { 3.568808822704872*^9, 3.5688088761037655`*^9}, {3.56880890727262*^9, 3.568808937365073*^9}, {3.568809097608555*^9, 3.5688091594002633`*^9}, 3.5688092418308077`*^9, {3.5688092753552666`*^9, 3.568809302577315*^9}, { 3.5688096257786827`*^9, 3.568809703201618*^9}, {3.5688097456336927`*^9, 3.568809788268568*^9}, 3.568809850575077*^9, {3.568809946562046*^9, 3.5688099535664577`*^9}, {3.5688100003197403`*^9, 3.568810008556555*^9}, { 3.568810075371472*^9, 3.56881007981748*^9}, {3.568813431920168*^9, 3.568813533086345*^9}, {3.5688135960324554`*^9, 3.5688136045344706`*^9}, { 3.5688136407889347`*^9, 3.568813707494652*^9}, {3.568813812529636*^9, 3.5688138275524626`*^9}, {3.5688140611940727`*^9, 3.568814062738476*^9}, { 3.568814121113778*^9, 3.5688141965555105`*^9}, {3.568814326675339*^9, 3.568814331074547*^9}, 3.568814589473401*^9, {3.5688146762563534`*^9, 3.5688146829331646`*^9}, 3.5688147231032352`*^9, 3.5688216336189733`*^9}], Cell[BoxData[ RowBox[{ RowBox[{"um", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"un", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"num", "[", RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.568807971988578*^9, 3.5688080897687845`*^9}, { 3.5688081769573374`*^9, 3.5688081865669546`*^9}, {3.568808274317109*^9, 3.5688082853151283`*^9}, {3.5688085365847692`*^9, 3.568808578330443*^9}, { 3.568808822704872*^9, 3.5688088761037655`*^9}, {3.56880890727262*^9, 3.568808937365073*^9}, {3.568809097608555*^9, 3.5688091594002633`*^9}, 3.5688092418308077`*^9, {3.5688092753552666`*^9, 3.568809302577315*^9}, { 3.5688096257786827`*^9, 3.568809703201618*^9}, {3.5688097456336927`*^9, 3.568809788268568*^9}, 3.568809850575077*^9, {3.568809946562046*^9, 3.5688099535664577`*^9}, {3.5688100003197403`*^9, 3.568810008556555*^9}, { 3.568810075371472*^9, 3.56881007981748*^9}, {3.568813431920168*^9, 3.568813533086345*^9}, {3.5688135960324554`*^9, 3.5688136045344706`*^9}, { 3.5688136407889347`*^9, 3.568813707494652*^9}, {3.568813812529636*^9, 3.5688138275524626`*^9}, {3.5688140611940727`*^9, 3.568814062738476*^9}, { 3.568814121113778*^9, 3.5688141965555105`*^9}, {3.568814326675339*^9, 3.568814331074547*^9}, 3.568814589473401*^9, {3.5688146762563534`*^9, 3.5688146885023746`*^9}, 3.5688147203576307`*^9, 3.56882164316619*^9}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"us", "=", RowBox[{"Sort", "[", "um", "]"}]}]], "Input", CellChangeTimes->{{3.568807971988578*^9, 3.5688080897687845`*^9}, { 3.5688081769573374`*^9, 3.5688081865669546`*^9}, {3.568808274317109*^9, 3.5688082853151283`*^9}, {3.5688085365847692`*^9, 3.568808578330443*^9}, { 3.568808822704872*^9, 3.5688088761037655`*^9}, {3.56880890727262*^9, 3.568808937365073*^9}, {3.568809097608555*^9, 3.5688091594002633`*^9}, 3.5688092418308077`*^9, {3.5688092753552666`*^9, 3.568809302577315*^9}, { 3.5688096257786827`*^9, 3.568809703201618*^9}, {3.5688097456336927`*^9, 3.568809788268568*^9}, 3.568809850575077*^9, {3.568809946562046*^9, 3.5688099535664577`*^9}, {3.5688100003197403`*^9, 3.568810008556555*^9}, { 3.568810075371472*^9, 3.56881007981748*^9}, {3.568813431920168*^9, 3.568813533086345*^9}, {3.5688135960324554`*^9, 3.5688136045344706`*^9}, { 3.5688136407889347`*^9, 3.568813707494652*^9}, {3.568813812529636*^9, 3.5688138275524626`*^9}, {3.5688140611940727`*^9, 3.568814062738476*^9}, { 3.568814121113778*^9, 3.5688141965555105`*^9}, {3.568814326675339*^9, 3.568814331074547*^9}, 3.568814589473401*^9, {3.5688146762563534`*^9, 3.5688146945707855`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"\<\"C\"\>", ",", "1"}], "}"}], ",", RowBox[{"{", RowBox[{"\<\"Z\"\>", ",", "2"}], "}"}], ",", RowBox[{"{", RowBox[{"\<\"Z\"\>", ",", "3"}], "}"}], ",", RowBox[{"{", RowBox[{"\<\"Z\"\>", ",", "4"}], "}"}], ",", RowBox[{"{", RowBox[{"\<\"Z\"\>", ",", "5"}], "}"}]}], "}"}]], "Output", CellChangeTimes->{ 3.5689997295421886`*^9, 3.5690600924296207`*^9, 3.5690652301720853`*^9, 3.5690652767381673`*^9, {3.5690653769371433`*^9, 3.5690654061559944`*^9}, 3.569065561859868*^9, 3.5690656436508117`*^9, 3.569065837434352*^9, 3.5690663620164733`*^9, 3.569066404479748*^9, {3.5690665341939754`*^9, 3.569066555472413*^9}, {3.56906661625012*^9, 3.5690666341277514`*^9}, { 3.5690666758734245`*^9, 3.5690667053886766`*^9}, {3.5690667505195556`*^9, 3.569066770581191*^9}, {3.5690668203140783`*^9, 3.5690668471461253`*^9}, 3.569066880545784*^9, 3.5690694995283566`*^9, 3.5690695661560736`*^9, 3.569069721579147*^9, 3.569069809688102*^9, 3.569069855021781*^9, 3.5690699454707403`*^9, 3.5690701756335444`*^9, 3.5690702504044757`*^9, 3.569070440178809*^9, 3.56907076533018*^9, 3.569071525028715*^9, 3.5690715824524155`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"num", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"us", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.56881016518083*^9, 3.5688102745058217`*^9}, { 3.5688103430679426`*^9, 3.5688103543155622`*^9}, {3.5688124958093233`*^9, 3.568812496386524*^9}, {3.5688126345871673`*^9, 3.568812635491969*^9}, { 3.568812673025635*^9, 3.5688126732128353`*^9}, {3.5688127492161684`*^9, 3.568812848541543*^9}, 3.5688129280080824`*^9, {3.568813062152718*^9, 3.5688131360032477`*^9}, 3.568813172678912*^9, 3.568821701869093*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{"1", ",", "2", ",", "3", ",", "4", ",", "5"}], "}"}]], "Output", CellChangeTimes->{ 3.568999729557789*^9, 3.569060092445221*^9, 3.5690652301876855`*^9, 3.5690652767381673`*^9, {3.5690653769527435`*^9, 3.5690654061715946`*^9}, 3.5690655618754683`*^9, 3.5690656436664114`*^9, 3.569065837449952*^9, 3.5690663620320735`*^9, 3.569066404479748*^9, {3.5690665342095757`*^9, 3.5690665554880133`*^9}, {3.56906661625012*^9, 3.5690666341433516`*^9}, { 3.5690666758734245`*^9, 3.5690667054042764`*^9}, {3.5690667505351553`*^9, 3.5690667705967913`*^9}, {3.5690668203296785`*^9, 3.5690668471617255`*^9}, 3.569066880561384*^9, 3.5690694995439568`*^9, 3.569069566171674*^9, 3.5690697215947466`*^9, 3.5690698097037015`*^9, 3.569069855037381*^9, 3.56906994548634*^9, 3.569070175649144*^9, 3.5690702504200754`*^9, 3.569070440194409*^9, 3.56907076534578*^9, 3.569071525044315*^9, 3.5690715824680157`*^9}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "<", RowBox[{"n", "+", "1"}]}], ",", RowBox[{"i", "++"}], ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"us", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "\[Equal]", "\"\\""}], ",", RowBox[{ RowBox[{"us", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "=", "\"\\""}]}], "]"}]}], "]"}], ";"}]], "Input", CellChangeTimes->{{3.568277663590458*^9, 3.568277712652544*^9}, 3.5682785838808746`*^9, {3.56880682872497*^9, 3.568806862015428*^9}, 3.5688069650848093`*^9, 3.568807273668951*^9, 3.5688073087846127`*^9, { 3.568807384553946*^9, 3.5688073893119545`*^9}, {3.56881489654434*^9, 3.5688149048591547`*^9}, {3.568821714349115*^9, 3.568821718951123*^9}}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{"u1", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"us", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "\[Equal]", "\"\\""}], ",", "0.255", ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"us", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "\[Equal]", "\"\\""}], ",", "0.1", ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"us", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "\[Equal]", "\"\\""}], ",", "0.24", ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"us", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "\[Equal]", "\"\\""}], ",", "0.178"}], "]"}]}], "]"}]}], "]"}]}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], "\[IndentingNewLine]", RowBox[{"u2", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"us", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "\[Equal]", "\"\\""}], ",", "0.25", ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"us", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "\[Equal]", "\"\\""}], ",", "0.53", ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"us", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "\[Equal]", "\"\\""}], ",", "0.22522", ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"us", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "\[Equal]", "\"\\""}], ",", "0.178"}], "]"}]}], "]"}]}], "]"}]}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}]}], "Input", CellChangeTimes->{{3.5679480369450207`*^9, 3.567948069112277*^9}, { 3.5679486805553513`*^9, 3.5679487948723516`*^9}, 3.567950206581231*^9, { 3.5681216077716007`*^9, 3.5681216839309344`*^9}, {3.568266380090639*^9, 3.5682664108850937`*^9}, {3.568266481350417*^9, 3.56826649971165*^9}, 3.5682795170433135`*^9, {3.5682795636873956`*^9, 3.5682797843809834`*^9}, { 3.568279837187076*^9, 3.5682798459854913`*^9}, {3.568279877201146*^9, 3.568279880711152*^9}, {3.5682799570264864`*^9, 3.5682799599124913`*^9}, { 3.5683770044498715`*^9, 3.5683771523693314`*^9}, 3.5688133425788107`*^9, { 3.568814931192001*^9, 3.5688150075073347`*^9}, {3.5688765350957685`*^9, 3.5688765981198792`*^9}, {3.568876634452343*^9, 3.5688766352167444`*^9}, { 3.568894358551657*^9, 3.5688944054453397`*^9}, {3.568894522929146*^9, 3.5688945259087515`*^9}, {3.5688946101644993`*^9, 3.5688946171221113`*^9}, {3.568894651067771*^9, 3.568894651956972*^9}, 3.5688946827982264`*^9, {3.568894734527918*^9, 3.568894735495119*^9}, { 3.5688947726699843`*^9, 3.568894773496786*^9}, {3.5688948078948464`*^9, 3.568895097914956*^9}, {3.5689558039020233`*^9, 3.568955804619624*^9}, { 3.5689559169086213`*^9, 3.5689559179070234`*^9}, {3.5689559496998787`*^9, 3.568955950589081*^9}, {3.568956090817727*^9, 3.568956118273775*^9}, 3.5689561556982408`*^9, {3.5690652215764704`*^9, 3.5690652227464724`*^9}, 3.569065370681532*^9, 3.569065829525138*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{"0.255`", ",", "0.1`", ",", "0.1`", ",", "0.1`", ",", "0.1`"}], "}"}]], "Output", CellChangeTimes->{ 3.568999729573389*^9, 3.569060092460821*^9, 3.569065230203285*^9, 3.569065276769367*^9, {3.569065376983943*^9, 3.5690654061871943`*^9}, 3.569065561891068*^9, 3.5690656436820116`*^9, 3.569065837481152*^9, 3.569066362047673*^9, 3.569066404510948*^9, {3.569066534225176*^9, 3.5690665555036135`*^9}, {3.56906661628132*^9, 3.5690666341589513`*^9}, { 3.5690666758890247`*^9, 3.569066705419876*^9}, {3.5690667505507555`*^9, 3.569066770612391*^9}, {3.5690668203452787`*^9, 3.5690668471773252`*^9}, 3.569066880592584*^9, 3.569069499559557*^9, 3.569069566187274*^9, 3.569069721610347*^9, 3.5690698097193017`*^9, 3.5690698550529814`*^9, 3.56906994550194*^9, 3.5690701756647444`*^9, 3.5690702504356756`*^9, 3.569070440210009*^9, 3.56907076536138*^9, 3.5690715250599146`*^9, 3.569071582483616*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{"0.25`", ",", "0.53`", ",", "0.53`", ",", "0.53`", ",", "0.53`"}], "}"}]], "Output", CellChangeTimes->{ 3.568999729573389*^9, 3.569060092460821*^9, 3.569065230203285*^9, 3.569065276769367*^9, {3.569065376983943*^9, 3.5690654061871943`*^9}, 3.569065561891068*^9, 3.5690656436820116`*^9, 3.569065837481152*^9, 3.569066362047673*^9, 3.569066404510948*^9, {3.569066534225176*^9, 3.5690665555036135`*^9}, {3.56906661628132*^9, 3.5690666341589513`*^9}, { 3.5690666758890247`*^9, 3.569066705419876*^9}, {3.5690667505507555`*^9, 3.569066770612391*^9}, {3.5690668203452787`*^9, 3.5690668471773252`*^9}, 3.569066880592584*^9, 3.569069499559557*^9, 3.569069566187274*^9, 3.569069721610347*^9, 3.5690698097193017`*^9, 3.5690698550529814`*^9, 3.56906994550194*^9, 3.5690701756647444`*^9, 3.5690702504356756`*^9, 3.569070440210009*^9, 3.56907076536138*^9, 3.5690715250599146`*^9, 3.569071582483616*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"u", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"us", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.5679472243083935`*^9, 3.567947394894693*^9}, { 3.567947453425996*^9, 3.567947471241227*^9}, {3.5679475415037503`*^9, 3.5679475448733563`*^9}, {3.5679488138107853`*^9, 3.567948820877598*^9}, { 3.568820296244224*^9, 3.5688203477867146`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{"\<\"C\"\>", ",", "\<\"H\"\>", ",", "\<\"H\"\>", ",", "\<\"H\"\>", ",", "\<\"H\"\>"}], "}"}]], "Output", CellChangeTimes->{ 3.5689997295889893`*^9, 3.5690600924920206`*^9, 3.5690652302188854`*^9, 3.5690652767849674`*^9, {3.569065376999543*^9, 3.569065406202794*^9}, 3.5690655619066677`*^9, 3.5690656437132115`*^9, 3.569065837496752*^9, 3.5690663620632734`*^9, 3.569066404526548*^9, {3.569066534240776*^9, 3.569066555519213*^9}, {3.56906661629692*^9, 3.5690666341745515`*^9}, { 3.569066675920225*^9, 3.5690667054354763`*^9}, {3.5690667505663557`*^9, 3.5690667706279907`*^9}, {3.5690668203608785`*^9, 3.5690668471929255`*^9}, 3.569066880608184*^9, 3.5690694995751567`*^9, 3.569069566202874*^9, 3.569069721625947*^9, 3.5690698097349014`*^9, 3.569069855068581*^9, 3.56906994551754*^9, 3.5690701756803446`*^9, 3.569070250451276*^9, 3.5690704402256093`*^9, 3.5690707653769803`*^9, 3.569071525075515*^9, 3.5690715824992156`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ 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[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.5688224161787477`*^9, 3.5688224173799496`*^9}, { 3.5688224905908785`*^9, 3.568822500434496*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{"6.`", ",", "1.`", ",", "1.`", ",", "1.`", ",", "1.`"}], "}"}]], "Output", CellChangeTimes->{ 3.568999729604589*^9, 3.569060092507621*^9, 3.5690652302344856`*^9, 3.569065276800567*^9, {3.5690653770151434`*^9, 3.5690654062183943`*^9}, 3.569065561922268*^9, 3.5690656437288113`*^9, 3.569065837527952*^9, 3.5690663620788736`*^9, 3.569066404542148*^9, {3.5690665342563763`*^9, 3.569066555534813*^9}, {3.56906661631252*^9, 3.569066634190151*^9}, { 3.569066675920225*^9, 3.5690667054510765`*^9}, {3.569066750581956*^9, 3.569066770643591*^9}, {3.569066820376478*^9, 3.5690668472085257`*^9}, 3.569066880623784*^9, 3.569069499590757*^9, 3.569069566218474*^9, 3.5690697216415467`*^9, 3.5690698097505016`*^9, 3.5690698550841813`*^9, 3.56906994553314*^9, 3.5690701756959443`*^9, 3.569070250451276*^9, 3.569070440241209*^9, 3.56907076539258*^9, 3.569071525091115*^9, 3.5690715825148153`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"bc", "=", "0"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", " ", RowBox[{"i", "<", RowBox[{"n", "+", "1"}]}], ",", RowBox[{"i", "++"}], ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"[", "i", "]"}], "]"}], "\[NotEqual]", "\"\\""}], ",", " ", RowBox[{"bc", "=", RowBox[{"bc", "+", "1"}]}]}], "]"}]}], "]"}], ";"}], "\[IndentingNewLine]", "bc", "\[IndentingNewLine]", RowBox[{"bh", "=", RowBox[{"n", "-", "bc"}]}]}], "Input", CellChangeTimes->{{3.5688281182607965`*^9, 3.5688281434392405`*^9}}], Cell[BoxData["1"], "Output", CellChangeTimes->{ 3.5689997296201887`*^9, 3.569060092523221*^9, 3.569065230250086*^9, 3.569065276800567*^9, {3.569065377030743*^9, 3.5690654062339945`*^9}, 3.569065561937868*^9, 3.5690656437444115`*^9, 3.5690658375435524`*^9, 3.569066362094474*^9, 3.569066404557748*^9, {3.569066534271976*^9, 3.569066555550413*^9}, {3.56906661631252*^9, 3.5690666342057514`*^9}, { 3.569066675935825*^9, 3.5690667054666767`*^9}, {3.569066750597556*^9, 3.569066770643591*^9}, {3.5690668203920784`*^9, 3.569066847224126*^9}, 3.5690668806393843`*^9, 3.569069499606357*^9, 3.569069566234074*^9, 3.569069721657147*^9, 3.569069809766102*^9, 3.569069855099781*^9, 3.5690699455487404`*^9, 3.5690701757115445`*^9, 3.569070250466876*^9, 3.569070440256809*^9, 3.56907076540818*^9, 3.5690715251067147`*^9, 3.5690715825304155`*^9}], Cell[BoxData["4"], "Output", CellChangeTimes->{ 3.5689997296201887`*^9, 3.569060092523221*^9, 3.569065230250086*^9, 3.569065276800567*^9, {3.569065377030743*^9, 3.5690654062339945`*^9}, 3.569065561937868*^9, 3.5690656437444115`*^9, 3.5690658375435524`*^9, 3.569066362094474*^9, 3.569066404557748*^9, {3.569066534271976*^9, 3.569066555550413*^9}, {3.56906661631252*^9, 3.5690666342057514`*^9}, { 3.569066675935825*^9, 3.5690667054666767`*^9}, {3.569066750597556*^9, 3.569066770643591*^9}, {3.5690668203920784`*^9, 3.569066847224126*^9}, 3.5690668806393843`*^9, 3.569069499606357*^9, 3.569069566234074*^9, 3.569069721657147*^9, 3.569069809766102*^9, 3.569069855099781*^9, 3.5690699455487404`*^9, 3.5690701757115445`*^9, 3.569070250466876*^9, 3.569070440256809*^9, 3.56907076540818*^9, 3.5690715251067147`*^9, 3.5690715825304155`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ RowBox[{ "DeleteFile", "[", " ", "\"\\"", "]"}], ";", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "<", RowBox[{"n", "+", "1"}]}], ",", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"u", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"Extract", "[", RowBox[{"t1", ",", RowBox[{"{", "i", "}"}]}], "]"}]}], "}"}], ">>>", " ", "\"\\""}]}], "]"}]}], "*)"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"ChemicalData", "\[NoBreak]", "[", RowBox[{"\"\\"", ",", "\"\\""}], "]"}]}]], "Input", CellChangeTimes->{{3.567869377256353*^9, 3.5678694101568108`*^9}, { 3.5678695244894114`*^9, 3.56786953511303*^9}, 3.56786988243204*^9, 3.5678699252853155`*^9, {3.567870069179968*^9, 3.5678700766835814`*^9}, { 3.5678703216664114`*^9, 3.567870325738019*^9}, 3.5679271251072817`*^9, { 3.567942840466694*^9, 3.567942854381918*^9}, {3.567942932272855*^9, 3.567942950509287*^9}, 3.5679430681802936`*^9, {3.568278272060727*^9, 3.56827827391713*^9}, 3.5688732159559383`*^9, 3.568910867631465*^9, 3.5689997098549547`*^9, {3.5690600808232*^9, 3.569060082227203*^9}, 3.5690651556507545`*^9}], Cell[BoxData["\<\"CH4\"\>"], "Output", CellChangeTimes->{ 3.5689997296201887`*^9, 3.569060092538821*^9, 3.5690652302812853`*^9, 3.569065276831767*^9, {3.569065377046343*^9, 3.569065406265195*^9}, 3.5690655619534683`*^9, 3.5690656437600117`*^9, 3.569065837559152*^9, 3.5690663621100736`*^9, 3.569066404573348*^9, {3.5690665342875757`*^9, 3.5690665555660133`*^9}, {3.56906661632812*^9, 3.5690666342213516`*^9}, { 3.569066675967025*^9, 3.569066705482277*^9}, {3.5690667506131563`*^9, 3.569066770659191*^9}, {3.5690668204076786`*^9, 3.5690668472553253`*^9}, 3.5690668806549845`*^9, 3.569069499621957*^9, 3.569069566249674*^9, 3.5690697216727467`*^9, 3.569069809781702*^9, 3.5690698551153812`*^9, 3.56906994556434*^9, 3.5690701757271442`*^9, 3.5690702504980755`*^9, 3.569070440288009*^9, 3.5690707654237804`*^9, 3.5690715251223145`*^9, 3.5690715825460157`*^9}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"t", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"t1", "[", RowBox[{"[", RowBox[{"num", "[", RowBox[{"[", "i", "]"}], "]"}], "]"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.567869377256353*^9, 3.567869394073182*^9}, { 3.5678694248832364`*^9, 3.5678694526200852`*^9}, {3.5678694891241493`*^9, 3.5678694955045605`*^9}, {3.567869596046737*^9, 3.567869602645549*^9}, 3.56786971175214*^9, {3.5678697908910794`*^9, 3.567869798800293*^9}, { 3.5679434548425727`*^9, 3.5679434621121855`*^9}, 3.568011875875622*^9, 3.56811435956207*^9, {3.5682782821071444`*^9, 3.5682782839479475`*^9}, { 3.5688205455482616`*^9, 3.5688206262784038`*^9}, 3.5688208026055136`*^9, 3.5688732159559383`*^9, 3.568910867631465*^9, {3.5689830694736805`*^9, 3.568983072687286*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"Plot1", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"Graphics3D", "[", RowBox[{"{", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\\""}], ",", RowBox[{"Lighter", "[", "Blue", "]"}], ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\\""}], ",", RowBox[{"Lighter", "[", "Red", "]"}], ",", RowBox[{"Lighter", "[", "Brown", "]"}]}], "]"}]}], "]"}], ",", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], ",", RowBox[{"t", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], ",", RowBox[{"t", "[", RowBox[{"[", RowBox[{"i", ",", "3"}], "]"}], "]"}]}], "}"}], ",", RowBox[{"u1", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}], "}"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.5679440895764875`*^9, 3.5679441000441055`*^9}, { 3.5679441319773617`*^9, 3.5679441523041973`*^9}, {3.567944240210352*^9, 3.5679442813476243`*^9}, {3.567944323966899*^9, 3.567944345931738*^9}, { 3.567944397567828*^9, 3.567944456941533*^9}, {3.567944491557994*^9, 3.5679445468444905`*^9}, {3.5679446505690727`*^9, 3.567944683875131*^9}, { 3.567944736244423*^9, 3.5679448177233667`*^9}, {3.567945549707852*^9, 3.5679455499730525`*^9}, {3.5679456039959474`*^9, 3.567945606195551*^9}, { 3.567945704460124*^9, 3.5679457164877453`*^9}, {3.5679457857362666`*^9, 3.5679458050023003`*^9}, {3.567946353670864*^9, 3.5679464127169676`*^9}, { 3.5679466013368993`*^9, 3.5679466311173515`*^9}, {3.5679467590999765`*^9, 3.567946915739851*^9}, {3.5679469466123056`*^9, 3.567946946815106*^9}, { 3.5679471385082426`*^9, 3.5679471497870626`*^9}, {3.5679483907224417`*^9, 3.5679484036080647`*^9}, {3.567948441937332*^9, 3.5679484421869326`*^9}, { 3.567948496740228*^9, 3.56794849772303*^9}, {3.567948847912445*^9, 3.5679488949933276`*^9}, {3.568011728877214*^9, 3.5680117308213253`*^9}, { 3.568011766325356*^9, 3.568011773436763*^9}, {3.5681173314452896`*^9, 3.5681173370612993`*^9}, {3.568117919800323*^9, 3.5681179678952074`*^9}, { 3.56811809459863*^9, 3.568118105409449*^9}, {3.5681181818963833`*^9, 3.5681182389144835`*^9}, {3.568118274888147*^9, 3.5681182883977704`*^9}, { 3.568118518513775*^9, 3.56811852739019*^9}, {3.5681188250543127`*^9, 3.568118862634779*^9}, {3.568118972396572*^9, 3.5681189862493963`*^9}, { 3.568119020694257*^9, 3.568119047136303*^9}, {3.5681190903015785`*^9, 3.568119136930061*^9}, {3.5681191671005135`*^9, 3.568119202808976*^9}, { 3.5681194788046613`*^9, 3.568119487603077*^9}, {3.5681195353235607`*^9, 3.5681195410331707`*^9}, 3.56811959201406*^9, {3.5681196426985493`*^9, 3.568119659702579*^9}, {3.568119726813897*^9, 3.56811974027672*^9}, { 3.5682659184078283`*^9, 3.5682659267694435`*^9}, {3.568265972477524*^9, 3.568266002725977*^9}, 3.568277802649102*^9, {3.568820898608082*^9, 3.568820922429324*^9}, 3.5688217528499823`*^9, 3.568821789307246*^9}], Cell[BoxData[ RowBox[{ RowBox[{"Plot2", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"Graphics3D", "[", RowBox[{"Text", "[", RowBox[{"i", ",", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "-", "0.3"}], ",", RowBox[{"t", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], ",", RowBox[{"t", "[", RowBox[{"[", RowBox[{"i", ",", "3"}], "]"}], "]"}]}], "}"}], " ", ",", RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}]}], "]"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.5679440895764875`*^9, 3.5679441000441055`*^9}, { 3.5679441319773617`*^9, 3.5679441523041973`*^9}, {3.567944240210352*^9, 3.5679442813476243`*^9}, {3.567944323966899*^9, 3.567944345931738*^9}, { 3.567944397567828*^9, 3.567944456941533*^9}, {3.567944491557994*^9, 3.5679445468444905`*^9}, {3.5679446505690727`*^9, 3.567944683875131*^9}, { 3.567944736244423*^9, 3.5679448177233667`*^9}, {3.567945549707852*^9, 3.5679455499730525`*^9}, {3.5679456039959474`*^9, 3.567945606195551*^9}, { 3.567945704460124*^9, 3.5679457164877453`*^9}, {3.5679457857362666`*^9, 3.5679458050023003`*^9}, {3.567946353670864*^9, 3.5679464127169676`*^9}, { 3.5679466013368993`*^9, 3.5679466311173515`*^9}, {3.5679467590999765`*^9, 3.567946915739851*^9}, {3.5679469466123056`*^9, 3.567946946815106*^9}, { 3.5679471385082426`*^9, 3.5679471497870626`*^9}, {3.5679483907224417`*^9, 3.5679484036080647`*^9}, {3.567948441937332*^9, 3.5679484421869326`*^9}, { 3.567948496740228*^9, 3.56794849772303*^9}, {3.567948847912445*^9, 3.5679488949933276`*^9}, {3.568011728877214*^9, 3.5680117308213253`*^9}, { 3.568011766325356*^9, 3.568011773436763*^9}, {3.5681173314452896`*^9, 3.5681173370612993`*^9}, {3.568117919800323*^9, 3.5681179678952074`*^9}, { 3.56811809459863*^9, 3.568118105409449*^9}, {3.5681181818963833`*^9, 3.5681182389144835`*^9}, {3.568118274888147*^9, 3.5681182883977704`*^9}, { 3.568118518513775*^9, 3.56811852739019*^9}, {3.5681188250543127`*^9, 3.568118862634779*^9}, {3.568118972396572*^9, 3.5681189862493963`*^9}, { 3.568119020694257*^9, 3.568119047136303*^9}, {3.5681190903015785`*^9, 3.568119136930061*^9}, {3.5681191671005135`*^9, 3.568119202808976*^9}, { 3.5681194788046613`*^9, 3.568119487603077*^9}, {3.5681195353235607`*^9, 3.5681195410331707`*^9}, 3.56811959201406*^9, {3.5681196426985493`*^9, 3.568119659702579*^9}, {3.568119726813897*^9, 3.56811974027672*^9}, { 3.5681199008322024`*^9, 3.568119920691037*^9}, {3.5681199536070952`*^9, 3.568119958521104*^9}, {3.568120000641178*^9, 3.568120067206495*^9}, { 3.5681211528280015`*^9, 3.5681211643408213`*^9}, 3.5681213201694956`*^9, { 3.568121503313817*^9, 3.568121508961027*^9}, {3.568266015439999*^9, 3.5682660461096525`*^9}, 3.568277806081108*^9, {3.568381589141924*^9, 3.568381589391525*^9}, 3.5688217602599955`*^9, {3.568821800305266*^9, 3.568821801397268*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"ChemicalData", "\[NoBreak]", "[", RowBox[{"\"\\"", ",", "\"\\""}], "]"}]], "Input", CellChangeTimes->{{3.5678630319943895`*^9, 3.5678630956269007`*^9}, { 3.5678632280087337`*^9, 3.567863228336334*^9}, {3.5678632621727934`*^9, 3.5678633509213495`*^9}, {3.5678644293824434`*^9, 3.567864451908883*^9}, { 3.567864923404111*^9, 3.56786494544695*^9}, {3.5678650143210707`*^9, 3.5678650845367947`*^9}, {3.567869231084096*^9, 3.5678692412397137`*^9}, 3.5678693168686466`*^9, {3.5678704450470285`*^9, 3.5678704621290584`*^9}, 3.5682783281584253`*^9, 3.5688732159715385`*^9, 3.568910867647065*^9, 3.5689997098549547`*^9}], Cell[BoxData[ GraphicsBox[{{ {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[ StyleBox["\<\"C\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.26354303745649343`], FontWeight->Plain], {253.69000000000003`, 0.}], LineBox[{{271.01102614831365`, 9.999102618207111}, {280.54, 15.5}}], {GrayLevel[0], InsetBox[ StyleBox["\<\"H\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.26354303745649343`], FontWeight->Plain], {307.39, 31.}], LineBox[{{280.54, 15.5}, {290.0689738516864, 21.00089738179288}}]}}, {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[ StyleBox["\<\"C\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.26354303745649343`], FontWeight->Plain], {253.69000000000003`, 0.}], LineBox[{{236.36978025515313`, -10.000499387041417`}, { 226.84500000000003`, -15.5}}], {GrayLevel[0], InsetBox[ StyleBox["\<\"H\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.26354303745649343`], FontWeight->Plain], {200., -31.}], LineBox[{{226.84500000000003`, -15.5}, { 217.3202197448469, -20.99950061295859}}]}}, {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[ StyleBox["\<\"C\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.26354303745649343`], FontWeight->Plain], {253.69000000000003`, 0.}], LineBox[{{243.6895006129586, 17.320219744846884`}, {238.19, 26.845000000000002`}}], {GrayLevel[0], InsetBox[ StyleBox["\<\"H\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.26354303745649343`], FontWeight->Plain], {222.69, 53.690000000000005`}], LineBox[{{238.19, 26.845000000000002`}, {232.69049938704143`, 36.36978025515312}}]}}, {GrayLevel[0], AbsoluteThickness[0.55], InsetBox[ StyleBox["\<\"C\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.26354303745649343`], FontWeight->Plain], {253.69000000000003`, 0.}], LineBox[{{263.6904993870414, -17.320219744846902`}, { 269.19, -26.845000000000002`}}], {GrayLevel[0], InsetBox[ StyleBox["\<\"H\"\>", StripOnInput->False, ScriptSizeMultipliers->0.8, ScriptMinSize->8, FontSize->Scaled[0.26354303745649343`], FontWeight->Plain], {284.69, -53.690000000000005`}], LineBox[{{269.19, -26.845000000000002`}, { 274.68950061295857`, -36.36978025515312}}]}}}, {}, {}}, ImageSize->{{0., 45.53336000000001}, {0., 49.224921203597276`}}, PlotRange->{{177.5, 325.39}, {-83.69, 76.19}}]], "Output", CellChangeTimes->{ 3.5689997296669893`*^9, 3.569060092585621*^9, 3.5690652303124857`*^9, 3.5690652768629675`*^9, {3.5690653770931435`*^9, 3.569065406296395*^9}, 3.569065562000268*^9, 3.569065643791212*^9, 3.569065837590352*^9, 3.5690663621412735`*^9, 3.569066404604548*^9, {3.5690665343343763`*^9, 3.5690665556128135`*^9}, {3.5690666163593197`*^9, 3.5690666342681513`*^9}, {3.5690666759982247`*^9, 3.569066705513477*^9}, { 3.569066750644356*^9, 3.5690667706903915`*^9}, {3.569066820438878*^9, 3.5690668472865257`*^9}, 3.5690668806861844`*^9, 3.569069499653157*^9, 3.5690695662808743`*^9, 3.569069721703947*^9, 3.5690698098129015`*^9, 3.569069855162182*^9, 3.5690699456111403`*^9, 3.5690701757583447`*^9, 3.569070250529276*^9, 3.569070440319209*^9, 3.5690707654549804`*^9, 3.569071525169115*^9, 3.569071582577216*^9}] }, Open ]], Cell[BoxData[{ RowBox[{ RowBox[{"x", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"y", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"z", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{"i", ",", "3"}], "]"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}]}], "Input", CellChangeTimes->{{3.5678630319943895`*^9, 3.5678630956269007`*^9}, { 3.5678632280087337`*^9, 3.567863228336334*^9}, {3.5678632621727934`*^9, 3.5678633509213495`*^9}, {3.5678644293824434`*^9, 3.567864451908883*^9}, { 3.567864923404111*^9, 3.56786494544695*^9}, {3.5678650143210707`*^9, 3.5678650845367947`*^9}, {3.567869231084096*^9, 3.567869278118179*^9}, { 3.5682778102931156`*^9, 3.5682778156907253`*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"(*", " ", RowBox[{"Distance", " ", "matrix", " ", "Bohr"}], " ", "*)"}], "Text"], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Dij", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{"i", "\[GreaterEqual]", "j"}], ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"sq", "=", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"x", "[", RowBox[{"[", "i", "]"}], "]"}], "-", RowBox[{"x", "[", RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"y", "[", RowBox[{"[", "i", "]"}], "]"}], "-", RowBox[{"y", "[", RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"z", "[", RowBox[{"[", "i", "]"}], "]"}], "-", RowBox[{"z", "[", RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}]}], "]"}]}], ";", " ", RowBox[{"sq", "<", "3.0"}]}], ",", "sq", ",", "0"}], "]"}], ",", "0"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"MatrixForm", "[", "Dij", "]"}]}]}]], "Input", CellChangeTimes->{{3.567863535157673*^9, 3.5678635412104836`*^9}, { 3.5678638182825704`*^9, 3.567863827283786*^9}, {3.5679503910983553`*^9, 3.5679504266352177`*^9}, {3.5679504653700857`*^9, 3.5679505048069553`*^9}, 3.567950534915008*^9, {3.567951107623214*^9, 3.567951136046464*^9}, 3.567951302545556*^9, {3.5679562986474605`*^9, 3.567956305355472*^9}, { 3.5679563590975666`*^9, 3.567956370688387*^9}, {3.567956414930065*^9, 3.567956430467692*^9}, {3.567956484740187*^9, 3.567956500667815*^9}, { 3.5679565456270943`*^9, 3.5679565868267665`*^9}, {3.5679569267045636`*^9, 3.5679569448941956`*^9}, {3.567956976562251*^9, 3.567957000711094*^9}, { 3.5679570348595533`*^9, 3.5679570632672033`*^9}, {3.5679573263928657`*^9, 3.5679573286548696`*^9}, {3.5679573611029263`*^9, 3.5679573637393312`*^9}, {3.5679592820903006`*^9, 3.567959317127962*^9}, { 3.567959397655304*^9, 3.5679594004945087`*^9}, {3.5679594704762316`*^9, 3.5679594724574347`*^9}, {3.567959655898157*^9, 3.567959673073787*^9}, 3.5681143121223865`*^9, 3.568121866388855*^9, 3.568220668221574*^9, { 3.568220718968463*^9, 3.568220723742071*^9}, {3.568220862051914*^9, 3.5682208641111183`*^9}, {3.568220906761593*^9, 3.5682209083527956`*^9}, { 3.5682209451844606`*^9, 3.5682209651680956`*^9}, {3.568221147782016*^9, 3.5682211695908546`*^9}, {3.568221237279373*^9, 3.568221253612602*^9}, { 3.568221305420293*^9, 3.568221306247094*^9}, {3.5682215161142635`*^9, 3.568221546206716*^9}, 3.56822158247678*^9, {3.568277819185131*^9, 3.5682778284359474`*^9}, 3.56882183684053*^9, {3.5688219650415554`*^9, 3.568821965915157*^9}, {3.5688769979329815`*^9, 3.5688771619984694`*^9}, { 3.568877204212144*^9, 3.5688772091729527`*^9}, {3.568877282493081*^9, 3.5688772843338842`*^9}, {3.568877374112042*^9, 3.568877378074449*^9}, { 3.56887744125456*^9, 3.5688774664954042`*^9}, {3.56887758366721*^9, 3.568877604867647*^9}, {3.568877667220957*^9, 3.5688776710429635`*^9}, { 3.5688778111156096`*^9, 3.56887782262843*^9}, 3.569065552687052*^9}], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"0.`", "0", "0", "0", "0"}, {"2.060314249067387`", "0.`", "0", "0", "0"}, {"2.0603203349966233`", "0", "0.`", "0", "0"}, {"2.060354568301202`", "0", "0", "0.`", "0"}, {"2.0602326842173277`", "0", "0", "0", "0.`"} }, GridBoxAlignment->{ "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, "RowsIndexed" -> {}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.7]}, Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Output", CellChangeTimes->{ 3.569065562015868*^9, 3.5690656438068123`*^9, 3.5690658376059523`*^9, 3.5690663621568737`*^9, 3.569066404620148*^9, {3.569066534349976*^9, 3.569066555628413*^9}, {3.56906661637492*^9, 3.5690666342837515`*^9}, { 3.5690666760138245`*^9, 3.5690667055290766`*^9}, {3.569066750659956*^9, 3.569066770721591*^9}, {3.5690668204700785`*^9, 3.569066847302126*^9}, 3.569066880701784*^9, 3.5690694996687574`*^9, 3.569069566296474*^9, 3.5690697217195473`*^9, 3.5690698098285017`*^9, 3.5690698551777816`*^9, 3.5690699456267405`*^9, 3.569070175773945*^9, 3.569070250544876*^9, 3.5690704403348093`*^9, 3.56907076547058*^9, 3.5690715251847153`*^9, 3.569071582608416*^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", "<", "i"}], ",", 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.6"}], ",", 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", "]"}], "]"}], "=", "0.5241"}]}], "]"}]}], "]"}]}], "]"}]}], "]"}]}], "]"}]}], "]"}]}], "]"}], "\[IndentingNewLine]", "u2"}], "Input", CellChangeTimes->{{3.5688800293771057`*^9, 3.5688801774057655`*^9}, { 3.5688802290574565`*^9, 3.568880330208034*^9}, {3.5688803807053227`*^9, 3.5688804136525807`*^9}, 3.568880491434317*^9, 3.5688805648948464`*^9, { 3.5688806280437574`*^9, 3.568880655234605*^9}, {3.568880712424306*^9, 3.568880716324312*^9}, {3.5688808087544746`*^9, 3.568880827786508*^9}, 3.5688808733385878`*^9, {3.568881358936241*^9, 3.5688813959239063`*^9}, { 3.568881446561595*^9, 3.5688815849182377`*^9}, {3.5688816223739033`*^9, 3.568881642607139*^9}, 3.568881728953291*^9, {3.568881797624612*^9, 3.568881809870633*^9}, 3.568881863612727*^9, {3.568883098916497*^9, 3.5688831014749017`*^9}, {3.56888372124679*^9, 3.568883843239004*^9}, 3.568888211750622*^9, {3.5688882521546926`*^9, 3.5688885232519693`*^9}, { 3.5688885890840845`*^9, 3.5688885895052853`*^9}, {3.568888675773437*^9, 3.5688886914514647`*^9}, {3.568888736285943*^9, 3.5688887717760057`*^9}, { 3.5688888256429005`*^9, 3.5688888541441507`*^9}, {3.568892442649892*^9, 3.5688924499507046`*^9}, {3.568892758441247*^9, 3.5688927597360487`*^9}, { 3.5688928234309607`*^9, 3.5688928379389863`*^9}, 3.5690652691877537`*^9, 3.5690715187575035`*^9, {3.5690715736228*^9, 3.5690715745432014`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{"0.25`", ",", "0.6`", ",", "0.6`", ",", "0.6`", ",", "0.6`"}], "}"}]], "Output", CellChangeTimes->{ 3.5689997296981893`*^9, 3.5690600926168213`*^9, 3.569065230328086*^9, 3.5690652768941674`*^9, {3.5690653771087437`*^9, 3.5690654063119946`*^9}, 3.5690655620314684`*^9, 3.569065643822412*^9, 3.5690658376215525`*^9, 3.569066362172474*^9, 3.569066404635748*^9, {3.569066534349976*^9, 3.5690665556440134`*^9}, {3.56906661639052*^9, 3.5690666342993517`*^9}, { 3.5690666760294247`*^9, 3.569066705544677*^9}, {3.569066750675556*^9, 3.569066770737191*^9}, {3.5690668204856787`*^9, 3.5690668473177257`*^9}, 3.5690668807173843`*^9, 3.569069499684357*^9, 3.5690695663120737`*^9, 3.5690697217351475`*^9, 3.569069809844102*^9, 3.5690698551933813`*^9, 3.5690699456423407`*^9, 3.5690701757895446`*^9, 3.5690702505604763`*^9, 3.5690704403504095`*^9, 3.5690707654861803`*^9, 3.5690715252003155`*^9, 3.569071582624016*^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.5683793122399254`*^9, 3.568379313019926*^9}, { 3.568379456758579*^9, 3.56837947986222*^9}, 3.5683796000760307`*^9, { 3.5683796473285136`*^9, 3.5683796831617765`*^9}, {3.568381744939398*^9, 3.5683817743298492`*^9}, {3.56846523720809*^9, 3.568465240218895*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"Nij", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{"i", "\[LessEqual]", "j"}], ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"Lij", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], ">", "0"}], ",", RowBox[{"2.25", "/", RowBox[{ RowBox[{"Lij", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], "^", "2"}]}], ",", "0"}], "]"}], ",", "\"\< \>\""}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"(*", RowBox[{"MatrixForm", "[", "Nij", "]"}], "*)"}]}]}], "Input", CellChangeTimes->{{3.5679628944794455`*^9, 3.5679629710755796`*^9}, { 3.567963889963994*^9, 3.567963933800071*^9}, {3.5680478391738477`*^9, 3.5680478414358516`*^9}, 3.568048231171336*^9, {3.5680482823550262`*^9, 3.56804828449223*^9}, 3.5680483349271183`*^9, 3.568114273059918*^9, 3.5681219175101447`*^9, {3.568277851305588*^9, 3.568277857295998*^9}, { 3.5688722635742655`*^9, 3.568872266616271*^9}, 3.5688723032763357`*^9, 3.5688725225503206`*^9}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Tval", "=", RowBox[{"Sum", "[", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{"i", "\[LessEqual]", "j"}], ",", RowBox[{"Nij", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], ",", "0"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", "n"}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.567963992456174*^9, 3.567964003890994*^9}, { 3.5679640974287586`*^9, 3.5679642038209453`*^9}, {3.56796433180357*^9, 3.567964353487608*^9}, {3.567964403641696*^9, 3.567964453046983*^9}, { 3.5679646160984693`*^9, 3.5679646168316708`*^9}, {3.567964688747797*^9, 3.567964762255126*^9}, {3.568277859480002*^9, 3.568277863598409*^9}, { 3.568379511358675*^9, 3.5683795129030776`*^9}, 3.5688214977895346`*^9, { 3.56882774467134*^9, 3.56882781892747*^9}}], Cell[BoxData["6.144018671283403`"], "Output", CellChangeTimes->{ 3.568999729729389*^9, 3.5690600926324215`*^9, 3.5690652303592854`*^9, 3.5690652769253674`*^9, {3.5690653771399436`*^9, 3.569065406343195*^9}, 3.5690655620626683`*^9, 3.569065643853612*^9, 3.569065837668352*^9, 3.5690663622036734`*^9, 3.5690664046669483`*^9, {3.569066534381176*^9, 3.569066555675214*^9}, {3.5690666164217205`*^9, 3.5690666343305516`*^9}, { 3.569066676060625*^9, 3.5690667055758767`*^9}, {3.569066750706756*^9, 3.5690667707527914`*^9}, {3.569066820516879*^9, 3.5690668473489256`*^9}, 3.5690668807485847`*^9, 3.569069499715557*^9, 3.569069566343274*^9, 3.569069721766347*^9, 3.5690698098753023`*^9, 3.5690698552245817`*^9, 3.56906994567354*^9, 3.5690701758207445`*^9, 3.5690702505916758`*^9, 3.569070440381609*^9, 3.5690707655173807`*^9, 3.569071525247115*^9, 3.569071582655216*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ 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"}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{ 3.5688724058309155`*^9, {3.568872452646598*^9, 3.568872478480243*^9}, { 3.5688764461756124`*^9, 3.5688764647552447`*^9}}], Cell[BoxData["36.`"], "Output", CellChangeTimes->{ 3.5689997297449894`*^9, 3.569060092648021*^9, 3.5690652303748856`*^9, 3.5690652769409676`*^9, {3.5690653771399436`*^9, 3.5690654063587947`*^9}, 3.569065562078268*^9, 3.569065643869212*^9, 3.5690658376839523`*^9, 3.5690663622192736`*^9, 3.5690664046669483`*^9, {3.569066534396776*^9, 3.569066555675214*^9}, {3.56906661643732*^9, 3.5690666343461514`*^9}, { 3.5690666760762253`*^9, 3.569066705591477*^9}, {3.569066750722356*^9, 3.569066770783991*^9}, {3.569066820516879*^9, 3.569066847364526*^9}, 3.5690668807641845`*^9, 3.5690694997311573`*^9, 3.5690695663588743`*^9, 3.569069721781947*^9, 3.569069809890902*^9, 3.569069855240182*^9, 3.56906994567354*^9, 3.5690701758363447`*^9, 3.569070250607276*^9, 3.569070440397209*^9, 3.5690707655329804`*^9, 3.5690715252627153`*^9, 3.569071582670816*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Ttot", "=", RowBox[{"Tval", "+", "Tcor", " ", RowBox[{"(*", " ", RowBox[{"plus", " ", "\[Pi]", " ", "and", " ", "lone", " ", RowBox[{"pairs", ":", " ", RowBox[{"+", " ", "8.7616"}]}]}], " ", "*)"}]}]}]], "Input", CellChangeTimes->{ 3.5688724058309155`*^9, {3.568872452646598*^9, 3.568872478480243*^9}, { 3.5688764461756124`*^9, 3.568876459248435*^9}, {3.568876733122516*^9, 3.5688767685033784`*^9}, {3.568895188410715*^9, 3.5688952349923964`*^9}}], Cell[BoxData["42.144018671283405`"], "Output", CellChangeTimes->{ 3.5689997297449894`*^9, 3.569060092679221*^9, 3.569065230390486*^9, 3.5690652769409676`*^9, {3.5690653771555433`*^9, 3.569065406374395*^9}, 3.5690655620938683`*^9, 3.5690656438848124`*^9, 3.5690658376995525`*^9, 3.569066362234874*^9, 3.5690664046825485`*^9, {3.569066534412376*^9, 3.5690665556908135`*^9}, {3.56906661645292*^9, 3.5690666343617516`*^9}, { 3.569066676091825*^9, 3.5690667056070766`*^9}, {3.569066750737956*^9, 3.569066770799591*^9}, {3.569066820532479*^9, 3.569066847380126*^9}, 3.569066880779784*^9, 3.5690694997467575`*^9, 3.5690695663744745`*^9, 3.5690697217975473`*^9, 3.569069809906502*^9, 3.5690698552557817`*^9, 3.5690699456891403`*^9, 3.569070175851945*^9, 3.569070250622876*^9, 3.5690704404128094`*^9, 3.56907076554858*^9, 3.5690715252783155`*^9, 3.5690715826864157`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"TvalH", "=", RowBox[{"Sum", "[", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"i", "\[LessEqual]", "j"}], " ", "&&", " ", RowBox[{ RowBox[{"u", "[", RowBox[{"[", "j", "]"}], "]"}], "\[Equal]", " ", "\"\\""}]}], ")"}], ",", RowBox[{"Nij", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], ",", "0"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "n"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", "n"}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{ 3.5688724058309155`*^9, {3.568872452646598*^9, 3.568872478480243*^9}, { 3.5688764461756124`*^9, 3.568876459248435*^9}, {3.568876733122516*^9, 3.5688767642445707`*^9}}], Cell[BoxData["6.144018671283403`"], "Output", CellChangeTimes->{ 3.568999729760589*^9, 3.5690600926948214`*^9, 3.569065230406086*^9, 3.569065276956568*^9, {3.5690653771711435`*^9, 3.5690654063899946`*^9}, 3.5690655621094685`*^9, 3.569065643900412*^9, 3.5690658377151527`*^9, 3.569066362234874*^9, 3.5690664046981487`*^9, {3.569066534427976*^9, 3.5690665557064133`*^9}, {3.56906661646852*^9, 3.5690666343773518`*^9}, { 3.569066676091825*^9, 3.569066705622677*^9}, {3.5690667507535563`*^9, 3.569066770799591*^9}, {3.5690668205480785`*^9, 3.569066847395726*^9}, 3.5690668807953844`*^9, 3.569069499762357*^9, 3.5690695663900747`*^9, 3.5690697218131475`*^9, 3.569069809922102*^9, 3.5690698552713814`*^9, 3.5690699457047405`*^9, 3.569070175851945*^9, 3.569070250638476*^9, 3.5690704404128094`*^9, 3.5690707655641804`*^9, 3.5690715252939153`*^9, 3.569071582702016*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"TvalC", "=", RowBox[{"Tval", "-", "TvalH"}]}]], "Input", CellChangeTimes->{{3.5688726088808727`*^9, 3.568872610581275*^9}, { 3.5688726952426243`*^9, 3.5688727146490583`*^9}}], Cell[BoxData["0.`"], "Output", CellChangeTimes->{ 3.568999729776189*^9, 3.569060092710421*^9, 3.5690652304216857`*^9, 3.5690652769721675`*^9, {3.5690653771867437`*^9, 3.569065406405595*^9}, 3.5690655621250687`*^9, 3.569065643916012*^9, 3.5690658377307525`*^9, 3.569066362250474*^9, 3.5690664047137485`*^9, {3.569066534443576*^9, 3.5690665557220135`*^9}, {3.56906661646852*^9, 3.5690666343929515`*^9}, { 3.5690666761074247`*^9, 3.569066705638277*^9}, {3.5690667507535563`*^9, 3.5690667708151913`*^9}, {3.5690668205636787`*^9, 3.569066847411326*^9}, 3.5690668808109846`*^9, 3.569069499762357*^9, 3.5690695664056745`*^9, 3.5690697218287473`*^9, 3.569069809937702*^9, 3.5690698552869816`*^9, 3.5690699457203407`*^9, 3.569070175867545*^9, 3.569070250654076*^9, 3.5690704404284096`*^9, 3.5690707655797806`*^9, 3.569071525309515*^9, 3.569071582717616*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"Vnn", "=", "0.0"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"For", "[", " ", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "<", "n"}], ",", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"j", "=", RowBox[{"i", "+", "1"}]}], ",", RowBox[{"j", "<", RowBox[{"n", "+", "1"}]}], ",", RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", RowBox[{"Vnn", "=", RowBox[{"Vnn", "+", RowBox[{ RowBox[{"cn", "[", RowBox[{"[", "i", "]"}], "]"}], "*", RowBox[{ RowBox[{"cn", "[", RowBox[{"[", "j", "]"}], "]"}], "/", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"x", "[", RowBox[{"[", "i", "]"}], "]"}], "-", RowBox[{"x", "[", RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"y", "[", RowBox[{"[", "i", "]"}], "]"}], "-", RowBox[{"y", "[", RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"z", "[", RowBox[{"[", "i", "]"}], "]"}], "-", RowBox[{"z", "[", RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}]}], "]"}]}]}]}]}]}], "]"}]}], "]"}], ";"}], "\[IndentingNewLine]", "Vnn"}], "Input"], Cell[BoxData["13.432102855050447`"], "Output", CellChangeTimes->{ 3.568999729791789*^9, 3.5690600927260213`*^9, 3.569065230437286*^9, 3.5690652769877677`*^9, {3.569065377202344*^9, 3.569065406405595*^9}, 3.5690655621406684`*^9, 3.569065643947212*^9, 3.569065837746352*^9, 3.5690663622660737`*^9, 3.5690664047449484`*^9, {3.569066534459176*^9, 3.5690665557376137`*^9}, {3.5690666164841204`*^9, 3.5690666344085517`*^9}, {3.569066676123025*^9, 3.569066705653877*^9}, { 3.569066750769156*^9, 3.5690667708307915`*^9}, {3.569066820579279*^9, 3.5690668474269257`*^9}, 3.5690668808109846`*^9, 3.569069499777957*^9, 3.569069566421274*^9, 3.569069721844347*^9, 3.569069809953302*^9, 3.5690698552869816`*^9, 3.5690699457359405`*^9, 3.569070175883145*^9, 3.569070250669676*^9, 3.56907044044401*^9, 3.5690707655953803`*^9, 3.569071525325115*^9, 3.569071582717616*^9}] }, Open ]], Cell[BoxData[ RowBox[{"\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Mij", "=", RowBox[{"Prepend", "[", RowBox[{"Qij", ",", "u"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Oij", "=", RowBox[{"Transpose", "[", "Mij", "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Pij", "=", RowBox[{"Prepend", "[", RowBox[{"Oij", ",", RowBox[{"Prepend", "[", RowBox[{"u", ",", "\"\< \>\""}], "]"}]}], "]"}]}], ";"}]}]}]], "Input", CellChangeTimes->{{3.568045828002715*^9, 3.56804588761042*^9}, 3.5680459828797874`*^9, {3.568046033938677*^9, 3.5680460518787084`*^9}, { 3.568046989190755*^9, 3.5680470284404235`*^9}, {3.568047061028881*^9, 3.568047085567724*^9}, {3.5680471562982483`*^9, 3.568047169199471*^9}, { 3.568047223393966*^9, 3.568047364059413*^9}, {3.568047417286707*^9, 3.5680474225127163`*^9}, {3.5680474583459787`*^9, 3.568047487861231*^9}, { 3.5680475345521126`*^9, 3.5680475468605347`*^9}, {3.5680479168775845`*^9, 3.5680479311516094`*^9}, {3.568048066918648*^9, 3.568048067667449*^9}, { 3.5680483816336*^9, 3.5680483930216208`*^9}, {3.568048424159275*^9, 3.568048450211321*^9}, {3.5680486326224413`*^9, 3.5680486370996494`*^9}, { 3.5680487358166227`*^9, 3.5680487967191296`*^9}, 3.568048834377596*^9, 3.568048984418659*^9, {3.5680490255559316`*^9, 3.5680490790796256`*^9}, { 3.5680491274553103`*^9, 3.5680491280949116`*^9}, {3.5680491782333994`*^9, 3.568049231944294*^9}, {3.568049361533722*^9, 3.5680493691933355`*^9}, { 3.568049422171028*^9, 3.568049427490638*^9}, {3.568049476427923*^9, 3.5680495489368505`*^9}, {3.568049620634577*^9, 3.5680496274361887`*^9}, { 3.5680496655314555`*^9, 3.568049693876705*^9}, 3.568049746573598*^9, { 3.568049846460573*^9, 3.568049932759925*^9}, {3.568050009059659*^9, 3.5680500167348723`*^9}, {3.5680501572755194`*^9, 3.568050179911159*^9}, 3.5681087021497335`*^9}], Cell[BoxData[ RowBox[{"(*", RowBox[{ RowBox[{"MatrixForm", "[", "Pij", "]"}], ";"}], "*)"}]], "Input", CellChangeTimes->{{3.568045828002715*^9, 3.56804588761042*^9}, 3.5680459828797874`*^9, {3.568046033938677*^9, 3.5680460518787084`*^9}, { 3.568046989190755*^9, 3.5680470284404235`*^9}, {3.568047061028881*^9, 3.568047085567724*^9}, {3.5680471562982483`*^9, 3.568047169199471*^9}, { 3.568047223393966*^9, 3.568047364059413*^9}, {3.568047417286707*^9, 3.5680474225127163`*^9}, {3.5680474583459787`*^9, 3.568047487861231*^9}, { 3.5680475345521126`*^9, 3.5680475468605347`*^9}, {3.5680479168775845`*^9, 3.5680479311516094`*^9}, {3.568048066918648*^9, 3.568048067667449*^9}, { 3.5680483816336*^9, 3.5680483930216208`*^9}, {3.568048424159275*^9, 3.568048450211321*^9}, {3.5680486326224413`*^9, 3.5680486370996494`*^9}, { 3.5680487358166227`*^9, 3.5680487967191296`*^9}, 3.568048834377596*^9, 3.568048984418659*^9, {3.5680490255559316`*^9, 3.5680490790796256`*^9}, { 3.5680491274553103`*^9, 3.5680491280949116`*^9}, {3.5680491782333994`*^9, 3.568049231944294*^9}, {3.568049361533722*^9, 3.5680493691933355`*^9}, { 3.568049422171028*^9, 3.568049427490638*^9}, {3.568049476427923*^9, 3.5680495489368505`*^9}, {3.568049620634577*^9, 3.5680496274361887`*^9}, { 3.5680496655314555`*^9, 3.568049693876705*^9}, 3.568049746573598*^9, { 3.568049846460573*^9, 3.568049932759925*^9}, {3.568050009059659*^9, 3.5680500167348723`*^9}, {3.5680501572755194`*^9, 3.568050179911159*^9}, { 3.5681087021497335`*^9, 3.5681087223049684`*^9}, 3.568114284401138*^9, 3.5681220092695055`*^9}], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{"sigma", " ", "bond", " ", "cloud", " ", "radii"}], " ", "*)"}]], "Input", CellChangeTimes->{{3.568045828002715*^9, 3.56804588761042*^9}, 3.5680459828797874`*^9, {3.568046033938677*^9, 3.5680460518787084`*^9}, { 3.568046989190755*^9, 3.5680470284404235`*^9}, {3.568047061028881*^9, 3.568047085567724*^9}, {3.5680471562982483`*^9, 3.568047169199471*^9}, { 3.568047223393966*^9, 3.568047364059413*^9}, {3.568047417286707*^9, 3.5680474225127163`*^9}, {3.5680474583459787`*^9, 3.568047487861231*^9}, { 3.5680475345521126`*^9, 3.5680475468605347`*^9}, {3.5680479168775845`*^9, 3.5680479311516094`*^9}, {3.568048066918648*^9, 3.568048067667449*^9}, { 3.5680483816336*^9, 3.5680483930216208`*^9}, {3.568048424159275*^9, 3.568048450211321*^9}, {3.5680486326224413`*^9, 3.5680486370996494`*^9}, { 3.5680487358166227`*^9, 3.5680487967191296`*^9}, 3.568048834377596*^9, 3.568048984418659*^9, {3.5680490255559316`*^9, 3.5680490790796256`*^9}, { 3.5680491274553103`*^9, 3.5680491280949116`*^9}, {3.5680491782333994`*^9, 3.568049231944294*^9}, {3.568049361533722*^9, 3.5680493691933355`*^9}, { 3.568049422171028*^9, 3.568049427490638*^9}, {3.568049476427923*^9, 3.5680495489368505`*^9}, {3.568049620634577*^9, 3.5680496274361887`*^9}, { 3.5680496655314555`*^9, 3.568049693876705*^9}, 3.568049746573598*^9, { 3.568049846460573*^9, 3.568049932759925*^9}, {3.568050009059659*^9, 3.5680500167348723`*^9}, {3.5680501572755194`*^9, 3.568050179911159*^9}, { 3.5681087021497335`*^9, 3.5681087223049684`*^9}, {3.5681133837179556`*^9, 3.5681133888191648`*^9}}], Cell[BoxData[ RowBox[{"(*", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "2"}], ",", RowBox[{"i", "<", RowBox[{"n", "+", "2"}]}], ",", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"j", "=", "2"}], ",", RowBox[{"j", "\[LessEqual]", "i"}], ",", RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"Pij", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], ">", "0"}], ",", RowBox[{"Print", "[", RowBox[{ RowBox[{"u", "[", RowBox[{"[", RowBox[{"i", "-", "1"}], "]"}], "]"}], ",", RowBox[{"i", "-", "1"}], ",", "\"\< \>\"", ",", RowBox[{"u", "[", RowBox[{"[", RowBox[{"j", "-", "1"}], "]"}], "]"}], ",", RowBox[{"j", "-", "1"}], ",", "\"\< \>\"", ",", RowBox[{"Extract", "[", RowBox[{"Pij", ",", RowBox[{"{", RowBox[{"i", ",", "j"}], "}"}]}], "]"}]}], "]"}]}], "]"}]}], "]"}]}], "]"}], "*)"}]], "Input", CellChangeTimes->{{3.568108824079547*^9, 3.568109012824279*^9}, { 3.568109558076436*^9, 3.568109588168889*^9}, {3.5681097523435774`*^9, 3.5681097637783976`*^9}, {3.5681100381828794`*^9, 3.5681101616726966`*^9}, {3.5681102015307665`*^9, 3.5681102288464146`*^9}, {3.568110263743676*^9, 3.568110267331682*^9}, { 3.5681103163469687`*^9, 3.5681103268925867`*^9}, {3.5681106443375444`*^9, 3.5681106916368275`*^9}, 3.568110873002746*^9, {3.5681109388660617`*^9, 3.5681110268970165`*^9}, {3.568111092526332*^9, 3.568111123461186*^9}, { 3.5681111884196997`*^9, 3.568111192756508*^9}, {3.5681112268113675`*^9, 3.568111366213212*^9}, {3.568111434946933*^9, 3.568111463354583*^9}, { 3.568111544396725*^9, 3.568111547625931*^9}, {3.5681116856081734`*^9, 3.5681117131422215`*^9}, {3.568111864696488*^9, 3.5681119835998964`*^9}, 3.5681120218355637`*^9, {3.56811206543764*^9, 3.568112066389242*^9}, { 3.568112332369709*^9, 3.5681123404505234`*^9}, {3.568112505935614*^9, 3.56811251481203*^9}, {3.5681125623609133`*^9, 3.568112569006525*^9}, { 3.568112648067464*^9, 3.568112734210815*^9}, {3.568112764443668*^9, 3.5681127697320776`*^9}, {3.5681128018681335`*^9, 3.568112843426607*^9}, { 3.568112884641879*^9, 3.5681128887602863`*^9}, {3.5681129663080225`*^9, 3.5681129699584293`*^9}, {3.568113133789917*^9, 3.568113227936082*^9}, { 3.5681132946885996`*^9, 3.5681133044698167`*^9}, {3.568113502200164*^9, 3.568113507348173*^9}, 3.56811354556824*^9, {3.5681222886191964`*^9, 3.5681222896019983`*^9}, {3.568207424138049*^9, 3.5682074415164795`*^9}, 3.568277923134914*^9, 3.568883508322016*^9, 3.568883588350157*^9, 3.5688836441982546`*^9}], Cell[BoxData[{ RowBox[{ RowBox[{"DeleteFile", "[", " ", "\"\\"", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "2"}], ",", RowBox[{"i", "<", RowBox[{"n", "+", "2"}]}], ",", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"j", "=", "2"}], ",", RowBox[{"j", "\[LessEqual]", "i"}], ",", RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"Pij", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], ">", "0"}], ",", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"i", "-", "1"}], ",", RowBox[{"j", "-", "1"}], ",", RowBox[{"Extract", "[", RowBox[{"Pij", ",", RowBox[{"{", RowBox[{"i", ",", "j"}], "}"}]}], "]"}]}], "}"}], ">>>", " ", "\"\\""}]}], "]"}]}], "]"}]}], "]"}]}], "Input", CellChangeTimes->{{3.5681154444347754`*^9, 3.568115497147268*^9}, { 3.568127214924899*^9, 3.568127265796589*^9}, {3.568128215276656*^9, 3.568128238536297*^9}, {3.568130718566253*^9, 3.5681307485339055`*^9}, 3.5681308175952272`*^9, {3.5681308916497574`*^9, 3.5681309458598523`*^9}, { 3.5681309983383446`*^9, 3.5681310118947687`*^9}, {3.5681310723448744`*^9, 3.568131074357278*^9}, {3.5681315215476637`*^9, 3.568131584509374*^9}, { 3.56813162789305*^9, 3.5681316427754765`*^9}, {3.5681317113063965`*^9, 3.568131715268804*^9}, {3.5681323203938665`*^9, 3.5681323260098763`*^9}, { 3.568132387567585*^9, 3.568132440014877*^9}, {3.5681325121026034`*^9, 3.568132590945142*^9}, {3.56813274298981*^9, 3.5681327548302307`*^9}, { 3.568132877867647*^9, 3.568132888335265*^9}, {3.568133052345154*^9, 3.568133072344389*^9}, {3.5681331740533676`*^9, 3.5681331951446047`*^9}, { 3.568133267793932*^9, 3.5681332684803333`*^9}, {3.5681335079111547`*^9, 3.568133536833606*^9}, {3.5681336272201643`*^9, 3.5681336314945717`*^9}, { 3.568133737994359*^9, 3.568133744000369*^9}, {3.568133978515581*^9, 3.5681339917912045`*^9}, {3.5681340324916763`*^9, 3.5681340418204927`*^9}, {3.568134244262048*^9, 3.5681342454944506`*^9}, { 3.5681343354754086`*^9, 3.568134336099409*^9}, {3.5681345828762426`*^9, 3.568134583952645*^9}, {3.5681348406666956`*^9, 3.5681348474683075`*^9}, 3.568134889151581*^9, {3.5681352914606876`*^9, 3.568135296265496*^9}, { 3.568135518877887*^9, 3.568135519595488*^9}, 3.5681357593679094`*^9, { 3.5681787035360527`*^9, 3.568178728449296*^9}, 3.5682779128076954`*^9}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{"a", "=", RowBox[{ "ReadList", "[", "\"\\"", "]"}]}], "\[IndentingNewLine]", RowBox[{"nb", "=", RowBox[{"Length", "[", "a", "]"}]}], "\[IndentingNewLine]", RowBox[{"cw", "=", RowBox[{"nb", "-", "bh"}]}]}], "Input", CellChangeTimes->{{3.5681343270045934`*^9, 3.568134331481801*^9}, 3.5681343823222904`*^9, {3.5681345934218616`*^9, 3.5681345986322703`*^9}, { 3.568134635479535*^9, 3.5681346521559644`*^9}, {3.5681351451792307`*^9, 3.568135145428831*^9}, {3.5681351781888885`*^9, 3.5681352349417877`*^9}, { 3.568135340163973*^9, 3.568135340335573*^9}, {3.568135523495495*^9, 3.5681355244314966`*^9}, {3.5681357674799232`*^9, 3.5681357771675406`*^9}, {3.568136284433632*^9, 3.5681362987700567`*^9}, { 3.568136395053426*^9, 3.5681363965510283`*^9}, {3.568381170343589*^9, 3.568381179844005*^9}, {3.5689533222924633`*^9, 3.568953376658559*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"2", ",", "1", ",", "1.210314249067387`"}], "}"}], ",", RowBox[{"{", RowBox[{"3", ",", "1", ",", "1.2103203349966232`"}], "}"}], ",", RowBox[{"{", RowBox[{"4", ",", "1", ",", "1.2103545683012018`"}], "}"}], ",", RowBox[{"{", RowBox[{"5", ",", "1", ",", "1.2102326842173277`"}], "}"}]}], "}"}]], "Output", CellChangeTimes->{ 3.5689997298229895`*^9, 3.569060092757221*^9, 3.569065230468486*^9, 3.5690652770189676`*^9, {3.5690653772491436`*^9, 3.569065406436795*^9}, 3.5690655621718683`*^9, 3.569065643978412*^9, 3.5690658377775526`*^9, 3.5690663622972736`*^9, 3.5690664047761483`*^9, {3.5690665344903765`*^9, 3.5690665557688136`*^9}, {3.5690666165309205`*^9, 3.5690666344397516`*^9}, {3.5690666761542253`*^9, 3.5690667056850767`*^9}, {3.569066750815956*^9, 3.5690667708619914`*^9}, { 3.569066820610479*^9, 3.569066847442526*^9}, 3.5690668808421845`*^9, 3.5690694998091574`*^9, 3.5690695664524746`*^9, 3.5690697218755474`*^9, 3.569069809984502*^9, 3.5690698553181815`*^9, 3.5690699457671404`*^9, 3.569070175914345*^9, 3.569070250700876*^9, 3.5690704404752097`*^9, 3.5690707656265807`*^9, 3.5690715253563156`*^9, 3.5690715827488165`*^9}], Cell[BoxData["4"], "Output", CellChangeTimes->{ 3.5689997298229895`*^9, 3.569060092757221*^9, 3.569065230468486*^9, 3.5690652770189676`*^9, {3.5690653772491436`*^9, 3.569065406436795*^9}, 3.5690655621718683`*^9, 3.569065643978412*^9, 3.5690658377775526`*^9, 3.5690663622972736`*^9, 3.5690664047761483`*^9, {3.5690665344903765`*^9, 3.5690665557688136`*^9}, {3.5690666165309205`*^9, 3.5690666344397516`*^9}, {3.5690666761542253`*^9, 3.5690667056850767`*^9}, {3.569066750815956*^9, 3.5690667708619914`*^9}, { 3.569066820610479*^9, 3.569066847442526*^9}, 3.5690668808421845`*^9, 3.5690694998091574`*^9, 3.5690695664524746`*^9, 3.5690697218755474`*^9, 3.569069809984502*^9, 3.5690698553181815`*^9, 3.5690699457671404`*^9, 3.569070175914345*^9, 3.569070250700876*^9, 3.5690704404752097`*^9, 3.5690707656265807`*^9, 3.5690715253563156`*^9, 3.569071582764416*^9}], Cell[BoxData["0"], "Output", CellChangeTimes->{ 3.5689997298229895`*^9, 3.569060092757221*^9, 3.569065230468486*^9, 3.5690652770189676`*^9, {3.5690653772491436`*^9, 3.569065406436795*^9}, 3.5690655621718683`*^9, 3.569065643978412*^9, 3.5690658377775526`*^9, 3.5690663622972736`*^9, 3.5690664047761483`*^9, {3.5690665344903765`*^9, 3.5690665557688136`*^9}, {3.5690666165309205`*^9, 3.5690666344397516`*^9}, {3.5690666761542253`*^9, 3.5690667056850767`*^9}, {3.569066750815956*^9, 3.5690667708619914`*^9}, { 3.569066820610479*^9, 3.569066847442526*^9}, 3.5690668808421845`*^9, 3.5690694998091574`*^9, 3.5690695664524746`*^9, 3.5690697218755474`*^9, 3.569069809984502*^9, 3.5690698553181815`*^9, 3.5690699457671404`*^9, 3.569070175914345*^9, 3.569070250700876*^9, 3.5690704404752097`*^9, 3.5690707656265807`*^9, 3.5690715253563156`*^9, 3.569071582764416*^9}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"wxc", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "\[NotEqual]", "\"\\""}], ",", RowBox[{ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "3"}], "]"}], "]"}], "+", RowBox[{"u2", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}]}], ")"}], "/", RowBox[{"(", RowBox[{ RowBox[{"u2", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "+", RowBox[{"2", "*", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{"u2", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ")"}]}], "*", RowBox[{"(", RowBox[{ RowBox[{"x", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}], "-", RowBox[{"x", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}]}], ")"}]}], "+", RowBox[{"x", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}]}], ",", "0"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "nb"}], "}"}]}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.5681365021632137`*^9, 3.5681365171392403`*^9}, { 3.568136560351316*^9, 3.5681365863097615`*^9}, {3.5681369135895367`*^9, 3.5681369265999594`*^9}, {3.568137821916732*^9, 3.5681378378443604`*^9}, { 3.5682072492149415`*^9, 3.5682072529901485`*^9}, 3.568209802830227*^9, { 3.5682216503836985`*^9, 3.5682217583046885`*^9}, {3.5682217912831464`*^9, 3.568221843465238*^9}, {3.5682219760654707`*^9, 3.5682221133457117`*^9}, { 3.568222159615393*^9, 3.5682221693498106`*^9}, {3.5682222565227633`*^9, 3.5682222578487663`*^9}, 3.5682224736283445`*^9, 3.5682227227607822`*^9, 3.5682227558640404`*^9, {3.568222818810151*^9, 3.5682228486062036`*^9}, { 3.5682230057764797`*^9, 3.5682230079916835`*^9}, {3.568223041952943*^9, 3.5682230548229656`*^9}, 3.5682232619601297`*^9, 3.5682233000709963`*^9, { 3.568223536645412*^9, 3.5682235381118145`*^9}, {3.568223591151908*^9, 3.568223591526308*^9}, {3.5682236491528096`*^9, 3.568223657483224*^9}, { 3.568223713222122*^9, 3.5682237532361927`*^9}, {3.568223800379475*^9, 3.568223827648323*^9}, {3.5682238887380304`*^9, 3.5682238893152313`*^9}, { 3.5682240287170763`*^9, 3.5682240612431335`*^9}, 3.5682241088076167`*^9, { 3.568224198772975*^9, 3.5682242396450467`*^9}, {3.568224282170721*^9, 3.5682243185031853`*^9}, {3.568224437890195*^9, 3.568224439590598*^9}, { 3.568224572253231*^9, 3.5682245867144566`*^9}, {3.5682247081762695`*^9, 3.568224733557514*^9}, {3.56822480520844*^9, 3.5682248433505073`*^9}, { 3.568224884206979*^9, 3.5682248959069996`*^9}, {3.568224982830352*^9, 3.5682250047639904`*^9}, {3.56822503870965*^9, 3.5682250803305235`*^9}, { 3.5682251223413973`*^9, 3.568225141451431*^9}, {3.568225234177994*^9, 3.5682252929588966`*^9}, {3.5682253775890455`*^9, 3.568225482140429*^9}, { 3.56822551694409*^9, 3.5682255297205124`*^9}, {3.568225573197789*^9, 3.5682255809666023`*^9}, 3.568225660698343*^9, 3.5682257301964645`*^9, { 3.568225765078126*^9, 3.5682257828933573`*^9}, {3.568225927026843*^9, 3.568225927416844*^9}, {3.568267062420238*^9, 3.568267081358671*^9}, { 3.568267188328059*^9, 3.5682672830358257`*^9}, {3.5682676367040462`*^9, 3.568267747979042*^9}, 3.5682678404716043`*^9, 3.568267882030078*^9, { 3.5682679252577534`*^9, 3.568267996612279*^9}, {3.5682680488879704`*^9, 3.568268061133992*^9}, {3.568268100914062*^9, 3.5682681150944867`*^9}, { 3.5682681672609787`*^9, 3.5682681709113846`*^9}, {3.568268272904364*^9, 3.5682683182536435`*^9}, {3.5682685202739983`*^9, 3.568268599194537*^9}, { 3.5682686348873997`*^9, 3.568268637773405*^9}, {3.5682686944951043`*^9, 3.568268872273017*^9}, {3.5682689420207396`*^9, 3.568269016947671*^9}, { 3.5682695040741262`*^9, 3.568269601153097*^9}, {3.5682696652380095`*^9, 3.5682696728352227`*^9}, {3.5682697489633565`*^9, 3.5682698308635006`*^9}, {3.568271156491429*^9, 3.568271160063835*^9}, { 3.568271244927984*^9, 3.568271369010602*^9}, {3.5682716250070515`*^9, 3.5682717312276382`*^9}, {3.5682717843301315`*^9, 3.5682718576814604`*^9}, {3.568271940673606*^9, 3.5682720536490045`*^9}, { 3.5682722259669075`*^9, 3.5682723223126764`*^9}, {3.568272370813162*^9, 3.568272411466833*^9}, {3.5682725893539457`*^9, 3.568272622238803*^9}, { 3.5682726936713285`*^9, 3.5682726960113325`*^9}, {3.568272741079812*^9, 3.56827275143823*^9}, {3.5682729218373294`*^9, 3.568272991070251*^9}, { 3.5682734285574193`*^9, 3.5682734302890224`*^9}, {3.5682734898811274`*^9, 3.5682735111595645`*^9}, {3.5682736115925407`*^9, 3.5682736313577757`*^9}, {3.5682738601945777`*^9, 3.568273987272401*^9}, { 3.5682763944346285`*^9, 3.5682764313910933`*^9}, 3.5682765128232365`*^9, { 3.5682766527398825`*^9, 3.568276660976697*^9}, {3.5682774496048822`*^9, 3.568277461663703*^9}, {3.5683784032731285`*^9, 3.568378422991563*^9}, { 3.568378459308427*^9, 3.5683784812420654`*^9}, {3.568378513440522*^9, 3.568378867046343*^9}, {3.568378915468828*^9, 3.56837896203491*^9}, { 3.5683790764299107`*^9, 3.568379270041851*^9}, {3.568380233530743*^9, 3.568380234388745*^9}, 3.5683804214954734`*^9, 3.568380471867962*^9, 3.5683805399464817`*^9, {3.5683812396077104`*^9, 3.568381293537005*^9}, { 3.568381339541486*^9, 3.5683813935175805`*^9}, 3.5684646396646404`*^9, 3.5684646988823442`*^9, 3.568465086387025*^9, {3.5684655859935026`*^9, 3.568465605087936*^9}, 3.568821235038273*^9, 3.5688233741452303`*^9, 3.568823449742963*^9, 3.5688274593312387`*^9, {3.5688284679978104`*^9, 3.5688284689026117`*^9}, {3.5688285131598897`*^9, 3.5688285201487017`*^9}, 3.5688285866360188`*^9, {3.5688292397623663`*^9, 3.568829246688778*^9}, { 3.5688780328231993`*^9, 3.568878075458074*^9}, 3.5688953020257144`*^9}], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"wyc", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "\[NotEqual]", "\"\\""}], ",", RowBox[{ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "3"}], "]"}], "]"}], "+", RowBox[{"u2", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}]}], ")"}], "/", RowBox[{"(", RowBox[{ RowBox[{"u2", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "+", RowBox[{"2", "*", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{"u2", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ")"}]}], "*", RowBox[{"(", RowBox[{ RowBox[{"y", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}], "-", RowBox[{"y", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}]}], ")"}]}], "+", RowBox[{"y", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}]}], ",", "0"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "nb"}], "}"}]}], "]"}]}], ";"}], " "}]], "Input", CellChangeTimes->{{3.5681365021632137`*^9, 3.5681365171392403`*^9}, { 3.568136560351316*^9, 3.5681365863097615`*^9}, {3.5681369135895367`*^9, 3.5681369265999594`*^9}, {3.568137821916732*^9, 3.5681378378443604`*^9}, { 3.5682072492149415`*^9, 3.5682072529901485`*^9}, 3.568209802830227*^9, { 3.5682216503836985`*^9, 3.5682217583046885`*^9}, {3.5682217912831464`*^9, 3.568221843465238*^9}, {3.5682219760654707`*^9, 3.5682221133457117`*^9}, { 3.568222159615393*^9, 3.5682221693498106`*^9}, {3.5682222565227633`*^9, 3.5682222578487663`*^9}, 3.5682224736283445`*^9, 3.5682227227607822`*^9, 3.5682227558640404`*^9, {3.568222818810151*^9, 3.5682228486062036`*^9}, { 3.5682230057764797`*^9, 3.5682230079916835`*^9}, {3.568223041952943*^9, 3.5682230548229656`*^9}, 3.5682232619601297`*^9, 3.5682233000709963`*^9, { 3.568223536645412*^9, 3.5682235381118145`*^9}, {3.568223591151908*^9, 3.568223591526308*^9}, {3.5682236491528096`*^9, 3.568223657483224*^9}, { 3.568223713222122*^9, 3.5682237532361927`*^9}, {3.568223800379475*^9, 3.568223827648323*^9}, {3.5682238887380304`*^9, 3.5682238893152313`*^9}, { 3.5682240287170763`*^9, 3.5682240612431335`*^9}, 3.5682241088076167`*^9, { 3.568224198772975*^9, 3.5682242396450467`*^9}, {3.568224282170721*^9, 3.5682243185031853`*^9}, {3.568224437890195*^9, 3.568224439590598*^9}, { 3.568224572253231*^9, 3.5682245867144566`*^9}, {3.5682247081762695`*^9, 3.568224733557514*^9}, {3.56822480520844*^9, 3.5682248433505073`*^9}, { 3.568224884206979*^9, 3.5682248959069996`*^9}, {3.568224982830352*^9, 3.5682250047639904`*^9}, {3.56822503870965*^9, 3.5682250803305235`*^9}, { 3.5682251223413973`*^9, 3.568225141451431*^9}, {3.568225234177994*^9, 3.5682252929588966`*^9}, {3.5682253775890455`*^9, 3.568225482140429*^9}, { 3.56822551694409*^9, 3.5682255297205124`*^9}, {3.568225573197789*^9, 3.5682255809666023`*^9}, 3.568225660698343*^9, 3.5682257301964645`*^9, { 3.568225765078126*^9, 3.5682257828933573`*^9}, {3.568225927026843*^9, 3.568225927416844*^9}, {3.568267062420238*^9, 3.568267081358671*^9}, { 3.568267188328059*^9, 3.5682672830358257`*^9}, {3.5682676367040462`*^9, 3.568267747979042*^9}, 3.5682678404716043`*^9, 3.568267882030078*^9, { 3.5682679252577534`*^9, 3.568267996612279*^9}, {3.5682680488879704`*^9, 3.568268061133992*^9}, {3.568268100914062*^9, 3.5682681150944867`*^9}, { 3.5682681672609787`*^9, 3.5682681709113846`*^9}, {3.568268272904364*^9, 3.5682683182536435`*^9}, {3.5682685202739983`*^9, 3.568268599194537*^9}, { 3.5682686348873997`*^9, 3.568268637773405*^9}, {3.5682686944951043`*^9, 3.568268872273017*^9}, {3.5682689420207396`*^9, 3.568269016947671*^9}, { 3.5682695040741262`*^9, 3.568269601153097*^9}, {3.5682696652380095`*^9, 3.5682696728352227`*^9}, {3.5682697489633565`*^9, 3.5682698308635006`*^9}, {3.568271156491429*^9, 3.568271160063835*^9}, { 3.568271244927984*^9, 3.568271369010602*^9}, {3.5682716250070515`*^9, 3.5682717312276382`*^9}, {3.5682717843301315`*^9, 3.5682718576814604`*^9}, {3.568271940673606*^9, 3.5682720536490045`*^9}, { 3.5682722259669075`*^9, 3.5682723223126764`*^9}, {3.568272370813162*^9, 3.568272411466833*^9}, {3.5682725893539457`*^9, 3.568272622238803*^9}, { 3.5682726936713285`*^9, 3.5682726960113325`*^9}, {3.568272741079812*^9, 3.56827275143823*^9}, {3.5682729218373294`*^9, 3.568272991070251*^9}, { 3.5682734285574193`*^9, 3.5682734302890224`*^9}, {3.5682734898811274`*^9, 3.5682735111595645`*^9}, {3.5682736115925407`*^9, 3.5682736313577757`*^9}, {3.5682738601945777`*^9, 3.568273987272401*^9}, { 3.5682763944346285`*^9, 3.5682764313910933`*^9}, 3.5682765128232365`*^9, { 3.5682766527398825`*^9, 3.568276660976697*^9}, {3.5682774496048822`*^9, 3.568277461663703*^9}, {3.5683784032731285`*^9, 3.568378422991563*^9}, { 3.568378459308427*^9, 3.5683784812420654`*^9}, {3.568378513440522*^9, 3.568378867046343*^9}, {3.568378915468828*^9, 3.56837896203491*^9}, { 3.5683790764299107`*^9, 3.568379270041851*^9}, {3.568380233530743*^9, 3.568380234388745*^9}, 3.5683804214954734`*^9, 3.568380471867962*^9, 3.5683805399464817`*^9, {3.5683812396077104`*^9, 3.568381293537005*^9}, { 3.568381339541486*^9, 3.5683813935175805`*^9}, 3.5684646396646404`*^9, 3.5684646988823442`*^9, 3.568465086387025*^9, {3.5684655859935026`*^9, 3.568465605087936*^9}, 3.568821235038273*^9, 3.5688233741452303`*^9, 3.568823449742963*^9, 3.5688274593312387`*^9, {3.5688284679978104`*^9, 3.5688284689026117`*^9}, {3.5688285131598897`*^9, 3.5688285201487017`*^9}, 3.5688285866360188`*^9, {3.5688292397623663`*^9, 3.568829246688778*^9}, { 3.5688780328231993`*^9, 3.5688780904809003`*^9}, 3.5688953042565184`*^9}], Cell[BoxData[ RowBox[{ RowBox[{"wzc", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "\[NotEqual]", "\"\\""}], ",", RowBox[{ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "3"}], "]"}], "]"}], "+", RowBox[{"u2", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}]}], ")"}], "/", RowBox[{"(", RowBox[{ RowBox[{"u2", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "+", RowBox[{"2", "*", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "3"}], "]"}], "]"}]}], "+", RowBox[{"u2", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ")"}]}], "*", RowBox[{"(", RowBox[{ RowBox[{"z", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}], "-", RowBox[{"z", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}]}], ")"}]}], "+", RowBox[{"z", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}]}], ",", "0"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "nb"}], "}"}]}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.5681365021632137`*^9, 3.5681365171392403`*^9}, { 3.568136560351316*^9, 3.5681365863097615`*^9}, {3.5681369135895367`*^9, 3.5681369265999594`*^9}, {3.568137821916732*^9, 3.5681378378443604`*^9}, { 3.5682072492149415`*^9, 3.5682072529901485`*^9}, 3.568209802830227*^9, { 3.5682216503836985`*^9, 3.5682217583046885`*^9}, {3.5682217912831464`*^9, 3.568221843465238*^9}, {3.5682219760654707`*^9, 3.5682221133457117`*^9}, { 3.568222159615393*^9, 3.5682221693498106`*^9}, {3.5682222565227633`*^9, 3.5682222578487663`*^9}, 3.5682224736283445`*^9, 3.5682227227607822`*^9, 3.5682227558640404`*^9, {3.568222818810151*^9, 3.5682228486062036`*^9}, { 3.5682230057764797`*^9, 3.5682230079916835`*^9}, {3.568223041952943*^9, 3.5682230548229656`*^9}, 3.5682232619601297`*^9, 3.5682233000709963`*^9, { 3.568223536645412*^9, 3.5682235381118145`*^9}, {3.568223591151908*^9, 3.568223591526308*^9}, {3.5682236491528096`*^9, 3.568223657483224*^9}, { 3.568223713222122*^9, 3.5682237532361927`*^9}, {3.568223800379475*^9, 3.568223827648323*^9}, {3.5682238887380304`*^9, 3.5682238893152313`*^9}, { 3.5682240287170763`*^9, 3.5682240612431335`*^9}, 3.5682241088076167`*^9, { 3.568224198772975*^9, 3.5682242396450467`*^9}, {3.568224282170721*^9, 3.5682243185031853`*^9}, {3.568224437890195*^9, 3.568224439590598*^9}, { 3.568224572253231*^9, 3.5682245867144566`*^9}, {3.5682247081762695`*^9, 3.568224733557514*^9}, {3.56822480520844*^9, 3.5682248433505073`*^9}, { 3.568224884206979*^9, 3.5682248959069996`*^9}, {3.568224982830352*^9, 3.5682250047639904`*^9}, {3.56822503870965*^9, 3.5682250803305235`*^9}, { 3.5682251223413973`*^9, 3.568225141451431*^9}, {3.568225234177994*^9, 3.5682252929588966`*^9}, {3.5682253775890455`*^9, 3.568225482140429*^9}, { 3.56822551694409*^9, 3.5682255297205124`*^9}, {3.568225573197789*^9, 3.5682255809666023`*^9}, 3.568225660698343*^9, 3.5682257301964645`*^9, { 3.568225765078126*^9, 3.5682257828933573`*^9}, {3.568225927026843*^9, 3.568225927416844*^9}, {3.568267062420238*^9, 3.568267081358671*^9}, { 3.568267188328059*^9, 3.5682672830358257`*^9}, {3.5682676367040462`*^9, 3.568267747979042*^9}, 3.5682678404716043`*^9, 3.568267882030078*^9, { 3.5682679252577534`*^9, 3.568267996612279*^9}, {3.5682680488879704`*^9, 3.568268061133992*^9}, {3.568268100914062*^9, 3.5682681150944867`*^9}, { 3.5682681672609787`*^9, 3.5682681709113846`*^9}, {3.568268272904364*^9, 3.5682683182536435`*^9}, {3.5682685202739983`*^9, 3.568268599194537*^9}, { 3.5682686348873997`*^9, 3.568268637773405*^9}, {3.5682686944951043`*^9, 3.568268872273017*^9}, {3.5682689420207396`*^9, 3.568269016947671*^9}, { 3.5682695040741262`*^9, 3.568269601153097*^9}, {3.5682696652380095`*^9, 3.5682696728352227`*^9}, {3.5682697489633565`*^9, 3.5682698308635006`*^9}, {3.568271156491429*^9, 3.568271160063835*^9}, { 3.568271244927984*^9, 3.568271369010602*^9}, {3.5682716250070515`*^9, 3.5682717312276382`*^9}, {3.5682717843301315`*^9, 3.5682718576814604`*^9}, {3.568271940673606*^9, 3.5682720536490045`*^9}, { 3.5682722259669075`*^9, 3.5682723223126764`*^9}, {3.568272370813162*^9, 3.568272411466833*^9}, {3.5682725893539457`*^9, 3.568272622238803*^9}, { 3.5682726936713285`*^9, 3.5682726960113325`*^9}, {3.568272741079812*^9, 3.56827275143823*^9}, {3.5682729218373294`*^9, 3.568272991070251*^9}, { 3.5682734285574193`*^9, 3.5682734302890224`*^9}, {3.5682734898811274`*^9, 3.5682735111595645`*^9}, {3.5682736115925407`*^9, 3.5682736313577757`*^9}, {3.5682738601945777`*^9, 3.568273987272401*^9}, { 3.5682763944346285`*^9, 3.5682764313910933`*^9}, 3.5682765128232365`*^9, { 3.5682766527398825`*^9, 3.568276660976697*^9}, {3.5682774496048822`*^9, 3.568277461663703*^9}, {3.5683784032731285`*^9, 3.568378422991563*^9}, { 3.568378459308427*^9, 3.5683784812420654`*^9}, {3.568378513440522*^9, 3.568378867046343*^9}, {3.568378915468828*^9, 3.56837896203491*^9}, { 3.5683790764299107`*^9, 3.568379270041851*^9}, {3.568380233530743*^9, 3.568380234388745*^9}, 3.5683804214954734`*^9, 3.568380471867962*^9, 3.5683805399464817`*^9, {3.5683812396077104`*^9, 3.568381293537005*^9}, { 3.568381339541486*^9, 3.5683813935175805`*^9}, 3.5684646396646404`*^9, 3.5684646988823442`*^9, 3.568465086387025*^9, {3.5684655859935026`*^9, 3.568465605087936*^9}, 3.568821235038273*^9, 3.5688233741452303`*^9, 3.568823449742963*^9, 3.5688274593312387`*^9, {3.5688284679978104`*^9, 3.5688284689026117`*^9}, {3.5688285131598897`*^9, 3.5688285201487017`*^9}, 3.5688285866360188`*^9, {3.5688292397623663`*^9, 3.568829246688778*^9}, { 3.5688780328231993`*^9, 3.568878099435316*^9}, 3.5688953066901226`*^9}], Cell[BoxData[ RowBox[{ RowBox[{"Plot7", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"Graphics3D", "[", RowBox[{"{", RowBox[{ RowBox[{"Darker", "[", "Gray", "]"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Opacity", "[", "0.5", "]"}], ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"wxc", "[", RowBox[{"[", "i", "]"}], "]"}], "\[NotEqual]", "0"}], ",", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"wxc", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"wyc", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"wzc", "[", RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "3"}], "]"}], "]"}]}], "]"}]}], "]"}]}], "}"}]}], "}"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "nb"}], "}"}]}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.5681365021632137`*^9, 3.5681365171392403`*^9}, { 3.568136560351316*^9, 3.5681365863097615`*^9}, {3.5681369135895367`*^9, 3.5681369265999594`*^9}, {3.568137821916732*^9, 3.5681378378443604`*^9}, { 3.5682072492149415`*^9, 3.5682072529901485`*^9}, 3.568209802830227*^9, { 3.5682216503836985`*^9, 3.5682217583046885`*^9}, {3.5682217912831464`*^9, 3.568221843465238*^9}, {3.5682219760654707`*^9, 3.5682221133457117`*^9}, { 3.568222159615393*^9, 3.5682221693498106`*^9}, {3.5682222565227633`*^9, 3.5682222578487663`*^9}, 3.5682224736283445`*^9, 3.5682227227607822`*^9, 3.5682227558640404`*^9, {3.568222818810151*^9, 3.5682228486062036`*^9}, { 3.5682230057764797`*^9, 3.5682230079916835`*^9}, {3.568223041952943*^9, 3.5682230548229656`*^9}, 3.5682232619601297`*^9, 3.5682233000709963`*^9, { 3.568223536645412*^9, 3.5682235381118145`*^9}, {3.568223591151908*^9, 3.568223591526308*^9}, {3.5682236491528096`*^9, 3.568223657483224*^9}, { 3.568223713222122*^9, 3.5682237532361927`*^9}, {3.568223800379475*^9, 3.568223827648323*^9}, {3.5682238887380304`*^9, 3.5682238893152313`*^9}, { 3.5682240287170763`*^9, 3.5682240612431335`*^9}, 3.5682241088076167`*^9, { 3.568224198772975*^9, 3.5682242396450467`*^9}, {3.568224282170721*^9, 3.5682243185031853`*^9}, {3.568224437890195*^9, 3.568224439590598*^9}, { 3.568224572253231*^9, 3.5682245867144566`*^9}, {3.5682247081762695`*^9, 3.568224733557514*^9}, {3.56822480520844*^9, 3.5682248433505073`*^9}, { 3.568224884206979*^9, 3.5682248959069996`*^9}, {3.568224982830352*^9, 3.5682250047639904`*^9}, {3.56822503870965*^9, 3.5682250803305235`*^9}, { 3.5682251223413973`*^9, 3.568225141451431*^9}, {3.568225234177994*^9, 3.5682252929588966`*^9}, {3.5682253775890455`*^9, 3.568225482140429*^9}, { 3.56822551694409*^9, 3.5682255297205124`*^9}, {3.568225573197789*^9, 3.5682255809666023`*^9}, 3.568225660698343*^9, 3.5682257301964645`*^9, { 3.568225765078126*^9, 3.5682257828933573`*^9}, {3.568225927026843*^9, 3.568225927416844*^9}, {3.568267062420238*^9, 3.568267081358671*^9}, { 3.568267188328059*^9, 3.5682672830358257`*^9}, {3.5682676367040462`*^9, 3.568267747979042*^9}, 3.5682678404716043`*^9, 3.568267882030078*^9, { 3.5682679252577534`*^9, 3.568267996612279*^9}, {3.5682680488879704`*^9, 3.568268061133992*^9}, {3.568268100914062*^9, 3.5682681150944867`*^9}, { 3.5682681672609787`*^9, 3.5682681709113846`*^9}, {3.568268272904364*^9, 3.5682683182536435`*^9}, {3.5682685202739983`*^9, 3.568268599194537*^9}, { 3.5682686348873997`*^9, 3.568268637773405*^9}, {3.5682686944951043`*^9, 3.568268872273017*^9}, {3.5682689420207396`*^9, 3.568269016947671*^9}, { 3.5682695040741262`*^9, 3.568269601153097*^9}, {3.5682696652380095`*^9, 3.5682696728352227`*^9}, {3.5682697489633565`*^9, 3.5682698308635006`*^9}, {3.568271156491429*^9, 3.568271160063835*^9}, { 3.568271244927984*^9, 3.568271369010602*^9}, {3.5682716250070515`*^9, 3.5682717312276382`*^9}, {3.5682717843301315`*^9, 3.5682718576814604`*^9}, {3.568271940673606*^9, 3.5682720536490045`*^9}, { 3.5682722259669075`*^9, 3.5682723223126764`*^9}, {3.568272370813162*^9, 3.568272411466833*^9}, {3.5682725893539457`*^9, 3.568272622238803*^9}, { 3.5682726936713285`*^9, 3.5682726960113325`*^9}, {3.568272741079812*^9, 3.56827275143823*^9}, {3.5682729218373294`*^9, 3.568272991070251*^9}, { 3.5682734285574193`*^9, 3.5682734302890224`*^9}, {3.5682734898811274`*^9, 3.5682735111595645`*^9}, {3.5682736115925407`*^9, 3.5682736313577757`*^9}, {3.5682738601945777`*^9, 3.568273987272401*^9}, { 3.5682763944346285`*^9, 3.5682764313910933`*^9}, 3.5682765128232365`*^9, { 3.5682766527398825`*^9, 3.568276660976697*^9}, {3.5682774496048822`*^9, 3.568277461663703*^9}, {3.5683784032731285`*^9, 3.568378422991563*^9}, { 3.568378459308427*^9, 3.5683784812420654`*^9}, {3.568378513440522*^9, 3.568378867046343*^9}, {3.568378915468828*^9, 3.56837896203491*^9}, { 3.5683790764299107`*^9, 3.568379270041851*^9}, {3.568380233530743*^9, 3.568380234388745*^9}, 3.5683804214954734`*^9, 3.568380471867962*^9, 3.5683805399464817`*^9, {3.5683812396077104`*^9, 3.568381293537005*^9}, { 3.568381339541486*^9, 3.5683813935175805`*^9}, 3.5684646396646404`*^9, 3.5684646988823442`*^9, 3.568465086387025*^9, {3.5684655859935026`*^9, 3.568465605087936*^9}, 3.568821235038273*^9, 3.5688233741452303`*^9, 3.568823449742963*^9, 3.5688274593312387`*^9, {3.5688284679978104`*^9, 3.5688284689026117`*^9}, {3.5688285131598897`*^9, 3.5688285201487017`*^9}, 3.5688285866360188`*^9, {3.5688292397623663`*^9, 3.568829246688778*^9}, { 3.5688780328231993`*^9, 3.56887811871695*^9}, {3.568895331556566*^9, 3.5688953367201753`*^9}, 3.568895995276532*^9}], Cell[BoxData[ RowBox[{ RowBox[{"wxh", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "\[Equal]", "\"\\""}], ",", RowBox[{ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "3"}], "]"}], "]"}], "+", RowBox[{"u2", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ")"}], "/", RowBox[{"(", RowBox[{ RowBox[{"u2", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "+", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "3"}], "]"}], "]"}], "+", RowBox[{"u2", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ")"}]}], "*", RowBox[{"(", RowBox[{ RowBox[{"x", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "-", RowBox[{"x", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ")"}]}], "+", RowBox[{"x", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ",", "0"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "nb"}], "}"}]}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.5681365021632137`*^9, 3.5681365171392403`*^9}, { 3.568136560351316*^9, 3.5681365863097615`*^9}, {3.5681369135895367`*^9, 3.5681369265999594`*^9}, {3.568137821916732*^9, 3.5681378378443604`*^9}, { 3.5682072492149415`*^9, 3.5682072529901485`*^9}, 3.568209802830227*^9, { 3.5682216503836985`*^9, 3.5682217583046885`*^9}, {3.5682217912831464`*^9, 3.568221843465238*^9}, {3.5682219760654707`*^9, 3.5682221133457117`*^9}, { 3.568222159615393*^9, 3.5682221693498106`*^9}, {3.5682222565227633`*^9, 3.5682222578487663`*^9}, 3.5682224736283445`*^9, 3.5682227227607822`*^9, 3.5682227558640404`*^9, {3.568222818810151*^9, 3.5682228486062036`*^9}, { 3.5682230057764797`*^9, 3.5682230079916835`*^9}, {3.568223041952943*^9, 3.5682230548229656`*^9}, 3.5682232619601297`*^9, 3.5682233000709963`*^9, { 3.568223536645412*^9, 3.5682235381118145`*^9}, {3.568223591151908*^9, 3.568223591526308*^9}, {3.5682236491528096`*^9, 3.568223657483224*^9}, { 3.568223713222122*^9, 3.5682237532361927`*^9}, {3.568223800379475*^9, 3.568223827648323*^9}, {3.5682238887380304`*^9, 3.5682238893152313`*^9}, { 3.5682240287170763`*^9, 3.5682240612431335`*^9}, 3.5682241088076167`*^9, { 3.568224198772975*^9, 3.5682242396450467`*^9}, {3.568224282170721*^9, 3.5682243185031853`*^9}, {3.568224437890195*^9, 3.568224439590598*^9}, { 3.568224572253231*^9, 3.5682245867144566`*^9}, {3.5682247081762695`*^9, 3.568224733557514*^9}, {3.56822480520844*^9, 3.5682248433505073`*^9}, { 3.568224884206979*^9, 3.5682248959069996`*^9}, {3.568224982830352*^9, 3.5682250047639904`*^9}, {3.56822503870965*^9, 3.5682250803305235`*^9}, { 3.5682251223413973`*^9, 3.568225141451431*^9}, {3.568225234177994*^9, 3.5682252929588966`*^9}, {3.5682253775890455`*^9, 3.568225482140429*^9}, { 3.56822551694409*^9, 3.5682255297205124`*^9}, {3.568225573197789*^9, 3.5682255809666023`*^9}, 3.568225660698343*^9, 3.5682257301964645`*^9, { 3.568225765078126*^9, 3.5682257828933573`*^9}, {3.568225927026843*^9, 3.568225927416844*^9}, {3.568267062420238*^9, 3.568267081358671*^9}, { 3.568267188328059*^9, 3.5682672830358257`*^9}, {3.5682676367040462`*^9, 3.568267747979042*^9}, 3.5682678404716043`*^9, 3.568267882030078*^9, { 3.5682679252577534`*^9, 3.568267996612279*^9}, {3.5682680488879704`*^9, 3.568268061133992*^9}, {3.568268100914062*^9, 3.5682681150944867`*^9}, { 3.5682681672609787`*^9, 3.5682681709113846`*^9}, {3.568268272904364*^9, 3.5682683182536435`*^9}, {3.5682685202739983`*^9, 3.568268599194537*^9}, { 3.5682686348873997`*^9, 3.568268637773405*^9}, {3.5682686944951043`*^9, 3.568268872273017*^9}, {3.5682689420207396`*^9, 3.568269016947671*^9}, { 3.5682695040741262`*^9, 3.568269601153097*^9}, {3.5682696652380095`*^9, 3.5682696728352227`*^9}, {3.5682697489633565`*^9, 3.5682698308635006`*^9}, {3.568271156491429*^9, 3.568271160063835*^9}, { 3.568271244927984*^9, 3.568271369010602*^9}, {3.5682716250070515`*^9, 3.5682717312276382`*^9}, {3.5682717843301315`*^9, 3.5682718576814604`*^9}, {3.568271940673606*^9, 3.5682720536490045`*^9}, { 3.5682722259669075`*^9, 3.5682723223126764`*^9}, {3.568272370813162*^9, 3.568272411466833*^9}, {3.5682725893539457`*^9, 3.568272622238803*^9}, { 3.5682726936713285`*^9, 3.5682726960113325`*^9}, {3.568272741079812*^9, 3.56827275143823*^9}, {3.5682729218373294`*^9, 3.568272991070251*^9}, { 3.5682734285574193`*^9, 3.5682734302890224`*^9}, {3.5682734898811274`*^9, 3.5682735111595645`*^9}, {3.5682736115925407`*^9, 3.5682736313577757`*^9}, {3.5682738601945777`*^9, 3.568273987272401*^9}, { 3.5682763944346285`*^9, 3.5682764313910933`*^9}, 3.5682765128232365`*^9, { 3.5682766527398825`*^9, 3.568276660976697*^9}, {3.5682774496048822`*^9, 3.568277461663703*^9}, {3.5683784032731285`*^9, 3.568378422991563*^9}, { 3.568378459308427*^9, 3.5683784812420654`*^9}, {3.568378513440522*^9, 3.568378867046343*^9}, {3.568378915468828*^9, 3.56837896203491*^9}, { 3.5683790764299107`*^9, 3.568379270041851*^9}, {3.568380233530743*^9, 3.568380234388745*^9}, 3.5683804214954734`*^9, 3.568380471867962*^9, 3.5683805399464817`*^9, {3.5683812396077104`*^9, 3.568381293537005*^9}, { 3.568381339541486*^9, 3.5683813935175805`*^9}, 3.5684646396646404`*^9, 3.5684646988823442`*^9, 3.568465086387025*^9, {3.5684655859935026`*^9, 3.568465605087936*^9}, 3.568821235038273*^9, 3.5688233741452303`*^9, 3.568823449742963*^9, 3.5688274593312387`*^9, {3.5688284679978104`*^9, 3.5688284689026117`*^9}, {3.5688285131598897`*^9, 3.5688285201487017`*^9}, 3.5688285866360188`*^9, {3.5688292397623663`*^9, 3.568829246688778*^9}, { 3.5688780328231993`*^9, 3.5688781280925665`*^9}, 3.568895312524533*^9}], Cell[BoxData[ RowBox[{ RowBox[{"wyh", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "\[Equal]", "\"\\""}], ",", RowBox[{ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "3"}], "]"}], "]"}], "+", RowBox[{"u2", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ")"}], "/", RowBox[{"(", RowBox[{ RowBox[{"u2", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "+", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "3"}], "]"}], "]"}], "+", RowBox[{"u2", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ")"}]}], "*", RowBox[{"(", RowBox[{ RowBox[{"y", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "-", RowBox[{"y", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ")"}]}], "+", RowBox[{"y", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ",", "0"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "nb"}], "}"}]}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.5681365021632137`*^9, 3.5681365171392403`*^9}, { 3.568136560351316*^9, 3.5681365863097615`*^9}, {3.5681369135895367`*^9, 3.5681369265999594`*^9}, {3.568137821916732*^9, 3.5681378378443604`*^9}, { 3.5682072492149415`*^9, 3.5682072529901485`*^9}, 3.568209802830227*^9, { 3.5682216503836985`*^9, 3.5682217583046885`*^9}, {3.5682217912831464`*^9, 3.568221843465238*^9}, {3.5682219760654707`*^9, 3.5682221133457117`*^9}, { 3.568222159615393*^9, 3.5682221693498106`*^9}, {3.5682222565227633`*^9, 3.5682222578487663`*^9}, 3.5682224736283445`*^9, 3.5682227227607822`*^9, 3.5682227558640404`*^9, {3.568222818810151*^9, 3.5682228486062036`*^9}, { 3.5682230057764797`*^9, 3.5682230079916835`*^9}, {3.568223041952943*^9, 3.5682230548229656`*^9}, 3.5682232619601297`*^9, 3.5682233000709963`*^9, { 3.568223536645412*^9, 3.5682235381118145`*^9}, {3.568223591151908*^9, 3.568223591526308*^9}, {3.5682236491528096`*^9, 3.568223657483224*^9}, { 3.568223713222122*^9, 3.5682237532361927`*^9}, {3.568223800379475*^9, 3.568223827648323*^9}, {3.5682238887380304`*^9, 3.5682238893152313`*^9}, { 3.5682240287170763`*^9, 3.5682240612431335`*^9}, 3.5682241088076167`*^9, { 3.568224198772975*^9, 3.5682242396450467`*^9}, {3.568224282170721*^9, 3.5682243185031853`*^9}, {3.568224437890195*^9, 3.568224439590598*^9}, { 3.568224572253231*^9, 3.5682245867144566`*^9}, {3.5682247081762695`*^9, 3.568224733557514*^9}, {3.56822480520844*^9, 3.5682248433505073`*^9}, { 3.568224884206979*^9, 3.5682248959069996`*^9}, {3.568224982830352*^9, 3.5682250047639904`*^9}, {3.56822503870965*^9, 3.5682250803305235`*^9}, { 3.5682251223413973`*^9, 3.568225141451431*^9}, {3.568225234177994*^9, 3.5682252929588966`*^9}, {3.5682253775890455`*^9, 3.568225482140429*^9}, { 3.56822551694409*^9, 3.5682255297205124`*^9}, {3.568225573197789*^9, 3.5682255809666023`*^9}, 3.568225660698343*^9, 3.5682257301964645`*^9, { 3.568225765078126*^9, 3.5682257828933573`*^9}, {3.568225927026843*^9, 3.568225927416844*^9}, {3.568267062420238*^9, 3.568267081358671*^9}, { 3.568267188328059*^9, 3.5682672830358257`*^9}, {3.5682676367040462`*^9, 3.568267747979042*^9}, 3.5682678404716043`*^9, 3.568267882030078*^9, { 3.5682679252577534`*^9, 3.568267996612279*^9}, {3.5682680488879704`*^9, 3.568268061133992*^9}, {3.568268100914062*^9, 3.5682681150944867`*^9}, { 3.5682681672609787`*^9, 3.5682681709113846`*^9}, {3.568268272904364*^9, 3.5682683182536435`*^9}, {3.5682685202739983`*^9, 3.568268599194537*^9}, { 3.5682686348873997`*^9, 3.568268637773405*^9}, {3.5682686944951043`*^9, 3.568268872273017*^9}, {3.5682689420207396`*^9, 3.568269016947671*^9}, { 3.5682695040741262`*^9, 3.568269601153097*^9}, {3.5682696652380095`*^9, 3.5682696728352227`*^9}, {3.5682697489633565`*^9, 3.5682698308635006`*^9}, {3.568271156491429*^9, 3.568271160063835*^9}, { 3.568271244927984*^9, 3.568271369010602*^9}, {3.5682716250070515`*^9, 3.5682717312276382`*^9}, {3.5682717843301315`*^9, 3.5682718576814604`*^9}, {3.568271940673606*^9, 3.5682720536490045`*^9}, { 3.5682722259669075`*^9, 3.5682723223126764`*^9}, {3.568272370813162*^9, 3.568272411466833*^9}, {3.5682725893539457`*^9, 3.568272622238803*^9}, { 3.5682726936713285`*^9, 3.5682726960113325`*^9}, {3.568272741079812*^9, 3.56827275143823*^9}, {3.5682729218373294`*^9, 3.568272991070251*^9}, { 3.5682734285574193`*^9, 3.5682734302890224`*^9}, {3.5682734898811274`*^9, 3.5682735111595645`*^9}, {3.5682736115925407`*^9, 3.5682736313577757`*^9}, {3.5682738601945777`*^9, 3.568273987272401*^9}, { 3.5682763944346285`*^9, 3.5682764313910933`*^9}, 3.5682765128232365`*^9, { 3.5682766527398825`*^9, 3.568276660976697*^9}, {3.5682774496048822`*^9, 3.568277461663703*^9}, {3.5683784032731285`*^9, 3.568378422991563*^9}, { 3.568378459308427*^9, 3.5683784812420654`*^9}, {3.568378513440522*^9, 3.568378867046343*^9}, {3.568378915468828*^9, 3.56837896203491*^9}, { 3.5683790764299107`*^9, 3.568379270041851*^9}, {3.568380233530743*^9, 3.568380234388745*^9}, 3.5683804214954734`*^9, 3.568380471867962*^9, 3.5683805399464817`*^9, {3.5683812396077104`*^9, 3.568381293537005*^9}, { 3.568381339541486*^9, 3.5683813935175805`*^9}, 3.5684646396646404`*^9, 3.5684646988823442`*^9, 3.568465086387025*^9, {3.5684655859935026`*^9, 3.568465605087936*^9}, 3.568821235038273*^9, 3.5688233741452303`*^9, 3.568823449742963*^9, 3.5688274593312387`*^9, {3.5688284679978104`*^9, 3.5688284689026117`*^9}, {3.5688285131598897`*^9, 3.5688285201487017`*^9}, 3.5688285866360188`*^9, {3.5688292397623663`*^9, 3.568829246688778*^9}, { 3.5688780328231993`*^9, 3.5688781383261843`*^9}, 3.5688953148957367`*^9}], Cell[BoxData[ RowBox[{ RowBox[{"wzh", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"u", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "\[Equal]", "\"\\""}], ",", RowBox[{ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "3"}], "]"}], "]"}], "+", RowBox[{"u2", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ")"}], "/", RowBox[{"(", RowBox[{ RowBox[{"u2", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "+", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "3"}], "]"}], "]"}], "+", RowBox[{"u2", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ")"}]}], "*", RowBox[{"(", RowBox[{ RowBox[{"z", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "-", RowBox[{"z", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ")"}]}], "+", RowBox[{"z", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ",", "0"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "nb"}], "}"}]}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.5681365021632137`*^9, 3.5681365171392403`*^9}, { 3.568136560351316*^9, 3.5681365863097615`*^9}, {3.5681369135895367`*^9, 3.5681369265999594`*^9}, {3.568137821916732*^9, 3.5681378378443604`*^9}, { 3.5682072492149415`*^9, 3.5682072529901485`*^9}, 3.568209802830227*^9, { 3.5682216503836985`*^9, 3.5682217583046885`*^9}, {3.5682217912831464`*^9, 3.568221843465238*^9}, {3.5682219760654707`*^9, 3.5682221133457117`*^9}, { 3.568222159615393*^9, 3.5682221693498106`*^9}, {3.5682222565227633`*^9, 3.5682222578487663`*^9}, 3.5682224736283445`*^9, 3.5682227227607822`*^9, 3.5682227558640404`*^9, {3.568222818810151*^9, 3.5682228486062036`*^9}, { 3.5682230057764797`*^9, 3.5682230079916835`*^9}, {3.568223041952943*^9, 3.5682230548229656`*^9}, 3.5682232619601297`*^9, 3.5682233000709963`*^9, { 3.568223536645412*^9, 3.5682235381118145`*^9}, {3.568223591151908*^9, 3.568223591526308*^9}, {3.5682236491528096`*^9, 3.568223657483224*^9}, { 3.568223713222122*^9, 3.5682237532361927`*^9}, {3.568223800379475*^9, 3.568223827648323*^9}, {3.5682238887380304`*^9, 3.5682238893152313`*^9}, { 3.5682240287170763`*^9, 3.5682240612431335`*^9}, 3.5682241088076167`*^9, { 3.568224198772975*^9, 3.5682242396450467`*^9}, {3.568224282170721*^9, 3.5682243185031853`*^9}, {3.568224437890195*^9, 3.568224439590598*^9}, { 3.568224572253231*^9, 3.5682245867144566`*^9}, {3.5682247081762695`*^9, 3.568224733557514*^9}, {3.56822480520844*^9, 3.5682248433505073`*^9}, { 3.568224884206979*^9, 3.5682248959069996`*^9}, {3.568224982830352*^9, 3.5682250047639904`*^9}, {3.56822503870965*^9, 3.5682250803305235`*^9}, { 3.5682251223413973`*^9, 3.568225141451431*^9}, {3.568225234177994*^9, 3.5682252929588966`*^9}, {3.5682253775890455`*^9, 3.568225482140429*^9}, { 3.56822551694409*^9, 3.5682255297205124`*^9}, {3.568225573197789*^9, 3.5682255809666023`*^9}, 3.568225660698343*^9, 3.5682257301964645`*^9, { 3.568225765078126*^9, 3.5682257828933573`*^9}, {3.568225927026843*^9, 3.568225927416844*^9}, {3.568267062420238*^9, 3.568267081358671*^9}, { 3.568267188328059*^9, 3.5682672830358257`*^9}, {3.5682676367040462`*^9, 3.568267747979042*^9}, 3.5682678404716043`*^9, 3.568267882030078*^9, { 3.5682679252577534`*^9, 3.568267996612279*^9}, {3.5682680488879704`*^9, 3.568268061133992*^9}, {3.568268100914062*^9, 3.5682681150944867`*^9}, { 3.5682681672609787`*^9, 3.5682681709113846`*^9}, {3.568268272904364*^9, 3.5682683182536435`*^9}, {3.5682685202739983`*^9, 3.568268599194537*^9}, { 3.5682686348873997`*^9, 3.568268637773405*^9}, {3.5682686944951043`*^9, 3.568268872273017*^9}, {3.5682689420207396`*^9, 3.568269016947671*^9}, { 3.5682695040741262`*^9, 3.568269601153097*^9}, {3.5682696652380095`*^9, 3.5682696728352227`*^9}, {3.5682697489633565`*^9, 3.5682698308635006`*^9}, {3.568271156491429*^9, 3.568271160063835*^9}, { 3.568271244927984*^9, 3.568271369010602*^9}, {3.5682716250070515`*^9, 3.5682717312276382`*^9}, {3.5682717843301315`*^9, 3.5682718576814604`*^9}, {3.568271940673606*^9, 3.5682720536490045`*^9}, { 3.5682722259669075`*^9, 3.5682723223126764`*^9}, {3.568272370813162*^9, 3.568272411466833*^9}, {3.5682725893539457`*^9, 3.568272622238803*^9}, { 3.5682726936713285`*^9, 3.5682726960113325`*^9}, {3.568272741079812*^9, 3.56827275143823*^9}, {3.5682729218373294`*^9, 3.568272991070251*^9}, { 3.5682734285574193`*^9, 3.5682734302890224`*^9}, {3.5682734898811274`*^9, 3.5682735111595645`*^9}, {3.5682736115925407`*^9, 3.5682736313577757`*^9}, {3.5682738601945777`*^9, 3.568273987272401*^9}, { 3.5682763944346285`*^9, 3.5682764313910933`*^9}, 3.5682765128232365`*^9, { 3.5682766527398825`*^9, 3.568276660976697*^9}, {3.5682774496048822`*^9, 3.568277461663703*^9}, {3.5683784032731285`*^9, 3.568378422991563*^9}, { 3.568378459308427*^9, 3.5683784812420654`*^9}, {3.568378513440522*^9, 3.568378867046343*^9}, {3.568378915468828*^9, 3.56837896203491*^9}, { 3.5683790764299107`*^9, 3.568379270041851*^9}, {3.568380233530743*^9, 3.568380234388745*^9}, 3.5683804214954734`*^9, 3.568380471867962*^9, 3.5683805399464817`*^9, {3.5683812396077104`*^9, 3.568381293537005*^9}, { 3.568381339541486*^9, 3.5683813935175805`*^9}, 3.5684646396646404`*^9, 3.5684646988823442`*^9, 3.568465086387025*^9, {3.5684655859935026`*^9, 3.568465605087936*^9}, 3.568821235038273*^9, 3.5688233741452303`*^9, 3.568823449742963*^9, 3.5688274593312387`*^9, {3.5688284679978104`*^9, 3.5688284689026117`*^9}, {3.5688285131598897`*^9, 3.5688285201487017`*^9}, 3.5688285866360188`*^9, {3.5688292397623663`*^9, 3.568829246688778*^9}, { 3.5688780328231993`*^9, 3.5688781532398105`*^9}, 3.568895317251341*^9}], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"Plot8", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"Graphics3D", "[", RowBox[{"{", RowBox[{ RowBox[{"Opacity", "[", "0.5", "]"}], ",", RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"wxh", "[", RowBox[{"[", "i", "]"}], "]"}], "\[NotEqual]", "0"}], ",", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"wxh", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"wyh", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"wzh", "[", RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "3"}], "]"}], "]"}]}], "]"}]}], "]"}]}], "}"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "nb"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]"}]], "Input", CellChangeTimes->{{3.5681365021632137`*^9, 3.5681365171392403`*^9}, { 3.568136560351316*^9, 3.5681365863097615`*^9}, {3.5681369135895367`*^9, 3.5681369265999594`*^9}, {3.568137821916732*^9, 3.5681378378443604`*^9}, { 3.5682072492149415`*^9, 3.5682072529901485`*^9}, 3.568209802830227*^9, { 3.5682216503836985`*^9, 3.5682217583046885`*^9}, {3.5682217912831464`*^9, 3.568221843465238*^9}, {3.5682219760654707`*^9, 3.5682221133457117`*^9}, { 3.568222159615393*^9, 3.5682221693498106`*^9}, {3.5682222565227633`*^9, 3.5682222578487663`*^9}, 3.5682224736283445`*^9, 3.5682227227607822`*^9, 3.5682227558640404`*^9, {3.568222818810151*^9, 3.5682228486062036`*^9}, { 3.5682230057764797`*^9, 3.5682230079916835`*^9}, {3.568223041952943*^9, 3.5682230548229656`*^9}, 3.5682232619601297`*^9, 3.5682233000709963`*^9, { 3.568223536645412*^9, 3.5682235381118145`*^9}, {3.568223591151908*^9, 3.568223591526308*^9}, {3.5682236491528096`*^9, 3.568223657483224*^9}, { 3.568223713222122*^9, 3.5682237532361927`*^9}, {3.568223800379475*^9, 3.568223827648323*^9}, {3.5682238887380304`*^9, 3.5682238893152313`*^9}, { 3.5682240287170763`*^9, 3.5682240612431335`*^9}, 3.5682241088076167`*^9, { 3.568224198772975*^9, 3.5682242396450467`*^9}, {3.568224282170721*^9, 3.5682243185031853`*^9}, {3.568224437890195*^9, 3.568224439590598*^9}, { 3.568224572253231*^9, 3.5682245867144566`*^9}, {3.5682247081762695`*^9, 3.568224733557514*^9}, {3.56822480520844*^9, 3.5682248433505073`*^9}, { 3.568224884206979*^9, 3.5682248959069996`*^9}, {3.568224982830352*^9, 3.5682250047639904`*^9}, {3.56822503870965*^9, 3.5682250803305235`*^9}, { 3.5682251223413973`*^9, 3.568225141451431*^9}, {3.568225234177994*^9, 3.5682252929588966`*^9}, {3.5682253775890455`*^9, 3.568225482140429*^9}, { 3.56822551694409*^9, 3.5682255297205124`*^9}, {3.568225573197789*^9, 3.5682255809666023`*^9}, 3.568225660698343*^9, 3.5682257301964645`*^9, { 3.568225765078126*^9, 3.5682257828933573`*^9}, {3.568225927026843*^9, 3.568225927416844*^9}, {3.568267062420238*^9, 3.568267081358671*^9}, { 3.568267188328059*^9, 3.5682672830358257`*^9}, {3.5682676367040462`*^9, 3.568267747979042*^9}, 3.5682678404716043`*^9, 3.568267882030078*^9, { 3.5682679252577534`*^9, 3.568267996612279*^9}, {3.5682680488879704`*^9, 3.568268061133992*^9}, {3.568268100914062*^9, 3.5682681150944867`*^9}, { 3.5682681672609787`*^9, 3.5682681709113846`*^9}, {3.568268272904364*^9, 3.5682683182536435`*^9}, {3.5682685202739983`*^9, 3.568268599194537*^9}, { 3.5682686348873997`*^9, 3.568268637773405*^9}, {3.5682686944951043`*^9, 3.568268872273017*^9}, {3.5682689420207396`*^9, 3.568269016947671*^9}, { 3.5682695040741262`*^9, 3.568269601153097*^9}, {3.5682696652380095`*^9, 3.5682696728352227`*^9}, {3.5682697489633565`*^9, 3.5682698308635006`*^9}, {3.568271156491429*^9, 3.568271160063835*^9}, { 3.568271244927984*^9, 3.568271369010602*^9}, {3.5682716250070515`*^9, 3.5682717312276382`*^9}, {3.5682717843301315`*^9, 3.5682718576814604`*^9}, {3.568271940673606*^9, 3.5682720536490045`*^9}, { 3.5682722259669075`*^9, 3.5682723223126764`*^9}, {3.568272370813162*^9, 3.568272411466833*^9}, {3.5682725893539457`*^9, 3.568272622238803*^9}, { 3.5682726936713285`*^9, 3.5682726960113325`*^9}, {3.568272741079812*^9, 3.56827275143823*^9}, {3.5682729218373294`*^9, 3.568272991070251*^9}, { 3.5682734285574193`*^9, 3.5682734302890224`*^9}, {3.5682734898811274`*^9, 3.5682735111595645`*^9}, {3.5682736115925407`*^9, 3.5682736313577757`*^9}, {3.5682738601945777`*^9, 3.568273987272401*^9}, { 3.5682763944346285`*^9, 3.5682764313910933`*^9}, 3.5682765128232365`*^9, { 3.5682766527398825`*^9, 3.568276660976697*^9}, {3.5682774496048822`*^9, 3.568277461663703*^9}, {3.5683784032731285`*^9, 3.568378422991563*^9}, { 3.568378459308427*^9, 3.5683784812420654`*^9}, {3.568378513440522*^9, 3.568378867046343*^9}, {3.568378915468828*^9, 3.56837896203491*^9}, { 3.5683790764299107`*^9, 3.568379270041851*^9}, {3.568380233530743*^9, 3.568380234388745*^9}, 3.5683804214954734`*^9, 3.568380471867962*^9, 3.5683805399464817`*^9, {3.5683812396077104`*^9, 3.568381293537005*^9}, { 3.568381339541486*^9, 3.5683813935175805`*^9}, 3.5684646396646404`*^9, 3.5684646988823442`*^9, 3.568465086387025*^9, {3.5684655859935026`*^9, 3.568465605087936*^9}, 3.568821235038273*^9, 3.5688233741452303`*^9, 3.568823449742963*^9, 3.5688274593312387`*^9, {3.5688284679978104`*^9, 3.5688284689026117`*^9}, {3.5688285131598897`*^9, 3.5688285201487017`*^9}, 3.5688285866360188`*^9, {3.5688292397623663`*^9, 3.568829246688778*^9}, { 3.5688780328231993`*^9, 3.568878147967001*^9}, {3.56889532230575*^9, 3.5688953242557535`*^9}, 3.56889600529175*^9, 3.5688961595292206`*^9}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"Plot9", "=", " ", RowBox[{"Table", "[", RowBox[{ RowBox[{"Graphics3D", "[", RowBox[{"{", RowBox[{ RowBox[{"Darker", "[", RowBox[{"Red", ",", "0.8"}], "]"}], ",", RowBox[{"Cylinder", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"x", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], ",", RowBox[{"y", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], ",", RowBox[{"z", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"x", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}], ",", RowBox[{"y", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}], ",", RowBox[{"z", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], "}"}]}], "}"}], ",", "0.05"}], "]"}]}], "}"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "nb"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"Show", "[", RowBox[{"Plot1", ",", "Plot7", ",", "Plot8", ",", "Plot9", ",", RowBox[{"{", RowBox[{ RowBox[{"AspectRatio", "\[Rule]", "Automatic"}], ",", RowBox[{"Boxed", "->", "False"}]}], "}"}]}], "]"}], "\[IndentingNewLine]", RowBox[{"Show", "[", RowBox[{"Plot1", ",", "Plot2", ",", "Plot9", ",", RowBox[{"{", RowBox[{"AspectRatio", "\[Rule]", "Automatic"}], "}"}]}], "]"}]}], "Input", CellChangeTimes->{{3.5681365021632137`*^9, 3.5681365171392403`*^9}, { 3.568136560351316*^9, 3.5681365863097615`*^9}, {3.5681369135895367`*^9, 3.5681369265999594`*^9}, {3.568137821916732*^9, 3.5681378378443604`*^9}, { 3.5682072492149415`*^9, 3.5682072529901485`*^9}, 3.568209802830227*^9, { 3.5682216503836985`*^9, 3.5682217583046885`*^9}, {3.5682217912831464`*^9, 3.568221843465238*^9}, {3.5682219760654707`*^9, 3.5682221133457117`*^9}, { 3.568222159615393*^9, 3.5682221693498106`*^9}, {3.5682222565227633`*^9, 3.5682222578487663`*^9}, 3.5682224736283445`*^9, 3.5682227227607822`*^9, 3.5682227558640404`*^9, {3.568222818810151*^9, 3.5682228486062036`*^9}, { 3.5682230057764797`*^9, 3.5682230079916835`*^9}, {3.568223041952943*^9, 3.5682230548229656`*^9}, 3.5682232619601297`*^9, 3.5682233000709963`*^9, { 3.568223536645412*^9, 3.5682235381118145`*^9}, {3.568223591151908*^9, 3.568223591526308*^9}, {3.5682236491528096`*^9, 3.568223657483224*^9}, { 3.568223713222122*^9, 3.5682237532361927`*^9}, {3.568223800379475*^9, 3.568223827648323*^9}, {3.5682238887380304`*^9, 3.5682238893152313`*^9}, { 3.5682240287170763`*^9, 3.5682240612431335`*^9}, 3.5682241088076167`*^9, { 3.568224198772975*^9, 3.5682242396450467`*^9}, {3.568224282170721*^9, 3.5682243185031853`*^9}, {3.568224437890195*^9, 3.568224439590598*^9}, { 3.568224572253231*^9, 3.5682245867144566`*^9}, {3.5682247081762695`*^9, 3.568224733557514*^9}, {3.56822480520844*^9, 3.5682248433505073`*^9}, { 3.568224884206979*^9, 3.5682248959069996`*^9}, {3.568224982830352*^9, 3.5682250047639904`*^9}, {3.56822503870965*^9, 3.5682250803305235`*^9}, { 3.5682251223413973`*^9, 3.568225141451431*^9}, {3.568225234177994*^9, 3.5682252929588966`*^9}, {3.5682253775890455`*^9, 3.568225482140429*^9}, { 3.56822551694409*^9, 3.5682255297205124`*^9}, {3.568225573197789*^9, 3.5682255809666023`*^9}, 3.568225660698343*^9, 3.5682257301964645`*^9, { 3.568225765078126*^9, 3.5682257828933573`*^9}, {3.568225927026843*^9, 3.568225927416844*^9}, {3.568267062420238*^9, 3.568267081358671*^9}, { 3.568267188328059*^9, 3.5682672830358257`*^9}, {3.5682676367040462`*^9, 3.568267747979042*^9}, 3.5682678404716043`*^9, 3.568267882030078*^9, { 3.5682679252577534`*^9, 3.568267996612279*^9}, {3.5682680488879704`*^9, 3.568268061133992*^9}, {3.568268100914062*^9, 3.5682681150944867`*^9}, { 3.5682681672609787`*^9, 3.5682681709113846`*^9}, {3.568268272904364*^9, 3.5682683182536435`*^9}, {3.5682685202739983`*^9, 3.568268599194537*^9}, { 3.5682686348873997`*^9, 3.568268637773405*^9}, {3.5682686944951043`*^9, 3.568268872273017*^9}, {3.5682689420207396`*^9, 3.568269016947671*^9}, { 3.5682695040741262`*^9, 3.568269601153097*^9}, {3.5682696652380095`*^9, 3.5682696728352227`*^9}, {3.5682697489633565`*^9, 3.5682698308635006`*^9}, {3.568271156491429*^9, 3.568271160063835*^9}, { 3.568271244927984*^9, 3.568271369010602*^9}, {3.5682716250070515`*^9, 3.5682717312276382`*^9}, {3.5682717843301315`*^9, 3.5682718576814604`*^9}, {3.568271940673606*^9, 3.5682720536490045`*^9}, { 3.5682722259669075`*^9, 3.5682723223126764`*^9}, {3.568272370813162*^9, 3.568272411466833*^9}, {3.5682725893539457`*^9, 3.568272622238803*^9}, { 3.5682726936713285`*^9, 3.5682726960113325`*^9}, {3.568272741079812*^9, 3.56827275143823*^9}, {3.5682729218373294`*^9, 3.568272991070251*^9}, { 3.5682734285574193`*^9, 3.5682734302890224`*^9}, {3.5682734898811274`*^9, 3.5682735111595645`*^9}, {3.5682736115925407`*^9, 3.5682736313577757`*^9}, {3.5682738601945777`*^9, 3.568273987272401*^9}, { 3.5682763944346285`*^9, 3.5682764313910933`*^9}, 3.5682765128232365`*^9, { 3.5682766527398825`*^9, 3.568276660976697*^9}, {3.5682774496048822`*^9, 3.568277461663703*^9}, {3.5683784032731285`*^9, 3.568378422991563*^9}, { 3.568378459308427*^9, 3.5683784812420654`*^9}, {3.568378513440522*^9, 3.568378867046343*^9}, {3.568378915468828*^9, 3.56837896203491*^9}, { 3.5683790764299107`*^9, 3.568379270041851*^9}, {3.568380233530743*^9, 3.568380234388745*^9}, 3.5683804214954734`*^9, 3.568380471867962*^9, 3.5683805399464817`*^9, {3.5683812396077104`*^9, 3.568381293537005*^9}, { 3.568381339541486*^9, 3.568381420224828*^9}, {3.5684645223056345`*^9, 3.5684645572340956`*^9}, {3.5684651558539467`*^9, 3.5684651618911576`*^9}, {3.5688209923798466`*^9, 3.568820997683856*^9}, { 3.568821158754139*^9, 3.568821160750942*^9}, 3.56882222404881*^9, { 3.568822649274357*^9, 3.5688226508031597`*^9}, 3.568822723109287*^9}], Cell[BoxData[ Graphics3DBox[{ {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], SphereBox[{ 0.00021354623828926078`, -0.00011922259401907318`, \ -0.00004316425121513322}, 0.255]}, {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], SphereBox[{-0.12934170283061472`, 2.0235255842401028`, 0.3646100108646693}, 0.1]}, {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], SphereBox[{-0.9885757516053951, -0.4479637748223944, \ -1.7512287238963824`}, 0.1]}, {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], SphereBox[{-0.8630092084766477, -1.0344278187894953`, 1.5588388552416423`}, 0.1]}, {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], SphereBox[{ 1.9807093422551116`, -0.5410118171750784, -0.17217775172391747`}, 0.1]}, {RGBColor[ 0.33333333333333337`, 0.33333333333333337`, 0.33333333333333337`], {Opacity[0.5]}}, {RGBColor[ 0.33333333333333337`, 0.33333333333333337`, 0.33333333333333337`], {Opacity[0.5]}}, {RGBColor[ 0.33333333333333337`, 0.33333333333333337`, 0.33333333333333337`], {Opacity[0.5]}}, {RGBColor[ 0.33333333333333337`, 0.33333333333333337`, 0.33333333333333337`], {Opacity[0.5]}}, {Opacity[0.5], SphereBox[{-0.0916129197206941, 1.4342043752788571`, 0.25841654781961776`}, 1.210314249067387]}, {Opacity[0.5], SphereBox[{-0.7006236459891667, -0.31754389458689425`, \ -1.2412539797074782`}, 1.2103203349966232`]}, {Opacity[0.5], SphereBox[{-0.6116283729654048, -0.7332247310952404, 1.1048736853209575`}, 1.2103545683012018`]}, {Opacity[0.5], SphereBox[{ 1.4039310556503042`, -0.383488077631103, -0.12204712654025901`}, 1.2102326842173277`]}, {RGBColor[0.19999999999999996`, 0., 0.], CylinderBox[{{-0.12934170283061472`, 2.0235255842401028`, 0.3646100108646693}, { 0.00021354623828926078`, -0.00011922259401907318`, \ -0.00004316425121513322}}, 0.05]}, {RGBColor[0.19999999999999996`, 0., 0.], CylinderBox[{{-0.9885757516053951, -0.4479637748223944, \ -1.7512287238963824`}, { 0.00021354623828926078`, -0.00011922259401907318`, \ -0.00004316425121513322}}, 0.05]}, {RGBColor[0.19999999999999996`, 0., 0.], CylinderBox[{{-0.8630092084766477, -1.0344278187894953`, 1.5588388552416423`}, { 0.00021354623828926078`, -0.00011922259401907318`, \ -0.00004316425121513322}}, 0.05]}, {RGBColor[0.19999999999999996`, 0., 0.], CylinderBox[{{ 1.9807093422551116`, -0.5410118171750784, -0.17217775172391747`}, { 0.00021354623828926078`, -0.00011922259401907318`, \ -0.00004316425121513322}}, 0.05]}}, AspectRatio->Automatic, AutomaticImageSize->True, Boxed->False, ImageSize->{381.4681669931018, 404.1912114674104}, ViewPoint->{-2.7589253365402997`, -0.6679367093606237, 1.8417903082828118`}, ViewVertical->{-0.17062622309585398`, 0.2931727128124413, 0.9455719376914227}]], "Output", CellChangeTimes->{ 3.5689997298853893`*^9, 3.5690600928196216`*^9, 3.569065230515286*^9, 3.5690652770657673`*^9, {3.5690653772959437`*^9, 3.569065406483595*^9}, 3.569065562218669*^9, 3.5690656440252123`*^9, 3.5690658378399525`*^9, 3.569066362359674*^9, 3.5690664048229485`*^9, {3.569066534537176*^9, 3.5690665558156137`*^9}, {3.5690666165777206`*^9, 3.5690666344865522`*^9}, {3.569066676201025*^9, 3.569066705716277*^9}, { 3.5690667508627567`*^9, 3.5690667709087915`*^9}, {3.569066820657279*^9, 3.569066847489326*^9}, 3.5690668808889847`*^9, 3.5690694998559575`*^9, 3.5690695664992743`*^9, 3.569069721922347*^9, 3.5690698100313025`*^9, 3.5690698553649817`*^9, 3.569069945813941*^9, 3.5690701759767447`*^9, 3.5690702507476764`*^9, 3.5690704405376096`*^9, 3.569070765688981*^9, 3.5690715254031153`*^9, 3.5690715828112164`*^9}], Cell[BoxData[ Graphics3DBox[{ {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], SphereBox[{ 0.00021354623828926078`, -0.00011922259401907318`, \ -0.00004316425121513322}, 0.255]}, {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], SphereBox[{-0.12934170283061472`, 2.0235255842401028`, 0.3646100108646693}, 0.1]}, {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], SphereBox[{-0.9885757516053951, -0.4479637748223944, \ -1.7512287238963824`}, 0.1]}, {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], SphereBox[{-0.8630092084766477, -1.0344278187894953`, 1.5588388552416423`}, 0.1]}, {RGBColor[0.7333333333333333, 0.6, 0.4666666666666667], SphereBox[{ 1.9807093422551116`, -0.5410118171750784, -0.17217775172391747`}, 0.1]}, Text3DBox[ "1", {-0.2997864537617107, -0.00011922259401907318`, \ -0.00004316425121513322}, {1, 0}], Text3DBox["2", {-0.4293417028306147, 2.0235255842401028`, 0.3646100108646693}, {1, 0}], Text3DBox[ "3", {-1.288575751605395, -0.4479637748223944, -1.7512287238963824`}, {1, 0}], Text3DBox["4", {-1.1630092084766477`, -1.0344278187894953`, 1.5588388552416423`}, {1, 0}], Text3DBox["5", { 1.6807093422551116`, -0.5410118171750784, -0.17217775172391747`}, {1, 0}], {RGBColor[0.19999999999999996`, 0., 0.], CylinderBox[{{-0.12934170283061472`, 2.0235255842401028`, 0.3646100108646693}, { 0.00021354623828926078`, -0.00011922259401907318`, \ -0.00004316425121513322}}, 0.05]}, {RGBColor[0.19999999999999996`, 0., 0.], CylinderBox[{{-0.9885757516053951, -0.4479637748223944, \ -1.7512287238963824`}, { 0.00021354623828926078`, -0.00011922259401907318`, \ -0.00004316425121513322}}, 0.05]}, {RGBColor[0.19999999999999996`, 0., 0.], CylinderBox[{{-0.8630092084766477, -1.0344278187894953`, 1.5588388552416423`}, { 0.00021354623828926078`, -0.00011922259401907318`, \ -0.00004316425121513322}}, 0.05]}, {RGBColor[0.19999999999999996`, 0., 0.], CylinderBox[{{ 1.9807093422551116`, -0.5410118171750784, -0.17217775172391747`}, { 0.00021354623828926078`, -0.00011922259401907318`, \ -0.00004316425121513322}}, 0.05]}}, AspectRatio->Automatic, AutomaticImageSize->True, ImageSize->{381.4681669931018, 404.1912114674104}, ViewPoint->{-2.7589253365402997`, -0.6679367093606237, 1.8417903082828118`}, ViewVertical->{-0.17062622309585398`, 0.2931727128124413, 0.9455719376914227}]], "Output", CellChangeTimes->{ 3.5689997298853893`*^9, 3.5690600928196216`*^9, 3.569065230515286*^9, 3.5690652770657673`*^9, {3.5690653772959437`*^9, 3.569065406483595*^9}, 3.569065562218669*^9, 3.5690656440252123`*^9, 3.5690658378399525`*^9, 3.569066362359674*^9, 3.5690664048229485`*^9, {3.569066534537176*^9, 3.5690665558156137`*^9}, {3.5690666165777206`*^9, 3.5690666344865522`*^9}, {3.569066676201025*^9, 3.569066705716277*^9}, { 3.5690667508627567`*^9, 3.5690667709087915`*^9}, {3.569066820657279*^9, 3.569066847489326*^9}, 3.5690668808889847`*^9, 3.5690694998559575`*^9, 3.5690695664992743`*^9, 3.569069721922347*^9, 3.5690698100313025`*^9, 3.5690698553649817`*^9, 3.569069945813941*^9, 3.5690701759767447`*^9, 3.5690702507476764`*^9, 3.5690704405376096`*^9, 3.569070765688981*^9, 3.5690715254031153`*^9, 3.5690715828112164`*^9}] }, Open ]], Cell[BoxData[{ RowBox[{ RowBox[{"twc", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"wxc", "[", RowBox[{"[", "i", "]"}], "]"}], "\[NotEqual]", "0"}], ",", RowBox[{"{", RowBox[{ RowBox[{"wxc", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"wyc", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"wzc", "[", RowBox[{"[", "i", "]"}], "]"}]}], "}"}]}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "cw"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"twh", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"If", "[", RowBox[{ RowBox[{ RowBox[{"wxh", "[", RowBox[{"[", RowBox[{"i", "+", "cw"}], "]"}], "]"}], "\[NotEqual]", "0"}], ",", RowBox[{"{", RowBox[{ RowBox[{"wxh", "[", RowBox[{"[", RowBox[{"i", "+", "cw"}], "]"}], "]"}], ",", RowBox[{"wyh", "[", RowBox[{"[", RowBox[{"i", "+", "cw"}], "]"}], "]"}], ",", RowBox[{"wzh", "[", RowBox[{"[", RowBox[{"i", "+", "cw"}], "]"}], "]"}]}], "}"}]}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "bh"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"tw", "=", RowBox[{"Join", "[", RowBox[{"twc", ",", "twh"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"MatrixForm", "[", "tw", "]"}], "*)"}]}], "\[IndentingNewLine]"}], "Input", CellChangeTimes->{{3.568895397419882*^9, 3.5688955157148895`*^9}, { 3.568895549629349*^9, 3.5688955855718126`*^9}, {3.568895617317868*^9, 3.568895845234268*^9}, {3.5688960797194805`*^9, 3.5688961015127187`*^9}, { 3.568896219074525*^9, 3.5688962276233397`*^9}, {3.568896262910602*^9, 3.5688962903198504`*^9}, 3.5688964250728865`*^9, 3.568896456912543*^9, { 3.5688981540863237`*^9, 3.5688981869399815`*^9}, {3.5689021233572955`*^9, 3.5689021341525145`*^9}, {3.5689111169043026`*^9, 3.5689111218495116`*^9}, {3.568911235261711*^9, 3.568911242437723*^9}, 3.5689112948850155`*^9, {3.5689116510180407`*^9, 3.568911800419503*^9}, { 3.568911871820829*^9, 3.5689119012892804`*^9}, {3.5689119993042526`*^9, 3.5689120006770554`*^9}, {3.568912144306507*^9, 3.568912165085744*^9}, { 3.5689528310632005`*^9, 3.5689528569592457`*^9}, {3.5689534065794115`*^9, 3.5689534297298517`*^9}, {3.568954496214126*^9, 3.5689545054337425`*^9}, 3.5689569073699613`*^9}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"Vee", "=", "0.0"}], ";", " ", RowBox[{"Vne", "=", "0.0"}], ";"}], "\[IndentingNewLine]", RowBox[{"bw", "=", RowBox[{"Length", "[", "tw", "]"}]}]}], "Input", CellChangeTimes->{{3.5688982251756487`*^9, 3.5688986440675845`*^9}, { 3.5688986873108606`*^9, 3.56889871560931*^9}, {3.5688987970414534`*^9, 3.5688988526711507`*^9}, {3.568898887147211*^9, 3.568898888644814*^9}, { 3.5688989206248703`*^9, 3.568899041759083*^9}, {3.5688990876855636`*^9, 3.568899258349863*^9}, 3.568899300267137*^9, {3.568899351326027*^9, 3.5688993831032825`*^9}, {3.5688994370481772`*^9, 3.568899468217032*^9}, { 3.5688995813016305`*^9, 3.5688997165538683`*^9}, {3.5688998495129013`*^9, 3.568899932364647*^9}, {3.568899967667509*^9, 3.568900024529609*^9}, { 3.5689000794261055`*^9, 3.568900136257005*^9}, {3.5689001767858763`*^9, 3.56890020149632*^9}, {3.568902068726*^9, 3.568902074139209*^9}, 3.56890216040736*^9, 3.568902344113283*^9, {3.5689024237982235`*^9, 3.56890243907065*^9}, 3.5689025770684924`*^9, 3.568910195129884*^9, 3.568910620995032*^9, {3.568912393532545*^9, 3.5689124022841606`*^9}}], Cell[BoxData["4"], "Output", CellChangeTimes->{ 3.5689997299165897`*^9, 3.5690600928508215`*^9, 3.569065230546486*^9, 3.5690652770969677`*^9, {3.5690653773115435`*^9, 3.569065406561595*^9}, 3.5690655622498684`*^9, 3.5690656440564127`*^9, 3.569065837871153*^9, 3.5690663624220743`*^9, 3.5690664048697486`*^9, {3.5690665345995765`*^9, 3.5690665558780136`*^9}, {3.5690666166245203`*^9, 3.569066634548952*^9}, { 3.569066676263425*^9, 3.569066705778677*^9}, {3.5690667509095564`*^9, 3.5690667709711914`*^9}, {3.569066820719679*^9, 3.5690668475517263`*^9}, 3.5690668809513845`*^9, 3.5690694999495573`*^9, 3.5690695665616746`*^9, 3.569069721984748*^9, 3.5690698100937023`*^9, 3.569069855427382*^9, 3.569069945876341*^9, 3.569070176039145*^9, 3.5690702508100767`*^9, 3.56907044060001*^9, 3.5690707657357807`*^9, 3.5690715254343157`*^9, 3.5690715828268166`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "<", "bw"}], ",", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"j", "=", RowBox[{"i", "+", "1"}]}], ",", RowBox[{"j", "<", RowBox[{"bw", "+", "1"}]}], ",", RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", RowBox[{"Vee", "=", RowBox[{"Vee", "+", RowBox[{"4", "/", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"tw", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "-", RowBox[{"tw", "[", RowBox[{"[", RowBox[{"j", ",", "1"}], "]"}], "]"}]}], ")"}], "^", "2"}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"tw", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "-", RowBox[{"tw", "[", RowBox[{"[", RowBox[{"j", ",", "2"}], "]"}], "]"}]}], ")"}], "^", "2"}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"tw", "[", RowBox[{"[", RowBox[{"i", ",", "3"}], "]"}], "]"}], "-", RowBox[{"tw", "[", RowBox[{"[", RowBox[{"j", ",", "3"}], "]"}], "]"}]}], ")"}], "^", "2"}]}], "]"}]}]}]}]}], "]"}]}], "]"}], "\[IndentingNewLine]", "Vee"}], "Input", CellChangeTimes->{{3.5688982251756487`*^9, 3.5688986440675845`*^9}, { 3.5688986873108606`*^9, 3.56889871560931*^9}, {3.5688987970414534`*^9, 3.5688988526711507`*^9}, {3.568898887147211*^9, 3.568898888644814*^9}, { 3.5688989206248703`*^9, 3.568899041759083*^9}, {3.5688990876855636`*^9, 3.568899258349863*^9}, 3.568899300267137*^9, {3.568899351326027*^9, 3.5688993831032825`*^9}, {3.5688994370481772`*^9, 3.568899468217032*^9}, { 3.5688995813016305`*^9, 3.5688997165538683`*^9}, {3.5688998495129013`*^9, 3.568899932364647*^9}, {3.568899967667509*^9, 3.568900024529609*^9}, { 3.5689000794261055`*^9, 3.568900136257005*^9}, {3.5689001767858763`*^9, 3.56890020149632*^9}, {3.568902068726*^9, 3.568902074139209*^9}, 3.56890216040736*^9, 3.568902344113283*^9, {3.5689024237982235`*^9, 3.56890243907065*^9}, 3.5689025770684924`*^9, {3.568910195129884*^9, 3.568910199825492*^9}, 3.568910617609826*^9, {3.56891263539537*^9, 3.568912638842976*^9}, {3.568952611212014*^9, 3.5689526227248344`*^9}, { 3.5689527476498537`*^9, 3.5689527480086546`*^9}, {3.5689530637688093`*^9, 3.5689530642368097`*^9}}], Cell[BoxData["10.064291151956446`"], "Output", CellChangeTimes->{ 3.5689997299165897`*^9, 3.5690600928664217`*^9, 3.569065230562086*^9, 3.569065277112568*^9, {3.5690653773271437`*^9, 3.569065406561595*^9}, 3.5690655622654686`*^9, 3.5690656440720124`*^9, 3.569065837886753*^9, 3.569066362437674*^9, 3.5690664048853483`*^9, {3.5690665345995765`*^9, 3.569066555893614*^9}, {3.5690666166401205`*^9, 3.5690666345645523`*^9}, { 3.569066676279025*^9, 3.569066705794277*^9}, {3.5690667509251566`*^9, 3.5690667709867916`*^9}, {3.569066820735279*^9, 3.5690668475673265`*^9}, 3.5690668809669847`*^9, 3.5690694999651575`*^9, 3.5690695665772743`*^9, 3.5690697220003476`*^9, 3.5690698101093025`*^9, 3.5690698554429817`*^9, 3.569069945891941*^9, 3.569070176039145*^9, 3.5690702508256764`*^9, 3.56907044060001*^9, 3.569070765751381*^9, 3.569071525449916*^9, 3.5690715828424163`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "<", "bc"}], ",", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"j", "=", RowBox[{"i", "+", "1"}]}], ",", RowBox[{"j", "<", RowBox[{"bc", "+", "1"}]}], ",", RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", RowBox[{"Vee", "=", RowBox[{"Vee", "+", RowBox[{"4", "/", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "-", RowBox[{"t", "[", RowBox[{"[", RowBox[{"j", ",", "1"}], "]"}], "]"}]}], ")"}], "^", "2"}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "-", RowBox[{"t", "[", RowBox[{"[", RowBox[{"j", ",", "2"}], "]"}], "]"}]}], ")"}], "^", "2"}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{"i", ",", "3"}], "]"}], "]"}], "-", RowBox[{"t", "[", RowBox[{"[", RowBox[{"j", ",", "3"}], "]"}], "]"}]}], ")"}], "^", "2"}]}], "]"}]}]}]}]}], "]"}]}], "]"}], "\[IndentingNewLine]", "Vee"}], "Input", CellChangeTimes->{{3.5688982251756487`*^9, 3.5688986440675845`*^9}, { 3.5688986873108606`*^9, 3.56889871560931*^9}, {3.5688987970414534`*^9, 3.5688988526711507`*^9}, {3.568898887147211*^9, 3.568898888644814*^9}, { 3.5688989206248703`*^9, 3.568899041759083*^9}, {3.5688990876855636`*^9, 3.568899258349863*^9}, 3.568899300267137*^9, {3.568899351326027*^9, 3.5688993831032825`*^9}, {3.5688994370481772`*^9, 3.568899468217032*^9}, { 3.5688995813016305`*^9, 3.5688997165538683`*^9}, {3.5688998495129013`*^9, 3.568899932364647*^9}, {3.568899967667509*^9, 3.568900024529609*^9}, { 3.5689000794261055`*^9, 3.568900136257005*^9}, {3.5689001767858763`*^9, 3.56890020149632*^9}, {3.568902068726*^9, 3.568902074139209*^9}, 3.56890216040736*^9, 3.568902344113283*^9, {3.5689024237982235`*^9, 3.56890243907065*^9}, 3.5689025770684924`*^9, {3.568910195129884*^9, 3.5689102048487005`*^9}, 3.5689106144742203`*^9, {3.568912644739786*^9, 3.56891264654939*^9}, {3.568952467364841*^9, 3.5689524768652573`*^9}, { 3.5689525940994635`*^9, 3.568952604472802*^9}, {3.568954607114721*^9, 3.5689546170051384`*^9}}], Cell[BoxData["10.064291151956446`"], "Output", CellChangeTimes->{ 3.56899972993219*^9, 3.569060092882022*^9, 3.5690652305776863`*^9, 3.569065277128168*^9, {3.569065377342744*^9, 3.569065406577195*^9}, 3.569065562281069*^9, 3.569065644087612*^9, 3.569065837902353*^9, 3.569066362453274*^9, 3.5690664049009485`*^9, {3.569066534615176*^9, 3.569066555909214*^9}, {3.5690666166557207`*^9, 3.569066634580152*^9}, { 3.5690666762946253`*^9, 3.569066705794277*^9}, {3.5690667509407563`*^9, 3.5690667709867916`*^9}, {3.5690668207508793`*^9, 3.5690668475829263`*^9}, 3.569066880982585*^9, 3.5690694999807577`*^9, 3.5690695665928745`*^9, 3.5690697220159473`*^9, 3.5690698101249027`*^9, 3.569069855458582*^9, 3.569069945907541*^9, 3.569070176054745*^9, 3.5690702508256764`*^9, 3.56907044061561*^9, 3.569070765751381*^9, 3.5690715254655156`*^9, 3.5690715828580165`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "<", RowBox[{"bw", "+", "1"}]}], ",", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"j", "=", "1"}], ",", RowBox[{"j", "<", RowBox[{"bc", "+", "1"}]}], ",", RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", RowBox[{"Vee", "=", RowBox[{"Vee", "+", RowBox[{"4", "/", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"tw", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "-", RowBox[{"t", "[", RowBox[{"[", RowBox[{"j", ",", "1"}], "]"}], "]"}]}], ")"}], "^", "2"}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"tw", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "-", RowBox[{"t", "[", RowBox[{"[", RowBox[{"j", ",", "2"}], "]"}], "]"}]}], ")"}], "^", "2"}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"tw", "[", RowBox[{"[", RowBox[{"i", ",", "3"}], "]"}], "]"}], "-", RowBox[{"t", "[", RowBox[{"[", RowBox[{"j", ",", "3"}], "]"}], "]"}]}], ")"}], "^", "2"}]}], "]"}]}]}]}]}], "]"}]}], "]"}], "\[IndentingNewLine]", "Vee"}], "Input", CellChangeTimes->{{3.5688982251756487`*^9, 3.5688986440675845`*^9}, { 3.5688986873108606`*^9, 3.56889871560931*^9}, {3.5688987970414534`*^9, 3.5688988526711507`*^9}, {3.568898887147211*^9, 3.568898888644814*^9}, { 3.5688989206248703`*^9, 3.568899041759083*^9}, {3.5688990876855636`*^9, 3.568899258349863*^9}, 3.568899300267137*^9, {3.568899351326027*^9, 3.5688993831032825`*^9}, {3.5688994370481772`*^9, 3.568899468217032*^9}, { 3.5688995813016305`*^9, 3.5688997165538683`*^9}, {3.5688998495129013`*^9, 3.568899932364647*^9}, {3.568899967667509*^9, 3.568900024529609*^9}, { 3.5689000794261055`*^9, 3.568900136257005*^9}, {3.5689001767858763`*^9, 3.56890020149632*^9}, {3.568902068726*^9, 3.568902074139209*^9}, 3.56890216040736*^9, 3.568902344113283*^9, {3.5689024237982235`*^9, 3.56890243907065*^9}, 3.5689025770684924`*^9, {3.568910195129884*^9, 3.5689102112915125`*^9}, {3.5689106097162123`*^9, 3.568910610371413*^9}, { 3.568912654084203*^9, 3.5689126558314056`*^9}, 3.568952500358899*^9, 3.568952653737689*^9, {3.568953090039255*^9, 3.5689530906476564`*^9}, { 3.568954639047977*^9, 3.5689546392507772`*^9}}], Cell[BoxData["21.020902944929748`"], "Output", CellChangeTimes->{ 3.5689997299477897`*^9, 3.5690600928976216`*^9, 3.569065230593286*^9, 3.5690652771437683`*^9, {3.569065377358344*^9, 3.569065406577195*^9}, 3.569065562296669*^9, 3.5690656441032124`*^9, 3.5690658379179525`*^9, 3.569066362453274*^9, 3.5690664049165487`*^9, {3.5690665346307764`*^9, 3.569066555909214*^9}, {3.569066616671321*^9, 3.569066634595752*^9}, { 3.5690666763102255`*^9, 3.5690667058098774`*^9}, {3.569066750956356*^9, 3.569066771002392*^9}, {3.5690668207508793`*^9, 3.569066847598526*^9}, 3.5690668809981847`*^9, 3.569069499996358*^9, 3.5690695666084747`*^9, 3.5690697220315475`*^9, 3.5690698101405025`*^9, 3.569069855474182*^9, 3.5690699459231405`*^9, 3.569070176070345*^9, 3.569070250841276*^9, 3.56907044063121*^9, 3.569070765766981*^9, 3.5690715254811153`*^9, 3.569071582873616*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"sc", "=", "0.4"}], ";"}], "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], " ", ",", RowBox[{"i", "<", RowBox[{"bc", "+", "1"}]}], ",", RowBox[{"i", "++"}], ",", RowBox[{"Vee", "=", RowBox[{"Vee", "+", RowBox[{"3", "*", RowBox[{"sc", "/", RowBox[{"u2", "[", RowBox[{"[", "i", "]"}], "]"}]}]}]}]}]}], "]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"sw", "=", "0.4"}], ";"}], "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], " ", ",", RowBox[{"i", "<", RowBox[{"nb", "+", "1"}]}], ",", RowBox[{"i", "++"}], ",", RowBox[{"Vee", "=", RowBox[{"Vee", "+", RowBox[{"3", "*", RowBox[{"sw", "/", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "3"}], "]"}], "]"}]}]}]}]}]}], "]"}], "\[IndentingNewLine]", "Vee", "\[IndentingNewLine]"}], "Input", CellChangeTimes->{{3.5689548890696163`*^9, 3.568955116299615*^9}, { 3.5689552553894596`*^9, 3.5689552575266633`*^9}, {3.5689552936875267`*^9, 3.5689553203323736`*^9}, {3.568955404229321*^9, 3.568955433245372*^9}, 3.5689554882822685`*^9, 3.5689562157115464`*^9, {3.568956283961666*^9, 3.5689562858180695`*^9}, {3.5689566111878405`*^9, 3.5689566291278725`*^9}, {3.5689566695319433`*^9, 3.5689567220104356`*^9}, {3.5690653987303815`*^9, 3.569065401881587*^9}, { 3.5690666985090647`*^9, 3.5690667010674686`*^9}}], Cell[BoxData["29.78684391486262`"], "Output", CellChangeTimes->{ 3.5689997299633894`*^9, 3.5690600929132214`*^9, 3.569065230608886*^9, 3.569065277159368*^9, {3.5690653773739443`*^9, 3.5690654065927954`*^9}, 3.5690655623122687`*^9, 3.5690656441188126`*^9, 3.569065837949153*^9, 3.569066362468874*^9, 3.569066404932149*^9, {3.5690665346463766`*^9, 3.569066555924814*^9}, {3.5690666166869206`*^9, 3.569066634626952*^9}, { 3.5690666763258257`*^9, 3.5690667058098774`*^9}, {3.5690667509719563`*^9, 3.5690667710179915`*^9}, {3.5690668207664795`*^9, 3.569066847614126*^9}, 3.569066881013785*^9, 3.569069499996358*^9, 3.569069566624075*^9, 3.5690697220471478`*^9, 3.569069810156102*^9, 3.5690698554897823`*^9, 3.5690699459387407`*^9, 3.569070176085945*^9, 3.5690702508568764`*^9, 3.5690704406468096`*^9, 3.569070765782581*^9, 3.5690715254967155`*^9, 3.5690715828892164`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "<", RowBox[{"bw", "+", "1"}]}], ",", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"j", "=", "1"}], ",", RowBox[{"j", "<", RowBox[{"n", "+", "1"}]}], ",", RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", RowBox[{"Vne", "=", RowBox[{"Vne", "-", RowBox[{"2", "*", RowBox[{ RowBox[{"cn", "[", RowBox[{"[", "j", "]"}], "]"}], "/", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"tw", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "-", RowBox[{"t", "[", RowBox[{"[", RowBox[{"j", ",", "1"}], "]"}], "]"}]}], ")"}], "^", "2"}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"tw", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "-", RowBox[{"t", "[", RowBox[{"[", RowBox[{"j", ",", "2"}], "]"}], "]"}]}], ")"}], "^", "2"}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"tw", "[", RowBox[{"[", RowBox[{"i", ",", "3"}], "]"}], "]"}], "-", RowBox[{"t", "[", RowBox[{"[", RowBox[{"j", ",", "3"}], "]"}], "]"}]}], ")"}], "^", "2"}]}], "]"}]}]}]}]}]}], "]"}]}], "]"}], "\[IndentingNewLine]", "Vne"}], "Input", CellChangeTimes->{{3.5688982251756487`*^9, 3.5688986440675845`*^9}, { 3.5688986873108606`*^9, 3.56889871560931*^9}, {3.5688987970414534`*^9, 3.5688988526711507`*^9}, {3.568898887147211*^9, 3.568898888644814*^9}, { 3.5688989206248703`*^9, 3.568899041759083*^9}, {3.5688990876855636`*^9, 3.568899258349863*^9}, 3.568899300267137*^9, {3.568899351326027*^9, 3.5688993831032825`*^9}, {3.5688994370481772`*^9, 3.568899468217032*^9}, { 3.5688995813016305`*^9, 3.5688997165538683`*^9}, {3.5688998495129013`*^9, 3.568899932364647*^9}, {3.568899967667509*^9, 3.568900024529609*^9}, { 3.5689000794261055`*^9, 3.568900136257005*^9}, {3.5689001767858763`*^9, 3.56890020149632*^9}, {3.568902068726*^9, 3.568902074139209*^9}, 3.56890216040736*^9, 3.568902344113283*^9, {3.5689024237982235`*^9, 3.56890243907065*^9}, 3.5689025770684924`*^9, {3.568910195129884*^9, 3.568910216923122*^9}, {3.568910624395838*^9, 3.5689106253786397`*^9}, { 3.568912441143829*^9, 3.568912441424629*^9}, 3.5689124745746875`*^9, { 3.5689126628670177`*^9, 3.568912665066622*^9}, {3.568953170722597*^9, 3.568953194637439*^9}}], Cell[BoxData[ RowBox[{"-", "54.49228165073883`"}]], "Output", CellChangeTimes->{ 3.5689997299789896`*^9, 3.5690600929288216`*^9, 3.569065230624486*^9, 3.569065277174968*^9, {3.569065377389544*^9, 3.5690654065927954`*^9}, 3.5690655623278685`*^9, 3.5690656441344123`*^9, 3.569065837964753*^9, 3.569066362484474*^9, 3.569066404947749*^9, {3.569066534661977*^9, 3.569066555940414*^9}, {3.5690666166869206`*^9, 3.569066634626952*^9}, { 3.5690666763258257`*^9, 3.569066705825477*^9}, {3.5690667509875565`*^9, 3.5690667710335913`*^9}, {3.569066820782079*^9, 3.5690668476297264`*^9}, 3.5690668810293846`*^9, 3.5690695000119576`*^9, 3.569069566639675*^9, 3.569069722062748*^9, 3.5690698101717024`*^9, 3.569069855505382*^9, 3.569069945954341*^9, 3.5690701761015453`*^9, 3.5690702508724766`*^9, 3.56907044066241*^9, 3.569070765782581*^9, 3.5690715255123158`*^9, 3.5690715829048166`*^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", "=", "i"}], ",", RowBox[{"j", "<", RowBox[{"n", "+", "1"}]}], ",", RowBox[{"j", "++"}], ",", RowBox[{"If", "[", RowBox[{ RowBox[{"i", "\[NotEqual]", "j"}], ",", "\[IndentingNewLine]", RowBox[{"Vne", "=", RowBox[{"Vne", "-", RowBox[{"2", "*", RowBox[{ RowBox[{"cn", "[", RowBox[{"[", "j", "]"}], "]"}], "/", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "-", RowBox[{"t", "[", RowBox[{"[", RowBox[{"j", ",", "1"}], "]"}], "]"}]}], ")"}], "^", "2"}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{"i", ",", "2"}], "]"}], "]"}], "-", RowBox[{"t", "[", RowBox[{"[", RowBox[{"j", ",", "2"}], "]"}], "]"}]}], ")"}], "^", "2"}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"t", "[", RowBox[{"[", RowBox[{"i", ",", "3"}], "]"}], "]"}], "-", RowBox[{"t", "[", RowBox[{"[", RowBox[{"j", ",", "3"}], "]"}], "]"}]}], ")"}], "^", "2"}]}], "]"}]}]}]}]}], ",", "0"}], "]"}]}], "]"}]}], "]"}], "\[IndentingNewLine]", "Vne"}], "Input", CellChangeTimes->{{3.5688982251756487`*^9, 3.5688986440675845`*^9}, { 3.5688986873108606`*^9, 3.56889871560931*^9}, {3.5688987970414534`*^9, 3.5688988526711507`*^9}, {3.568898887147211*^9, 3.568898888644814*^9}, { 3.5688989206248703`*^9, 3.568899041759083*^9}, {3.5688990876855636`*^9, 3.568899258349863*^9}, 3.568899300267137*^9, {3.568899351326027*^9, 3.5688993831032825`*^9}, {3.5688994370481772`*^9, 3.568899468217032*^9}, { 3.5688995813016305`*^9, 3.5688997165538683`*^9}, {3.5688998495129013`*^9, 3.568899932364647*^9}, {3.568899967667509*^9, 3.568900024529609*^9}, { 3.5689000794261055`*^9, 3.568900136257005*^9}, {3.5689001767858763`*^9, 3.56890020149632*^9}, {3.568902068726*^9, 3.568902074139209*^9}, 3.56890216040736*^9, 3.568902344113283*^9, {3.5689024237982235`*^9, 3.56890243907065*^9}, 3.5689025770684924`*^9, {3.568910195129884*^9, 3.5689102216343307`*^9}, 3.568910629060246*^9, {3.5689126699494305`*^9, 3.5689126747230387`*^9}}], Cell[BoxData[ RowBox[{"-", "61.94189044206585`"}]], "Output", CellChangeTimes->{ 3.56899972999459*^9, 3.569060092944422*^9, 3.569065230640086*^9, 3.569065277190568*^9, {3.5690653774051437`*^9, 3.569065406608395*^9}, 3.5690655623434687`*^9, 3.5690656441500125`*^9, 3.569065837980353*^9, 3.5690663625000744`*^9, 3.569066404963349*^9, {3.5690665346775765`*^9, 3.5690665559560137`*^9}, {3.5690666167181206`*^9, 3.5690666346425524`*^9}, {3.5690666763414254`*^9, 3.569066705825477*^9}, { 3.5690667510031567`*^9, 3.5690667710491915`*^9}, {3.569066820797679*^9, 3.5690668476453266`*^9}, 3.569066881044985*^9, 3.5690695000275574`*^9, 3.569069566655275*^9, 3.5690697220783477`*^9, 3.5690698101873026`*^9, 3.569069855520982*^9, 3.569069945969941*^9, 3.5690701761327453`*^9, 3.5690702508880763`*^9, 3.56907044067801*^9, 3.5690707657981806`*^9, 3.569071525527916*^9, 3.5690715829204164`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "<", RowBox[{"bc", "+", "1"}]}], ",", RowBox[{"i", "++"}], ",", RowBox[{"Vne", "=", RowBox[{"Vne", "-", RowBox[{"3", "*", RowBox[{ RowBox[{"cn", "[", RowBox[{"[", "i", "]"}], "]"}], "/", RowBox[{"u2", "[", RowBox[{"[", "i", "]"}], "]"}]}]}]}]}]}], "]"}], "\[IndentingNewLine]", "Vne", "\[IndentingNewLine]"}], "Input", CellChangeTimes->{{3.5688982251756487`*^9, 3.5688986440675845`*^9}, { 3.5688986873108606`*^9, 3.56889871560931*^9}, {3.5688987970414534`*^9, 3.5688988526711507`*^9}, {3.568898887147211*^9, 3.568898888644814*^9}, { 3.5688989206248703`*^9, 3.568899041759083*^9}, {3.5688990876855636`*^9, 3.568899258349863*^9}, 3.568899300267137*^9, {3.568899351326027*^9, 3.5688993831032825`*^9}, {3.5688994370481772`*^9, 3.568899468217032*^9}, { 3.5688995813016305`*^9, 3.5688997165538683`*^9}, {3.5688998495129013`*^9, 3.568899932364647*^9}, {3.568899967667509*^9, 3.568900024529609*^9}, { 3.5689000794261055`*^9, 3.568900136257005*^9}, {3.5689001767858763`*^9, 3.56890020149632*^9}, {3.568902068726*^9, 3.568902074139209*^9}, 3.56890216040736*^9, 3.568902344113283*^9, {3.5689024237982235`*^9, 3.56890243907065*^9}, 3.5689025770684924`*^9, {3.568910195129884*^9, 3.5689102262831383`*^9}, {3.5689126851126575`*^9, 3.5689126875306616`*^9}, 3.569070277688924*^9, 3.5690703146453886`*^9}], Cell[BoxData[ RowBox[{"-", "133.94189044206584`"}]], "Output", CellChangeTimes->{ 3.5689997300101895`*^9, 3.569060092960022*^9, 3.5690652306556864`*^9, 3.569065277206168*^9, {3.569065377420744*^9, 3.5690654066239953`*^9}, 3.569065562359069*^9, 3.569065644165612*^9, 3.5690658379959526`*^9, 3.5690663625312743`*^9, 3.5690664049789486`*^9, {3.5690665346931763`*^9, 3.569066555971614*^9}, {3.569066616733721*^9, 3.569066634658152*^9}, { 3.5690666763570256`*^9, 3.569066705841077*^9}, {3.569066751018757*^9, 3.5690667710647917`*^9}, {3.569066820813279*^9, 3.5690668476609263`*^9}, 3.569066881060585*^9, 3.5690695000431576`*^9, 3.5690695666708746`*^9, 3.5690697220939474`*^9, 3.5690698102029023`*^9, 3.569069855536582*^9, 3.5690699459855413`*^9, 3.569070176148345*^9, 3.5690702509036765`*^9, 3.56907044069361*^9, 3.569070765813781*^9, 3.5690715255435157`*^9, 3.569071582936016*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", RowBox[{"cw", "+", "1"}]}], ",", RowBox[{"i", "<", RowBox[{"nb", "+", "1"}]}], ",", RowBox[{"i", "++"}], ",", " ", RowBox[{ RowBox[{"Print", "[", RowBox[{ RowBox[{"u2", "[", RowBox[{"[", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], ",", "\"\< \>\"", ",", RowBox[{"a", "[", RowBox[{"[", RowBox[{"i", ",", "3"}], "]"}], "]"}]}], "]"}], ";", 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", "\[IndentingNewLine]"}], "Input", CellChangeTimes->{{3.5688982251756487`*^9, 3.5688986440675845`*^9}, { 3.5688986873108606`*^9, 3.56889871560931*^9}, {3.5688987970414534`*^9, 3.5688988526711507`*^9}, {3.568898887147211*^9, 3.568898888644814*^9}, { 3.5688989206248703`*^9, 3.568899041759083*^9}, {3.5688990876855636`*^9, 3.568899258349863*^9}, 3.568899300267137*^9, {3.568899351326027*^9, 3.5688993831032825`*^9}, {3.5688994370481772`*^9, 3.568899468217032*^9}, { 3.5688995813016305`*^9, 3.5688997165538683`*^9}, {3.5688998495129013`*^9, 3.568899932364647*^9}, {3.568899967667509*^9, 3.568900024529609*^9}, { 3.5689000794261055`*^9, 3.568900136257005*^9}, {3.5689001767858763`*^9, 3.56890020149632*^9}, {3.568902068726*^9, 3.568902074139209*^9}, 3.56890216040736*^9, 3.568902344113283*^9, {3.5689024237982235`*^9, 3.56890243907065*^9}, 3.5689025770684924`*^9, {3.568910195129884*^9, 3.5689102484819775`*^9}, {3.568910416821873*^9, 3.568910421314681*^9}, { 3.5689104658527594`*^9, 3.568910491858005*^9}, {3.5689105421368933`*^9, 3.5689105444144974`*^9}, {3.568912692943871*^9, 3.56891270358309*^9}, { 3.569066249867876*^9, 3.569066267152707*^9}, {3.569066330504418*^9, 3.569066351580055*^9}, 3.569066397116535*^9, {3.569066515349143*^9, 3.5690665477347994`*^9}, {3.5690666015392942`*^9, 3.569066668775412*^9}, { 3.5690667401611376`*^9, 3.5690668426533175`*^9}, 3.5690668756161757`*^9, { 3.5690694931947455`*^9, 3.569069493413146*^9}, {3.569069542990033*^9, 3.569069557092458*^9}, {3.5690697105343275`*^9, 3.5690697132643323`*^9}, 3.5690698018412876`*^9, {3.569069843321761*^9, 3.569069845178164*^9}, { 3.5690699336771193`*^9, 3.5690699355803227`*^9}, 3.569070167958331*^9, { 3.5690702075356007`*^9, 3.569070240436058*^9}, {3.5690704108975573`*^9, 3.5690704125199604`*^9}}], Cell[CellGroupData[{ Cell[BoxData[ InterpretationBox[ RowBox[{ "0.6`", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "1.210314249067387`"}], SequenceForm[0.6, " ", 1.210314249067387], Editable->False]], "Print", CellChangeTimes->{ 3.569066404994549*^9, {3.5690665347087765`*^9, 3.569066555987214*^9}, { 3.569066616749321*^9, 3.5690666346737523`*^9}, {3.5690666763726254`*^9, 3.569066705856677*^9}, {3.5690667510343566`*^9, 3.569066771080392*^9}, { 3.5690668208288794`*^9, 3.569066847676526*^9}, 3.569066881076185*^9, 3.569069500058758*^9, 3.569069566686475*^9, 3.5690697221095476`*^9, 3.569069810218502*^9, 3.569069855552182*^9, 3.569069946001141*^9, 3.569070176163945*^9, 3.5690702509192767`*^9, 3.56907044070921*^9, 3.569070765829381*^9, 3.5690715255591154`*^9, 3.5690715829516163`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "0.6`", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "1.2103203349966232`"}], SequenceForm[0.6, " ", 1.2103203349966232`], Editable->False]], "Print", CellChangeTimes->{ 3.569066404994549*^9, {3.5690665347087765`*^9, 3.569066555987214*^9}, { 3.569066616749321*^9, 3.5690666346737523`*^9}, {3.5690666763726254`*^9, 3.569066705856677*^9}, {3.5690667510343566`*^9, 3.569066771080392*^9}, { 3.5690668208288794`*^9, 3.569066847676526*^9}, 3.569066881076185*^9, 3.569069500058758*^9, 3.569069566686475*^9, 3.5690697221095476`*^9, 3.569069810218502*^9, 3.569069855552182*^9, 3.569069946001141*^9, 3.569070176163945*^9, 3.5690702509192767`*^9, 3.56907044070921*^9, 3.569070765829381*^9, 3.5690715255591154`*^9, 3.5690715829672165`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "0.6`", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "1.2103545683012018`"}], SequenceForm[0.6, " ", 1.2103545683012018`], Editable->False]], "Print", CellChangeTimes->{ 3.569066404994549*^9, {3.5690665347087765`*^9, 3.569066555987214*^9}, { 3.569066616749321*^9, 3.5690666346737523`*^9}, {3.5690666763726254`*^9, 3.569066705856677*^9}, {3.5690667510343566`*^9, 3.569066771080392*^9}, { 3.5690668208288794`*^9, 3.569066847676526*^9}, 3.569066881076185*^9, 3.569069500058758*^9, 3.569069566686475*^9, 3.5690697221095476`*^9, 3.569069810218502*^9, 3.569069855552182*^9, 3.569069946001141*^9, 3.569070176163945*^9, 3.5690702509192767`*^9, 3.56907044070921*^9, 3.569070765829381*^9, 3.5690715255591154`*^9, 3.5690715829672165`*^9}], Cell[BoxData[ InterpretationBox[ RowBox[{ "0.6`", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", "1.2102326842173277`"}], SequenceForm[0.6, " ", 1.2102326842173277`], Editable->False]], "Print", CellChangeTimes->{ 3.569066404994549*^9, {3.5690665347087765`*^9, 3.569066555987214*^9}, { 3.569066616749321*^9, 3.5690666346737523`*^9}, {3.5690666763726254`*^9, 3.569066705856677*^9}, {3.5690667510343566`*^9, 3.569066771080392*^9}, { 3.5690668208288794`*^9, 3.569066847676526*^9}, 3.569066881076185*^9, 3.569069500058758*^9, 3.569069566686475*^9, 3.5690697221095476`*^9, 3.569069810218502*^9, 3.569069855552182*^9, 3.569069946001141*^9, 3.569070176163945*^9, 3.5690702509192767`*^9, 3.56907044070921*^9, 3.569070765829381*^9, 3.5690715255591154`*^9, 3.5690715829672165`*^9}] }, Open ]], Cell[BoxData[ RowBox[{"-", "143.04451568572898`"}]], "Output", CellChangeTimes->{ 3.569066404994549*^9, {3.5690665347087765`*^9, 3.5690665560028143`*^9}, { 3.569066617045721*^9, 3.5690666346893525`*^9}, {3.5690666763882256`*^9, 3.569066705856677*^9}, {3.5690667510343566`*^9, 3.569066771080392*^9}, { 3.56906682101608*^9, 3.569066847848127*^9}, 3.569066881076185*^9, 3.569069500058758*^9, 3.5690695668736753`*^9, 3.569069722125148*^9, 3.5690698102341022`*^9, 3.569069855552182*^9, 3.569069946001141*^9, 3.569070176163945*^9, 3.5690702509348764`*^9, 3.5690704407248096`*^9, 3.569070765829381*^9, 3.5690715255591154`*^9, 3.5690715829672165`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{"Vee", "\[IndentingNewLine]", "Vnn", "\[IndentingNewLine]", \ "Vne", "\[IndentingNewLine]", RowBox[{"Vtot", "=", RowBox[{ "Vne", "+", "Vee", "+", "Vnn"}]}], "\[IndentingNewLine]", "Ttot", "\[IndentingNewLine]", RowBox[{ RowBox[{"-", "Vtot"}], "/", "Ttot"}], "\[IndentingNewLine]", RowBox[{"Etot", "=", RowBox[{"Ttot", "+", "Vtot"}]}], "\[IndentingNewLine]"}], "Input", CellChangeTimes->{{3.5689074950860367`*^9, 3.5689075109980645`*^9}, { 3.568910298074465*^9, 3.5689103015064707`*^9}, {3.568910724064413*^9, 3.5689107693044925`*^9}, {3.568955368910859*^9, 3.568955378629676*^9}, { 3.568955601366867*^9, 3.568955616218093*^9}}], Cell[BoxData["29.78684391486262`"], "Output", CellChangeTimes->{ 3.5689997300257893`*^9, 3.5690600929756217`*^9, 3.569065230671286*^9, 3.569065277221768*^9, {3.5690653774519444`*^9, 3.5690654066239953`*^9}, 3.569065562374669*^9, 3.5690656441812124`*^9, 3.569065838027153*^9, 3.5690663626872745`*^9, 3.569066405010149*^9, {3.5690665347243767`*^9, 3.569066556018414*^9}, {3.5690666170613213`*^9, 3.5690666347049522`*^9}, { 3.5690666764038258`*^9, 3.5690667058722773`*^9}, {3.5690667510499563`*^9, 3.569066771095992*^9}, {3.5690668210316796`*^9, 3.5690668478637266`*^9}, 3.569066881091785*^9, 3.569069500074358*^9, 3.5690695668892756`*^9, 3.569069722140748*^9, 3.5690698102497025`*^9, 3.5690698555677824`*^9, 3.569069946016741*^9, 3.5690701761795454`*^9, 3.5690702509504766`*^9, 3.56907044074041*^9, 3.569070765844981*^9, 3.5690715255747156`*^9, 3.5690715829828167`*^9}], Cell[BoxData["13.432102855050447`"], "Output", CellChangeTimes->{ 3.5689997300257893`*^9, 3.5690600929756217`*^9, 3.569065230671286*^9, 3.569065277221768*^9, {3.5690653774519444`*^9, 3.5690654066239953`*^9}, 3.569065562374669*^9, 3.5690656441812124`*^9, 3.569065838027153*^9, 3.5690663626872745`*^9, 3.569066405010149*^9, {3.5690665347243767`*^9, 3.569066556018414*^9}, {3.5690666170613213`*^9, 3.5690666347049522`*^9}, { 3.5690666764038258`*^9, 3.5690667058722773`*^9}, {3.5690667510499563`*^9, 3.569066771095992*^9}, {3.5690668210316796`*^9, 3.5690668478637266`*^9}, 3.569066881091785*^9, 3.569069500074358*^9, 3.5690695668892756`*^9, 3.569069722140748*^9, 3.5690698102497025`*^9, 3.5690698555677824`*^9, 3.569069946016741*^9, 3.5690701761795454`*^9, 3.5690702509504766`*^9, 3.56907044074041*^9, 3.569070765844981*^9, 3.5690715255747156`*^9, 3.5690715829828167`*^9}], Cell[BoxData[ RowBox[{"-", "143.04451568572898`"}]], "Output", CellChangeTimes->{ 3.5689997300257893`*^9, 3.5690600929756217`*^9, 3.569065230671286*^9, 3.569065277221768*^9, {3.5690653774519444`*^9, 3.5690654066239953`*^9}, 3.569065562374669*^9, 3.5690656441812124`*^9, 3.569065838027153*^9, 3.5690663626872745`*^9, 3.569066405010149*^9, {3.5690665347243767`*^9, 3.569066556018414*^9}, {3.5690666170613213`*^9, 3.5690666347049522`*^9}, { 3.5690666764038258`*^9, 3.5690667058722773`*^9}, {3.5690667510499563`*^9, 3.569066771095992*^9}, {3.5690668210316796`*^9, 3.5690668478637266`*^9}, 3.569066881091785*^9, 3.569069500074358*^9, 3.5690695668892756`*^9, 3.569069722140748*^9, 3.5690698102497025`*^9, 3.5690698555677824`*^9, 3.569069946016741*^9, 3.5690701761795454`*^9, 3.5690702509504766`*^9, 3.56907044074041*^9, 3.569070765844981*^9, 3.5690715255747156`*^9, 3.5690715829828167`*^9}], Cell[BoxData[ RowBox[{"-", "99.82556891581592`"}]], "Output", CellChangeTimes->{ 3.5689997300257893`*^9, 3.5690600929756217`*^9, 3.569065230671286*^9, 3.569065277221768*^9, {3.5690653774519444`*^9, 3.5690654066239953`*^9}, 3.569065562374669*^9, 3.5690656441812124`*^9, 3.569065838027153*^9, 3.5690663626872745`*^9, 3.569066405010149*^9, {3.5690665347243767`*^9, 3.569066556018414*^9}, {3.5690666170613213`*^9, 3.5690666347049522`*^9}, { 3.5690666764038258`*^9, 3.5690667058722773`*^9}, {3.5690667510499563`*^9, 3.569066771095992*^9}, {3.5690668210316796`*^9, 3.5690668478637266`*^9}, 3.569066881091785*^9, 3.569069500074358*^9, 3.5690695668892756`*^9, 3.569069722140748*^9, 3.5690698102497025`*^9, 3.5690698555677824`*^9, 3.569069946016741*^9, 3.5690701761795454`*^9, 3.5690702509504766`*^9, 3.56907044074041*^9, 3.569070765844981*^9, 3.5690715255747156`*^9, 3.5690715829828167`*^9}], Cell[BoxData["42.144018671283405`"], "Output", CellChangeTimes->{ 3.5689997300257893`*^9, 3.5690600929756217`*^9, 3.569065230671286*^9, 3.569065277221768*^9, {3.5690653774519444`*^9, 3.5690654066239953`*^9}, 3.569065562374669*^9, 3.5690656441812124`*^9, 3.569065838027153*^9, 3.5690663626872745`*^9, 3.569066405010149*^9, {3.5690665347243767`*^9, 3.569066556018414*^9}, {3.5690666170613213`*^9, 3.5690666347049522`*^9}, { 3.5690666764038258`*^9, 3.5690667058722773`*^9}, {3.5690667510499563`*^9, 3.569066771095992*^9}, {3.5690668210316796`*^9, 3.5690668478637266`*^9}, 3.569066881091785*^9, 3.569069500074358*^9, 3.5690695668892756`*^9, 3.569069722140748*^9, 3.5690698102497025`*^9, 3.5690698555677824`*^9, 3.569069946016741*^9, 3.5690701761795454`*^9, 3.5690702509504766`*^9, 3.56907044074041*^9, 3.569070765844981*^9, 3.5690715255747156`*^9, 3.5690715829828167`*^9}], Cell[BoxData["2.3686770285111955`"], "Output", CellChangeTimes->{ 3.5689997300257893`*^9, 3.5690600929756217`*^9, 3.569065230671286*^9, 3.569065277221768*^9, {3.5690653774519444`*^9, 3.5690654066239953`*^9}, 3.569065562374669*^9, 3.5690656441812124`*^9, 3.569065838027153*^9, 3.5690663626872745`*^9, 3.569066405010149*^9, {3.5690665347243767`*^9, 3.569066556018414*^9}, {3.5690666170613213`*^9, 3.5690666347049522`*^9}, { 3.5690666764038258`*^9, 3.5690667058722773`*^9}, {3.5690667510499563`*^9, 3.569066771095992*^9}, {3.5690668210316796`*^9, 3.5690668478637266`*^9}, 3.569066881091785*^9, 3.569069500074358*^9, 3.5690695668892756`*^9, 3.569069722140748*^9, 3.5690698102497025`*^9, 3.5690698555677824`*^9, 3.569069946016741*^9, 3.5690701761795454`*^9, 3.5690702509504766`*^9, 3.56907044074041*^9, 3.569070765844981*^9, 3.5690715255747156`*^9, 3.569071582998417*^9}], Cell[BoxData[ RowBox[{"-", "57.68155024453252`"}]], "Output", CellChangeTimes->{ 3.5689997300257893`*^9, 3.5690600929756217`*^9, 3.569065230671286*^9, 3.569065277221768*^9, {3.5690653774519444`*^9, 3.5690654066239953`*^9}, 3.569065562374669*^9, 3.5690656441812124`*^9, 3.569065838027153*^9, 3.5690663626872745`*^9, 3.569066405010149*^9, {3.5690665347243767`*^9, 3.569066556018414*^9}, {3.5690666170613213`*^9, 3.5690666347049522`*^9}, { 3.5690666764038258`*^9, 3.5690667058722773`*^9}, {3.5690667510499563`*^9, 3.569066771095992*^9}, {3.5690668210316796`*^9, 3.5690668478637266`*^9}, 3.569066881091785*^9, 3.569069500074358*^9, 3.5690695668892756`*^9, 3.569069722140748*^9, 3.5690698102497025`*^9, 3.5690698555677824`*^9, 3.569069946016741*^9, 3.5690701761795454`*^9, 3.5690702509504766`*^9, 3.56907044074041*^9, 3.569070765844981*^9, 3.5690715255747156`*^9, 3.569071582998417*^9}] }, Open ]] }, WindowSize->{1050, 939}, WindowMargins->{{51, Automatic}, {Automatic, -2}}, 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, 943, 20, 50, "Input"], Cell[CellGroupData[{ Cell[1525, 44, 415, 6, 29, "Input"], Cell[1943, 52, 1968, 35, 148, "Output"] }, Open ]], Cell[3926, 90, 906, 17, 50, "Input"], Cell[4835, 109, 762, 12, 29, "Input"], Cell[CellGroupData[{ Cell[5622, 125, 852, 13, 29, "Input"], Cell[6477, 140, 1009, 16, 28, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[7523, 161, 251, 5, 29, "Input"], Cell[7777, 168, 891, 13, 28, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[8705, 186, 1560, 29, 29, "Input"], Cell[10268, 217, 999, 16, 28, "Output"] }, Open ]], Cell[11282, 236, 1390, 22, 29, "Input"], Cell[12675, 260, 1566, 28, 29, "Input"], Cell[CellGroupData[{ Cell[14266, 292, 1234, 18, 29, "Input"], Cell[15503, 312, 1255, 25, 28, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[16795, 342, 697, 14, 29, "Input"], Cell[17495, 358, 969, 15, 28, "Output"] }, Open ]], Cell[18479, 376, 908, 24, 29, "Input"], Cell[CellGroupData[{ Cell[19412, 404, 3667, 89, 50, "Input"], Cell[23082, 495, 973, 16, 28, "Output"], Cell[24058, 513, 976, 16, 28, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[25071, 534, 519, 12, 29, "Input"], Cell[25593, 548, 1010, 16, 28, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[26640, 569, 1070, 31, 29, "Input"], Cell[27713, 602, 966, 16, 28, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[28716, 623, 701, 21, 90, "Input"], Cell[29420, 646, 877, 13, 28, "Output"], Cell[30300, 661, 877, 13, 28, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[31214, 679, 1494, 34, 110, "Input"], Cell[32711, 715, 901, 13, 28, "Output"] }, Open ]], Cell[33627, 731, 935, 19, 29, "Input"], Cell[34565, 752, 3535, 66, 70, "Input"], Cell[38103, 820, 3216, 56, 29, "Input"], Cell[CellGroupData[{ Cell[41344, 880, 699, 11, 29, "Input"], Cell[42046, 893, 3874, 86, 57, "Output"] }, Open ]], Cell[45935, 982, 1254, 35, 70, "Input"], Cell[CellGroupData[{ Cell[47214, 1021, 3842, 77, 90, "Input"], Cell[51059, 1100, 1502, 31, 97, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[52598, 1136, 3387, 76, 130, "Input"], Cell[55988, 1214, 992, 16, 28, "Output"] }, Open ]], Cell[56995, 1233, 1708, 43, 90, "Input"], Cell[58706, 1278, 511, 11, 50, "Input"], Cell[59220, 1291, 2048, 48, 90, "Input"], Cell[61271, 1341, 2180, 51, 70, "Input"], Cell[63454, 1394, 1345, 34, 50, "Input"], Cell[CellGroupData[{ Cell[64824, 1432, 923, 20, 29, "Input"], Cell[65750, 1454, 907, 13, 28, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[66694, 1472, 611, 17, 29, "Input"], Cell[67308, 1491, 888, 13, 28, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[68233, 1509, 507, 10, 29, "Input"], Cell[68743, 1521, 907, 13, 28, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[69687, 1539, 823, 23, 29, "Input"], Cell[70513, 1564, 907, 13, 28, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[71457, 1582, 204, 4, 29, "Input"], Cell[71664, 1588, 894, 13, 28, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[72595, 1606, 1679, 48, 110, "Input"], Cell[74277, 1656, 900, 13, 28, "Output"] }, Open ]], Cell[75192, 1672, 1971, 35, 90, "Input"], Cell[77166, 1709, 1638, 24, 29, "Input"], Cell[78807, 1735, 1655, 24, 29, "Input"], Cell[80465, 1761, 2877, 57, 70, "Input"], Cell[83345, 1820, 2696, 53, 90, "Input"], Cell[CellGroupData[{ Cell[86066, 1877, 939, 17, 70, "Input"], Cell[87008, 1896, 1278, 24, 28, "Output"], Cell[88289, 1922, 901, 13, 28, "Output"], Cell[89193, 1937, 901, 13, 28, "Output"] }, Open ]], Cell[90109, 1953, 6921, 129, 90, "Input"], Cell[97033, 2084, 7013, 131, 90, "Input"], Cell[104049, 2217, 6921, 129, 90, "Input"], Cell[110973, 2348, 5982, 97, 29, "Input"], Cell[116958, 2447, 6882, 128, 90, "Input"], Cell[123843, 2577, 6884, 128, 90, "Input"], Cell[130730, 2707, 6882, 128, 90, "Input"], Cell[137615, 2837, 5931, 96, 50, "Input"], Cell[CellGroupData[{ Cell[143571, 2937, 6916, 125, 150, "Input"], Cell[150490, 3064, 3849, 82, 377, "Output"], Cell[154342, 3148, 3389, 67, 377, "Output"] }, Open ]], Cell[157746, 3218, 2645, 66, 110, "Input"], Cell[CellGroupData[{ Cell[160416, 3288, 1163, 18, 50, "Input"], Cell[161582, 3308, 898, 13, 28, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[162517, 3326, 2746, 64, 90, "Input"], Cell[165266, 3392, 912, 13, 28, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[166215, 3410, 2740, 64, 90, "Input"], Cell[168958, 3476, 901, 13, 28, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[169896, 3494, 2764, 64, 90, "Input"], Cell[172663, 3560, 904, 13, 28, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[173604, 3578, 1505, 39, 130, "Input"], Cell[175112, 3619, 912, 13, 28, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[176061, 3637, 2865, 67, 90, "Input"], Cell[178929, 3706, 922, 14, 28, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[179888, 3725, 2900, 68, 90, "Input"], Cell[182791, 3795, 925, 14, 28, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[183753, 3814, 1502, 29, 70, "Input"], Cell[185258, 3845, 923, 14, 28, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[186218, 3864, 3087, 65, 70, "Input"], Cell[CellGroupData[{ Cell[189330, 3933, 820, 15, 21, "Print"], Cell[190153, 3950, 823, 15, 21, "Print"], Cell[190979, 3967, 823, 15, 21, "Print"], Cell[191805, 3984, 823, 15, 21, "Print"] }, Open ]], Cell[192643, 4002, 673, 10, 28, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[193353, 4017, 669, 13, 170, "Input"], Cell[194025, 4032, 913, 13, 28, "Output"], Cell[194941, 4047, 914, 13, 28, "Output"], Cell[195858, 4062, 931, 14, 28, "Output"], Cell[196792, 4078, 930, 14, 28, "Output"], Cell[197725, 4094, 914, 13, 28, "Output"], Cell[198642, 4109, 912, 13, 28, "Output"], Cell[199557, 4124, 928, 14, 60, "Output"] }, Open ]] } ] *) (* End of internal cache information *)