(* 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[ 111779, 2776] NotebookOptionsPosition[ 109492, 2703] NotebookOutlinePosition[ 109839, 2718] CellTagsIndexPosition[ 109796, 2715] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[BoxData[ RowBox[{"(*", " ", RowBox[{"Ethene", ",", " ", RowBox[{"21.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.5676639574765873`*^9, 3.567663958334589*^9}, 3.5676703779330645`*^9, {3.567785379855941*^9, 3.5677853802771416`*^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", ",", "\[IndentingNewLine]", "k1", ",", "k2", ",", "k3", ",", "k4", ",", "s1", ",", "s2", ",", "s3", ",", "s4", ",", "pch"}], "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"nc", "=", "11"}], ";", " ", RowBox[{"nn", "=", "6"}], ";", " ", RowBox[{"Z", "=", "6.0"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"bohr", "=", "0.529177"}], ";", RowBox[{"rad", "=", "57.29578"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"k1", "=", "1.032"}], ";", " ", RowBox[{"k2", "=", "1.12"}], ";", " ", RowBox[{"k3", "=", "1.23"}], ";", " ", RowBox[{"k4", "=", "1.6"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"s1", "=", "0.3"}], ";", RowBox[{"s2", "=", "0.3"}], ";", " ", RowBox[{"s3", "=", "0.3"}], ";"}], " "}]}], "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.5676650828937645`*^9, 3.5676651249046383`*^9}, { 3.567669337395637*^9, 3.56766933904924*^9}, {3.5676694423058214`*^9, 3.5676694436474237`*^9}, 3.5676697706707983`*^9, {3.5677840604836235`*^9, 3.567784065054432*^9}, 3.567784100684894*^9, {3.5677841378753595`*^9, 3.5677841860482445`*^9}, {3.5677842222091074`*^9, 3.567784250211157*^9}, { 3.567784700490348*^9, 3.567784725325591*^9}}], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{ RowBox[{"6", " ", RowBox[{"Variable", ":", " ", "Cr"}]}], ",", "CC", ",", "CH", ",", "Cpi", ",", "pch", ",", "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}], Cell[BoxData[ RowBox[{"(*", RowBox[{ RowBox[{"Cr", "=", "0.26"}], ";", " ", RowBox[{"CC", "=", "1.0"}], ";", " ", RowBox[{"CH", "=", "1.23"}], ";", " ", RowBox[{"Cpi", "=", "1.4"}], ";", " ", RowBox[{"w", "=", "1.02"}], ";", " ", RowBox[{"pch", "=", "0.56"}], ";"}], "*)"}]], "Input", CellChangeTimes->{{3.5676156777655764`*^9, 3.567615748683301*^9}, { 3.567615781240558*^9, 3.5676157870125685`*^9}, {3.5676162471041765`*^9, 3.567616264030206*^9}, 3.567616341141142*^9, {3.567616615342823*^9, 3.567616620444032*^9}, {3.567616774650303*^9, 3.567616775617505*^9}, { 3.5676171308145285`*^9, 3.5676171359469376`*^9}, 3.567666136363615*^9, { 3.5676694799486876`*^9, 3.5676694860950985`*^9}, {3.5676697424503484`*^9, 3.5676697434331503`*^9}, {3.567669776786009*^9, 3.56766977748801*^9}, { 3.5677834632365746`*^9, 3.567783468446984*^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"}]}]}], ";"}]}]}]], "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.567666104570759*^9, {3.567669344322049*^9, 3.567669345367251*^9}}], Cell[BoxData[ RowBox[{"(*", RowBox[{"bonding", " ", "clouds", " ", "and", " ", "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}, { 3.5676669451938353`*^9, 3.5676669476430397`*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"Ekin", "=", RowBox[{"Ekin", "+", " ", RowBox[{"2.25", "*", RowBox[{"(", RowBox[{ RowBox[{"k2", "/", RowBox[{"CC", "^", "2"}]}], "+", RowBox[{"4", "*", 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", "*", RowBox[{"s3", "/", "CH"}]}]}], ")"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"vne", "=", RowBox[{"vne", "-", RowBox[{"4", "*", RowBox[{ RowBox[{"(", RowBox[{"3.0", "-", RowBox[{ RowBox[{"(", RowBox[{"pch", "/", "CH"}], ")"}], "^", "2"}]}], ")"}], "/", "CH"}]}]}]}], ";"}]}], "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.567666081389118*^9, { 3.5676693508428607`*^9, 3.567669351716462*^9}, {3.567684266942669*^9, 3.5676842693294725`*^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", "+", "pch"}], ")"}], "*", "csw"}]}]}], ";", " ", RowBox[{"chy", "=", " ", RowBox[{ RowBox[{"(", RowBox[{"Cr", "+", "CH", "+", "pch"}], ")"}], "*", "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", "=", "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.5676693609048786`*^9, 3.5676693619344797`*^9}, {3.567784653378265*^9, 3.5677846542830667`*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"cn", "=", RowBox[{"{", RowBox[{ "Z", ",", " ", "Z", ",", " ", "1", ",", " ", "1", ",", " ", "1", ",", " ", "1"}], "}"}]}], ";"}], "\[IndentingNewLine]", 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.567666020471011*^9, 3.5676660624974847`*^9, {3.5676661584220533`*^9, 3.56766620225813*^9}, { 3.5676693672540894`*^9, 3.5676693682368913`*^9}, {3.5676839196548586`*^9, 3.5676839242724667`*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"cc", "=", RowBox[{"{", RowBox[{ RowBox[{"-", "2"}], ",", " ", RowBox[{"-", "2"}], ",", " ", RowBox[{"-", "2"}], ",", " ", RowBox[{"-", "2"}], ",", " ", RowBox[{"-", "2"}], ",", " ", RowBox[{"-", "2"}], ",", RowBox[{"-", "2"}], ",", RowBox[{ RowBox[{"-", "1"}], "/", "2"}], ",", RowBox[{ RowBox[{"-", "1"}], "/", "2"}], ",", RowBox[{ RowBox[{"-", "1"}], "/", "2"}], ",", RowBox[{ RowBox[{"-", "1"}], "/", "2"}]}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"xc", "=", RowBox[{"{", RowBox[{"cx", ",", RowBox[{"-", "cx"}], ",", " ", "cchx", ",", " ", "cchx", ",", RowBox[{"-", "cchx"}], ",", RowBox[{"-", "cchx"}], ",", " ", "0", ",", " ", "cx", ",", " ", "cx", ",", " ", RowBox[{"-", "cx"}], ",", " ", RowBox[{"-", "cx"}]}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"yc", "=", RowBox[{"{", " ", RowBox[{"0", ",", " ", "0", ",", " ", "cchy", ",", RowBox[{"-", "cchy"}], ",", " ", "cchy", ",", RowBox[{"-", "cchy"}], ",", " ", "0", ",", " ", "0", ",", " ", "0", ",", " ", "0", ",", " ", "0"}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"zc", "=", RowBox[{"{", " ", RowBox[{ "0", ",", " ", "0", ",", " ", "0", ",", " ", "0", ",", " ", "0", ",", " ", "0", ",", " ", "0", ",", " ", "ccpi", ",", RowBox[{"-", "ccpi"}], ",", "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.5676660330758333`*^9, 3.5676660669122925`*^9, {3.567666211150146*^9, 3.5676662904606853`*^9}, { 3.5676663213799396`*^9, 3.5676664350261393`*^9}, {3.5676664797202177`*^9, 3.567666585800404*^9}, {3.567669373962101*^9, 3.5676693750229034`*^9}, { 3.5677833575151887`*^9, 3.567783390400047*^9}, {3.567783878150503*^9, 3.5677839827330875`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"(*", "ww", "*)"}], "\[IndentingNewLine]", RowBox[{ 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"}]}], "]"}]}]}]}]}]}], "]"}]}], "]"}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"vee", "=", RowBox[{"vee", "-", RowBox[{ RowBox[{"(", RowBox[{"1", "/", "4"}], ")"}], "/", "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.5676014051115084`*^9}, {3.5676177595268326`*^9, 3.567617780103269*^9}, { 3.5676178274649525`*^9, 3.5676178737346334`*^9}, {3.5676608539135365`*^9, 3.5676608549587383`*^9}, {3.56766792142795*^9, 3.567667928603963*^9}, { 3.5676679666368294`*^9, 3.5676679998648877`*^9}, {3.5676680373517537`*^9, 3.56766807504142*^9}, {3.567668660557248*^9, 3.567668705422927*^9}, { 3.567668743424594*^9, 3.5676688469307756`*^9}, {3.5676688805644345`*^9, 3.567668901437271*^9}, {3.5676689491889553`*^9, 3.5676689567705684`*^9}, { 3.567668992510231*^9, 3.567669029326296*^9}, {3.5676690617431526`*^9, 3.5676692209726324`*^9}, {3.56766938469492*^9, 3.5676693871285243`*^9}, { 3.567670317857359*^9, 3.567670330306181*^9}, {3.5677835407063107`*^9, 3.5677835537323337`*^9}, 3.5677836061640253`*^9, {3.567783659172919*^9, 3.567783668439335*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"vnn", "=", "0.0"}], ";"}]], "Input"], 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"}]}], "]"}]}]}]}]}]}], "]"}]}], "]"}], ";"}]}]], "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.5676672208775196`*^9, { 3.567669392167333*^9, 3.567669393025335*^9}, {3.567670066291317*^9, 3.567670127459025*^9}, {3.567670195911945*^9, 3.567670213227975*^9}, { 3.567670271681278*^9, 3.567670272242879*^9}, 3.5676832749993267`*^9, { 3.5677835115186596`*^9, 3.567783526619486*^9}, 3.5677835998460145`*^9, 3.567783693477379*^9}], Cell[BoxData[ RowBox[{ RowBox[{"(*", "nw", "*)"}], "\[IndentingNewLine]", RowBox[{ 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[{"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[{"zc", "[", RowBox[{"[", "i", "]"}], "]"}], ")"}], "^", "2"}]}], "]"}]}]}]}]}]}], "]"}]}], "]"}]}], "]"}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{ RowBox[{"{", RowBox[{"Method", "\[Rule]", "\"\\""}], "}"}], ",", RowBox[{"{", 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.567601421335537*^9}, { 3.5676156065414515`*^9, 3.5676156341635*^9}, 3.567663665537675*^9, { 3.5676674885115895`*^9, 3.567667493924799*^9}, {3.567669397830143*^9, 3.567669398703745*^9}, {3.5676701421386504`*^9, 3.5676701592518806`*^9}, 3.5677837050837994`*^9}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"Epot", "=", RowBox[{"vne", "+", "vee", "+", "vnn"}]}], ";"}], "\n", RowBox[{ RowBox[{"func", "=", RowBox[{"Ekin", "+", "Epot"}]}], ";"}], "\n", RowBox[{"t", "=", RowBox[{"FindMinimum", "[", RowBox[{"func", ",", RowBox[{"{", RowBox[{"Cr", ",", "0.26"}], "}"}], ",", RowBox[{"{", RowBox[{"CC", ",", "1.0"}], "}"}], ",", RowBox[{"{", RowBox[{"CH", ",", "1.25"}], "}"}], ",", RowBox[{"{", RowBox[{"Cpi", ",", "1.41"}], "}"}], ",", RowBox[{"{", RowBox[{"w", ",", "1.05"}], "}"}], ",", RowBox[{"{", RowBox[{"pch", ",", "0.55"}], "}"}], ",", 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.567663923998929*^9, 3.5676639414397593`*^9}, {3.5676654228807616`*^9, 3.567665457918423*^9}, {3.567667087403685*^9, 3.567667091381692*^9}, { 3.5676694073305597`*^9, 3.5676694083757615`*^9}, {3.567669500103923*^9, 3.5676695340807824`*^9}, {3.567783415500491*^9, 3.567783421256901*^9}, { 3.5677847533900404`*^9, 3.5677847859940977`*^9}, {3.5677851285863*^9, 3.5677851328763075`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"-", "78.01491117629921`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Cr", "\[Rule]", "0.26420244686116506`"}], ",", RowBox[{"CC", "\[Rule]", "0.991214880714466`"}], ",", RowBox[{"CH", "\[Rule]", "1.2429612744569003`"}], ",", RowBox[{"Cpi", "\[Rule]", "1.4323847188578354`"}], ",", RowBox[{"w", "\[Rule]", "1.053835233730336`"}], ",", RowBox[{"pch", "\[Rule]", "0.5454675488039527`"}]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{ 3.5677836114680347`*^9, 3.5677838254380107`*^9, 3.5677840748356485`*^9, 3.5677841065349045`*^9, {3.5677841440061703`*^9, 3.567784193052656*^9}, { 3.567784232349126*^9, 3.5677842556087666`*^9}, 3.5677846611626787`*^9, 3.567784731752803*^9, 3.567784792280909*^9, 3.567785042630149*^9, 3.5677851391163177`*^9}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ RowBox[{"t", "=", RowBox[{"{", RowBox[{"func", ",", RowBox[{"{", RowBox[{ RowBox[{"Cr", "\[Rule]", "0.26"}], ",", RowBox[{"CC", "\[Rule]", "1.0"}], ",", RowBox[{"CH", "\[Rule]", "1.23"}], ",", RowBox[{"Cpi", "\[Rule]", "1.4"}], ",", RowBox[{"w", "\[Rule]", "1.02"}], ",", RowBox[{"pch", "\[Rule]", "0.56"}]}], "}"}]}], "}"}]}], ";"}], "*)"}], "\[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.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.567663923998929*^9, 3.5676639414397593`*^9}, {3.5676654228807616`*^9, 3.567665457918423*^9}, {3.567667087403685*^9, 3.567667091381692*^9}, { 3.5676694073305597`*^9, 3.5676694083757615`*^9}, {3.567669500103923*^9, 3.5676695340807824`*^9}, {3.567669641486971*^9, 3.5676696483977833`*^9}, { 3.5676698093588657`*^9, 3.5676698227592897`*^9}, {3.5677834276685123`*^9, 3.5677834352813253`*^9}, 3.567783586897992*^9}], 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", "+", "pch"}], ")"}], "*", "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}}], Cell[BoxData["78.0149111775737`"], "Output", CellChangeTimes->{ 3.567783611483635*^9, 3.567783825453611*^9, 3.5677840748512487`*^9, 3.5677841065505047`*^9, {3.5677841440217705`*^9, 3.5677841930682564`*^9}, { 3.5677842323647256`*^9, 3.5677842556243663`*^9}, 3.567784661178279*^9, 3.567784731768403*^9, 3.567784792296509*^9, 3.567785042645749*^9, 3.567785139131918*^9}], Cell[BoxData[ RowBox[{"-", "249.22054196817865`"}]], "Output", CellChangeTimes->{ 3.567783611483635*^9, 3.567783825453611*^9, 3.5677840748512487`*^9, 3.5677841065505047`*^9, {3.5677841440217705`*^9, 3.5677841930682564`*^9}, { 3.5677842323647256`*^9, 3.5677842556243663`*^9}, 3.567784661178279*^9, 3.567784731768403*^9, 3.567784792296509*^9, 3.567785042645749*^9, 3.567785139131918*^9}], Cell[BoxData["59.739155892521524`"], "Output", CellChangeTimes->{ 3.567783611483635*^9, 3.567783825453611*^9, 3.5677840748512487`*^9, 3.5677841065505047`*^9, {3.5677841440217705`*^9, 3.5677841930682564`*^9}, { 3.5677842323647256`*^9, 3.5677842556243663`*^9}, 3.567784661178279*^9, 3.567784731768403*^9, 3.567784792296509*^9, 3.567785042645749*^9, 3.567785139131918*^9}], Cell[BoxData["33.45156372178421`"], "Output", CellChangeTimes->{ 3.567783611483635*^9, 3.567783825453611*^9, 3.5677840748512487`*^9, 3.5677841065505047`*^9, {3.5677841440217705`*^9, 3.5677841930682564`*^9}, { 3.5677842323647256`*^9, 3.5677842556243663`*^9}, 3.567784661178279*^9, 3.567784731768403*^9, 3.567784792296509*^9, 3.567785042645749*^9, 3.567785139147518*^9}], Cell[BoxData["1.9999999999836644`"], "Output", CellChangeTimes->{ 3.567783611483635*^9, 3.567783825453611*^9, 3.5677840748512487`*^9, 3.5677841065505047`*^9, {3.5677841440217705`*^9, 3.5677841930682564`*^9}, { 3.5677842323647256`*^9, 3.5677842556243663`*^9}, 3.567784661178279*^9, 3.567784731768403*^9, 3.567784792296509*^9, 3.567785042645749*^9, 3.567785139147518*^9}], Cell[BoxData["120.76062341612383`"], "Output", CellChangeTimes->{ 3.567783611483635*^9, 3.567783825453611*^9, 3.5677840748512487`*^9, 3.5677841065505047`*^9, {3.5677841440217705`*^9, 3.5677841930682564`*^9}, { 3.5677842323647256`*^9, 3.5677842556243663`*^9}, 3.567784661178279*^9, 3.567784731768403*^9, 3.567784792296509*^9, 3.567785042645749*^9, 3.567785139147518*^9}], Cell[BoxData["1.3286759503089796`"], "Output", CellChangeTimes->{ 3.567783611483635*^9, 3.567783825453611*^9, 3.5677840748512487`*^9, 3.5677841065505047`*^9, {3.5677841440217705`*^9, 3.5677841930682564`*^9}, { 3.5677842323647256`*^9, 3.5677842556243663`*^9}, 3.567784661178279*^9, 3.567784731768403*^9, 3.567784792296509*^9, 3.567785042645749*^9, 3.567785139147518*^9}], Cell[BoxData["1.086205257629359`"], "Output", CellChangeTimes->{ 3.567783611483635*^9, 3.567783825453611*^9, 3.5677840748512487`*^9, 3.5677841065505047`*^9, {3.5677841440217705`*^9, 3.5677841930682564`*^9}, { 3.5677842323647256`*^9, 3.5677842556243663`*^9}, 3.567784661178279*^9, 3.567784731768403*^9, 3.567784792296509*^9, 3.567785042645749*^9, 3.567785139147518*^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[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "10", "]"}], "]"}], ",", RowBox[{"yc", "[", RowBox[{"[", "10", "]"}], "]"}]}], "}"}], ",", "Cpi"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "11", "]"}], "]"}], ",", RowBox[{"yc", "[", RowBox[{"[", "11", "]"}], "]"}]}], "}"}], ",", "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", "]"}], ",", "\n", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "1", "]"}], "]"}]}], "}"}]}], "}"}], "]"}], ",", "\n", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "3", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "3", "]"}], "]"}]}], "}"}]}], "}"}], "]"}], ",", "\n", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "4", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "4", "]"}], "]"}]}], "}"}]}], "}"}], "]"}], ",", "\n", 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.5676637340373955`*^9, 3.5676637518214264`*^9}, {3.5676844549385986`*^9, 3.5676844556562*^9}, { 3.56778379094635*^9, 3.5677837984967637`*^9}, {3.5677848901086807`*^9, 3.5677848906546817`*^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.5676637340373955`*^9, 3.5676637518214264`*^9}}], Cell[BoxData[ GraphicsBox[{CircleBox[{1.255417327575631, 0}, 0.26420244686116506`], CircleBox[{-1.255417327575631, 0}, 0.26420244686116506`], CircleBox[{2.0003188684550133`, 1.3102153171341118`}, 1.2429612744569003`], CircleBox[{2.0003188684550133`, -1.3102153171341118`}, 1.2429612744569003`], CircleBox[{-2.0003188684550133`, 1.3102153171341118`}, 1.2429612744569003`], CircleBox[{-2.0003188684550133`, -1.3102153171341118`}, 1.2429612744569003`], CircleBox[{0, 0}, 0.991214880714466], CircleBox[{1.255417327575631, 0}, 1.4323847188578354`], CircleBox[{1.255417327575631, 0}, 1.4323847188578354`], CircleBox[{-1.255417327575631, 0}, 1.4323847188578354`], CircleBox[{-1.255417327575631, 0}, 1.4323847188578354`], DiskBox[{1.255417327575631, 0}, 0.08], DiskBox[{-1.255417327575631, 0}, 0.08], DiskBox[{2.26991109115556, 1.7844039718461064`}, 0.08], DiskBox[{2.26991109115556, -1.7844039718461064`}, 0.08], DiskBox[{-2.26991109115556, 1.7844039718461064`}, 0.08], DiskBox[{-2.26991109115556, -1.7844039718461064`}, 0.08], {Thickness[0.008], LineBox[{{-1.255417327575631, 0}, {1.255417327575631, 0}}], LineBox[{{1.255417327575631, 0}, {2.26991109115556, 1.7844039718461064`}}], LineBox[{{1.255417327575631, 0}, { 2.26991109115556, -1.7844039718461064`}}], LineBox[{{-1.255417327575631, 0}, {-2.26991109115556, 1.7844039718461064`}}], LineBox[{{-1.255417327575631, 0}, {-2.26991109115556, -1.7844039718461064`}}]}}, AspectRatio->Automatic, Axes->True, Frame->True, GridLines->Automatic, PlotRange->Automatic]], "Output", CellChangeTimes->{ 3.567783611514835*^9, 3.567783825484811*^9, 3.567784074882449*^9, 3.5677841065817046`*^9, {3.5677841440529704`*^9, 3.5677841930838566`*^9}, { 3.5677842323959255`*^9, 3.567784255655566*^9}, 3.567784661209479*^9, 3.5677847317996025`*^9, 3.567784792327709*^9, 3.5677850426769485`*^9, 3.5677851391631184`*^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[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "10", "]"}], "]"}], ",", RowBox[{"zc", "[", RowBox[{"[", "10", "]"}], "]"}]}], "}"}], ",", "Cpi"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "11", "]"}], "]"}], ",", RowBox[{"zc", "[", RowBox[{"[", "11", "]"}], "]"}]}], "}"}], ",", "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", "]"}], ",", "\n", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "1", "]"}], "]"}]}], "}"}]}], "}"}], "]"}], ",", "\n", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "3", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "3", "]"}], "]"}]}], "}"}]}], "}"}], "]"}], ",", "\n", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "4", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "4", "]"}], "]"}]}], "}"}]}], "}"}], "]"}], ",", "\n", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "5", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "5", "]"}], "]"}]}], "}"}]}], "}"}], "]"}], ",", 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.567664762672002*^9, 3.567664763623604*^9}, { 3.567664799082466*^9, 3.5676648803742085`*^9}, {3.5676844666386194`*^9, 3.567684487324256*^9}, 3.5676846491589403`*^9, {3.567684686037405*^9, 3.5676847139614534`*^9}, {3.5677838129579887`*^9, 3.5677838199000015`*^9}, {3.567785070210997*^9, 3.567785070678998*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{"xz", "-", "plane"}], " ", "*)"}], "\n", 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[{"-", "3.5"}], ",", "3.5"}], "}"}]}], "}"}]}], ",", " ", RowBox[{"Frame", " ", "->", " ", "True"}]}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.5676647867584443`*^9, 3.5676647880064464`*^9}, { 3.567664891278628*^9, 3.56766489241743*^9}, {3.5676821221261015`*^9, 3.567682182435807*^9}}], Cell[BoxData[ GraphicsBox[{CircleBox[{1.255417327575631, 0}, 0.26420244686116506`], CircleBox[{-1.255417327575631, 0}, 0.26420244686116506`], CircleBox[{2.0003188684550133`, 0}, 1.2429612744569003`], CircleBox[{2.0003188684550133`, 0}, 1.2429612744569003`], CircleBox[{-2.0003188684550133`, 0}, 1.2429612744569003`], CircleBox[{-2.0003188684550133`, 0}, 1.2429612744569003`], CircleBox[{0, 0}, 0.991214880714466], CircleBox[{1.255417327575631, 1.4323847188578354`}, 1.4323847188578354`], CircleBox[{1.255417327575631, -1.4323847188578354`}, 1.4323847188578354`], CircleBox[{-1.255417327575631, 1.4323847188578354`}, 1.4323847188578354`], CircleBox[{-1.255417327575631, -1.4323847188578354`}, 1.4323847188578354`], DiskBox[{1.255417327575631, 0}, 0.08], DiskBox[{-1.255417327575631, 0}, 0.08], DiskBox[{2.26991109115556, 0}, 0.08], DiskBox[{2.26991109115556, 0}, 0.08], DiskBox[{-2.26991109115556, 0}, 0.08], DiskBox[{-2.26991109115556, 0}, 0.08], {Thickness[0.008], LineBox[{{-1.255417327575631, 0}, {1.255417327575631, 0}}], LineBox[{{1.255417327575631, 0}, {2.26991109115556, 0}}], LineBox[{{1.255417327575631, 0}, {2.26991109115556, 0}}], LineBox[{{-1.255417327575631, 0}, {-2.26991109115556, 0}}], LineBox[{{-1.255417327575631, 0}, {-2.26991109115556, 0}}]}}, AspectRatio->Automatic, Axes->True, Frame->True, GridLines->Automatic, PlotRange->{{-3.5, 3.5}, {-3.5, 3.5}}]], "Output", CellChangeTimes->{ 3.5677836115460353`*^9, 3.5677838255160108`*^9, 3.567784074898049*^9, 3.5677841066129045`*^9, {3.56778414406857*^9, 3.5677841931150565`*^9}, { 3.5677842324115257`*^9, 3.5677842556711664`*^9}, 3.567784661256279*^9, 3.567784731830803*^9, 3.5677847923433094`*^9, 3.5677850426925488`*^9, 3.5677851391943183`*^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[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "10", "]"}], "]"}], ",", RowBox[{"yc", "[", RowBox[{"[", "10", "]"}], "]"}], ",", RowBox[{"zc", "[", RowBox[{"[", "10", "]"}], "]"}]}], "}"}], ",", "Cpi"}], "]"}], ",", "\n", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "11", "]"}], "]"}], ",", RowBox[{"yc", "[", RowBox[{"[", "11", "]"}], "]"}], ",", RowBox[{"zc", "[", RowBox[{"[", "11", "]"}], "]"}]}], "}"}], ",", "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.5677849413391705`*^9, 3.5677849413391705`*^9}}], Cell[BoxData[ Graphics3DBox[{SphereBox[{1.255417327575631, 0, 0}, 0.26420244686116506`], SphereBox[{-1.255417327575631, 0, 0}, 0.26420244686116506`], {Opacity[0.6], SphereBox[{2.0003188684550133`, 1.3102153171341118`, 0}, 1.2429612744569003`], SphereBox[{2.0003188684550133`, -1.3102153171341118`, 0}, 1.2429612744569003`], SphereBox[{-2.0003188684550133`, 1.3102153171341118`, 0}, 1.2429612744569003`], SphereBox[{-2.0003188684550133`, -1.3102153171341118`, 0}, 1.2429612744569003`], SphereBox[{0, 0, 0}, 0.991214880714466], {Opacity[0.4], {RGBColor[ NCache[ Rational[1, 3], 0.3333333333333333], 1, NCache[ Rational[1, 3], 0.3333333333333333]], SphereBox[{1.255417327575631, 0, 1.4323847188578354`}, 1.4323847188578354`], SphereBox[{1.255417327575631, 0, -1.4323847188578354`}, 1.4323847188578354`], SphereBox[{-1.255417327575631, 0, 1.4323847188578354`}, 1.4323847188578354`], SphereBox[{-1.255417327575631, 0, -1.4323847188578354`}, 1.4323847188578354`]}}, {RGBColor[0, 0, 0], SphereBox[{1.255417327575631, 0, 0}, 0.08], SphereBox[{-1.255417327575631, 0, 0}, 0.08], SphereBox[{2.26991109115556, 1.7844039718461064`, 0}, 0.08], SphereBox[{2.26991109115556, -1.7844039718461064`, 0}, 0.08], SphereBox[{-2.26991109115556, 1.7844039718461064`, 0}, 0.08], SphereBox[{-2.26991109115556, -1.7844039718461064`, 0}, 0.08], CylinderBox[{{-1.255417327575631, -0.1, 0}, { 1.255417327575631, -0.1, 0}}, 0.05], CylinderBox[{{-1.255417327575631, 0.1, 0}, {1.255417327575631, 0.1, 0}}, 0.05], CylinderBox[{{2.26991109115556, 1.7844039718461064`, 0}, { 1.255417327575631, 0, 0}}, 0.05], CylinderBox[{{1.255417327575631, 0, 0}, { 2.26991109115556, -1.7844039718461064`, 0}}, 0.05], CylinderBox[{{-1.255417327575631, 0, 0}, {-2.26991109115556, 1.7844039718461064`, 0}}, 0.05], CylinderBox[{{-2.26991109115556, -1.7844039718461064`, 0}, {-1.255417327575631, 0, 0}}, 0.05]}}}, AspectRatio->Automatic, AutomaticImageSize->True, Boxed->False, ImageSize->{385.45918281140575`, 410.5368322903669}, ViewPoint->{0.4491978220218571, -1.5871504529584444`, 2.9545176859116316`}, ViewVertical->{-0.34506603591729473`, 0.37933026621964583`, 1.0073743849008687`}]], "Output", CellChangeTimes->{3.5677850427237487`*^9, 3.5677851392099185`*^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]", StyleBox[ RowBox[{"(*", " ", RowBox[{ "Cartesian", " ", "Coordinates", " ", "of", " ", "stationary", " ", "structure"}], " ", "*)"}], "Text"]}], "\[IndentingNewLine]", RowBox[{"u", "=", RowBox[{"MatrixForm", "[", "s", "]"}]}]}], "Input"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"C", "0.6643379751544898`", "0.`", "0.`"}, {"C", RowBox[{"-", "0.6643379751544898`"}], "0.`", "0.`"}, {"H", "1.2011847414844259`", "0.9442655406096071`", "0.`"}, {"H", "1.2011847414844259`", RowBox[{"-", "0.9442655406096071`"}], "0.`"}, {"H", RowBox[{"-", "1.2011847414844259`"}], "0.9442655406096071`", "0.`"}, {"H", RowBox[{"-", "1.2011847414844259`"}], RowBox[{"-", "0.9442655406096071`"}], "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.5677850427237487`*^9, 3.567785139225518*^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.3286759503089796`", "1.0862052576293593`", "1.0862052576293593`", "2.090887949527328`", "2.090887949527328`"}, {"1.3286759503089796`", "0.`", "2.090887949527328`", "2.090887949527328`", "1.0862052576293593`", "1.0862052576293593`"}, {"1.0862052576293593`", "2.090887949527328`", "0.`", "1.8885310812192142`", "2.4023694829688518`", "3.05580247683502`"}, {"1.0862052576293593`", "2.090887949527328`", "1.8885310812192142`", "0.`", "3.05580247683502`", "2.4023694829688518`"}, {"2.090887949527328`", "1.0862052576293593`", "2.4023694829688518`", "3.05580247683502`", "0.`", "1.8885310812192142`"}, {"2.090887949527328`", "1.0862052576293593`", "3.05580247683502`", "2.4023694829688518`", "1.8885310812192142`", "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.567785042739349*^9, 3.5677851392411184`*^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[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "10", "]"}], "]"}], ",", RowBox[{"yc", "[", RowBox[{"[", "10", "]"}], "]"}], ",", RowBox[{"zc", "[", RowBox[{"[", "10", "]"}], "]"}]}], "}"}], ",", "Cpi"}], "]"}], ",", "\n", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xc", "[", RowBox[{"[", "11", "]"}], "]"}], ",", RowBox[{"yc", "[", RowBox[{"[", "11", "]"}], "]"}], ",", RowBox[{"zc", "[", RowBox[{"[", "11", "]"}], "]"}]}], "}"}], ",", "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"], 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.255417327575631, 0, 0}, 0.26420244686116506`], SphereBox[{-1.255417327575631, 0, 0}, 0.26420244686116506`], {Opacity[0.6], {Opacity[0.4], {RGBColor[ NCache[ Rational[1, 3], 0.3333333333333333], 1, NCache[ Rational[1, 3], 0.3333333333333333]], SphereBox[{1.255417327575631, 0, 1.4323847188578354`}, 1.4323847188578354`], SphereBox[{1.255417327575631, 0, -1.4323847188578354`}, 1.4323847188578354`], SphereBox[{-1.255417327575631, 0, 1.4323847188578354`}, 1.4323847188578354`], SphereBox[{-1.255417327575631, 0, -1.4323847188578354`}, 1.4323847188578354`]}}, {RGBColor[0, 0, 0], SphereBox[{1.255417327575631, 0, 0}, 0.08], SphereBox[{-1.255417327575631, 0, 0}, 0.08], SphereBox[{2.26991109115556, 1.7844039718461064`, 0}, 0.08], SphereBox[{2.26991109115556, -1.7844039718461064`, 0}, 0.08], SphereBox[{-2.26991109115556, 1.7844039718461064`, 0}, 0.08], SphereBox[{-2.26991109115556, -1.7844039718461064`, 0}, 0.08], CylinderBox[{{2.26991109115556, 1.7844039718461064`, 0}, { 1.255417327575631, 0, 0}}, 0.05], CylinderBox[{{1.255417327575631, 0, 0}, { 2.26991109115556, -1.7844039718461064`, 0}}, 0.05], CylinderBox[{{-1.255417327575631, 0, 0}, {-2.26991109115556, 1.7844039718461064`, 0}}, 0.05], CylinderBox[{{-1.255417327575631, -0.1, 0}, { 1.255417327575631, -0.1, 0}}, 0.05], CylinderBox[{{-1.255417327575631, 0.1, 0}, {1.255417327575631, 0.1, 0}}, 0.05], CylinderBox[{{-2.26991109115556, -1.7844039718461064`, 0}, {-1.255417327575631, 0, 0}}, 0.06]}}}, AspectRatio->Automatic, AutomaticImageSize->True, Boxed->False, ImageSize->{372.2343412308929, 421.7673971606976}, ViewPoint->{1.2988110619305804`, -2.066183131782444, 2.3439234397361295`}, ViewVertical->{-0.059890489600677506`, 0.1859116022268886, 0.9910618113563924}]], "Output", CellChangeTimes->{3.567785042770549*^9, 3.5677851392567186`*^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", "]"}], "]"}]}]}], ";"}], "\n", "\[IndentingNewLine]", StyleBox[ RowBox[{"(*", " ", RowBox[{ "the", " ", "sigma", " ", "skeleton", " ", "of", " ", "the", " ", "molecule"}], " ", "*)"}], "Text"]}], "\[IndentingNewLine]", RowBox[{"Show", "[", RowBox[{"plot6", ",", RowBox[{"{", RowBox[{ RowBox[{"AspectRatio", "\[Rule]", "Automatic"}], ",", RowBox[{"Boxed", "\[Rule]", "False"}]}], "}"}]}], "]"}]}], "Input"], Cell[BoxData[ Graphics3DBox[{SphereBox[{1.255417327575631, 0, 0}, 0.26420244686116506`], SphereBox[{-1.255417327575631, 0, 0}, 0.26420244686116506`], {Opacity[0.6], SphereBox[{2.0003188684550133`, 1.3102153171341118`, 0}, 1.2429612744569003`], SphereBox[{2.0003188684550133`, -1.3102153171341118`, 0}, 1.2429612744569003`], SphereBox[{-2.0003188684550133`, 1.3102153171341118`, 0}, 1.2429612744569003`], SphereBox[{-2.0003188684550133`, -1.3102153171341118`, 0}, 1.2429612744569003`], SphereBox[{0, 0, 0}, 0.991214880714466]}, {RGBColor[0, 0, 0], SphereBox[{1.255417327575631, 0, 0}, 0.08], SphereBox[{-1.255417327575631, 0, 0}, 0.08], SphereBox[{2.26991109115556, 1.7844039718461064`, 0}, 0.08], SphereBox[{2.26991109115556, -1.7844039718461064`, 0}, 0.08], SphereBox[{-2.26991109115556, 1.7844039718461064`, 0}, 0.08], SphereBox[{-2.26991109115556, -1.7844039718461064`, 0}, 0.08], CylinderBox[{{-1.255417327575631, -0.1, 0}, {1.255417327575631, -0.1, 0}}, 0.05], CylinderBox[{{-1.255417327575631, 0.1, 0}, { 1.255417327575631, 0.1, 0}}, 0.05], CylinderBox[{{1.255417327575631, 0, 0}, { 2.26991109115556, -1.7844039718461064`, 0}}, 0.05], CylinderBox[{{1.255417327575631, 0, 0}, { 2.26991109115556, 1.7844039718461064`, 0}}, 0.05], CylinderBox[{{-1.255417327575631, 0, 0}, {-2.26991109115556, 1.7844039718461064`, 0}}, 0.05], CylinderBox[{{-2.26991109115556, -1.7844039718461064`, 0}, {-1.255417327575631, 0, 0}}, 0.05]}}, AspectRatio->Automatic, Boxed->False]], "Output", CellChangeTimes->{3.567785042786149*^9, 3.5677851392723184`*^9}] }, Open ]] }, WindowSize->{833, 1002}, WindowMargins->{{468, Automatic}, {Automatic, 13}}, 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, 905, 14, 31, "Input"], Cell[1465, 36, 2611, 47, 172, "Input"], Cell[4079, 85, 907, 16, 31, "Input"], Cell[4989, 103, 882, 16, 31, "Input"], Cell[5874, 121, 1314, 29, 92, "Input"], Cell[7191, 152, 844, 13, 31, "Input"], Cell[8038, 167, 1993, 50, 72, "Input"], Cell[10034, 219, 2028, 49, 112, "Input"], Cell[12065, 270, 1856, 40, 92, "Input"], Cell[13924, 312, 2867, 63, 92, "Input"], Cell[16794, 377, 3485, 77, 192, "Input"], Cell[20282, 456, 71, 2, 31, "Input"], Cell[20356, 460, 2492, 56, 92, "Input"], Cell[22851, 518, 3006, 71, 192, "Input"], Cell[CellGroupData[{ Cell[25882, 593, 2899, 56, 92, "Input"], Cell[28784, 651, 852, 18, 52, "Output"] }, Open ]], Cell[29651, 672, 2617, 45, 52, "Input"], Cell[CellGroupData[{ Cell[32293, 721, 1741, 45, 172, "Input"], Cell[34037, 768, 385, 6, 31, "Output"], Cell[34425, 776, 404, 7, 31, "Output"], Cell[34832, 785, 387, 6, 31, "Output"], Cell[35222, 793, 386, 6, 31, "Output"], Cell[35611, 801, 387, 6, 31, "Output"], Cell[36001, 809, 387, 6, 31, "Output"], Cell[36391, 817, 387, 6, 31, "Output"], Cell[36781, 825, 386, 6, 31, "Output"] }, Open ]], Cell[37182, 834, 9392, 254, 492, "Input"], Cell[CellGroupData[{ Cell[46599, 1092, 995, 19, 72, "Input"], Cell[47597, 1113, 2003, 40, 304, "Output"] }, Open ]], Cell[49615, 1156, 9087, 249, 492, "Input"], Cell[CellGroupData[{ Cell[58727, 1409, 929, 23, 72, "Input"], Cell[59659, 1434, 1831, 33, 375, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[61527, 1472, 13654, 347, 612, "Input"], Cell[75184, 1821, 2495, 50, 426, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[77716, 1876, 1448, 45, 152, "Input"], Cell[79167, 1923, 1085, 27, 116, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[80289, 1955, 1053, 33, 72, "Input"], Cell[81345, 1990, 1417, 28, 116, "Output"] }, Open ]], Cell[82777, 2021, 11028, 284, 472, "Input"], Cell[CellGroupData[{ Cell[93830, 2309, 344, 9, 52, "Input"], Cell[94177, 2320, 2083, 42, 469, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[96297, 2367, 11485, 300, 512, "Input"], Cell[107785, 2669, 1691, 31, 296, "Output"] }, Open ]] } ] *) (* End of internal cache information *)