(* 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[ 110385, 2692] NotebookOptionsPosition[ 108165, 2621] NotebookOutlinePosition[ 108511, 2636] CellTagsIndexPosition[ 108468, 2633] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[BoxData[ RowBox[{"(*", " ", RowBox[{"Ethene", ",", " ", RowBox[{"central", " ", "pi", " ", "cloud", " ", "23.01", ".2013"}]}], " ", "*)"}]], "Input", CellChangeTimes->{{3.5675987354944196`*^9, 3.567599332289068*^9}, { 3.567599364347124*^9, 3.567599388854767*^9}, {3.567599446824469*^9, 3.567599671277663*^9}, {3.5675997743158445`*^9, 3.567599879038828*^9}, { 3.567599920722101*^9, 3.567600053883935*^9}, {3.5676001118380375`*^9, 3.5676001705877404`*^9}, {3.5676002126766143`*^9, 3.5676003454796476`*^9}, {3.5676004083477583`*^9, 3.567600442527418*^9}, { 3.5676004745698743`*^9, 3.567600484132691*^9}, {3.5676005145995445`*^9, 3.5676005320091753`*^9}, {3.5676005791056576`*^9, 3.5676012559440465`*^9}, {3.5676013203097596`*^9, 3.5676013219321623`*^9}, 3.5676013524770164`*^9, {3.5676850872857094`*^9, 3.5676850878785105`*^9}, 3.567685750380474*^9, {3.5678518143429227`*^9, 3.567851822579737*^9}, { 3.567852306882588*^9, 3.567852307568989*^9}, {3.567921625688327*^9, 3.5679216259379272`*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"Clear", "[", RowBox[{ "Cr", ",", "CC", ",", "CH", ",", "Cpi", ",", "nc", ",", "nn", ",", "xn", ",", "yn", ",", "zn", ",", "xc", ",", "yc", ",", "zc", ",", "cx", ",", "chx", ",", "chy", ",", "\[IndentingNewLine]", "cchx", ",", "cchy", ",", "Ekin", ",", "vne", ",", "vee", ",", "vnn", ",", "Epot", ",", "w", ",", "csw", ",", "ssw", ",", "p", ",", "sq", ",", "sn", ",", "sr", ",", "\[IndentingNewLine]", "k1", ",", "k2", ",", "k3", ",", "k4", ",", "s1", ",", "s2", ",", "s3", ",", "s4", ",", "Z"}], "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"nc", "=", "9"}], ";", " ", RowBox[{"nn", "=", "6"}], ";", " ", RowBox[{"Z", "=", "6.0"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"bohr", "=", "0.529177"}], ";", RowBox[{"rad", "=", "57.29578"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"k1", "=", "1.029"}], ";", " ", RowBox[{"k2", "=", "1.21"}], ";", " ", RowBox[{"k3", "=", "1.25"}], ";", " ", RowBox[{"k4", "=", "1.1"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"s1", "=", "0.3"}], ";", RowBox[{"s2", "=", "0.3"}], ";", " ", RowBox[{"s3", "=", "0.3"}], ";"}], " "}], "\[IndentingNewLine]", RowBox[{ RowBox[{"as", "=", RowBox[{"{", RowBox[{ "\"\< \>\"", ",", "C", ",", "C", ",", "H", ",", "H", ",", "H", ",", "H"}], "}"}]}], ";", RowBox[{"(*", RowBox[{"atom", " ", "string"}], "*)"}]}]}], "Input", CellChangeTimes->{{3.5675987354944196`*^9, 3.567599332289068*^9}, { 3.567599364347124*^9, 3.567599388854767*^9}, {3.567599446824469*^9, 3.567599671277663*^9}, {3.5675997743158445`*^9, 3.567599879038828*^9}, { 3.567599920722101*^9, 3.567600053883935*^9}, {3.5676001118380375`*^9, 3.5676001705877404`*^9}, {3.5676002126766143`*^9, 3.5676003454796476`*^9}, {3.5676004083477583`*^9, 3.567600442527418*^9}, { 3.5676004745698743`*^9, 3.567600484132691*^9}, {3.5676005145995445`*^9, 3.5676005320091753`*^9}, {3.5676005791056576`*^9, 3.5676012559440465`*^9}, {3.5676013203097596`*^9, 3.5676013219321623`*^9}, 3.5676013524770164`*^9, {3.5676015878502293`*^9, 3.56760173054368*^9}, { 3.5676125832717414`*^9, 3.5676126027561755`*^9}, {3.5676126533938646`*^9, 3.5676126797423105`*^9}, {3.5676127131263695`*^9, 3.5676127364484105`*^9}, {3.567613357173501*^9, 3.567613363585112*^9}, { 3.5676168815728908`*^9, 3.567616936672188*^9}, {3.5676611367108335`*^9, 3.567661137740435*^9}, {3.5676854112046785`*^9, 3.567685413638283*^9}, 3.5676857759021187`*^9, {3.5676860335109715`*^9, 3.5676860647422266`*^9}, { 3.567687447325855*^9, 3.5676874507266607`*^9}, {3.567851420645031*^9, 3.5678514213782325`*^9}, {3.5678514547154913`*^9, 3.5678514824835396`*^9}, {3.5678516717430725`*^9, 3.567851672663474*^9}, { 3.5678521501335125`*^9, 3.5678521703511477`*^9}, {3.567852210802019*^9, 3.567852211816021*^9}, {3.567852407736765*^9, 3.567852432260008*^9}, { 3.5678525757958603`*^9, 3.5678525766694613`*^9}, {3.5678531332784395`*^9, 3.567853135010042*^9}, {3.5678815904829903`*^9, 3.5678816586707096`*^9}, { 3.5679215656750216`*^9, 3.567921590494665*^9}}], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{ RowBox[{"6", " ", RowBox[{"Variable", ":", " ", "Cr"}]}], ",", "CC", ",", "CH", ",", "Cpi", ",", "p", ",", "w"}], " ", "*)"}]], "Input", CellChangeTimes->{{3.5675987354944196`*^9, 3.567599332289068*^9}, { 3.567599364347124*^9, 3.567599388854767*^9}, {3.567599446824469*^9, 3.567599671277663*^9}, {3.5675997743158445`*^9, 3.567599879038828*^9}, { 3.567599920722101*^9, 3.567600053883935*^9}, {3.5676001118380375`*^9, 3.5676001705877404`*^9}, {3.5676002126766143`*^9, 3.5676003454796476`*^9}, {3.5676004083477583`*^9, 3.567600442527418*^9}, { 3.5676004745698743`*^9, 3.567600484132691*^9}, {3.5676005145995445`*^9, 3.5676005320091753`*^9}, {3.5676005791056576`*^9, 3.5676012559440465`*^9}, {3.5676013203097596`*^9, 3.5676013219321623`*^9}, {3.5676013524770164`*^9, 3.5676013644734373`*^9}, 3.56761275877205*^9, 3.5676857814557285`*^9}], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{"C", " ", "He", " ", "shell"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Ekin", "=", RowBox[{"2.", "*", "2.25", "*", RowBox[{"k1", "/", RowBox[{"Cr", "^", "2"}]}]}]}], ";"}], "\n", RowBox[{ RowBox[{"vee", "=", RowBox[{"2.", "*", "3.0", "*", RowBox[{"s1", "/", "Cr"}]}]}], ";"}], "\n", RowBox[{ RowBox[{"vne", "=", RowBox[{ RowBox[{"-", "2."}], "*", "3.0", "*", RowBox[{"Z", "/", "Cr"}]}]}], ";"}], "\[IndentingNewLine]"}]}]], "Input", CellChangeTimes->{{3.5675987354944196`*^9, 3.567599332289068*^9}, { 3.567599364347124*^9, 3.567599388854767*^9}, {3.567599446824469*^9, 3.567599671277663*^9}, {3.5675997743158445`*^9, 3.567599879038828*^9}, { 3.567599920722101*^9, 3.567600053883935*^9}, {3.5676001118380375`*^9, 3.5676001705877404`*^9}, {3.5676002126766143`*^9, 3.5676003454796476`*^9}, { 3.5676004083477583`*^9, 3.567600442527418*^9}, {3.5676004745698743`*^9, 3.567600484132691*^9}, {3.5676005145995445`*^9, 3.5676005320091753`*^9}, { 3.5676005791056576`*^9, 3.5676012559440465`*^9}, {3.5676013203097596`*^9, 3.5676013219321623`*^9}, {3.5676013524770164`*^9, 3.5676013644734373`*^9}, { 3.5676917901504855`*^9, 3.567691799978503*^9}}], Cell[BoxData[ RowBox[{"(*", RowBox[{"bonding", " ", "pairs"}], " ", "*)"}]], "Input", CellChangeTimes->{{3.5675987354944196`*^9, 3.567599332289068*^9}, { 3.567599364347124*^9, 3.567599388854767*^9}, {3.567599446824469*^9, 3.567599671277663*^9}, {3.5675997743158445`*^9, 3.567599879038828*^9}, { 3.567599920722101*^9, 3.567600053883935*^9}, {3.5676001118380375`*^9, 3.5676001705877404`*^9}, {3.5676002126766143`*^9, 3.5676003454796476`*^9}, { 3.5676004083477583`*^9, 3.567600442527418*^9}, {3.5676004745698743`*^9, 3.567600484132691*^9}, {3.5676005145995445`*^9, 3.5676005320091753`*^9}, { 3.5676005791056576`*^9, 3.5676012559440465`*^9}, {3.5676013203097596`*^9, 3.5676013219321623`*^9}, {3.5676013524770164`*^9, 3.567601370448248*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"Ekin", "=", RowBox[{"Ekin", "+", " ", RowBox[{"2.25", "*", RowBox[{"(", RowBox[{ RowBox[{"k2", "/", RowBox[{"CC", "^", "2"}]}], "+", RowBox[{"4.0", "*", RowBox[{"k3", "/", RowBox[{"CH", "^", "2"}]}]}], "+", RowBox[{"k4", "/", RowBox[{"Cpi", "^", "2"}]}]}], ")"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"vee", "=", RowBox[{"vee", "+", RowBox[{"3.0", "*", RowBox[{"(", RowBox[{ RowBox[{"s2", "/", "CC"}], "+", RowBox[{"4.0", "*", RowBox[{"s3", "/", "CH"}]}]}], ")"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"vne", "=", RowBox[{"vne", "-", RowBox[{"4.0", "*", RowBox[{ RowBox[{"(", RowBox[{"3.0", "-", RowBox[{ RowBox[{"(", RowBox[{"p", "/", "CH"}], ")"}], "^", "2"}]}], ")"}], "/", "CH"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"vnn", "=", "0.0"}], ";"}]}], "Input", CellChangeTimes->{{3.5675987354944196`*^9, 3.567599332289068*^9}, { 3.567599364347124*^9, 3.567599388854767*^9}, {3.567599446824469*^9, 3.567599671277663*^9}, {3.5675997743158445`*^9, 3.567599879038828*^9}, { 3.567599920722101*^9, 3.567600053883935*^9}, {3.5676001118380375`*^9, 3.5676001705877404`*^9}, {3.5676002126766143`*^9, 3.5676003454796476`*^9}, {3.5676004083477583`*^9, 3.567600442527418*^9}, { 3.5676004745698743`*^9, 3.567600484132691*^9}, {3.5676005145995445`*^9, 3.5676005320091753`*^9}, {3.5676005791056576`*^9, 3.5676012559440465`*^9}, {3.5676013203097596`*^9, 3.5676013219321623`*^9}, {3.5676013524770164`*^9, 3.567601370448248*^9}, { 3.5676046193892126`*^9, 3.5676046713841043`*^9}, {3.567611708656205*^9, 3.567611709685807*^9}, {3.567611936042205*^9, 3.5676119371030064`*^9}, { 3.567612182787838*^9, 3.5676121931306562`*^9}, 3.5676857932649493`*^9, { 3.5676918030049086`*^9, 3.5676918121465244`*^9}, {3.5676919908760386`*^9, 3.567691999190853*^9}, 3.567852293419764*^9}], Cell[BoxData[{ RowBox[{ RowBox[{"cn", "=", RowBox[{"{", RowBox[{"Z", ",", "Z", ",", "1.", ",", "1.", ",", "1.", ",", "1."}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"cc", "=", RowBox[{"{", RowBox[{ RowBox[{"-", "2."}], ",", RowBox[{"-", "2."}], ",", RowBox[{"-", "2."}], ",", RowBox[{"-", "2."}], ",", RowBox[{"-", "2."}], ",", RowBox[{"-", "2."}], ",", RowBox[{"-", "2."}], ",", RowBox[{"-", "1."}], ",", RowBox[{"-", "1."}]}], "}"}]}], ";"}]}], "Input", CellChangeTimes->{{3.5675987354944196`*^9, 3.567599332289068*^9}, { 3.567599364347124*^9, 3.567599388854767*^9}, {3.567599446824469*^9, 3.567599671277663*^9}, {3.5675997743158445`*^9, 3.567599879038828*^9}, { 3.567599920722101*^9, 3.567600053883935*^9}, {3.5676001118380375`*^9, 3.5676001705877404`*^9}, {3.5676002126766143`*^9, 3.5676003454796476`*^9}, {3.5676004083477583`*^9, 3.567600442527418*^9}, { 3.5676004745698743`*^9, 3.567600484132691*^9}, {3.5676005145995445`*^9, 3.5676005320091753`*^9}, {3.5676005791056576`*^9, 3.5676012559440465`*^9}, {3.5676013203097596`*^9, 3.5676013219321623`*^9}, {3.5676013524770164`*^9, 3.567601375861457*^9}, { 3.5676918198373375`*^9, 3.5676918343609633`*^9}, {3.5678516559090443`*^9, 3.5678516638338585`*^9}, 3.5678522868677526`*^9, {3.5678523210006123`*^9, 3.5678523224826155`*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"csw", "=", RowBox[{"Cos", "[", "w", "]"}]}], ";", " ", RowBox[{"ssw", "=", RowBox[{"Sin", "[", "w", "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"cx", "=", RowBox[{"Cr", "+", "CC"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"chx", "=", " ", RowBox[{"cx", "+", RowBox[{ RowBox[{"(", RowBox[{"Cr", "+", "CH", "+", "p"}], ")"}], "*", "csw"}]}]}], ";", " ", RowBox[{"chy", "=", " ", RowBox[{ RowBox[{"(", RowBox[{"Cr", "+", "CH", "+", "p"}], ")"}], "*", "ssw"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"cchx", "=", RowBox[{"cx", "+", RowBox[{ RowBox[{"(", RowBox[{"Cr", "+", "CH"}], ")"}], "*", "csw"}]}]}], ";", " ", RowBox[{"cchy", "=", RowBox[{ RowBox[{"(", RowBox[{"Cr", "+", "CH"}], ")"}], "*", "ssw"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"ccpi", "=", RowBox[{"Cr", "+", "Cpi"}]}], ";"}]}], "Input", CellChangeTimes->{{3.5675987354944196`*^9, 3.567599332289068*^9}, { 3.567599364347124*^9, 3.567599388854767*^9}, {3.567599446824469*^9, 3.567599671277663*^9}, {3.5675997743158445`*^9, 3.567599879038828*^9}, { 3.567599920722101*^9, 3.567600053883935*^9}, {3.5676001118380375`*^9, 3.5676001705877404`*^9}, {3.5676002126766143`*^9, 3.5676003454796476`*^9}, {3.5676004083477583`*^9, 3.567600442527418*^9}, { 3.5676004745698743`*^9, 3.567600484132691*^9}, {3.5676005145995445`*^9, 3.5676005320091753`*^9}, {3.5676005791056576`*^9, 3.5676012559440465`*^9}, {3.5676013203097596`*^9, 3.5676013219321623`*^9}, {3.5676013524770164`*^9, 3.567601384847073*^9}, { 3.567601484843249*^9, 3.5676015176657066`*^9}, 3.5676019519548693`*^9, { 3.5676050017146845`*^9, 3.567605007845495*^9}, {3.5676052437491093`*^9, 3.5676052486163177`*^9}, {3.567660817518673*^9, 3.5676608349127035`*^9}, 3.5676855136344585`*^9, {3.567685802016565*^9, 3.567685806400172*^9}, { 3.5678513921125813`*^9, 3.567851395326187*^9}, {3.5678523287226257`*^9, 3.56785233117183*^9}, {3.5678523823555202`*^9, 3.5678523836347227`*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"xn", "=", RowBox[{"{", RowBox[{"cx", ",", RowBox[{"-", "cx"}], ",", "chx", ",", " ", "chx", ",", RowBox[{"-", "chx"}], ",", RowBox[{"-", "chx"}]}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"yn", "=", RowBox[{"{", RowBox[{"0", ",", " ", "0", ",", "chy", ",", RowBox[{"-", "chy"}], ",", " ", "chy", ",", RowBox[{"-", "chy"}]}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"zn", "=", RowBox[{"{", RowBox[{"0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0"}], "}"}]}], ";"}]}], "Input", CellChangeTimes->{{3.5675987354944196`*^9, 3.567599332289068*^9}, { 3.567599364347124*^9, 3.567599388854767*^9}, {3.567599446824469*^9, 3.567599671277663*^9}, {3.5675997743158445`*^9, 3.567599879038828*^9}, { 3.567599920722101*^9, 3.567600053883935*^9}, {3.5676001118380375`*^9, 3.5676001705877404`*^9}, {3.5676002126766143`*^9, 3.5676003454796476`*^9}, {3.5676004083477583`*^9, 3.567600442527418*^9}, { 3.5676004745698743`*^9, 3.567600484132691*^9}, {3.5676005145995445`*^9, 3.5676005320091753`*^9}, {3.5676005791056576`*^9, 3.5676012559440465`*^9}, {3.5676013203097596`*^9, 3.5676013219321623`*^9}, {3.5676013524770164`*^9, 3.5676013901354823`*^9}, {3.5676048225327697`*^9, 3.567604838756798*^9}, { 3.5676114185892954`*^9, 3.5676114417709365`*^9}, {3.5676854961000276`*^9, 3.5676855041184416`*^9}, 3.5678522801441407`*^9}], Cell[BoxData[{ RowBox[{ RowBox[{"xc", "=", RowBox[{"{", RowBox[{"cx", ",", RowBox[{"-", "cx"}], ",", "cchx", ",", " ", "cchx", ",", RowBox[{"-", "cchx"}], ",", RowBox[{"-", "cchx"}], ",", "0", ",", " ", "0", ",", " ", "0"}], " ", "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"yc", "=", RowBox[{"{", RowBox[{"0", ",", " ", "0", ",", "cchy", ",", RowBox[{"-", "cchy"}], ",", " ", "cchy", ",", RowBox[{"-", "cchy"}], ",", "0", ",", " ", "0", ",", " ", "0"}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"zc", "=", RowBox[{"{", RowBox[{ "0", ",", " ", "0", ",", " ", "0", ",", " ", "0", ",", " ", "0", ",", " ", "0", ",", "0", ",", "ccpi", ",", RowBox[{"-", "ccpi"}]}], "}"}]}], ";"}]}], "Input", CellChangeTimes->{{3.5675987354944196`*^9, 3.567599332289068*^9}, { 3.567599364347124*^9, 3.567599388854767*^9}, {3.567599446824469*^9, 3.567599671277663*^9}, {3.5675997743158445`*^9, 3.567599879038828*^9}, { 3.567599920722101*^9, 3.567600053883935*^9}, {3.5676001118380375`*^9, 3.5676001705877404`*^9}, {3.5676002126766143`*^9, 3.5676003454796476`*^9}, {3.5676004083477583`*^9, 3.567600442527418*^9}, { 3.5676004745698743`*^9, 3.567600484132691*^9}, {3.5676005145995445`*^9, 3.5676005320091753`*^9}, {3.5676005791056576`*^9, 3.5676012559440465`*^9}, {3.5676013203097596`*^9, 3.5676013219321623`*^9}, {3.5676013524770164`*^9, 3.5676013952366915`*^9}, {3.567604906585717*^9, 3.5676049391429744`*^9}, { 3.5676607907646255`*^9, 3.567660802589446*^9}, {3.567685532853692*^9, 3.5676856314926653`*^9}, {3.567850970958641*^9, 3.567850994327482*^9}, { 3.5678513525821114`*^9, 3.567851378400157*^9}, {3.567851597627342*^9, 3.5678516463150277`*^9}, {3.5678516918827076`*^9, 3.567851696125915*^9}, 3.567852268990121*^9}], Cell[BoxData[ RowBox[{ RowBox[{"(*", "ww", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "<", "nc"}], ",", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"j", "=", RowBox[{"i", "+", "1"}]}], ",", RowBox[{"j", "<", RowBox[{"nc", "+", "1"}]}], ",", RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", RowBox[{"vee", "=", RowBox[{"vee", "+", RowBox[{ RowBox[{"cc", "[", RowBox[{"[", "i", "]"}], "]"}], "*", RowBox[{ RowBox[{"cc", "[", RowBox[{"[", "j", "]"}], "]"}], "/", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "i", "]"}], "]"}], "-", RowBox[{"xc", "[", RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"yc", "[", RowBox[{"[", "i", "]"}], "]"}], "-", RowBox[{"yc", "[", RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"zc", "[", RowBox[{"[", "i", "]"}], "]"}], "-", RowBox[{"zc", "[", RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}]}], "]"}]}]}]}]}]}], "]"}]}], "]"}], ";"}]}]], "Input", CellChangeTimes->{{3.5675987354944196`*^9, 3.567599332289068*^9}, { 3.567599364347124*^9, 3.567599388854767*^9}, {3.567599446824469*^9, 3.567599671277663*^9}, {3.5675997743158445`*^9, 3.567599879038828*^9}, { 3.567599920722101*^9, 3.567600053883935*^9}, {3.5676001118380375`*^9, 3.5676001705877404`*^9}, {3.5676002126766143`*^9, 3.5676003454796476`*^9}, {3.5676004083477583`*^9, 3.567600442527418*^9}, { 3.5676004745698743`*^9, 3.567600484132691*^9}, {3.5676005145995445`*^9, 3.5676005320091753`*^9}, {3.5676005791056576`*^9, 3.5676012559440465`*^9}, {3.5676013203097596`*^9, 3.5676013219321623`*^9}, {3.5676013524770164`*^9, 3.5676014051115084`*^9}, {3.5676177595268326`*^9, 3.567617780103269*^9}, { 3.5676178274649525`*^9, 3.5676178737346334`*^9}, {3.5676608539135365`*^9, 3.5676608549587383`*^9}, {3.567686936378157*^9, 3.5676870257195144`*^9}, { 3.567851087802847*^9, 3.5678510997212677`*^9}, 3.5678511305469217`*^9, { 3.5678815497981186`*^9, 3.5678815515453215`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"(*", "nn", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "<", "nn"}], ",", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"j", "=", RowBox[{"i", "+", "1"}]}], ",", RowBox[{"j", "<", RowBox[{"nn", "+", "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[{"xn", "[", RowBox[{"[", "i", "]"}], "]"}], "-", RowBox[{"xn", "[", RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"yn", "[", RowBox[{"[", "i", "]"}], "]"}], "-", RowBox[{"yn", "[", RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"zn", "[", RowBox[{"[", "i", "]"}], "]"}], "-", RowBox[{"zn", "[", RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}]}], "]"}]}]}]}]}]}], "]"}]}], "]"}], ";"}]}]], "Input", CellChangeTimes->{{3.5675987354944196`*^9, 3.567599332289068*^9}, { 3.567599364347124*^9, 3.567599388854767*^9}, {3.567599446824469*^9, 3.567599671277663*^9}, {3.5675997743158445`*^9, 3.567599879038828*^9}, { 3.567599920722101*^9, 3.567600053883935*^9}, {3.5676001118380375`*^9, 3.5676001705877404`*^9}, {3.5676002126766143`*^9, 3.5676003454796476`*^9}, {3.5676004083477583`*^9, 3.567600442527418*^9}, { 3.5676004745698743`*^9, 3.567600484132691*^9}, {3.5676005145995445`*^9, 3.5676005320091753`*^9}, {3.5676005791056576`*^9, 3.5676012559440465`*^9}, {3.5676013203097596`*^9, 3.5676013219321623`*^9}, {3.5676013524770164`*^9, 3.56760141183512*^9}, { 3.5676636515600505`*^9, 3.5676636621992693`*^9}, {3.5676872211098576`*^9, 3.5676872966451902`*^9}, 3.5676928599223647`*^9, {3.5678511084728827`*^9, 3.5678511270213156`*^9}, 3.5678522555584974`*^9}], Cell[BoxData[ RowBox[{ RowBox[{"(*", "nw", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"i", "=", "1"}], ",", RowBox[{"i", "<", RowBox[{"nc", "+", "1"}]}], ",", RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", RowBox[{"For", "[", RowBox[{ RowBox[{"j", "=", "1"}], ",", RowBox[{"j", "<", RowBox[{"nn", "+", "1"}]}], ",", RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", RowBox[{"If", "[", RowBox[{ RowBox[{"(", RowBox[{"i", "\[NotEqual]", "j"}], ")"}], ",", RowBox[{"vne", "=", RowBox[{"vne", "+", RowBox[{ RowBox[{"cc", "[", RowBox[{"[", "i", "]"}], "]"}], "*", RowBox[{ RowBox[{"cn", "[", RowBox[{"[", "j", "]"}], "]"}], "/", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "i", "]"}], "]"}], "-", RowBox[{"xn", "[", RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"yc", "[", RowBox[{"[", "i", "]"}], "]"}], "-", RowBox[{"yn", "[", RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"zc", "[", RowBox[{"[", "i", "]"}], "]"}], "-", RowBox[{"zn", "[", RowBox[{"[", "j", "]"}], "]"}]}], ")"}], "^", "2"}]}], "]"}]}]}]}]}]}], "]"}]}], "]"}]}], "]"}], ";"}]}]], "Input", CellChangeTimes->{{3.5675987354944196`*^9, 3.567599332289068*^9}, { 3.567599364347124*^9, 3.567599388854767*^9}, {3.567599446824469*^9, 3.567599671277663*^9}, {3.5675997743158445`*^9, 3.567599879038828*^9}, { 3.567599920722101*^9, 3.567600053883935*^9}, {3.5676001118380375`*^9, 3.5676001705877404`*^9}, {3.5676002126766143`*^9, 3.5676003454796476`*^9}, {3.5676004083477583`*^9, 3.567600442527418*^9}, { 3.5676004745698743`*^9, 3.567600484132691*^9}, {3.5676005145995445`*^9, 3.5676005320091753`*^9}, {3.5676005791056576`*^9, 3.5676012559440465`*^9}, {3.5676013203097596`*^9, 3.5676013219321623`*^9}, {3.5676013524770164`*^9, 3.567601421335537*^9}, { 3.5676156065414515`*^9, 3.5676156341635*^9}, 3.567663665537675*^9, { 3.5676849746847115`*^9, 3.5676849977571526`*^9}, 3.5676858120005827`*^9, { 3.5676873370648613`*^9, 3.5676874232706127`*^9}, 3.567692864898773*^9, { 3.5678511442905455`*^9, 3.567851167238186*^9}, {3.567851531732826*^9, 3.5678515386280384`*^9}}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"Epot", "=", RowBox[{"vne", "+", "vee", "+", "vnn"}]}], ";"}], "\n", RowBox[{ RowBox[{"func", "=", RowBox[{"Ekin", "+", "Epot"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"t", "=", RowBox[{"FindMinimum", "[", RowBox[{"func", ",", RowBox[{"{", RowBox[{"Cr", ",", "0.276"}], "}"}], ",", RowBox[{"{", RowBox[{"CC", ",", "0.986"}], "}"}], ",", RowBox[{"{", RowBox[{"CH", ",", "1.23"}], "}"}], ",", RowBox[{"{", RowBox[{"Cpi", ",", "2.0"}], "}"}], ",", RowBox[{"{", RowBox[{"w", ",", "1.02"}], "}"}], ",", RowBox[{"{", RowBox[{"p", ",", "0.56"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Method", "\[Rule]", "\"\\""}], ",", RowBox[{"MaxIterations", "\[Rule]", "800"}]}], "}"}]}], "]"}]}]}], "Input", CellChangeTimes->{{3.5675987354944196`*^9, 3.567599332289068*^9}, { 3.567599364347124*^9, 3.567599388854767*^9}, {3.567599446824469*^9, 3.567599671277663*^9}, {3.5675997743158445`*^9, 3.567599879038828*^9}, { 3.567599920722101*^9, 3.567600053883935*^9}, {3.5676001118380375`*^9, 3.5676001705877404`*^9}, {3.5676002126766143`*^9, 3.5676003454796476`*^9}, {3.5676004083477583`*^9, 3.567600442527418*^9}, { 3.5676004745698743`*^9, 3.567600484132691*^9}, {3.5676005145995445`*^9, 3.5676005320091753`*^9}, {3.5676005791056576`*^9, 3.5676012559440465`*^9}, {3.5676013203097596`*^9, 3.5676013219321623`*^9}, {3.5676013524770164`*^9, 3.5676014283711495`*^9}, {3.567604737731021*^9, 3.5676047599610596`*^9}, { 3.567611568287159*^9, 3.5676115840743866`*^9}, {3.5676124995619946`*^9, 3.5676125004199963`*^9}, {3.5676147541247544`*^9, 3.567614755169956*^9}, { 3.567615587415818*^9, 3.5676155981798368`*^9}, {3.5676156384067073`*^9, 3.5676156583747425`*^9}, {3.5676158998007665`*^9, 3.5676159408288383`*^9}, {3.567615989126523*^9, 3.567616164252431*^9}, { 3.5676162788658323`*^9, 3.567616283686241*^9}, {3.567616477313781*^9, 3.5676165242074633`*^9}, {3.5676165629735312`*^9, 3.5676165647363343`*^9}, {3.5676166087284117`*^9, 3.5676166312860513`*^9}, {3.5676169704150467`*^9, 3.567617073671628*^9}, { 3.567617106774886*^9, 3.5676171121256957`*^9}, 3.567617182341419*^9, { 3.567617294022015*^9, 3.567617336079689*^9}, {3.567660899793217*^9, 3.5676609125384398`*^9}, {3.5676609439100943`*^9, 3.567660949073704*^9}, { 3.5676610469482756`*^9, 3.567661068710314*^9}, 3.567661166756486*^9, { 3.5676850168515854`*^9, 3.5676850349164176`*^9}, {3.5676852124135294`*^9, 3.567685248761593*^9}, {3.5676856839243574`*^9, 3.5676856880427647`*^9}, { 3.5676858185057936`*^9, 3.567685824714605*^9}, 3.567686114506314*^9, { 3.567686148249173*^9, 3.5676861890588446`*^9}, {3.567686447551299*^9, 3.5676864707953396`*^9}, {3.567686506925003*^9, 3.567686565409506*^9}, { 3.567686809066334*^9, 3.567686826366764*^9}, {3.5676874811935143`*^9, 3.5676874818643155`*^9}, {3.5676921395598993`*^9, 3.5676921509791193`*^9}, {3.567692753654978*^9, 3.5676927950730505`*^9}, { 3.567698206184089*^9, 3.5676982494429655`*^9}, 3.5678510438419695`*^9, { 3.567851238998312*^9, 3.5678512624919534`*^9}, {3.5678515219984093`*^9, 3.567851522840811*^9}, {3.5678524805576925`*^9, 3.567852491056511*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"-", "78.03749588279042`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Cr", "\[Rule]", "0.2625736948837923`"}], ",", RowBox[{"CC", "\[Rule]", "0.9949332760919413`"}], ",", RowBox[{"CH", "\[Rule]", "1.21900954961684`"}], ",", RowBox[{"Cpi", "\[Rule]", "2.11331870905448`"}], ",", RowBox[{"w", "\[Rule]", "1.0232027350128836`"}], ",", RowBox[{"p", "\[Rule]", "0.5739955948290986`"}]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{ 3.5678519660531893`*^9, 3.567852018906082*^9, 3.567852050184137*^9, { 3.567852087936203*^9, 3.5678521163126535`*^9}, {3.567852155141121*^9, 3.567852175733157*^9}, 3.5678522190700336`*^9, 3.5678523363042393`*^9, { 3.5678523915439367`*^9, 3.567852415942379*^9}, 3.5678524507148404`*^9, 3.567852497639723*^9, {3.5678525612254343`*^9, 3.5678525819422708`*^9}, 3.5678527263829246`*^9, 3.5678527963334475`*^9, 3.567852832899912*^9, 3.567852992597392*^9, 3.567853042205479*^9, 3.567853143933258*^9, 3.5678815573017316`*^9, {3.5678815961145997`*^9, 3.567881636674671*^9}, 3.567881672071133*^9, {3.5679215716186323`*^9, 3.567921596407076*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{"Ekin", "/.", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}], "\n", RowBox[{"vne", "/.", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}], "\n", RowBox[{"vee", "/.", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}], "\n", RowBox[{"vnn", "/.", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}], "\n", RowBox[{ RowBox[{ RowBox[{"-", "Epot"}], "/", "Ekin"}], "/.", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}], "\n", RowBox[{ RowBox[{"2", "*", "w", "*", "rad"}], "/.", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}], "\n", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"2", "*", "Cr"}], "+", RowBox[{"2", "*", "CC"}]}], ")"}], "*", "bohr"}], "/.", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}], "\n", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"Cr", "+", "CH", "+", "p"}], ")"}], "*", "bohr"}], "/.", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}]}], "Input", CellChangeTimes->{{3.5675987354944196`*^9, 3.567599332289068*^9}, { 3.567599364347124*^9, 3.567599388854767*^9}, {3.567599446824469*^9, 3.567599671277663*^9}, {3.5675997743158445`*^9, 3.567599879038828*^9}, { 3.567599920722101*^9, 3.567600053883935*^9}, {3.5676001118380375`*^9, 3.5676001705877404`*^9}, {3.5676002126766143`*^9, 3.5676003454796476`*^9}, {3.5676004083477583`*^9, 3.567600442527418*^9}, { 3.5676004745698743`*^9, 3.567600484132691*^9}, {3.5676005145995445`*^9, 3.5676005320091753`*^9}, {3.5676005791056576`*^9, 3.5676012559440465`*^9}, {3.5676013203097596`*^9, 3.5676013219321623`*^9}, {3.5676013524770164`*^9, 3.5676014283711495`*^9}, {3.5676166536252904`*^9, 3.567616654654892*^9}, 3.567685840720233*^9}], Cell[BoxData["78.03749588279943`"], "Output", CellChangeTimes->{ 3.5678519660531893`*^9, 3.567852018921682*^9, 3.5678520501997366`*^9, { 3.567852087951803*^9, 3.5678521163282533`*^9}, {3.567852155156721*^9, 3.5678521757487574`*^9}, 3.5678522190700336`*^9, 3.5678523363198395`*^9, { 3.5678523915595365`*^9, 3.567852415957979*^9}, 3.56785245073044*^9, 3.567852497655323*^9, {3.5678525612410345`*^9, 3.5678525819578705`*^9}, 3.5678527263985243`*^9, 3.5678527963490477`*^9, 3.567852832899912*^9, 3.5678529926129923`*^9, 3.567853042205479*^9, 3.5678531439488583`*^9, 3.5678815573173313`*^9, {3.5678815961302*^9, 3.5678816366902714`*^9}, 3.567881672086733*^9, {3.567921571634232*^9, 3.5679215964226756`*^9}}], Cell[BoxData[ RowBox[{"-", "246.54328090904312`"}]], "Output", CellChangeTimes->{ 3.5678519660531893`*^9, 3.567852018921682*^9, 3.5678520501997366`*^9, { 3.567852087951803*^9, 3.5678521163282533`*^9}, {3.567852155156721*^9, 3.5678521757487574`*^9}, 3.5678522190700336`*^9, 3.5678523363198395`*^9, { 3.5678523915595365`*^9, 3.567852415957979*^9}, 3.56785245073044*^9, 3.567852497655323*^9, {3.5678525612410345`*^9, 3.5678525819578705`*^9}, 3.5678527263985243`*^9, 3.5678527963490477`*^9, 3.567852832899912*^9, 3.5678529926129923`*^9, 3.567853042205479*^9, 3.5678531439488583`*^9, 3.5678815573173313`*^9, {3.5678815961302*^9, 3.5678816366902714`*^9}, 3.567881672086733*^9, {3.567921571634232*^9, 3.5679215964226756`*^9}}], Cell[BoxData["57.12383513233799`"], "Output", CellChangeTimes->{ 3.5678519660531893`*^9, 3.567852018921682*^9, 3.5678520501997366`*^9, { 3.567852087951803*^9, 3.5678521163282533`*^9}, {3.567852155156721*^9, 3.5678521757487574`*^9}, 3.5678522190700336`*^9, 3.5678523363198395`*^9, { 3.5678523915595365`*^9, 3.567852415957979*^9}, 3.56785245073044*^9, 3.567852497655323*^9, {3.5678525612410345`*^9, 3.5678525819578705`*^9}, 3.5678527263985243`*^9, 3.5678527963490477`*^9, 3.567852832899912*^9, 3.5678529926129923`*^9, 3.567853042205479*^9, 3.5678531439488583`*^9, 3.5678815573173313`*^9, {3.5678815961302*^9, 3.5678816366902714`*^9}, 3.567881672086733*^9, {3.567921571634232*^9, 3.5679215964226756`*^9}}], Cell[BoxData["33.344454011115246`"], "Output", CellChangeTimes->{ 3.5678519660531893`*^9, 3.567852018921682*^9, 3.5678520501997366`*^9, { 3.567852087951803*^9, 3.5678521163282533`*^9}, {3.567852155156721*^9, 3.5678521757487574`*^9}, 3.5678522190700336`*^9, 3.5678523363198395`*^9, { 3.5678523915595365`*^9, 3.567852415957979*^9}, 3.56785245073044*^9, 3.567852497655323*^9, {3.5678525612410345`*^9, 3.5678525819578705`*^9}, 3.5678527263985243`*^9, 3.5678527963490477`*^9, 3.567852832899912*^9, 3.5678529926129923`*^9, 3.567853042205479*^9, 3.5678531439488583`*^9, 3.5678815573173313`*^9, {3.5678815961302*^9, 3.5678816366902714`*^9}, 3.567881672086733*^9, {3.567921571634232*^9, 3.5679215964226756`*^9}}], Cell[BoxData["1.999999999999885`"], "Output", CellChangeTimes->{ 3.5678519660531893`*^9, 3.567852018921682*^9, 3.5678520501997366`*^9, { 3.567852087951803*^9, 3.5678521163282533`*^9}, {3.567852155156721*^9, 3.5678521757487574`*^9}, 3.5678522190700336`*^9, 3.5678523363198395`*^9, { 3.5678523915595365`*^9, 3.567852415957979*^9}, 3.56785245073044*^9, 3.567852497655323*^9, {3.5678525612410345`*^9, 3.5678525819578705`*^9}, 3.5678527263985243`*^9, 3.5678527963490477`*^9, 3.567852832899912*^9, 3.5678529926129923`*^9, 3.567853042205479*^9, 3.5678531439488583`*^9, 3.5678815573173313`*^9, {3.5678815961302*^9, 3.5678816366902714`*^9}, 3.567881672086733*^9, {3.567921571634232*^9, 3.5679215964226756`*^9}}], Cell[BoxData["117.25039760139296`"], "Output", CellChangeTimes->{ 3.5678519660531893`*^9, 3.567852018921682*^9, 3.5678520501997366`*^9, { 3.567852087951803*^9, 3.5678521163282533`*^9}, {3.567852155156721*^9, 3.5678521757487574`*^9}, 3.5678522190700336`*^9, 3.5678523363198395`*^9, { 3.5678523915595365`*^9, 3.567852415957979*^9}, 3.56785245073044*^9, 3.567852497655323*^9, {3.5678525612410345`*^9, 3.5678525819578705`*^9}, 3.5678527263985243`*^9, 3.5678527963490477`*^9, 3.567852832899912*^9, 3.5678529926129923`*^9, 3.567853042205479*^9, 3.5678531439488583`*^9, 3.5678815573173313`*^9, {3.5678815961302*^9, 3.5678816366902714`*^9}, 3.567881672086733*^9, {3.567921571634232*^9, 3.5679215964226756`*^9}}], Cell[BoxData["1.3308875327600516`"], "Output", CellChangeTimes->{ 3.5678519660531893`*^9, 3.567852018921682*^9, 3.5678520501997366`*^9, { 3.567852087951803*^9, 3.5678521163282533`*^9}, {3.567852155156721*^9, 3.5678521757487574`*^9}, 3.5678522190700336`*^9, 3.5678523363198395`*^9, { 3.5678523915595365`*^9, 3.567852415957979*^9}, 3.56785245073044*^9, 3.567852497655323*^9, {3.5678525612410345`*^9, 3.5678525819578705`*^9}, 3.5678527263985243`*^9, 3.5678527963490477`*^9, 3.567852832899912*^9, 3.5678529926129923`*^9, 3.567853042205479*^9, 3.5678531439488583`*^9, 3.5678815573173313`*^9, {3.5678815961302*^9, 3.5678816366902714`*^9}, 3.567881672086733*^9, {3.567921571634232*^9, 3.5679215964226756`*^9}}], Cell[BoxData["1.087765043459989`"], "Output", CellChangeTimes->{ 3.5678519660531893`*^9, 3.567852018921682*^9, 3.5678520501997366`*^9, { 3.567852087951803*^9, 3.5678521163282533`*^9}, {3.567852155156721*^9, 3.5678521757487574`*^9}, 3.5678522190700336`*^9, 3.5678523363198395`*^9, { 3.5678523915595365`*^9, 3.567852415957979*^9}, 3.56785245073044*^9, 3.567852497655323*^9, {3.5678525612410345`*^9, 3.5678525819578705`*^9}, 3.5678527263985243`*^9, 3.5678527963490477`*^9, 3.567852832899912*^9, 3.5678529926129923`*^9, 3.567853042205479*^9, 3.5678531439488583`*^9, 3.5678815573173313`*^9, {3.5678815961302*^9, 3.5678816366902714`*^9}, 3.567881672086733*^9, {3.567921571634232*^9, 3.5679215964226756`*^9}}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"plot1", "=", RowBox[{ RowBox[{"Graphics", "[", RowBox[{"{", RowBox[{ RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"yc", "[", RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", "Cr"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"yc", "[", RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", "Cr"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "3", "]"}], "]"}], ",", RowBox[{"yc", "[", RowBox[{"[", "3", "]"}], "]"}]}], "}"}], ",", "CH"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "4", "]"}], "]"}], ",", RowBox[{"yc", "[", RowBox[{"[", "4", "]"}], "]"}]}], "}"}], ",", "CH"}], "]"}], ",", "\n", " ", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "5", "]"}], "]"}], ",", RowBox[{"yc", "[", RowBox[{"[", "5", "]"}], "]"}]}], "}"}], ",", "CH"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "6", "]"}], "]"}], ",", RowBox[{"yc", "[", RowBox[{"[", "6", "]"}], "]"}]}], "}"}], ",", "CH"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "7", "]"}], "]"}], ",", RowBox[{"yc", "[", RowBox[{"[", "7", "]"}], "]"}]}], "}"}], ",", "CC"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "8", "]"}], "]"}], ",", RowBox[{"yc", "[", RowBox[{"[", "8", "]"}], "]"}]}], "}"}], ",", "Cpi"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "9", "]"}], "]"}], ",", RowBox[{"yc", "[", RowBox[{"[", "9", "]"}], "]"}]}], "}"}], ",", "Cpi"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", " ", "\[IndentingNewLine]", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", " ", "\[IndentingNewLine]", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "3", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "3", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "4", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "4", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", " ", "\[IndentingNewLine]", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "5", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "5", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "6", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "6", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", "\n", RowBox[{"{", RowBox[{ RowBox[{"Thickness", "[", "0.008", "]"}], ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{ RowBox[{"yn", "[", RowBox[{"[", "2", "]"}], "]"}], "+", "0.06"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{ RowBox[{"yn", "[", RowBox[{"[", "1", "]"}], "]"}], "+", "0.06"}]}], "}"}]}], "}"}], "]"}], ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{ RowBox[{"yn", "[", RowBox[{"[", "2", "]"}], "]"}], "-", "0.06"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{ RowBox[{"yn", "[", RowBox[{"[", "1", "]"}], "]"}], "-", "0.06"}]}], "}"}]}], "}"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "3", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "3", "]"}], "]"}]}], "}"}]}], "}"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "4", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "4", "]"}], "]"}]}], "}"}]}], "}"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "5", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "5", "]"}], "]"}]}], "}"}]}], "}"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "6", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "6", "]"}], "]"}]}], "}"}]}], "}"}], "]"}]}], "}"}]}], "}"}], "]"}], " ", "/.", " ", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], "\n"}]], "Input", CellChangeTimes->{{3.567661245973425*^9, 3.567661380742062*^9}, { 3.567661415998124*^9, 3.5676614516285863`*^9}, {3.5676614877114496`*^9, 3.567661525089115*^9}, {3.5676615670531893`*^9, 3.5676615719827976`*^9}, { 3.5676616025120516`*^9, 3.567661696112216*^9}, {3.567663489803366*^9, 3.5676635070881968`*^9}, {3.567663544450262*^9, 3.567663576336718*^9}, { 3.567663626366006*^9, 3.567663629704412*^9}, {3.567692320255017*^9, 3.567692403044362*^9}, 3.5676924445872355`*^9, {3.567692490482516*^9, 3.56769258937109*^9}, {3.5676926267487555`*^9, 3.567692699054882*^9}, 3.5678517700076447`*^9, {3.567852103598631*^9, 3.5678521098542414`*^9}, 3.567852718614111*^9}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{"xy", "-", "plane"}], " ", "*)"}], "\n", RowBox[{"Show", "[", RowBox[{"plot1", ",", RowBox[{"{", RowBox[{ RowBox[{"AspectRatio", " ", "\[Rule]", " ", "Automatic"}], ",", RowBox[{"Axes", " ", "->", " ", "True"}], ",", RowBox[{"GridLines", " ", "->", " ", "Automatic"}], ",", " ", RowBox[{"PlotRange", " ", "\[Rule]", " ", "Automatic"}], ",", " ", RowBox[{"Frame", " ", "->", " ", "True"}]}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.567661245973425*^9, 3.567661380742062*^9}, { 3.567661415998124*^9, 3.5676614516285863`*^9}, {3.5676614877114496`*^9, 3.567661525089115*^9}, {3.5676615670531893`*^9, 3.5676615719827976`*^9}, { 3.5676616025120516`*^9, 3.567661696112216*^9}, {3.567663489803366*^9, 3.5676635070881968`*^9}, {3.567663544450262*^9, 3.567663576336718*^9}, { 3.567663626366006*^9, 3.567663629704412*^9}, {3.567692320255017*^9, 3.567692403044362*^9}, 3.5676924445872355`*^9}], Cell[BoxData[ GraphicsBox[{CircleBox[{1.2575069709757336`, 0}, 0.2625736948837923], CircleBox[{-1.2575069709757336`, 0}, 0.2625736948837923], CircleBox[{2.0288698715525646`, 1.2649459221637571`}, 1.21900954961684], CircleBox[{2.0288698715525646`, -1.2649459221637571`}, 1.21900954961684], CircleBox[{-2.0288698715525646`, 1.2649459221637571`}, 1.21900954961684], CircleBox[{-2.0288698715525646`, -1.2649459221637571`}, 1.21900954961684], CircleBox[{0, 0}, 0.9949332760919413], CircleBox[{0, 0}, 2.11331870905448], CircleBox[{0, 0}, 2.11331870905448], DiskBox[{1.2575069709757336`, 0}, 0.08], DiskBox[{-1.2575069709757336`, 0}, 0.08], DiskBox[{2.3277116062908165`, 1.7550117957581575`}, 0.08], DiskBox[{2.3277116062908165`, -1.7550117957581575`}, 0.08], DiskBox[{-2.3277116062908165`, 1.7550117957581575`}, 0.08], DiskBox[{-2.3277116062908165`, -1.7550117957581575`}, 0.08], {Thickness[0.008], LineBox[{{-1.2575069709757336`, 0.06}, {1.2575069709757336`, 0.06}}], LineBox[{{-1.2575069709757336`, -0.06}, {1.2575069709757336`, -0.06}}], LineBox[{{1.2575069709757336`, 0}, {2.3277116062908165`, 1.7550117957581575`}}], LineBox[{{1.2575069709757336`, 0}, { 2.3277116062908165`, -1.7550117957581575`}}], LineBox[{{-1.2575069709757336`, 0}, {-2.3277116062908165`, 1.7550117957581575`}}], LineBox[{{-1.2575069709757336`, 0}, {-2.3277116062908165`, -1.7550117957581575`}}]}}, AspectRatio->Automatic, Axes->True, Frame->True, GridLines->Automatic, PlotRange->Automatic]], "Output", CellChangeTimes->{ 3.567851966084389*^9, 3.567852018952882*^9, 3.567852050215337*^9, { 3.567852087983003*^9, 3.567852116359453*^9}, {3.5678521551723213`*^9, 3.5678521757643576`*^9}, 3.5678522191012335`*^9, 3.5678523363354397`*^9, { 3.5678523915907364`*^9, 3.5678524159891796`*^9}, 3.5678524507616405`*^9, 3.5678524976865225`*^9, {3.567852561272235*^9, 3.567852581989071*^9}, 3.5678527264297247`*^9, 3.5678527963802476`*^9, 3.5678528329311113`*^9, 3.5678529926441927`*^9, 3.567853042236679*^9, 3.5678531439800577`*^9, 3.5678815573485317`*^9, {3.5678815961614*^9, 3.567881636705871*^9}, 3.5678816721179333`*^9, {3.567921571681032*^9, 3.5679215964538755`*^9}}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"plot2", "=", RowBox[{ RowBox[{"Graphics", "[", RowBox[{"{", RowBox[{ RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"zc", "[", RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", "Cr"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"zc", "[", RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", "Cr"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "3", "]"}], "]"}], ",", RowBox[{"zc", "[", RowBox[{"[", "3", "]"}], "]"}]}], "}"}], ",", "CH"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "4", "]"}], "]"}], ",", RowBox[{"zc", "[", RowBox[{"[", "4", "]"}], "]"}]}], "}"}], ",", "CH"}], "]"}], ",", "\n", " ", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "5", "]"}], "]"}], ",", RowBox[{"zc", "[", RowBox[{"[", "5", "]"}], "]"}]}], "}"}], ",", "CH"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "6", "]"}], "]"}], ",", RowBox[{"zc", "[", RowBox[{"[", "6", "]"}], "]"}]}], "}"}], ",", "CH"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "7", "]"}], "]"}], ",", RowBox[{"zc", "[", RowBox[{"[", "7", "]"}], "]"}]}], "}"}], ",", "CC"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "8", "]"}], "]"}], ",", RowBox[{"zc", "[", RowBox[{"[", "8", "]"}], "]"}]}], "}"}], ",", "Cpi"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "9", "]"}], "]"}], ",", RowBox[{"zc", "[", RowBox[{"[", "9", "]"}], "]"}]}], "}"}], ",", "Cpi"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", " ", "\[IndentingNewLine]", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", " ", "\[IndentingNewLine]", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "3", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "3", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "4", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "4", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", " ", "\[IndentingNewLine]", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "5", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "5", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "6", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "6", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", "\n", RowBox[{"{", RowBox[{ RowBox[{"Thickness", "[", "0.008", "]"}], ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "1", "]"}], "]"}]}], "}"}]}], "}"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "3", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "3", "]"}], "]"}]}], "}"}]}], "}"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "4", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "4", "]"}], "]"}]}], "}"}]}], "}"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "5", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "5", "]"}], "]"}]}], "}"}]}], "}"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "6", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "6", "]"}], "]"}]}], "}"}]}], "}"}], "]"}]}], "}"}]}], "}"}], "]"}], " ", "/.", " ", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], "\n"}]], "Input", CellChangeTimes->{{3.567851861423805*^9, 3.56785193220113*^9}, { 3.567852075705782*^9, 3.5678520816805925`*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Show", "[", RowBox[{"plot2", ",", RowBox[{"{", RowBox[{ RowBox[{"AspectRatio", " ", "\[Rule]", " ", "Automatic"}], ",", RowBox[{"Axes", " ", "->", " ", "True"}], ",", RowBox[{"GridLines", " ", "->", " ", "Automatic"}], ",", " ", RowBox[{"PlotRange", " ", "\[Rule]", " ", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "3.5"}], ",", "3.5"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "4.5"}], ",", "4.5"}], "}"}]}], "}"}]}], ",", " ", RowBox[{"Frame", " ", "->", " ", "True"}]}], "}"}]}], "]"}]], "Input", CellChangeTimes->{{3.5678519457107534`*^9, 3.5678519464595547`*^9}, { 3.5678519871444263`*^9, 3.567852035878912*^9}}], Cell[BoxData[ GraphicsBox[{CircleBox[{1.2575069709757336`, 0}, 0.2625736948837923], CircleBox[{-1.2575069709757336`, 0}, 0.2625736948837923], CircleBox[{2.0288698715525646`, 0}, 1.21900954961684], CircleBox[{2.0288698715525646`, 0}, 1.21900954961684], CircleBox[{-2.0288698715525646`, 0}, 1.21900954961684], CircleBox[{-2.0288698715525646`, 0}, 1.21900954961684], CircleBox[{0, 0}, 0.9949332760919413], CircleBox[{0, 2.3758924039382725`}, 2.11331870905448], CircleBox[{0, -2.3758924039382725`}, 2.11331870905448], DiskBox[{1.2575069709757336`, 0}, 0.08], DiskBox[{-1.2575069709757336`, 0}, 0.08], DiskBox[{2.3277116062908165`, 0}, 0.08], DiskBox[{2.3277116062908165`, 0}, 0.08], DiskBox[{-2.3277116062908165`, 0}, 0.08], DiskBox[{-2.3277116062908165`, 0}, 0.08], {Thickness[0.008], LineBox[{{-1.2575069709757336`, 0}, {1.2575069709757336`, 0}}], LineBox[{{1.2575069709757336`, 0}, {2.3277116062908165`, 0}}], LineBox[{{1.2575069709757336`, 0}, {2.3277116062908165`, 0}}], LineBox[{{-1.2575069709757336`, 0}, {-2.3277116062908165`, 0}}], LineBox[{{-1.2575069709757336`, 0}, {-2.3277116062908165`, 0}}]}}, AspectRatio->Automatic, Axes->True, Frame->True, GridLines->Automatic, PlotRange->{{-3.5, 3.5}, {-4.5, 4.5}}]], "Output", CellChangeTimes->{ 3.567851966099989*^9, 3.567852018968482*^9, 3.567852050246537*^9, { 3.567852087998603*^9, 3.5678521163750534`*^9}, {3.5678521552035217`*^9, 3.5678521757955575`*^9}, 3.5678522191168337`*^9, 3.5678523363666396`*^9, { 3.5678523916063366`*^9, 3.5678524160047793`*^9}, 3.5678524507772408`*^9, 3.567852497717723*^9, {3.5678525612878346`*^9, 3.567852582004671*^9}, 3.567852726460925*^9, 3.5678527964114475`*^9, 3.5678528329623117`*^9, 3.567852992675392*^9, 3.5678530422678795`*^9, 3.567853143995658*^9, 3.567881557364132*^9, {3.567881596177*^9, 3.567881636737071*^9}, 3.567881672149133*^9, {3.5679215716966324`*^9, 3.5679215964694757`*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"plot3", "=", RowBox[{ RowBox[{"Graphics3D", "[", RowBox[{"{", "\n", RowBox[{ RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"yc", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"zc", "[", RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", "Cr"}], "]"}], ",", " ", "\n", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"yc", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"zc", "[", RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", "Cr"}], "]"}], ",", "\n", RowBox[{"{", RowBox[{ RowBox[{"Opacity", "[", "0.6", "]"}], ",", "\n", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "3", "]"}], "]"}], ",", RowBox[{"yc", "[", RowBox[{"[", "3", "]"}], "]"}], ",", RowBox[{"zc", "[", RowBox[{"[", "3", "]"}], "]"}]}], "}"}], ",", "CH"}], "]"}], ",", "\n", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "4", "]"}], "]"}], ",", RowBox[{"yc", "[", RowBox[{"[", "4", "]"}], "]"}], ",", RowBox[{"zc", "[", RowBox[{"[", "4", "]"}], "]"}]}], "}"}], ",", "CH"}], "]"}], ",", "\n", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "5", "]"}], "]"}], ",", RowBox[{"yc", "[", RowBox[{"[", "5", "]"}], "]"}], ",", RowBox[{"zc", "[", RowBox[{"[", "5", "]"}], "]"}]}], "}"}], ",", "CH"}], "]"}], ",", "\n", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "6", "]"}], "]"}], ",", RowBox[{"yc", "[", RowBox[{"[", "6", "]"}], "]"}], ",", RowBox[{"zc", "[", RowBox[{"[", "6", "]"}], "]"}]}], "}"}], ",", "CH"}], "]"}], ",", "\n", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "7", "]"}], "]"}], ",", RowBox[{"yc", "[", RowBox[{"[", "7", "]"}], "]"}], ",", RowBox[{"zc", "[", RowBox[{"[", "7", "]"}], "]"}]}], "}"}], ",", "CC"}], "]"}], ",", "\n", RowBox[{"{", RowBox[{ RowBox[{"Opacity", "[", "0.4", "]"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Lighter", "[", "Green", "]"}], ",", "\n", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "8", "]"}], "]"}], ",", RowBox[{"yc", "[", RowBox[{"[", "8", "]"}], "]"}], ",", RowBox[{"zc", "[", RowBox[{"[", "8", "]"}], "]"}]}], "}"}], ",", "Cpi"}], "]"}], ",", "\n", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "9", "]"}], "]"}], ",", RowBox[{"yc", "[", RowBox[{"[", "9", "]"}], "]"}], ",", RowBox[{"zc", "[", RowBox[{"[", "9", "]"}], "]"}]}], "}"}], ",", "Cpi"}], "]"}]}], "}"}]}], "}"}], ",", "\n", RowBox[{"{", RowBox[{ RowBox[{"Darker", "[", RowBox[{"Red", ",", "1"}], "]"}], ",", "\n", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", "\n", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", "\n", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "3", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "3", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "3", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", "\n", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "4", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "4", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "4", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", "\n", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "5", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "5", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "5", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", "\n", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "6", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "6", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "6", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", "\n", RowBox[{"Cylinder", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{ RowBox[{"yn", "[", RowBox[{"[", "2", "]"}], "]"}], "-", "0.1"}], ",", RowBox[{"zn", "[", RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{ RowBox[{"yn", "[", RowBox[{"[", "1", "]"}], "]"}], "-", "0.1"}], ",", RowBox[{"zn", "[", RowBox[{"[", "1", "]"}], "]"}]}], "}"}]}], "}"}], ",", "0.05"}], "]"}], ",", "\n", RowBox[{"Cylinder", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{ RowBox[{"yn", "[", RowBox[{"[", "2", "]"}], "]"}], "+", "0.1"}], ",", RowBox[{"zn", "[", RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{ RowBox[{"yn", "[", RowBox[{"[", "1", "]"}], "]"}], "+", "0.1"}], ",", RowBox[{"zn", "[", RowBox[{"[", "1", "]"}], "]"}]}], "}"}]}], "}"}], ",", "0.05"}], "]"}], ",", "\n", RowBox[{"Cylinder", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "3", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "3", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "3", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "1", "]"}], "]"}]}], "}"}]}], "}"}], ",", "0.05"}], "]"}], ",", "\n", RowBox[{"Cylinder", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "4", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "4", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "4", "]"}], "]"}]}], "}"}]}], "}"}], ",", "0.05"}], "]"}], ",", "\n", RowBox[{"Cylinder", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "5", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "5", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "5", "]"}], "]"}]}], "}"}]}], "}"}], ",", "0.05"}], "]"}], ",", "\n", RowBox[{"Cylinder", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "6", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "6", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "6", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "2", "]"}], "]"}]}], "}"}]}], "}"}], ",", "0.05"}], "]"}]}], "}"}]}], "}"}]}], "}"}], "]"}], " ", "/.", " ", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], "\n"}], "\n", RowBox[{"Show", "[", RowBox[{"plot3", ",", RowBox[{"{", RowBox[{ RowBox[{"AspectRatio", " ", "\[Rule]", " ", "Automatic"}], ",", RowBox[{"Boxed", "->", "False"}]}], "}"}]}], "]"}]}], "Input", CellChangeTimes->{{3.5678527653517933`*^9, 3.5678527897346363`*^9}, { 3.5678528275803022`*^9, 3.5678528279547033`*^9}}], Cell[BoxData[ Graphics3DBox[{SphereBox[{1.2575069709757336`, 0, 0}, 0.2625736948837923], SphereBox[{-1.2575069709757336`, 0, 0}, 0.2625736948837923], {Opacity[0.6], SphereBox[{2.0288698715525646`, 1.2649459221637571`, 0}, 1.21900954961684], SphereBox[{2.0288698715525646`, -1.2649459221637571`, 0}, 1.21900954961684], SphereBox[{-2.0288698715525646`, 1.2649459221637571`, 0}, 1.21900954961684], SphereBox[{-2.0288698715525646`, -1.2649459221637571`, 0}, 1.21900954961684], SphereBox[{0, 0, 0}, 0.9949332760919413], {Opacity[0.4], {RGBColor[ NCache[ Rational[1, 3], 0.3333333333333333], 1, NCache[ Rational[1, 3], 0.3333333333333333]], SphereBox[{0, 0, 2.3758924039382725`}, 2.11331870905448], SphereBox[{0, 0, -2.3758924039382725`}, 2.11331870905448]}}, {RGBColor[0, 0, 0], SphereBox[{1.2575069709757336`, 0, 0}, 0.08], SphereBox[{-1.2575069709757336`, 0, 0}, 0.08], SphereBox[{2.3277116062908165`, 1.7550117957581575`, 0}, 0.08], SphereBox[{2.3277116062908165`, -1.7550117957581575`, 0}, 0.08], SphereBox[{-2.3277116062908165`, 1.7550117957581575`, 0}, 0.08], SphereBox[{-2.3277116062908165`, -1.7550117957581575`, 0}, 0.08], CylinderBox[{{-1.2575069709757336`, -0.1, 0}, { 1.2575069709757336`, -0.1, 0}}, 0.05], CylinderBox[{{-1.2575069709757336`, 0.1, 0}, { 1.2575069709757336`, 0.1, 0}}, 0.05], CylinderBox[{{2.3277116062908165`, 1.7550117957581575`, 0}, { 1.2575069709757336`, 0, 0}}, 0.05], CylinderBox[{{1.2575069709757336`, 0, 0}, { 2.3277116062908165`, -1.7550117957581575`, 0}}, 0.05], CylinderBox[{{-1.2575069709757336`, 0, 0}, {-2.3277116062908165`, 1.7550117957581575`, 0}}, 0.05], CylinderBox[{{-2.3277116062908165`, -1.7550117957581575`, 0}, {-1.2575069709757336`, 0, 0}}, 0.05]}}}, AspectRatio->Automatic, Boxed->False]], "Output", CellChangeTimes->{ 3.5678528332275124`*^9, 3.5678529926909924`*^9, 3.5678530422834797`*^9, 3.5678531440268583`*^9, 3.567881557395332*^9, {3.5678815962082*^9, 3.5678816367526712`*^9}, 3.5678816721647334`*^9, {3.567921572024233*^9, 3.567921596485076*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{"as", ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"x", "=", RowBox[{ RowBox[{"bohr", "*", "xn"}], " ", "/.", " ", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"y", "=", RowBox[{ RowBox[{"bohr", "*", "yn"}], " ", "/.", " ", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"z", "=", RowBox[{ RowBox[{"bohr", "*", "zn"}], " ", "/.", " ", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"s", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"as", "[", RowBox[{"[", RowBox[{"l", "+", "1"}], "]"}], "]"}], ",", RowBox[{"x", "[", RowBox[{"[", "l", "]"}], "]"}], ",", RowBox[{"y", "[", RowBox[{"[", "l", "]"}], "]"}], ",", RowBox[{"z", "[", RowBox[{"[", "l", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"l", ",", "6"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", StyleBox[ RowBox[{"(*", " ", RowBox[{ "Cartesian", " ", "Coordinates", " ", "of", " ", "stationary", " ", "structure"}], " ", "*)"}], "Text"]}], "\[IndentingNewLine]", RowBox[{"u", "=", RowBox[{"MatrixForm", "[", "s", "]"}]}]}], "Input", CellChangeTimes->{3.567852914316455*^9}], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"C", "0.6654437663800258`", "0.`", "0.`"}, {"C", RowBox[{"-", "0.6654437663800258`"}], "0.`", "0.`"}, {"H", "1.2317714446821555`", "0.9287118770439146`", "0.`"}, {"H", "1.2317714446821555`", RowBox[{"-", "0.9287118770439146`"}], "0.`"}, {"H", RowBox[{"-", "1.2317714446821555`"}], "0.9287118770439146`", "0.`"}, {"H", RowBox[{"-", "1.2317714446821555`"}], RowBox[{"-", "0.9287118770439146`"}], "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.567852992893793*^9, 3.5678530424706798`*^9, 3.5678531440268583`*^9, 3.567881557395332*^9, {3.5678815962082*^9, 3.5678816367526712`*^9}, 3.5678816721647334`*^9, {3.567921572055433*^9, 3.567921596500676*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"(*", " ", RowBox[{"Distance", " ", "matrix", " ", "\[CapitalARing]"}], " ", "*)"}], "Text"], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Dij", "=", RowBox[{"Table", "[", RowBox[{ 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[{"i", ",", "6"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "6"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"MatrixForm", "[", "Dij", "]"}]}]}]], "Input"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"0.`", "1.3308875327600516`", "1.0877650434599888`", "1.0877650434599888`", "2.1123284090425303`", "2.1123284090425303`"}, {"1.3308875327600516`", "0.`", "2.1123284090425303`", "2.1123284090425303`", "1.0877650434599888`", "1.0877650434599888`"}, {"1.0877650434599888`", "2.1123284090425303`", "0.`", "1.8574237540878291`", "2.463542889364311`", "3.0852984571978093`"}, {"1.0877650434599888`", "2.1123284090425303`", "1.8574237540878291`", "0.`", "3.0852984571978093`", "2.463542889364311`"}, {"2.1123284090425303`", "1.0877650434599888`", "2.463542889364311`", "3.0852984571978093`", "0.`", "1.8574237540878291`"}, {"2.1123284090425303`", "1.0877650434599888`", "3.0852984571978093`", "2.463542889364311`", "1.8574237540878291`", "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.567852992909393*^9, 3.5678530424862795`*^9, 3.5678531440424585`*^9, 3.5678815574109316`*^9, {3.5678815962237997`*^9, 3.5678816367682714`*^9}, 3.567881672180333*^9, {3.5679215720710325`*^9, 3.567921596516276*^9}}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"plot4", "=", RowBox[{ RowBox[{"Graphics3D", "[", RowBox[{"{", "\n", RowBox[{ RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"yc", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"zc", "[", RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", "Cr"}], "]"}], ",", " ", "\n", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"yc", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"zc", "[", RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", "Cr"}], "]"}], ",", "\n", RowBox[{"{", RowBox[{ RowBox[{"Opacity", "[", "0.6", "]"}], ",", "\n", RowBox[{"{", RowBox[{ RowBox[{"Opacity", "[", "0.4", "]"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Lighter", "[", "Green", "]"}], ",", "\n", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "8", "]"}], "]"}], ",", RowBox[{"yc", "[", RowBox[{"[", "8", "]"}], "]"}], ",", RowBox[{"zc", "[", RowBox[{"[", "8", "]"}], "]"}]}], "}"}], ",", "Cpi"}], "]"}], ",", "\n", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "9", "]"}], "]"}], ",", RowBox[{"yc", "[", RowBox[{"[", "9", "]"}], "]"}], ",", RowBox[{"zc", "[", RowBox[{"[", "9", "]"}], "]"}]}], "}"}], ",", "Cpi"}], "]"}]}], "}"}]}], "}"}], ",", "\n", RowBox[{"{", RowBox[{ RowBox[{"Darker", "[", RowBox[{"Red", ",", "1"}], "]"}], ",", "\n", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", "\n", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", "\n", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "3", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "3", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "3", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", "\n", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "4", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "4", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "4", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", "\n", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "5", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "5", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "5", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", "\n", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "6", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "6", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "6", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", "\n", RowBox[{"Cylinder", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "3", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "3", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "3", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "1", "]"}], "]"}]}], "}"}]}], "}"}], ",", "0.05"}], "]"}], ",", "\n", RowBox[{"Cylinder", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "4", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "4", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "4", "]"}], "]"}]}], "}"}]}], "}"}], ",", "0.05"}], "]"}], ",", "\n", RowBox[{"Cylinder", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "5", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "5", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "5", "]"}], "]"}]}], "}"}]}], "}"}], ",", "0.05"}], "]"}], ",", "\n", RowBox[{"Cylinder", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{ RowBox[{"yn", "[", RowBox[{"[", "2", "]"}], "]"}], "-", "0.1"}], ",", RowBox[{"zn", "[", RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{ RowBox[{"yn", "[", RowBox[{"[", "1", "]"}], "]"}], "-", "0.1"}], ",", RowBox[{"zn", "[", RowBox[{"[", "1", "]"}], "]"}]}], "}"}]}], "}"}], ",", "0.05"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Cylinder", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{ RowBox[{"yn", "[", RowBox[{"[", "2", "]"}], "]"}], "+", "0.1"}], ",", RowBox[{"zn", "[", RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{ RowBox[{"yn", "[", RowBox[{"[", "1", "]"}], "]"}], "+", "0.1"}], ",", RowBox[{"zn", "[", RowBox[{"[", "1", "]"}], "]"}]}], "}"}]}], "}"}], ",", "0.05"}], "]"}], ",", RowBox[{"Cylinder", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "6", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "6", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "6", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "2", "]"}], "]"}]}], "}"}]}], "}"}], ",", "0.06"}], "]"}]}], "}"}]}], "}"}]}], "}"}], "]"}], " ", "/.", " ", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}]], "Input", CellChangeTimes->{{3.56785297434536*^9, 3.5678529770129647`*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{"overlapped", " ", "\[Pi]", " ", "system"}], " ", "*)"}], "\n", RowBox[{"Show", "[", RowBox[{"plot4", ",", RowBox[{"{", RowBox[{ RowBox[{"AspectRatio", " ", "\[Rule]", " ", "Automatic"}], ",", RowBox[{"Boxed", "->", "False"}]}], "}"}]}], "]"}]}]], "Input"], Cell[BoxData[ Graphics3DBox[{SphereBox[{1.2575069709757336`, 0, 0}, 0.2625736948837923], SphereBox[{-1.2575069709757336`, 0, 0}, 0.2625736948837923], {Opacity[0.6], {Opacity[0.4], {RGBColor[ NCache[ Rational[1, 3], 0.3333333333333333], 1, NCache[ Rational[1, 3], 0.3333333333333333]], SphereBox[{0, 0, 2.3758924039382725`}, 2.11331870905448], SphereBox[{0, 0, -2.3758924039382725`}, 2.11331870905448]}}, {RGBColor[0, 0, 0], SphereBox[{1.2575069709757336`, 0, 0}, 0.08], SphereBox[{-1.2575069709757336`, 0, 0}, 0.08], SphereBox[{2.3277116062908165`, 1.7550117957581575`, 0}, 0.08], SphereBox[{2.3277116062908165`, -1.7550117957581575`, 0}, 0.08], SphereBox[{-2.3277116062908165`, 1.7550117957581575`, 0}, 0.08], SphereBox[{-2.3277116062908165`, -1.7550117957581575`, 0}, 0.08], CylinderBox[{{2.3277116062908165`, 1.7550117957581575`, 0}, { 1.2575069709757336`, 0, 0}}, 0.05], CylinderBox[{{1.2575069709757336`, 0, 0}, { 2.3277116062908165`, -1.7550117957581575`, 0}}, 0.05], CylinderBox[{{-1.2575069709757336`, 0, 0}, {-2.3277116062908165`, 1.7550117957581575`, 0}}, 0.05], CylinderBox[{{-1.2575069709757336`, -0.1, 0}, { 1.2575069709757336`, -0.1, 0}}, 0.05], CylinderBox[{{-1.2575069709757336`, 0.1, 0}, { 1.2575069709757336`, 0.1, 0}}, 0.05], CylinderBox[{{-2.3277116062908165`, -1.7550117957581575`, 0}, {-1.2575069709757336`, 0, 0}}, 0.06]}}}, AspectRatio->Automatic, AutomaticImageSize->True, Boxed->False, ImageSize->{299.53443694076105`, 453.16208730374865`}, ViewPoint->{0.9203107977619439, -2.297577814780321, 2.307415051643735}, ViewVertical->{-0.09379359490747588, 0.22502713901592988`, 0.9930033672250635}]], "Output", CellChangeTimes->{ 3.567852992940593*^9, 3.5678530425018797`*^9, 3.5678531440580583`*^9, 3.567881557442132*^9, {3.5678815962394*^9, 3.567881636799471*^9}, 3.5678816722115335`*^9, {3.567921572102233*^9, 3.5679215965318756`*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"plot6", "=", RowBox[{ RowBox[{"Graphics3D", "[", RowBox[{"{", RowBox[{ RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"yc", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"zc", "[", RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", "Cr"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"yc", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"zc", "[", RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", "Cr"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"{", RowBox[{ RowBox[{"Opacity", "[", "0.6", "]"}], ",", "\[IndentingNewLine]", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "3", "]"}], "]"}], ",", RowBox[{"yc", "[", RowBox[{"[", "3", "]"}], "]"}], ",", RowBox[{"zc", "[", RowBox[{"[", "3", "]"}], "]"}]}], "}"}], ",", "CH"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "4", "]"}], "]"}], ",", RowBox[{"yc", "[", RowBox[{"[", "4", "]"}], "]"}], ",", RowBox[{"zc", "[", RowBox[{"[", "4", "]"}], "]"}]}], "}"}], ",", "CH"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "5", "]"}], "]"}], ",", RowBox[{"yc", "[", RowBox[{"[", "5", "]"}], "]"}], ",", RowBox[{"zc", "[", RowBox[{"[", "5", "]"}], "]"}]}], "}"}], ",", "CH"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "6", "]"}], "]"}], ",", RowBox[{"yc", "[", RowBox[{"[", "6", "]"}], "]"}], ",", RowBox[{"zc", "[", RowBox[{"[", "6", "]"}], "]"}]}], "}"}], ",", "CH"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "7", "]"}], "]"}], ",", RowBox[{"yc", "[", RowBox[{"[", "7", "]"}], "]"}], ",", RowBox[{"zc", "[", RowBox[{"[", "7", "]"}], "]"}]}], "}"}], ",", "CC"}], "]"}]}], "}"}], ",", "\[IndentingNewLine]", RowBox[{"{", RowBox[{ RowBox[{"Darker", "[", RowBox[{"Red", ",", "1"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "3", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "3", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "3", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "4", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "4", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "4", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "5", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "5", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "5", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "6", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "6", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "6", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], ",", RowBox[{"Cylinder", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{ RowBox[{"yn", "[", RowBox[{"[", "2", "]"}], "]"}], "-", "0.1"}], ",", RowBox[{"zn", "[", RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{ RowBox[{"yn", "[", RowBox[{"[", "1", "]"}], "]"}], "-", "0.1"}], ",", RowBox[{"zn", "[", RowBox[{"[", "1", "]"}], "]"}]}], "}"}]}], "}"}], ",", "0.05"}], "]"}], ",", RowBox[{"Cylinder", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{ RowBox[{"yn", "[", RowBox[{"[", "2", "]"}], "]"}], "+", "0.1"}], ",", RowBox[{"zn", "[", RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{ RowBox[{"yn", "[", RowBox[{"[", "1", "]"}], "]"}], "+", "0.1"}], ",", RowBox[{"zn", "[", RowBox[{"[", "1", "]"}], "]"}]}], "}"}]}], "}"}], ",", "0.05"}], "]"}], ",", RowBox[{"Cylinder", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "4", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "4", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "4", "]"}], "]"}]}], "}"}]}], "}"}], ",", "0.05"}], "]"}], ",", RowBox[{"Cylinder", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "3", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "3", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "3", "]"}], "]"}]}], "}"}]}], "}"}], ",", "0.05"}], "]"}], ",", RowBox[{"Cylinder", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "5", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "5", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "5", "]"}], "]"}]}], "}"}]}], "}"}], ",", "0.05"}], "]"}], ",", RowBox[{"Cylinder", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "6", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "6", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "6", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "2", "]"}], "]"}]}], "}"}]}], "}"}], ",", "0.05"}], "]"}]}], "}"}]}], "}"}], "]"}], "/.", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", StyleBox[ RowBox[{"(*", " ", RowBox[{ "the", " ", "sigma", " ", "skeleton", " ", "of", " ", "the", " ", "molecule"}], " ", "*)"}], "Text"]}], "\n", RowBox[{"Show", "[", RowBox[{"plot6", ",", RowBox[{"{", RowBox[{ RowBox[{"AspectRatio", "\[Rule]", "Automatic"}], ",", RowBox[{"Boxed", "\[Rule]", "False"}]}], "}"}]}], "]"}]}], "Input"], Cell[BoxData[ Graphics3DBox[{SphereBox[{1.2575069709757336`, 0, 0}, 0.2625736948837923], SphereBox[{-1.2575069709757336`, 0, 0}, 0.2625736948837923], {Opacity[0.6], SphereBox[{2.0288698715525646`, 1.2649459221637571`, 0}, 1.21900954961684], SphereBox[{2.0288698715525646`, -1.2649459221637571`, 0}, 1.21900954961684], SphereBox[{-2.0288698715525646`, 1.2649459221637571`, 0}, 1.21900954961684], SphereBox[{-2.0288698715525646`, -1.2649459221637571`, 0}, 1.21900954961684], SphereBox[{0, 0, 0}, 0.9949332760919413]}, {RGBColor[0, 0, 0], SphereBox[{1.2575069709757336`, 0, 0}, 0.08], SphereBox[{-1.2575069709757336`, 0, 0}, 0.08], SphereBox[{2.3277116062908165`, 1.7550117957581575`, 0}, 0.08], SphereBox[{2.3277116062908165`, -1.7550117957581575`, 0}, 0.08], SphereBox[{-2.3277116062908165`, 1.7550117957581575`, 0}, 0.08], SphereBox[{-2.3277116062908165`, -1.7550117957581575`, 0}, 0.08], CylinderBox[{{-1.2575069709757336`, -0.1, 0}, { 1.2575069709757336`, -0.1, 0}}, 0.05], CylinderBox[{{-1.2575069709757336`, 0.1, 0}, { 1.2575069709757336`, 0.1, 0}}, 0.05], CylinderBox[{{1.2575069709757336`, 0, 0}, { 2.3277116062908165`, -1.7550117957581575`, 0}}, 0.05], CylinderBox[{{1.2575069709757336`, 0, 0}, { 2.3277116062908165`, 1.7550117957581575`, 0}}, 0.05], CylinderBox[{{-1.2575069709757336`, 0, 0}, {-2.3277116062908165`, 1.7550117957581575`, 0}}, 0.05], CylinderBox[{{-2.3277116062908165`, -1.7550117957581575`, 0}, {-1.2575069709757336`, 0, 0}}, 0.05]}}, AspectRatio->Automatic, Boxed->False]], "Output", CellChangeTimes->{ 3.56785304253308*^9, 3.567853144073658*^9, 3.567881557457732*^9, { 3.5678815962706003`*^9, 3.567881636815071*^9}, 3.5678816722271338`*^9, { 3.567921572117833*^9, 3.567921596547476*^9}}] }, Open ]] }, WindowSize->{982, 1002}, WindowMargins->{{Automatic, 417}, {Automatic, 8}}, FrontEndVersion->"9.0 for Microsoft Windows (64-bit) (November 20, 2012)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[557, 20, 1051, 17, 31, "Input"], Cell[1611, 39, 3190, 60, 192, "Input"], Cell[4804, 101, 929, 16, 31, "Input"], Cell[5736, 119, 1307, 28, 112, "Input"], Cell[7046, 149, 761, 11, 31, "Input"], Cell[7810, 162, 2095, 52, 92, "Input"], Cell[9908, 216, 1424, 31, 52, "Input"], Cell[11335, 249, 2161, 50, 112, "Input"], Cell[13499, 301, 1496, 31, 72, "Input"], Cell[14998, 334, 1961, 39, 72, "Input"], Cell[16962, 375, 2710, 62, 92, "Input"], Cell[19675, 439, 2579, 60, 92, "Input"], Cell[22257, 501, 2896, 66, 112, "Input"], Cell[CellGroupData[{ Cell[25178, 571, 3324, 61, 92, "Input"], Cell[28505, 634, 1189, 22, 31, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[29731, 661, 1778, 47, 172, "Input"], Cell[31512, 710, 735, 10, 31, "Output"], Cell[32250, 722, 753, 11, 31, "Output"], Cell[33006, 735, 735, 10, 31, "Output"], Cell[33744, 747, 736, 10, 31, "Output"], Cell[34483, 759, 735, 10, 31, "Output"], Cell[35221, 771, 736, 10, 31, "Output"], Cell[35960, 783, 736, 10, 31, "Output"], Cell[36699, 795, 735, 10, 31, "Output"] }, Open ]], Cell[37449, 808, 9605, 257, 452, "Input"], Cell[CellGroupData[{ Cell[47079, 1069, 1021, 19, 72, "Input"], Cell[48103, 1090, 2288, 40, 297, "Output"] }, Open ]], Cell[50406, 1133, 8281, 229, 432, "Input"], Cell[CellGroupData[{ Cell[58712, 1366, 777, 19, 52, "Input"], Cell[59492, 1387, 2007, 36, 447, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[61536, 1428, 12791, 326, 552, "Input"], Cell[74330, 1756, 2246, 44, 447, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[76613, 1805, 1513, 46, 172, "Input"], Cell[78129, 1853, 1258, 30, 116, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[79424, 1888, 1053, 33, 72, "Input"], Cell[80480, 1923, 1608, 31, 116, "Output"] }, Open ]], Cell[82103, 1957, 10201, 263, 412, "Input"], Cell[CellGroupData[{ Cell[92329, 2224, 344, 9, 52, "Input"], Cell[92676, 2235, 2067, 40, 468, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[94780, 2280, 11488, 301, 492, "Input"], Cell[106271, 2583, 1878, 35, 325, "Output"] }, Open ]] } ] *) (* End of internal cache information *)