(* 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[ 88267, 1955] NotebookOptionsPosition[ 86010, 1879] NotebookOutlinePosition[ 86551, 1899] CellTagsIndexPosition[ 86508, 1896] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"(*", " ", RowBox[{ RowBox[{"B2H6", " ", "not", " ", "TSM"}], ",", " ", RowBox[{ "overlaps", " ", "exchange", " ", "corrected", " ", "ES", " ", "08.06", ".2013"}]}], " ", "*)"}], "Text"], "\n", RowBox[{ RowBox[{ RowBox[{"Clear", "[", RowBox[{ "z", ",", "nc", ",", "R1", ",", "R2", ",", "R3", ",", "R4", ",", "wi", ",", "p", ",", "vee", ",", "vne", ",", "vnn", ",", "\[Sigma]", ",", "a1", ",", "a2", ",", "a3", ",", "b1", ",", "b2", ",", "b3", ",", "cs", ",", "ss", ",", "d1", ",", "d2", ",", "d3", ",", "d4", ",", "d5", ",", "pi", ",", "pb", ",", "i", ",", "j", ",", "d", ",", "t", ",", "q", ",", "r", ",", "r2", ",", "ow", ",", "cn", ",", "c", ",", "w", ",", "n", ",", "dwr", ",", "dwnr", ",", "dnr", ",", "c2"}], "]"}], ";"}], "\n"}]}]], "Input", CellChangeTimes->{{3.5793386253714504`*^9, 3.5793386509182944`*^9}, { 3.579338894319522*^9, 3.5793389424768066`*^9}, {3.5793389747532635`*^9, 3.5793389925684943`*^9}, {3.579339023487749*^9, 3.5793390381361747`*^9}, { 3.5793390697886305`*^9, 3.5793390710990324`*^9}, {3.579339124373126*^9, 3.5793392373953247`*^9}, {3.579339355503132*^9, 3.579339448806896*^9}, { 3.5793399816258316`*^9, 3.579339981657032*^9}, {3.579340107627253*^9, 3.5793401136020637`*^9}, 3.5793409787327833`*^9, {3.57934431334344*^9, 3.5793443261510625`*^9}, {3.5793443573823175`*^9, 3.579344360112322*^9}, 3.5793443911407766`*^9, 3.5793444273796406`*^9, 3.5793444951773596`*^9, { 3.5796661192565126`*^9, 3.579666120005314*^9}, {3.5796719837948136`*^9, 3.579672026024088*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"z", "=", "5.0"}], ";", " ", RowBox[{"nc", "=", "8"}], ";", " ", RowBox[{"pi", "=", RowBox[{"-", "0.6"}]}], ";"}], "\n", RowBox[{ RowBox[{"bohr", "=", "0.529177"}], ";", " ", RowBox[{"rad", "=", "57.29578"}], ";"}], "\n", RowBox[{ RowBox[{"k", "=", RowBox[{"{", RowBox[{ "1.0135", ",", "1.0135", ",", "1.3", ",", "1.3", ",", "1.3", ",", "1.3", ",", "1.6", ",", "1.6"}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Sigma]", "=", RowBox[{"{", RowBox[{ "0.3", ",", "0.3", ",", "0.3", ",", "0.3", ",", "0.3", ",", "0.3", ",", "0.3", ",", "0.3"}], "}"}]}], ";"}]}], "Input", CellChangeTimes->{{3.5793386253714504`*^9, 3.5793386509182944`*^9}, { 3.579338894319522*^9, 3.5793389424768066`*^9}, {3.5793389747532635`*^9, 3.5793389925684943`*^9}, {3.579339023487749*^9, 3.5793390381361747`*^9}, { 3.5793390697886305`*^9, 3.5793390710990324`*^9}, {3.579339124373126*^9, 3.5793392373953247`*^9}, {3.579339355503132*^9, 3.579339448806896*^9}, { 3.5793399816258316`*^9, 3.579339981657032*^9}, {3.579340107627253*^9, 3.5793401136020637`*^9}, {3.5793402017266183`*^9, 3.579340301675994*^9}, { 3.5793406867626705`*^9, 3.5793406954674854`*^9}, 3.579349891502223*^9, 3.579350146796671*^9, 3.579350269615687*^9, 3.579350853930313*^9, 3.5793508955043864`*^9, {3.57935093144685*^9, 3.579350935253256*^9}, { 3.5793509795105343`*^9, 3.579351010772989*^9}, {3.579351043299046*^9, 3.5793510457950506`*^9}, {3.579364450837943*^9, 3.5793644593711576`*^9}, { 3.5793651957548513`*^9, 3.57936521787569*^9}, 3.579365250464147*^9, 3.579365312505456*^9, {3.579365342831909*^9, 3.5793653574179354`*^9}, 3.5795764865231023`*^9, {3.5795768362601166`*^9, 3.579576863872165*^9}, { 3.5795769265998755`*^9, 3.5795769273174767`*^9}, 3.5795769796867685`*^9, { 3.5795773238545732`*^9, 3.579577325055775*^9}, 3.5795773626830416`*^9, 3.579577395583499*^9, 3.57957744712599*^9, {3.57957755000817*^9, 3.579577550725772*^9}, 3.579577624435901*^9, {3.5795779205712214`*^9, 3.579577944501663*^9}, {3.5795779836109324`*^9, 3.579577984032133*^9}, { 3.579578023250602*^9, 3.579578023562602*^9}, {3.5796716344166*^9, 3.5796716368034043`*^9}, {3.5796716693138614`*^9, 3.5796716747894707`*^9}, {3.5796718550945873`*^9, 3.57967189630986*^9}}], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{ RowBox[{"C", " ", "He", " ", "shell"}], ",", " ", RowBox[{"i", "=", RowBox[{"j", " ", "each"}]}]}], " ", "*)"}]], "Input", CellChangeTimes->{{3.5793386253714504`*^9, 3.5793386509182944`*^9}, { 3.579338894319522*^9, 3.5793389424768066`*^9}, {3.5793389747532635`*^9, 3.5793389925684943`*^9}, {3.579339023487749*^9, 3.5793390381361747`*^9}, { 3.5793390697886305`*^9, 3.5793390710990324`*^9}, {3.579339124373126*^9, 3.5793392373953247`*^9}, {3.579339355503132*^9, 3.579339448806896*^9}, { 3.5793399816258316`*^9, 3.579339981657032*^9}, {3.579340107627253*^9, 3.579340120341275*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"dwr", "=", RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "nc"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", "nc"}], "}"}]}], "]"}]}], ";"}], "\n", RowBox[{ RowBox[{"dnr", "=", "dwr"}], ";", " ", RowBox[{"dwnr", "=", "dwr"}], ";"}], "\n", RowBox[{ RowBox[{"c2", "=", RowBox[{"Table", "[", RowBox[{"0", ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "nc"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"(*", RowBox[{ RowBox[{"R4", "=", "R2"}], ";"}], "*)"}]}]}], "Input", CellChangeTimes->{{3.5793386253714504`*^9, 3.5793386509182944`*^9}, { 3.579338894319522*^9, 3.5793389424768066`*^9}, {3.5793389747532635`*^9, 3.5793389925684943`*^9}, {3.579339023487749*^9, 3.5793390381361747`*^9}, { 3.5793390697886305`*^9, 3.5793390710990324`*^9}, {3.579339124373126*^9, 3.5793392373953247`*^9}, {3.579339355503132*^9, 3.579339448806896*^9}, { 3.5793399816258316`*^9, 3.579339981657032*^9}, {3.579340107627253*^9, 3.579340120341275*^9}, 3.5793407052955027`*^9, 3.5793437072199755`*^9, { 3.5793454193073826`*^9, 3.579345425344593*^9}, 3.5793503939479055`*^9}], Cell[BoxData[{ RowBox[{ RowBox[{"ow", "=", RowBox[{"{", RowBox[{ RowBox[{"-", "2"}], ",", RowBox[{"-", "2"}], ",", RowBox[{"-", "2"}], ",", RowBox[{"-", "2"}], ",", RowBox[{"-", "2"}], ",", RowBox[{"-", "2"}], ",", RowBox[{"-", "2"}], ",", RowBox[{"-", "2"}]}], "}"}]}], ";"}], "\n", RowBox[{ RowBox[{"cn", "=", RowBox[{"{", RowBox[{ "z", ",", "z", ",", "1", ",", "1", ",", "1", ",", "1", ",", "1", ",", "1"}], "}"}]}], ";"}], "\n", RowBox[{ RowBox[{"c", "=", RowBox[{"{", RowBox[{ "R1", ",", "R1", ",", "R2", ",", "R2", ",", "R2", ",", "R2", ",", "R4", ",", "R4"}], "}"}]}], ";"}], "\n", RowBox[{ RowBox[{"ec", "=", RowBox[{"{", RowBox[{ "0", ",", "0", ",", "p", ",", "p", ",", "p", ",", "p", ",", "pb", ",", "pb"}], "}"}]}], ";"}], "\n", RowBox[{ RowBox[{"cs", "=", RowBox[{"Cos", "[", "wi", "]"}]}], ";", " ", RowBox[{"ss", "=", RowBox[{"Sin", "[", "wi", "]"}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"d1", "=", RowBox[{"R1", "+", "R3"}]}], ";", " ", RowBox[{"d2", "=", RowBox[{"d1", "+", RowBox[{ RowBox[{"(", RowBox[{"R1", "+", "R2", "+", "p"}], ")"}], "*", "cs"}]}]}], ";", " ", RowBox[{"d3", "=", RowBox[{ RowBox[{"(", RowBox[{"R1", "+", "R2", "+", "p"}], ")"}], "*", "ss"}]}], ";"}], " "}], "\n", RowBox[{ RowBox[{"d4", "=", RowBox[{"d1", "+", RowBox[{ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], "*", "cs"}]}]}], ";", " ", RowBox[{"d5", "=", RowBox[{ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], "*", "ss"}]}], ";", " ", RowBox[{"d6", "=", RowBox[{"pi", "+", "R4", "+", "pb"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"r", "=", RowBox[{"pi", "+", "R4"}]}], ";", " ", RowBox[{"r2", "=", RowBox[{"2", "*", "r"}]}], ";"}], "\n", RowBox[{"(*", " ", RowBox[{ RowBox[{ RowBox[{"Nuclear", " ", "positions"}], ";", " ", "1"}], ",", "2", ",", RowBox[{"3", " ", "=", " ", "x"}], ",", "y", ",", "z"}], " ", "*)"}]}], "\n", RowBox[{ RowBox[{"b1", "=", RowBox[{"{", RowBox[{ RowBox[{"-", "d1"}], ",", "d1", ",", RowBox[{"-", "d2"}], ",", RowBox[{"-", "d2"}], ",", "d2", ",", "d2", ",", "0", ",", "0"}], "}"}]}], ";"}], "\n", RowBox[{ RowBox[{"b2", "=", RowBox[{"{", RowBox[{"0", ",", "0", ",", "d3", ",", RowBox[{"-", "d3"}], ",", RowBox[{"-", "d3"}], ",", "d3", ",", "0", ",", "0"}], "}"}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"b3", "=", RowBox[{"{", RowBox[{ "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "d6", ",", RowBox[{"-", "d6"}]}], "}"}]}], ";"}], " "}], "\n", RowBox[{ RowBox[{ RowBox[{"n", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"b1", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"b2", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"b3", "[", RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "nc"}], "}"}]}], "]"}]}], ";"}], "\n", RowBox[{"(*", " ", RowBox[{"Cloud", " ", "positions"}], " ", "*)"}]}], "\n", RowBox[{ RowBox[{"a1", "=", RowBox[{"{", RowBox[{ RowBox[{"-", "d1"}], ",", "d1", ",", RowBox[{"-", "d4"}], ",", RowBox[{"-", "d4"}], ",", "d4", ",", "d4", ",", "0", ",", "0"}], "}"}]}], ";"}], "\n", RowBox[{ RowBox[{"a2", "=", RowBox[{"{", RowBox[{"0", ",", "0", ",", "d5", ",", RowBox[{"-", "d5"}], ",", RowBox[{"-", "d5"}], ",", "d5", ",", "0", ",", "0"}], "}"}]}], ";"}], "\n", RowBox[{ RowBox[{"a3", "=", RowBox[{"{", RowBox[{ "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "r", ",", RowBox[{"-", "r"}]}], "}"}]}], ";"}], "\n", RowBox[{ RowBox[{"w", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"a1", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"a2", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"a3", "[", RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "nc"}], "}"}]}], "]"}]}], ";"}], "\n"}], "Input", CellChangeTimes->{{3.5793386253714504`*^9, 3.5793386509182944`*^9}, { 3.579338894319522*^9, 3.5793389424768066`*^9}, {3.5793389747532635`*^9, 3.5793389925684943`*^9}, {3.579339023487749*^9, 3.5793390381361747`*^9}, { 3.5793390697886305`*^9, 3.5793390710990324`*^9}, {3.579339124373126*^9, 3.5793392373953247`*^9}, {3.579339355503132*^9, 3.579339448806896*^9}, { 3.5793399816258316`*^9, 3.579339981657032*^9}, {3.579340107627253*^9, 3.579340120341275*^9}, {3.5793404189258*^9, 3.57934045884627*^9}, { 3.5793407200843287`*^9, 3.579340724171536*^9}, 3.5793435809845533`*^9, 3.579343698998761*^9, {3.5793440842570376`*^9, 3.5793441144118905`*^9}, { 3.5793445016201706`*^9, 3.5793445051301765`*^9}, {3.579345216756627*^9, 3.5793452263194437`*^9}, 3.579365706499748*^9, {3.5795771074665933`*^9, 3.579577115750208*^9}, {3.5795792011306705`*^9, 3.5795792047654767`*^9}, { 3.579579999789673*^9, 3.579580007496087*^9}, {3.5796657404098473`*^9, 3.579665809705169*^9}, {3.5796660974944744`*^9, 3.5796661080400934`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"For", "[", RowBox[{ RowBox[{"i", " ", "=", " ", "1"}], ",", " ", RowBox[{"i", " ", "\[LessEqual]", " ", "nc"}], ",", " ", RowBox[{"i", "++"}], ",", "\n", " ", RowBox[{ RowBox[{ RowBox[{"c2", "[", RowBox[{"[", "i", "]"}], "]"}], " ", "=", " ", RowBox[{ RowBox[{"k", "[", RowBox[{"[", "i", "]"}], "]"}], "/", RowBox[{ RowBox[{"c", "[", RowBox[{"[", "i", "]"}], "]"}], "^", "2"}]}]}], ";", "\n", " ", RowBox[{"For", "[", RowBox[{ RowBox[{"j", " ", "=", " ", "1"}], ",", " ", RowBox[{"j", " ", "\[LessEqual]", " ", "nc"}], ",", " ", RowBox[{"j", "++"}], ",", "\n", " ", RowBox[{ RowBox[{"d", " ", "=", " ", RowBox[{"Norm", "[", RowBox[{ RowBox[{"w", "[", RowBox[{"[", "i", "]"}], "]"}], " ", "-", " ", RowBox[{"w", "[", RowBox[{"[", "j", "]"}], "]"}]}], "]"}]}], ";", "\n", " ", RowBox[{ RowBox[{"dwr", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], " ", "=", " ", RowBox[{"If", "[", RowBox[{ RowBox[{"i", " ", "\[NotEqual]", " ", "j"}], ",", " ", RowBox[{"1", "/", "d"}], ",", RowBox[{"1.5", "*", RowBox[{"(", RowBox[{ RowBox[{"\[Sigma]", "[", RowBox[{"[", "i", "]"}], "]"}], "/", RowBox[{"c", "[", RowBox[{"[", "i", "]"}], "]"}]}], ")"}]}]}], "]"}]}], ";", "\n", " ", RowBox[{"d", "=", " ", RowBox[{"Norm", "[", RowBox[{ RowBox[{"n", "[", RowBox[{"[", "i", "]"}], "]"}], " ", "-", " ", RowBox[{"n", "[", RowBox[{"[", "j", "]"}], "]"}]}], "]"}]}], ";", "\n", " ", RowBox[{ RowBox[{"dnr", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], " ", "=", " ", RowBox[{"If", "[", RowBox[{ RowBox[{"i", " ", "\[NotEqual]", " ", "j"}], ",", " ", RowBox[{"1", "/", "d"}], ",", " ", "0"}], "]"}]}], ";", "\n", " ", RowBox[{"d", "=", " ", RowBox[{"Norm", "[", RowBox[{ RowBox[{"w", "[", RowBox[{"[", "i", "]"}], "]"}], " ", "-", " ", RowBox[{"n", "[", RowBox[{"[", "j", "]"}], "]"}]}], "]"}]}], ";", "\n", " ", RowBox[{ RowBox[{"dwnr", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], " ", "=", " ", RowBox[{"If", "[", RowBox[{ RowBox[{"i", " ", "\[NotEqual]", " ", "j"}], ",", " ", RowBox[{"1", "/", "d"}], ",", RowBox[{ RowBox[{"(", RowBox[{"1.5", " ", "-", " ", RowBox[{"0.5", "*", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"ec", "[", RowBox[{"[", "i", "]"}], "]"}], "/", RowBox[{"c", "[", RowBox[{"[", "i", "]"}], "]"}]}], ")"}], "^", "2"}]}]}], ")"}], "/", RowBox[{"c", "[", RowBox[{"[", "i", "]"}], "]"}]}]}], "]"}]}]}]}], "]"}]}]}], "]"}], ";"}], " "}]], "Input", CellChangeTimes->{{3.5793386253714504`*^9, 3.5793386509182944`*^9}, { 3.579338894319522*^9, 3.5793389424768066`*^9}, {3.5793389747532635`*^9, 3.5793389925684943`*^9}, {3.579339023487749*^9, 3.5793390381361747`*^9}, { 3.5793390697886305`*^9, 3.5793390710990324`*^9}, {3.579339124373126*^9, 3.5793392373953247`*^9}, {3.579339355503132*^9, 3.579339448806896*^9}, { 3.579339981657032*^9, 3.5793399816726317`*^9}, {3.579341222280411*^9, 3.5793412545412674`*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"Ekin", " ", "=", " ", RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"9", "/", "8"}], ")"}]}], "*", RowBox[{"ow", " ", ".", " ", "c2"}]}]}], ";"}], " "}], "\n", RowBox[{ RowBox[{ RowBox[{"vee", " ", "=", " ", RowBox[{"0.5", "*", RowBox[{"ow", " ", ".", " ", RowBox[{"(", RowBox[{"ow", " ", ".", " ", "dwr"}], ")"}]}]}]}], ";"}], " "}], "\n", RowBox[{ RowBox[{"delvee", "=", RowBox[{ RowBox[{"(", RowBox[{"4", "/", "R4"}], ")"}], "*", RowBox[{"(", RowBox[{"1.2", "-", RowBox[{"0.5", "*", RowBox[{ RowBox[{"(", RowBox[{"r2", "/", "R4"}], ")"}], "^", "2"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"3", "/", "16"}], ")"}], "*", RowBox[{ RowBox[{"(", RowBox[{"r2", "/", "R4"}], ")"}], "^", "3"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{"1", "/", "160"}], ")"}], "*", RowBox[{ RowBox[{"(", RowBox[{"r2", "/", "R4"}], ")"}], "^", "5"}]}]}], ")"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"vee", "=", RowBox[{"vee", "-", RowBox[{"4", "/", "r2"}], "+", "delvee", "-", RowBox[{ RowBox[{"(", RowBox[{"2", "/", "R4"}], ")"}], "*", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "0.5"}], "*", RowBox[{ RowBox[{"(", RowBox[{"r2", "/", "R4"}], ")"}], "^", "2"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"3", "/", "16"}], ")"}], "*", RowBox[{ RowBox[{"(", RowBox[{"r2", "/", "R4"}], ")"}], "^", "3"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{"1", "/", "160"}], ")"}], "*", RowBox[{ RowBox[{"(", RowBox[{"r2", "/", "R4"}], ")"}], "^", "5"}]}]}], ")"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"vnn", " ", "=", " ", RowBox[{"0.5", "*", RowBox[{"cn", " ", ".", " ", RowBox[{"(", RowBox[{"cn", " ", ".", " ", "dnr"}], ")"}]}]}]}], ";"}], " "}], "\n", RowBox[{ RowBox[{ RowBox[{"vne", " ", "=", " ", RowBox[{"cn", " ", ".", " ", RowBox[{"(", RowBox[{"ow", " ", ".", " ", "dwnr"}], ")"}]}]}], ";"}], " "}], "\n", RowBox[{ RowBox[{"ratio", "=", RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"vne", "+", RowBox[{"2", "*", "vnn"}]}], ")"}]}], "/", RowBox[{"(", RowBox[{"vee", "-", "vnn"}], ")"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Epot", " ", "=", " ", RowBox[{"vne", " ", "+", " ", "vee", " ", "+", " ", "vnn"}]}], ";"}], " "}], "\n", RowBox[{ RowBox[{"vir", "=", RowBox[{ RowBox[{"-", "Epot"}], "/", "Ekin"}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"func", " ", "=", " ", RowBox[{"Ekin", " ", "+", " ", "Epot"}]}], ";"}], " "}]}], "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.5793386253714504`*^9, 3.5793386509182944`*^9}, { 3.579338894319522*^9, 3.5793389424768066`*^9}, {3.5793389747532635`*^9, 3.5793389925684943`*^9}, {3.579339023487749*^9, 3.5793390381361747`*^9}, { 3.5793390697886305`*^9, 3.5793390710990324`*^9}, {3.579339124373126*^9, 3.5793392373953247`*^9}, {3.579339355503132*^9, 3.579339506199397*^9}, { 3.579339540051456*^9, 3.5793395413774586`*^9}, {3.5793395728427134`*^9, 3.579339683743308*^9}, {3.5793399816726317`*^9, 3.5793400106886826`*^9}, 3.579340927377493*^9, 3.5793453744885035`*^9, {3.579666680935499*^9, 3.5796668175293393`*^9}, {3.5796668656242237`*^9, 3.5796669598951893`*^9}, {3.5796671065042467`*^9, 3.5796671286562853`*^9}, {3.5800573498934617`*^9, 3.580057375604932*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"t", " ", "=", " ", RowBox[{"NMinimize", "[", RowBox[{ RowBox[{"{", RowBox[{"func", ",", RowBox[{ RowBox[{ RowBox[{ RowBox[{"2", "*", "d1"}], "-", "3.3"}], " ", "\[Equal]", "0"}], " ", "&&", " ", RowBox[{ RowBox[{"R1", "+", "R2", "+", "p"}], " ", "\[Equal]", " ", "2.249"}], " ", "&&", " ", RowBox[{ RowBox[{"d1", "/", "d6"}], "\[Equal]", "0.89946"}], " ", "&&", " ", RowBox[{"wi", "\[Equal]", "1.062"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"R1", ",", "0.31", ",", "0.332"}], "}"}], ",", RowBox[{"{", RowBox[{"R2", ",", "1.42", ",", "1.45"}], "}"}], ",", RowBox[{"{", RowBox[{"R3", ",", "1.3", ",", "1.35"}], "}"}], ",", RowBox[{"{", RowBox[{"R4", ",", "1.7", ",", "1.90"}], "}"}], ",", RowBox[{"{", RowBox[{"p", ",", "0.49", ",", "0.61"}], "}"}], ",", RowBox[{"{", RowBox[{"pb", ",", "0.58", ",", "0.65"}], "}"}], ",", RowBox[{"{", RowBox[{"wi", ",", "1.05", ",", "1.1"}], "}"}]}], "}"}]}], "]"}]}]], "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.5793386253714504`*^9, 3.5793386509182944`*^9}, { 3.579338894319522*^9, 3.5793389424768066`*^9}, {3.5793389747532635`*^9, 3.5793389925684943`*^9}, {3.579339023487749*^9, 3.5793390381361747`*^9}, { 3.5793390697886305`*^9, 3.5793390710990324`*^9}, {3.579339124373126*^9, 3.5793392373953247`*^9}, {3.579339355503132*^9, 3.579339506199397*^9}, { 3.579339540051456*^9, 3.5793395413774586`*^9}, {3.5793395728427134`*^9, 3.579339683743308*^9}, {3.5793399816726317`*^9, 3.5793400106886826`*^9}, 3.579340927377493*^9, {3.5793433833366065`*^9, 3.5793433837422075`*^9}, { 3.579344512243789*^9, 3.57934455821707*^9}, 3.579344675232876*^9, { 3.579344794885086*^9, 3.5793447997678947`*^9}, {3.5793452439006743`*^9, 3.579345253806692*^9}, {3.5793453515408635`*^9, 3.579345356361272*^9}, 3.579345443300225*^9, {3.579350403417122*^9, 3.579350409579133*^9}, { 3.579365415216037*^9, 3.5793654535453043`*^9}, {3.5794458283765016`*^9, 3.579445830700906*^9}, {3.5794458648805656`*^9, 3.5794458940682173`*^9}, { 3.579446370758254*^9, 3.5794464147659316`*^9}, {3.5794466902156153`*^9, 3.579446734426093*^9}, {3.5794468710979333`*^9, 3.579446890207967*^9}, 3.579576260743906*^9, {3.5795763810981174`*^9, 3.5795764136085744`*^9}, 3.579576825995299*^9, {3.579580198721223*^9, 3.5795802058504353`*^9}, { 3.579662269871752*^9, 3.5796623591663084`*^9}, {3.5796659620706367`*^9, 3.579665967827047*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"-", "52.99381989748658`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"R1", "\[Rule]", "0.3183387910180002`"}], ",", RowBox[{"R2", "\[Rule]", "1.438899692083698`"}], ",", RowBox[{"R3", "\[Rule]", "1.3316612089819997`"}], ",", RowBox[{"R4", "\[Rule]", "1.8271408028181948`"}], ",", RowBox[{"p", "\[Rule]", "0.49176151689830194`"}], ",", RowBox[{"pb", "\[Rule]", "0.6072931955188099`"}], ",", RowBox[{"wi", "\[Rule]", "1.062`"}]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{ 3.57957650798874*^9, 3.579576869394575*^9, 3.579576933510688*^9, 3.5795769848035774`*^9, 3.5795771258122253`*^9, 3.5795773313425865`*^9, 3.5795773667390485`*^9, 3.579577400606708*^9, 3.5795774517903976`*^9, 3.5795775558269806`*^9, 3.5795776291003094`*^9, 3.5795776737007875`*^9, { 3.57957792559443*^9, 3.5795779492908716`*^9}, 3.5795779886497407`*^9, 3.57957803343742*^9, 3.579579210085086*^9, 3.579579697383142*^9, 3.5795797319372025`*^9, 3.5795798842246704`*^9, 3.5795800145784993`*^9, 3.579580212121646*^9, 3.5795802717605505`*^9, 3.5795803440822783`*^9, 3.579662366701122*^9, 3.5796626466904135`*^9, {3.5796626963453007`*^9, 3.5796627140357323`*^9}, 3.5796658674564705`*^9, 3.579665987327081*^9, 3.579666028261553*^9, 3.579666999768859*^9, {3.579667113290259*^9, 3.5796671400755057`*^9}, 3.579671642981015*^9, 3.5796716811386824`*^9, 3.579671866560608*^9, 3.57967190203507*^9, {3.5796834754780064`*^9, 3.579683524680493*^9}, 3.5796835857390003`*^9, 3.5796837451556797`*^9, { 3.579683784233749*^9, 3.579683837913443*^9}, 3.579683922153591*^9, 3.579683968251672*^9, 3.5800574322331715`*^9, 3.5859073859597178`*^9, { 3.585907426003008*^9, 3.5859074504384055`*^9}, {3.585907570364265*^9, 3.5859075933605804`*^9}, 3.585907754221781*^9, 3.585907802952568*^9, 3.5859078753377085`*^9, 3.585907914804966*^9, {3.585907987503124*^9, 3.5859080145406704`*^9}, 3.5859081623461246`*^9, 3.5859092595369782`*^9, 3.5859092927788796`*^9, 3.5859094245634174`*^9, {3.5859095253401814`*^9, 3.585909554296838*^9}, 3.5859096401097455`*^9, 3.5859101671068883`*^9, 3.5859102090952897`*^9, 3.5859102595811777`*^9, 3.5859103240788665`*^9, 3.5859104487289963`*^9, {3.5859105536849995`*^9, 3.5859105729731026`*^9}, { 3.5859106320314803`*^9, 3.5859106523516426`*^9}, 3.585910817337079*^9, { 3.5859109613063135`*^9, 3.585910984821659*^9}, {3.585911016018443*^9, 3.585911032328376*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{"vee", " ", "/.", " ", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}], "\n", RowBox[{"vne", " ", "/.", " ", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}], "\n", RowBox[{"vnn", " ", "/.", " ", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}]}], "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.5793386253714504`*^9, 3.5793386509182944`*^9}, { 3.579338894319522*^9, 3.5793389424768066`*^9}, {3.5793389747532635`*^9, 3.5793389925684943`*^9}, {3.579339023487749*^9, 3.5793390381361747`*^9}, { 3.5793390697886305`*^9, 3.5793390710990324`*^9}, {3.579339124373126*^9, 3.5793392373953247`*^9}, {3.579339355503132*^9, 3.579339506199397*^9}, { 3.579339540051456*^9, 3.5793395413774586`*^9}, {3.5793395728427134`*^9, 3.579339683743308*^9}, {3.5793399816726317`*^9, 3.5793400106886826`*^9}, 3.579340927377493*^9, {3.5793433833366065`*^9, 3.5793433837422075`*^9}, { 3.579344512243789*^9, 3.57934455821707*^9}, 3.579344675232876*^9, { 3.579344794885086*^9, 3.5793447997678947`*^9}, {3.5793452439006743`*^9, 3.579345253806692*^9}, 3.5793453515408635`*^9, {3.5793504919472775`*^9, 3.5793505047393*^9}, {3.5793505769830265`*^9, 3.5793506236271086`*^9}, { 3.579350654125162*^9, 3.579350657697569*^9}, {3.579350819657053*^9, 3.5793508278626676`*^9}, {3.5793645321452856`*^9, 3.579364610285823*^9}, { 3.579364652530697*^9, 3.5793647063663917`*^9}, {3.580057397236169*^9, 3.580057418627393*^9}}], Cell[BoxData["48.57075260983609`"], "Output", CellChangeTimes->{ 3.5800574322501717`*^9, 3.5859073859717183`*^9, {3.585907426015009*^9, 3.585907450449406*^9}, {3.5859075703752656`*^9, 3.585907593371581*^9}, 3.585907754233782*^9, 3.5859078029635687`*^9, 3.585907875348709*^9, 3.585907914816966*^9, {3.5859079875151243`*^9, 3.5859080145526714`*^9}, 3.585908162358125*^9, 3.585909259548979*^9, 3.58590929279088*^9, 3.585909424574418*^9, {3.585909525352182*^9, 3.585909554307838*^9}, 3.585909640121746*^9, 3.585910167118889*^9, 3.5859102091072903`*^9, 3.585910259593178*^9, 3.585910324090867*^9, 3.585910448741997*^9, { 3.585910553696*^9, 3.585910572985103*^9}, {3.585910632041481*^9, 3.5859106523626432`*^9}, 3.58591081735108*^9, {3.5859109613153143`*^9, 3.5859109848356595`*^9}, {3.585911016032444*^9, 3.5859110323423767`*^9}}], Cell[BoxData[ RowBox[{"-", "186.6300480447529`"}]], "Output", CellChangeTimes->{ 3.5800574322501717`*^9, 3.5859073859717183`*^9, {3.585907426015009*^9, 3.585907450449406*^9}, {3.5859075703752656`*^9, 3.585907593371581*^9}, 3.585907754233782*^9, 3.5859078029635687`*^9, 3.585907875348709*^9, 3.585907914816966*^9, {3.5859079875151243`*^9, 3.5859080145526714`*^9}, 3.585908162358125*^9, 3.585909259548979*^9, 3.58590929279088*^9, 3.585909424574418*^9, {3.585909525352182*^9, 3.585909554307838*^9}, 3.585909640121746*^9, 3.585910167118889*^9, 3.5859102091072903`*^9, 3.585910259593178*^9, 3.585910324090867*^9, 3.585910448741997*^9, { 3.585910553696*^9, 3.585910572985103*^9}, {3.585910632041481*^9, 3.5859106523626432`*^9}, 3.58591081735108*^9, {3.5859109613153143`*^9, 3.5859109848356595`*^9}, {3.585911016032444*^9, 3.585911032344377*^9}}], Cell[BoxData["32.253163762256705`"], "Output", CellChangeTimes->{ 3.5800574322501717`*^9, 3.5859073859717183`*^9, {3.585907426015009*^9, 3.585907450449406*^9}, {3.5859075703752656`*^9, 3.585907593371581*^9}, 3.585907754233782*^9, 3.5859078029635687`*^9, 3.585907875348709*^9, 3.585907914816966*^9, {3.5859079875151243`*^9, 3.5859080145526714`*^9}, 3.585908162358125*^9, 3.585909259548979*^9, 3.58590929279088*^9, 3.585909424574418*^9, {3.585909525352182*^9, 3.585909554307838*^9}, 3.585909640121746*^9, 3.585910167118889*^9, 3.5859102091072903`*^9, 3.585910259593178*^9, 3.585910324090867*^9, 3.585910448741997*^9, { 3.585910553696*^9, 3.585910572985103*^9}, {3.585910632041481*^9, 3.5859106523626432`*^9}, 3.58591081735108*^9, {3.5859109613153143`*^9, 3.5859109848356595`*^9}, {3.585911016032444*^9, 3.585911032345377*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"ratio", " ", "/.", " ", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}]], "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.5793386253714504`*^9, 3.5793386509182944`*^9}, { 3.579338894319522*^9, 3.5793389424768066`*^9}, {3.5793389747532635`*^9, 3.5793389925684943`*^9}, {3.579339023487749*^9, 3.5793390381361747`*^9}, { 3.5793390697886305`*^9, 3.5793390710990324`*^9}, {3.579339124373126*^9, 3.5793392373953247`*^9}, {3.579339355503132*^9, 3.579339506199397*^9}, { 3.579339540051456*^9, 3.5793395413774586`*^9}, {3.5793395728427134`*^9, 3.579339683743308*^9}, {3.5793399816726317`*^9, 3.5793400106886826`*^9}, 3.579340927377493*^9, {3.5793433833366065`*^9, 3.5793433837422075`*^9}, { 3.579344512243789*^9, 3.57934455821707*^9}, 3.579344675232876*^9, { 3.579344794885086*^9, 3.5793447997678947`*^9}, {3.5793452439006743`*^9, 3.579345253806692*^9}, 3.5793453515408635`*^9, {3.5793504919472775`*^9, 3.5793505047393*^9}, {3.5793505769830265`*^9, 3.5793506236271086`*^9}, { 3.579350654125162*^9, 3.579350657697569*^9}, {3.579350819657053*^9, 3.5793508278626676`*^9}, {3.5793645321452856`*^9, 3.579364610285823*^9}, { 3.579364652530697*^9, 3.5793647063663917`*^9}, {3.580057397236169*^9, 3.580057423139651*^9}}], Cell[BoxData["7.484176838930208`"], "Output", CellChangeTimes->{ 3.580057432268173*^9, 3.585907385983719*^9, {3.58590742602901*^9, 3.585907450462407*^9}, {3.5859075703882666`*^9, 3.585907593386582*^9}, 3.5859077542467823`*^9, 3.5859078029765697`*^9, 3.5859078753657103`*^9, 3.5859079148329678`*^9, {3.5859079875281253`*^9, 3.585908014567672*^9}, 3.585908162373126*^9, 3.5859092595619793`*^9, 3.5859092928038807`*^9, 3.5859094245884185`*^9, {3.5859095253671827`*^9, 3.585909554322839*^9}, 3.585909640135747*^9, 3.5859101671328897`*^9, 3.585910209125292*^9, 3.5859102596111794`*^9, 3.585910324104868*^9, 3.5859104487559977`*^9, { 3.585910553708*^9, 3.585910572999104*^9}, {3.5859106320544815`*^9, 3.585910652376644*^9}, 3.585910817370081*^9, {3.5859109613283153`*^9, 3.585910984856661*^9}, {3.5859110160484447`*^9, 3.585911032362378*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"vir", " ", "/.", " ", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}]], "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.5793386253714504`*^9, 3.5793386509182944`*^9}, { 3.579338894319522*^9, 3.5793389424768066`*^9}, {3.5793389747532635`*^9, 3.5793389925684943`*^9}, {3.579339023487749*^9, 3.5793390381361747`*^9}, { 3.5793390697886305`*^9, 3.5793390710990324`*^9}, {3.579339124373126*^9, 3.5793392373953247`*^9}, {3.579339355503132*^9, 3.579339506199397*^9}, { 3.579339540051456*^9, 3.5793395413774586`*^9}, {3.5793395728427134`*^9, 3.579339683743308*^9}, {3.5793399816726317`*^9, 3.5793400106886826`*^9}, 3.579340927377493*^9, {3.5793433833366065`*^9, 3.5793433837422075`*^9}, { 3.579344512243789*^9, 3.57934455821707*^9}, 3.579344675232876*^9, { 3.579344794885086*^9, 3.5793447997678947`*^9}, {3.5793452439006743`*^9, 3.579345253806692*^9}, 3.5793453515408635`*^9, {3.5793504919472775`*^9, 3.5793505047393*^9}, {3.5793505769830265`*^9, 3.5793506236271086`*^9}, { 3.579350654125162*^9, 3.579350657697569*^9}, {3.579350819657053*^9, 3.5793508278626676`*^9}, {3.5793645321452856`*^9, 3.579364610285823*^9}, { 3.579364652530697*^9, 3.5793647063663917`*^9}, {3.580057397236169*^9, 3.580057418627393*^9}}], Cell[BoxData["2.003436852434821`"], "Output", CellChangeTimes->{ 3.57957650798874*^9, 3.5795768694101753`*^9, 3.5795769335262876`*^9, 3.5795769848191776`*^9, 3.5795771258278255`*^9, 3.5795773313581867`*^9, 3.5795773667546487`*^9, 3.579577400622308*^9, 3.579577451805998*^9, 3.5795775558269806`*^9, 3.5795776291159096`*^9, 3.5795776737007875`*^9, { 3.57957792561003*^9, 3.5795779492908716`*^9}, 3.579577988665341*^9, 3.57957803343742*^9, 3.579579210100686*^9, 3.579579697383142*^9, 3.5795797319528027`*^9, 3.5795798842402706`*^9, 3.579580014594099*^9, 3.579580212137246*^9, 3.5795802717761507`*^9, 3.579580344097878*^9, 3.5796623667167215`*^9, 3.5796626466904135`*^9, {3.579662696360901*^9, 3.5796627140357323`*^9}, 3.5796658674720707`*^9, 3.579665987327081*^9, 3.579666028277153*^9, 3.579666999784459*^9, {3.5796671133058586`*^9, 3.5796671400911055`*^9}, 3.579671642981015*^9, 3.579671681154282*^9, 3.579671866576208*^9, 3.57967190205067*^9, {3.5796834754936066`*^9, 3.5796835246960926`*^9}, 3.5796835857390003`*^9, 3.57968374517128*^9, { 3.5796837842493486`*^9, 3.579683837913443*^9}, 3.579683922153591*^9, 3.579683968251672*^9, 3.580057432283174*^9, 3.5859073859947195`*^9, { 3.58590742604001*^9, 3.5859074504734077`*^9}, {3.585907570403267*^9, 3.585907593397582*^9}, 3.585907754257783*^9, 3.5859078029915705`*^9, 3.585907875376711*^9, 3.585907914843968*^9, {3.585907987543126*^9, 3.5859080145786724`*^9}, 3.585908162387127*^9, 3.5859092595769806`*^9, 3.5859092928178816`*^9, 3.585909424599419*^9, {3.5859095253781834`*^9, 3.5859095543328395`*^9}, 3.5859096401507483`*^9, 3.5859101671438904`*^9, 3.5859102091412926`*^9, 3.58591025962718*^9, 3.585910324117869*^9, 3.5859104487699986`*^9, {3.585910553722001*^9, 3.585910573014105*^9}, { 3.5859106320644827`*^9, 3.5859106523866444`*^9}, 3.5859108173860817`*^9, { 3.5859109613373156`*^9, 3.585910984878662*^9}, {3.5859110160624456`*^9, 3.5859110323793793`*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"2", "*", "d1", "*", "bohr"}], " ", "/.", " ", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}], " ", RowBox[{"(*", " ", RowBox[{"B", "-", RowBox[{"B", " ", "1.753"}]}], " ", "*)"}]}]], "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.5793386253714504`*^9, 3.5793386509182944`*^9}, { 3.579338894319522*^9, 3.5793389424768066`*^9}, {3.5793389747532635`*^9, 3.5793389925684943`*^9}, {3.579339023487749*^9, 3.5793390381361747`*^9}, { 3.5793390697886305`*^9, 3.5793390710990324`*^9}, {3.579339124373126*^9, 3.5793392373953247`*^9}, {3.579339355503132*^9, 3.579339506199397*^9}, { 3.579339540051456*^9, 3.5793395413774586`*^9}, {3.5793395728427134`*^9, 3.579339683743308*^9}, {3.5793399816726317`*^9, 3.5793400106886826`*^9}, 3.579340927377493*^9, {3.5793433833366065`*^9, 3.5793433837422075`*^9}, { 3.579344512243789*^9, 3.57934455821707*^9}, 3.579344675232876*^9, { 3.579344794885086*^9, 3.5793447997678947`*^9}, {3.5793452439006743`*^9, 3.579345253806692*^9}, 3.5793453515408635`*^9, {3.5793504919472775`*^9, 3.5793505047393*^9}, {3.5793505769830265`*^9, 3.5793506236271086`*^9}, { 3.579350654125162*^9, 3.579350657697569*^9}, {3.579350819657053*^9, 3.5793508278626676`*^9}, {3.5793645321452856`*^9, 3.579364610285823*^9}, { 3.579364652530697*^9, 3.5793647063663917`*^9}, 3.579364778485318*^9, { 3.579364848888242*^9, 3.579364851321846*^9}}], Cell[BoxData["1.7462841`"], "Output", CellChangeTimes->{ 3.5795765080199404`*^9, 3.5795768694257755`*^9, 3.5795769335418873`*^9, 3.579576984834778*^9, 3.5795771258434258`*^9, 3.5795773313737864`*^9, 3.5795773667702484`*^9, 3.579577400637908*^9, 3.579577451821598*^9, 3.579577555842581*^9, 3.5795776291159096`*^9, 3.5795776737163877`*^9, { 3.5795779256256304`*^9, 3.579577949306472*^9}, 3.579577988665341*^9, 3.5795780334530196`*^9, 3.5795792101162863`*^9, 3.5795796973987417`*^9, 3.579579731968403*^9, 3.5795798842558703`*^9, 3.5795800146096992`*^9, 3.5795802121528463`*^9, 3.579580271791751*^9, 3.5795803441134777`*^9, 3.5796623667323217`*^9, 3.579662646721614*^9, {3.579662696376501*^9, 3.5796627140669317`*^9}, 3.579665867487671*^9, 3.579665987358281*^9, 3.579666028292753*^9, 3.5796669998000593`*^9, {3.579667113337059*^9, 3.5796671401067057`*^9}, 3.5796716429966154`*^9, 3.5796716811698823`*^9, 3.5796718665918083`*^9, 3.57967190206627*^9, {3.5796834755092063`*^9, 3.579683524711693*^9}, 3.5796835857546*^9, 3.57968374518688*^9, { 3.579683784264949*^9, 3.579683837929043*^9}, 3.579683922169191*^9, 3.579683968267272*^9, 3.5800574322951746`*^9, 3.5859073860087204`*^9, { 3.585907426051011*^9, 3.5859074504844084`*^9}, {3.5859075704142675`*^9, 3.585907593411583*^9}, 3.585907754272784*^9, 3.5859078030065713`*^9, 3.5859078753917117`*^9, 3.585907914859969*^9, {3.585907987553127*^9, 3.5859080145936737`*^9}, 3.5859081623981276`*^9, 3.5859092595939817`*^9, 3.585909292828882*^9, 3.58590942461342*^9, {3.585909525390184*^9, 3.5859095543478403`*^9}, 3.585909640166749*^9, 3.585910167158891*^9, 3.585910209153293*^9, 3.585910259643181*^9, 3.5859103241298695`*^9, 3.5859104487869997`*^9, {3.5859105537320023`*^9, 3.5859105730251055`*^9}, { 3.585910632075483*^9, 3.585910652397645*^9}, 3.585910817403083*^9, { 3.5859109613513165`*^9, 3.585910984894663*^9}, {3.5859110160784464`*^9, 3.58591103239538*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"R1", "+", "R2", "+", "p"}], ")"}], "*", "bohr"}], " ", "/.", " ", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}], " ", RowBox[{"(*", " ", RowBox[{"B", "-", RowBox[{"H2", " ", "1.19"}]}], " ", "*)"}]}]], "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.5793386253714504`*^9, 3.5793386509182944`*^9}, { 3.579338894319522*^9, 3.5793389424768066`*^9}, {3.5793389747532635`*^9, 3.5793389925684943`*^9}, {3.579339023487749*^9, 3.5793390381361747`*^9}, { 3.5793390697886305`*^9, 3.5793390710990324`*^9}, {3.579339124373126*^9, 3.5793392373953247`*^9}, {3.579339355503132*^9, 3.579339506199397*^9}, { 3.579339540051456*^9, 3.5793395413774586`*^9}, {3.5793395728427134`*^9, 3.579339683743308*^9}, {3.5793399816726317`*^9, 3.5793400106886826`*^9}, 3.579340927377493*^9, {3.5793433833366065`*^9, 3.5793433837422075`*^9}, { 3.579344512243789*^9, 3.57934455821707*^9}, 3.579344675232876*^9, { 3.579344794885086*^9, 3.5793447997678947`*^9}, {3.5793452439006743`*^9, 3.579345253806692*^9}, 3.5793453515408635`*^9, {3.5793504919472775`*^9, 3.5793505047393*^9}, {3.5793505769830265`*^9, 3.5793506236271086`*^9}, { 3.579350654125162*^9, 3.579350657697569*^9}, {3.579350819657053*^9, 3.5793508278626676`*^9}, {3.5793645321452856`*^9, 3.579364610285823*^9}, { 3.579364652530697*^9, 3.579364711623601*^9}, 3.5793647825569253`*^9, { 3.579364856142255*^9, 3.5793648764846907`*^9}, {3.579365039161776*^9, 3.57936504123658*^9}}], Cell[BoxData["1.190119073`"], "Output", CellChangeTimes->{ 3.5795765080199404`*^9, 3.5795768694257755`*^9, 3.5795769335418873`*^9, 3.579576984850378*^9, 3.5795771258434258`*^9, 3.5795773313737864`*^9, 3.5795773667858486`*^9, 3.579577400653508*^9, 3.579577451821598*^9, 3.579577555858181*^9, 3.5795776291315093`*^9, 3.579577673731988*^9, { 3.57957792564123*^9, 3.579577949322072*^9}, 3.5795779886809406`*^9, 3.5795780334686193`*^9, 3.5795792101162863`*^9, 3.579579697414342*^9, 3.579579731984003*^9, 3.57957988427147*^9, 3.5795800146252995`*^9, 3.5795802121684465`*^9, 3.579580271807351*^9, 3.579580344129078*^9, 3.579662366747922*^9, 3.579662646721614*^9, {3.579662696392101*^9, 3.5796627140669317`*^9}, 3.5796658675032706`*^9, 3.5796659873738813`*^9, 3.5796660283083534`*^9, 3.5796669998156595`*^9, {3.579667113337059*^9, 3.579667140122306*^9}, 3.579671643012215*^9, 3.579671681185482*^9, 3.579671866607408*^9, 3.57967190206627*^9, {3.5796834755248065`*^9, 3.579683524711693*^9}, 3.5796835857702*^9, 3.57968374518688*^9, { 3.579683784264949*^9, 3.579683837944643*^9}, 3.579683922184791*^9, 3.579683968282872*^9, 3.5800574323121758`*^9, 3.585907386018721*^9, { 3.585907426063011*^9, 3.585907450495409*^9}, {3.5859075704272685`*^9, 3.5859075934225836`*^9}, 3.585907754283785*^9, 3.585907803017572*^9, 3.585907875402712*^9, 3.58590791487497*^9, {3.5859079875681276`*^9, 3.585908014603674*^9}, 3.585908162413128*^9, 3.585909259604982*^9, 3.5859092928438835`*^9, 3.5859094246294208`*^9, {3.585909525401185*^9, 3.585909554358841*^9}, 3.5859096401777496`*^9, 3.585910167173892*^9, 3.585910209166294*^9, 3.585910259660182*^9, 3.5859103241448703`*^9, 3.5859104488030005`*^9, {3.585910553746003*^9, 3.5859105730391064`*^9}, { 3.5859106320854836`*^9, 3.5859106524076457`*^9}, 3.585910817420084*^9, { 3.5859109613613167`*^9, 3.585910984911664*^9}, {3.5859110160954475`*^9, 3.5859110324123807`*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"Sqrt", "[", RowBox[{ RowBox[{"d1", "^", "2"}], "+", RowBox[{"d6", "^", "2"}]}], "]"}], "*", "bohr"}], " ", "/.", " ", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}], " ", RowBox[{"(*", " ", RowBox[{"B", "-", RowBox[{"H1", " ", "1.31"}]}], " ", "*)"}], " "}]], "Input", CellChangeTimes->{{3.5793649976657033`*^9, 3.5793650337017665`*^9}, { 3.579446074950535*^9, 3.5794460764637375`*^9}, {3.579580323677442*^9, 3.5795803312434554`*^9}, {3.5796660029427085`*^9, 3.5796660088551188`*^9}}], Cell[BoxData["1.3056468628891476`"], "Output", CellChangeTimes->{ 3.5795765080355406`*^9, 3.5795768694413753`*^9, 3.5795769335574875`*^9, 3.579576984850378*^9, 3.579577125859026*^9, 3.579577331389386*^9, 3.5795773667858486`*^9, 3.579577400653508*^9, 3.5795774518371983`*^9, 3.579577555858181*^9, 3.5795776291471095`*^9, 3.579577673731988*^9, { 3.57957792564123*^9, 3.5795779493376718`*^9}, 3.579577988696541*^9, 3.5795780334842196`*^9, 3.5795792101318865`*^9, 3.579579697429942*^9, 3.579579731984003*^9, 3.57957988427147*^9, 3.579580014640899*^9, 3.5795802121840467`*^9, 3.579580271807351*^9, 3.579580344144678*^9, 3.579662366763522*^9, 3.5796626467372136`*^9, {3.5796626964077005`*^9, 3.579662714082532*^9}, 3.5796658675032706`*^9, 3.5796659873894815`*^9, 3.579666028323953*^9, 3.5796669998312597`*^9, {3.5796671133526587`*^9, 3.579667140137906*^9}, 3.579671643027815*^9, 3.579671681185482*^9, 3.579671866607408*^9, 3.57967190208187*^9, {3.579683475540406*^9, 3.579683524727293*^9}, 3.5796835857702*^9, 3.5796837452024803`*^9, { 3.579683784280549*^9, 3.579683837960243*^9}, 3.579683922200391*^9, 3.579683968298472*^9, 3.5800574323291764`*^9, 3.585907386033722*^9, { 3.585907426074012*^9, 3.5859074505064096`*^9}, {3.585907570438269*^9, 3.585907593437585*^9}, 3.5859077542967854`*^9, 3.585907803031573*^9, 3.585907875415713*^9, 3.58590791488597*^9, {3.5859079875791283`*^9, 3.585908014617675*^9}, 3.585908162429129*^9, 3.585909259617983*^9, 3.5859092928598843`*^9, 3.585909424640422*^9, {3.5859095254121857`*^9, 3.585909554371842*^9}, 3.58590964019075*^9, 3.5859101671848927`*^9, 3.5859102091782947`*^9, 3.585910259677183*^9, 3.5859103241608715`*^9, 3.585910448820001*^9, {3.5859105537560034`*^9, 3.585910573054107*^9}, { 3.585910632096484*^9, 3.5859106524186463`*^9}, 3.585910817437085*^9, { 3.5859109613703175`*^9, 3.5859109849276648`*^9}, {3.5859110161124487`*^9, 3.585911032429382*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"2", "*", "wi", "*", "rad"}], " ", "/.", " ", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}], " ", RowBox[{"(*", " ", RowBox[{"H2BH2", " ", "121.7"}], " ", "*)"}]}]], "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.5793386253714504`*^9, 3.5793386509182944`*^9}, { 3.579338894319522*^9, 3.5793389424768066`*^9}, {3.5793389747532635`*^9, 3.5793389925684943`*^9}, {3.579339023487749*^9, 3.5793390381361747`*^9}, { 3.5793390697886305`*^9, 3.5793390710990324`*^9}, {3.579339124373126*^9, 3.5793392373953247`*^9}, {3.579339355503132*^9, 3.579339506199397*^9}, { 3.579339540051456*^9, 3.5793395413774586`*^9}, {3.5793395728427134`*^9, 3.579339683743308*^9}, {3.5793399816726317`*^9, 3.5793400106886826`*^9}, 3.579340927377493*^9, {3.5793433833366065`*^9, 3.5793433837422075`*^9}, { 3.579344512243789*^9, 3.57934455821707*^9}, 3.579344675232876*^9, { 3.579344794885086*^9, 3.5793447997678947`*^9}, {3.5793452439006743`*^9, 3.579345253806692*^9}, 3.5793453515408635`*^9, {3.5793504919472775`*^9, 3.5793505047393*^9}, {3.5793505769830265`*^9, 3.5793506236271086`*^9}, { 3.579350654125162*^9, 3.579350657697569*^9}, {3.579350819657053*^9, 3.5793508278626676`*^9}, {3.5793645321452856`*^9, 3.579364610285823*^9}, { 3.579364652530697*^9, 3.579364717567211*^9}, {3.579364755553278*^9, 3.57936478491253*^9}, {3.579364865486671*^9, 3.5793648811646986`*^9}, { 3.579365054418603*^9, 3.5793650555730047`*^9}}], Cell[BoxData["121.69623672`"], "Output", CellChangeTimes->{ 3.5795765080511403`*^9, 3.579576869456975*^9, 3.5795769335730877`*^9, 3.579576984865978*^9, 3.5795771258746257`*^9, 3.5795773314049864`*^9, 3.579577366801449*^9, 3.579577400669108*^9, 3.579577451852798*^9, 3.5795775558737807`*^9, 3.5795776291627097`*^9, 3.5795776737475877`*^9, { 3.5795779256568303`*^9, 3.5795779493376718`*^9}, 3.579577988712141*^9, 3.5795780334842196`*^9, 3.5795792101474867`*^9, 3.5795796974455423`*^9, 3.579579731999603*^9, 3.5795798842870703`*^9, 3.579580014640899*^9, 3.5795802121996465`*^9, 3.5795802718229513`*^9, 3.5795803441602783`*^9, 3.5796623667791224`*^9, 3.5796626467528133`*^9, {3.5796626964233007`*^9, 3.579662714098132*^9}, 3.579665867518871*^9, 3.5796659873894815`*^9, 3.579666028339553*^9, 3.57966699984686*^9, {3.5796671133682585`*^9, 3.5796671401535063`*^9}, 3.579671643027815*^9, 3.579671681201082*^9, 3.579671866623008*^9, 3.5796719020974703`*^9, {3.579683475540406*^9, 3.579683524742893*^9}, 3.5796835857858*^9, 3.5796837452180805`*^9, { 3.579683784296149*^9, 3.5796838379758434`*^9}, 3.579683922215991*^9, 3.5796839683140717`*^9, 3.5800574323461776`*^9, 3.585907386043722*^9, { 3.585907426085013*^9, 3.58590745051641*^9}, {3.58590757045327*^9, 3.5859075934545856`*^9}, 3.585907754307786*^9, 3.5859078030415735`*^9, 3.585907875426714*^9, 3.585907914899971*^9, {3.585907987593129*^9, 3.5859080146286755`*^9}, 3.5859081624401298`*^9, 3.5859092596289835`*^9, 3.585909292870885*^9, 3.5859094246534224`*^9, {3.5859095254231863`*^9, 3.585909554382842*^9}, 3.5859096402017508`*^9, 3.5859101671988935`*^9, 3.5859102091932955`*^9, 3.585910259693184*^9, 3.5859103241768723`*^9, 3.5859104488370028`*^9, {3.5859105537660036`*^9, 3.5859105730651073`*^9}, { 3.585910632106485*^9, 3.585910652428647*^9}, 3.5859108174540863`*^9, { 3.585910961378318*^9, 3.585910984944666*^9}, {3.5859110161284494`*^9, 3.5859110324463825`*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"2", "*", RowBox[{"ArcTan", "[", RowBox[{"d1", "/", "d6"}], "]"}], "*", "rad"}], " ", "/.", " ", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}], " ", RowBox[{"(*", " ", RowBox[{"BH1B", " ", "83.94"}], " ", "*)"}], " "}]], "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.5793386253714504`*^9, 3.5793386509182944`*^9}, { 3.579338894319522*^9, 3.5793389424768066`*^9}, {3.5793389747532635`*^9, 3.5793389925684943`*^9}, {3.579339023487749*^9, 3.5793390381361747`*^9}, { 3.5793390697886305`*^9, 3.5793390710990324`*^9}, {3.579339124373126*^9, 3.5793392373953247`*^9}, {3.579339355503132*^9, 3.579339506199397*^9}, { 3.579339540051456*^9, 3.5793395413774586`*^9}, {3.5793395728427134`*^9, 3.579339683743308*^9}, {3.5793399816726317`*^9, 3.5793400106886826`*^9}, 3.579340927377493*^9, {3.5793433833366065`*^9, 3.5793433837422075`*^9}, { 3.579344512243789*^9, 3.57934455821707*^9}, 3.579344675232876*^9, { 3.579344794885086*^9, 3.5793447997678947`*^9}, {3.5793452439006743`*^9, 3.579345253806692*^9}, 3.5793453515408635`*^9, {3.5793504919472775`*^9, 3.5793505047393*^9}, {3.5793505769830265`*^9, 3.5793506236271086`*^9}, { 3.579350654125162*^9, 3.579350657697569*^9}, {3.579350819657053*^9, 3.5793508278626676`*^9}, {3.5793645321452856`*^9, 3.579364610285823*^9}, { 3.579364652530697*^9, 3.579364717567211*^9}, {3.579364755553278*^9, 3.5793647690629015`*^9}, {3.5793650698938303`*^9, 3.579365073965437*^9}, { 3.579446086369755*^9, 3.579446107585792*^9}, {3.579580259296129*^9, 3.5795802657857404`*^9}, {3.5796660160779314`*^9, 3.5796660212259407`*^9}}], Cell[BoxData["83.94022884384746`"], "Output", CellChangeTimes->{ 3.5795765080511403`*^9, 3.579576869456975*^9, 3.5795769335730877`*^9, 3.579576984881578*^9, 3.5795771258746257`*^9, 3.5795773314049864`*^9, 3.579577366817049*^9, 3.579577400684708*^9, 3.579577451852798*^9, 3.5795775558893805`*^9, 3.5795776291627097`*^9, 3.579577673763188*^9, { 3.5795779256724305`*^9, 3.579577949353272*^9}, 3.579577988712141*^9, 3.5795780334998198`*^9, 3.5795792101630864`*^9, 3.5795796974611425`*^9, 3.5795797320152025`*^9, 3.5795798843026705`*^9, 3.5795800146564994`*^9, 3.579580212215246*^9, 3.5795802718385515`*^9, 3.5795803441602783`*^9, 3.579662366794722*^9, 3.5796626467528133`*^9, {3.5796626964233007`*^9, 3.5796627141137323`*^9}, 3.5796658675344706`*^9, 3.5796659874050813`*^9, 3.579666028355153*^9, 3.5796669998624597`*^9, {3.5796671133838587`*^9, 3.5796671401535063`*^9}, 3.579671643043415*^9, 3.5796716812166824`*^9, 3.579671866638608*^9, 3.5796719021130705`*^9, {3.5796834755560064`*^9, 3.579683524758493*^9}, 3.5796835858014*^9, 3.5796837452336802`*^9, { 3.5796837843117485`*^9, 3.579683837991443*^9}, 3.579683922231591*^9, 3.5796839683140717`*^9, 3.5800574323641787`*^9, 3.5859073860537233`*^9, { 3.5859074260970135`*^9, 3.5859074505274105`*^9}, {3.5859075704642706`*^9, 3.585907593470587*^9}, 3.585907754322787*^9, 3.5859078030575743`*^9, 3.5859078754407144`*^9, 3.585907914910972*^9, {3.58590798760913*^9, 3.5859080146426764`*^9}, 3.5859081624531307`*^9, 3.5859092596439843`*^9, 3.585909292884886*^9, 3.585909424664423*^9, {3.585909525435187*^9, 3.5859095543978434`*^9}, 3.585909640217752*^9, 3.585910167209894*^9, 3.5859102092082963`*^9, 3.5859102597091846`*^9, 3.5859103241938734`*^9, 3.5859104488530035`*^9, {3.585910553781005*^9, 3.5859105730811086`*^9}, { 3.585910632117485*^9, 3.5859106524396477`*^9}, 3.5859108174710865`*^9, { 3.5859109613873186`*^9, 3.585910984961667*^9}, {3.5859110161454506`*^9, 3.5859110324623833`*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"plot1", "=", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"Graphics", "[", RowBox[{"{", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"a1", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"a2", "[", RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", RowBox[{"c", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}], "}"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "nc"}], "}"}]}], "]"}], " ", "/.", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], "\n", RowBox[{ RowBox[{"plot1a", "=", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"Graphics", "[", RowBox[{"{", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"b1", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"b2", "[", RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], "}"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "6"}], "}"}]}], "]"}], " ", "/.", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], "\n", RowBox[{ RowBox[{"plot1b", "=", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"Graphics", "[", RowBox[{"{", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"a1", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"a2", "[", RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], "}"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "nc"}], "}"}]}], "]"}], " ", "/.", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{"Show", "[", RowBox[{"plot1", ",", "plot1a", ",", "plot1b", ",", RowBox[{"{", RowBox[{ RowBox[{"AspectRatio", " ", "\[Rule]", " ", "Automatic"}], ",", RowBox[{"Axes", " ", "->", " ", "True"}], ",", RowBox[{"GridLines", " ", "->", " ", "Automatic"}], ",", " ", "\n", RowBox[{"PlotRange", " ", "\[Rule]", " ", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "4"}], ",", "4"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "4"}], ",", "4"}], "}"}]}], "}"}]}], ",", " ", RowBox[{"Frame", " ", "->", " ", "True"}]}], "}"}]}], "]"}], "\[IndentingNewLine]", StyleBox[ RowBox[{"xy", " ", "projection"}], "Text"]}], "Input", CellChangeTimes->{{3.579339709701754*^9, 3.5793398596336174`*^9}, { 3.5793399817038317`*^9, 3.579339981735032*^9}, {3.579579406271031*^9, 3.5795795916617565`*^9}, {3.5795796337662306`*^9, 3.579579644499049*^9}, { 3.5795797243243895`*^9, 3.5795797255879917`*^9}, {3.5795798997466974`*^9, 3.5795799348935595`*^9}, {3.5796626781556687`*^9, 3.579662707670921*^9}, { 3.5796717667828326`*^9, 3.5796717759712486`*^9}}], Cell[BoxData[ GraphicsBox[{CircleBox[{-1.65, 0}, 0.3183387910180002], CircleBox[{1.65, 0}, 0.3183387910180002], CircleBox[{-2.5059972462251365`, 1.5346516871748261`}, 1.438899692083698], CircleBox[{-2.5059972462251365`, -1.5346516871748261`}, 1.438899692083698], CircleBox[{2.5059972462251365`, -1.5346516871748261`}, 1.438899692083698], CircleBox[{2.5059972462251365`, 1.5346516871748261`}, 1.438899692083698], CircleBox[{0, 0}, 1.8271408028181948`], CircleBox[{0, 0}, 1.8271408028181948`], DiskBox[{-1.65, 0}, 0.08], DiskBox[{1.65, 0}, 0.08], DiskBox[{-2.7455472608147504`, 1.9641225010984673`}, 0.08], DiskBox[{-2.7455472608147504`, -1.9641225010984673`}, 0.08], DiskBox[{2.7455472608147504`, -1.9641225010984673`}, 0.08], DiskBox[{2.7455472608147504`, 1.9641225010984673`}, 0.08], CircleBox[{-1.65, 0}, 0.08], CircleBox[{1.65, 0}, 0.08], CircleBox[{-2.5059972462251365`, 1.5346516871748261`}, 0.08], CircleBox[{-2.5059972462251365`, -1.5346516871748261`}, 0.08], CircleBox[{2.5059972462251365`, -1.5346516871748261`}, 0.08], CircleBox[{2.5059972462251365`, 1.5346516871748261`}, 0.08], CircleBox[{0, 0}, 0.08], CircleBox[{0, 0}, 0.08]}, AspectRatio->Automatic, Axes->True, Frame->True, GridLines->Automatic, PlotRange->{{-4, 4}, {-4, 4}}]], "Output", CellChangeTimes->{ 3.5795765080667405`*^9, 3.579576869472575*^9, 3.579576933588688*^9, 3.5795769848971777`*^9, 3.5795771259058256`*^9, 3.5795773314205866`*^9, 3.5795773668326488`*^9, 3.579577400700308*^9, 3.5795774518683977`*^9, 3.5795775559049807`*^9, 3.5795776291783094`*^9, 3.5795776737787876`*^9, { 3.5795779256880302`*^9, 3.5795779493688717`*^9}, 3.5795779887277412`*^9, 3.5795780335154195`*^9, 3.579579210178686*^9, 3.5795796974767427`*^9, 3.5795797320308027`*^9, 3.5795798843026705`*^9, 3.579580014672099*^9, 3.5795802122308464`*^9, 3.5795802718541512`*^9, 3.579580344175878*^9, 3.579662366825922*^9, 3.5796626467684135`*^9, {3.579662696438901*^9, 3.579662714129332*^9}, 3.579665867550071*^9, 3.579665987420681*^9, 3.5796660283707533`*^9, 3.5796669998780594`*^9, {3.579667113399459*^9, 3.579667140169106*^9}, 3.5796716430746155`*^9, 3.5796716812322826`*^9, 3.579671866654208*^9, 3.5796719021286707`*^9, {3.579683475587207*^9, 3.5796835247740927`*^9}, 3.5796835858170004`*^9, 3.5796837452492805`*^9, { 3.5796837843273487`*^9, 3.579683838007043*^9}, 3.579683922231591*^9, 3.579683968329672*^9, 3.5800574323841796`*^9, 3.585907386071724*^9, { 3.5859074261120143`*^9, 3.5859074505414114`*^9}, {3.5859075704782715`*^9, 3.585907593485587*^9}, 3.5859077543377876`*^9, 3.5859078030765753`*^9, 3.585907875455715*^9, 3.5859079149299726`*^9, {3.5859079876241307`*^9, 3.585908014657677*^9}, 3.5859081624691315`*^9, 3.5859092596649857`*^9, 3.585909292899886*^9, 3.585909424683424*^9, {3.5859095254501877`*^9, 3.5859095544178443`*^9}, 3.585909640236753*^9, 3.585910167228895*^9, 3.585910209228297*^9, 3.585910259730186*^9, 3.585910324213874*^9, 3.5859104488700037`*^9, {3.5859105538010063`*^9, 3.5859105731011095`*^9}, { 3.5859106321324863`*^9, 3.5859106524546485`*^9}, 3.585910817491088*^9, { 3.585910961400319*^9, 3.5859109849796677`*^9}, {3.5859110161634517`*^9, 3.585911032480385*^9}}], Cell[BoxData[ RowBox[{"projection", " ", "xy"}]], "Output", CellChangeTimes->{ 3.5795765080667405`*^9, 3.579576869472575*^9, 3.579576933588688*^9, 3.5795769848971777`*^9, 3.5795771259058256`*^9, 3.5795773314205866`*^9, 3.5795773668326488`*^9, 3.579577400700308*^9, 3.5795774518683977`*^9, 3.5795775559049807`*^9, 3.5795776291783094`*^9, 3.5795776737787876`*^9, { 3.5795779256880302`*^9, 3.5795779493688717`*^9}, 3.5795779887277412`*^9, 3.5795780335154195`*^9, 3.579579210178686*^9, 3.5795796974767427`*^9, 3.5795797320308027`*^9, 3.5795798843026705`*^9, 3.579580014672099*^9, 3.5795802122308464`*^9, 3.5795802718541512`*^9, 3.579580344175878*^9, 3.579662366825922*^9, 3.5796626467684135`*^9, {3.579662696438901*^9, 3.579662714129332*^9}, 3.579665867550071*^9, 3.579665987420681*^9, 3.5796660283707533`*^9, 3.5796669998780594`*^9, {3.579667113399459*^9, 3.579667140169106*^9}, 3.5796716430746155`*^9, 3.5796716812322826`*^9, 3.579671866654208*^9, 3.5796719021286707`*^9, {3.579683475587207*^9, 3.5796835247740927`*^9}, 3.5796835858170004`*^9, 3.5796837452492805`*^9, { 3.5796837843273487`*^9, 3.579683838007043*^9}, 3.579683922231591*^9, 3.579683968329672*^9, 3.5800574323841796`*^9, 3.585907386071724*^9, { 3.5859074261120143`*^9, 3.5859074505414114`*^9}, {3.5859075704782715`*^9, 3.585907593485587*^9}, 3.5859077543377876`*^9, 3.5859078030765753`*^9, 3.585907875455715*^9, 3.5859079149299726`*^9, {3.5859079876241307`*^9, 3.585908014657677*^9}, 3.5859081624691315`*^9, 3.5859092596649857`*^9, 3.585909292899886*^9, 3.585909424683424*^9, {3.5859095254501877`*^9, 3.5859095544178443`*^9}, 3.585909640236753*^9, 3.585910167228895*^9, 3.585910209228297*^9, 3.585910259730186*^9, 3.585910324213874*^9, 3.5859104488700037`*^9, {3.5859105538010063`*^9, 3.5859105731011095`*^9}, { 3.5859106321324863`*^9, 3.5859106524546485`*^9}, 3.585910817491088*^9, { 3.585910961400319*^9, 3.5859109849796677`*^9}, {3.5859110161634517`*^9, 3.585911032484385*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"plot2", "=", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"Graphics", "[", RowBox[{"{", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"a1", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"a3", "[", RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", RowBox[{"c", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}], "}"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "nc"}], "}"}]}], "]"}], " ", "/.", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"plot2a", "=", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"Graphics", "[", RowBox[{"{", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"b1", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"b3", "[", RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], "}"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "nc"}], "}"}]}], "]"}], " ", "/.", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"plot2b", "=", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"Graphics", "[", RowBox[{"{", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"a1", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"a3", "[", RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}], "}"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "nc"}], "}"}]}], "]"}], " ", "/.", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{"Show", "[", RowBox[{"plot2", ",", "plot2a", ",", "plot2b", ",", RowBox[{"{", RowBox[{ RowBox[{"AspectRatio", " ", "\[Rule]", " ", "Automatic"}], ",", RowBox[{"Axes", " ", "->", " ", "True"}], ",", RowBox[{"GridLines", " ", "->", " ", "Automatic"}], ",", " ", "\n", RowBox[{"PlotRange", " ", "\[Rule]", " ", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "4"}], ",", "4"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "4"}], ",", "4"}], "}"}]}], "}"}]}], ",", " ", RowBox[{"Frame", " ", "->", " ", "True"}]}], "}"}]}], "]"}], "\[IndentingNewLine]", StyleBox[ RowBox[{"xz", " ", "projection"}], "Text"]}], "Input", CellChangeTimes->{{3.579339709701754*^9, 3.5793398772928486`*^9}, { 3.579339981735032*^9, 3.5793399817506323`*^9}, {3.5793507174300737`*^9, 3.579350731251698*^9}, {3.5795776651831727`*^9, 3.579577668599579*^9}, { 3.579579759954852*^9, 3.5795797795172863`*^9}, {3.579579823166163*^9, 3.57957986711144*^9}, {3.57957991244512*^9, 3.57957991831073*^9}, { 3.5796624612996883`*^9, 3.579662624101574*^9}, {3.5796717987004886`*^9, 3.5796718037080975`*^9}}], Cell[BoxData[ GraphicsBox[{CircleBox[{-1.65, 0}, 0.3183387910180002], CircleBox[{1.65, 0}, 0.3183387910180002], CircleBox[{-2.5059972462251365`, 0}, 1.438899692083698], CircleBox[{-2.5059972462251365`, 0}, 1.438899692083698], CircleBox[{2.5059972462251365`, 0}, 1.438899692083698], CircleBox[{2.5059972462251365`, 0}, 1.438899692083698], CircleBox[{0, 1.227140802818195}, 1.8271408028181948`], CircleBox[{0, -1.227140802818195}, 1.8271408028181948`], DiskBox[{-1.65, 0}, 0.08], DiskBox[{1.65, 0}, 0.08], DiskBox[{-2.7455472608147504`, 0}, 0.08], DiskBox[{-2.7455472608147504`, 0}, 0.08], DiskBox[{2.7455472608147504`, 0}, 0.08], DiskBox[{2.7455472608147504`, 0}, 0.08], DiskBox[{0, 1.8344339983370048`}, 0.08], DiskBox[{0, -1.8344339983370048`}, 0.08], CircleBox[{-1.65, 0}, 0.08], CircleBox[{1.65, 0}, 0.08], CircleBox[{-2.5059972462251365`, 0}, 0.08], CircleBox[{-2.5059972462251365`, 0}, 0.08], CircleBox[{2.5059972462251365`, 0}, 0.08], CircleBox[{2.5059972462251365`, 0}, 0.08], CircleBox[{0, 1.227140802818195}, 0.08], CircleBox[{0, -1.227140802818195}, 0.08]}, AspectRatio->Automatic, Axes->True, Frame->True, GridLines->Automatic, ImageSize->{364., Automatic}, PlotRange->{{-4, 4}, {-4, 4}}]], "Output", CellChangeTimes->{ 3.5795765080823402`*^9, 3.5795768694881754`*^9, 3.579576933604288*^9, 3.579576984912778*^9, 3.579577125921426*^9, 3.5795773314361863`*^9, 3.5795773668482485`*^9, 3.579577400715908*^9, 3.579577451883998*^9, 3.579577555920581*^9, 3.5795776291939096`*^9, 3.579577673794388*^9, { 3.5795779257036304`*^9, 3.579577949384472*^9}, 3.579577988743341*^9, 3.5795780335310197`*^9, 3.5795792101942863`*^9, 3.5795796974923425`*^9, 3.579579732046403*^9, 3.5795798843338704`*^9, 3.5795800146876993`*^9, 3.5795802122464466`*^9, 3.579580271869751*^9, 3.579580344207078*^9, 3.579662366841522*^9, 3.5796626467840137`*^9, {3.579662696454501*^9, 3.579662714144932*^9}, 3.579665867550071*^9, 3.5796659874518814`*^9, 3.5796660283863535`*^9, 3.5796669998936596`*^9, {3.5796671134306593`*^9, 3.579667140200306*^9}, 3.5796716430746155`*^9, 3.5796716812478824`*^9, 3.579671866669808*^9, 3.5796719021442704`*^9, {3.5796834756028066`*^9, 3.579683524789693*^9}, 3.5796835858326006`*^9, 3.57968374526488*^9, { 3.579683784342949*^9, 3.579683838022643*^9}, 3.579683922262791*^9, 3.5796839683608723`*^9, 3.580057432412181*^9, 3.5859073860927253`*^9, { 3.585907426132015*^9, 3.585907450561413*^9}, {3.585907570501273*^9, 3.5859075935095887`*^9}, 3.585907754362789*^9, 3.585907803099577*^9, 3.585907875481717*^9, 3.5859079149539747`*^9, {3.585907987650132*^9, 3.585908014682678*^9}, 3.5859081624901323`*^9, 3.585909259686987*^9, 3.585909292920888*^9, 3.5859094247084255`*^9, {3.5859095254711885`*^9, 3.5859095544388456`*^9}, 3.585909640259754*^9, 3.5859101672508965`*^9, 3.585910209255299*^9, 3.585910259754187*^9, 3.5859103242358756`*^9, 3.5859104488920054`*^9, {3.585910553823007*^9, 3.585910573124111*^9}, { 3.5859106321544876`*^9, 3.5859106524776497`*^9}, 3.5859108175160894`*^9, { 3.5859109614253206`*^9, 3.585910985004669*^9}, {3.5859110161884527`*^9, 3.5859110325053864`*^9}}], Cell[BoxData[ RowBox[{"projection", " ", "xz"}]], "Output", CellChangeTimes->{ 3.5795765080823402`*^9, 3.5795768694881754`*^9, 3.579576933604288*^9, 3.579576984912778*^9, 3.579577125921426*^9, 3.5795773314361863`*^9, 3.5795773668482485`*^9, 3.579577400715908*^9, 3.579577451883998*^9, 3.579577555920581*^9, 3.5795776291939096`*^9, 3.579577673794388*^9, { 3.5795779257036304`*^9, 3.579577949384472*^9}, 3.579577988743341*^9, 3.5795780335310197`*^9, 3.5795792101942863`*^9, 3.5795796974923425`*^9, 3.579579732046403*^9, 3.5795798843338704`*^9, 3.5795800146876993`*^9, 3.5795802122464466`*^9, 3.579580271869751*^9, 3.579580344207078*^9, 3.579662366841522*^9, 3.5796626467840137`*^9, {3.579662696454501*^9, 3.579662714144932*^9}, 3.579665867550071*^9, 3.5796659874518814`*^9, 3.5796660283863535`*^9, 3.5796669998936596`*^9, {3.5796671134306593`*^9, 3.579667140200306*^9}, 3.5796716430746155`*^9, 3.5796716812478824`*^9, 3.579671866669808*^9, 3.5796719021442704`*^9, {3.5796834756028066`*^9, 3.579683524789693*^9}, 3.5796835858326006`*^9, 3.57968374526488*^9, { 3.579683784342949*^9, 3.579683838022643*^9}, 3.579683922262791*^9, 3.5796839683608723`*^9, 3.580057432412181*^9, 3.5859073860927253`*^9, { 3.585907426132015*^9, 3.585907450561413*^9}, {3.585907570501273*^9, 3.5859075935095887`*^9}, 3.585907754362789*^9, 3.585907803099577*^9, 3.585907875481717*^9, 3.5859079149539747`*^9, {3.585907987650132*^9, 3.585908014682678*^9}, 3.5859081624901323`*^9, 3.585909259686987*^9, 3.585909292920888*^9, 3.5859094247084255`*^9, {3.5859095254711885`*^9, 3.5859095544388456`*^9}, 3.585909640259754*^9, 3.5859101672508965`*^9, 3.585910209255299*^9, 3.585910259754187*^9, 3.5859103242358756`*^9, 3.5859104488920054`*^9, {3.585910553823007*^9, 3.585910573124111*^9}, { 3.5859106321544876`*^9, 3.5859106524776497`*^9}, 3.5859108175160894`*^9, { 3.5859109614253206`*^9, 3.585910985004669*^9}, {3.5859110161884527`*^9, 3.5859110325103865`*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"plot3", "=", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"Graphics3D", "[", RowBox[{"{", RowBox[{ RowBox[{"Opacity", "[", "0.4", "]"}], ",", RowBox[{"Lighter", "[", RowBox[{"Gray", ",", "0.4"}], "]"}], ",", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"a1", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"a2", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"a3", "[", RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", RowBox[{"c", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}], "}"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", RowBox[{"nc", "-", "2"}]}], "}"}]}], "]"}], " ", "/.", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"plot3a", "=", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"Graphics3D", "[", RowBox[{"{", RowBox[{ RowBox[{"Darker", "[", RowBox[{"Red", ",", "1.0"}], "]"}], ",", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"b1", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"b2", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"b3", "[", RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}]}], "}"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "nc"}], "}"}]}], "]"}], " ", "/.", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"plot3aa", "=", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"Graphics3D", "[", RowBox[{"{", RowBox[{ RowBox[{"Opacity", "[", "0.8", "]"}], ",", RowBox[{"Lighter", "[", RowBox[{"White", ",", "0.2"}], "]"}], ",", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"b1", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"b2", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"b3", "[", RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", "0.3"}], "]"}]}], "}"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "3", ",", "nc"}], "}"}]}], "]"}], " ", "/.", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"plot3b", "=", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"Graphics3D", "[", RowBox[{"{", RowBox[{ RowBox[{"Opacity", "[", "0.6", "]"}], ",", RowBox[{"Lighter", "[", RowBox[{"Blue", ",", "0.7"}], "]"}], ",", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"a1", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"a2", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"a3", "[", RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", RowBox[{"0.9", "*", RowBox[{"c", "[", RowBox[{"[", "i", "]"}], "]"}]}]}], "]"}]}], "}"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", RowBox[{"nc", "-", "1"}], ",", "nc"}], "}"}]}], "]"}], " ", "/.", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"plot3c", "=", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"Graphics3D", "[", RowBox[{"{", RowBox[{ RowBox[{"Opacity", "[", "0.6", "]"}], ",", RowBox[{"Lighter", "[", RowBox[{"Green", ",", "0.2"}], "]"}], ",", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"a1", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"a2", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"a3", "[", RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", RowBox[{"c", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}], "}"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "2"}], "}"}]}], "]"}], " ", "/.", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Plot9", "=", " ", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"Graphics3D", "[", RowBox[{"{", RowBox[{ RowBox[{"Lighter", "[", RowBox[{"Orange", ",", "0.5"}], "]"}], ",", RowBox[{"Tube", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"b1", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"b2", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"b3", "[", RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"b1", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"b2", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"b3", "[", RowBox[{"[", "i", "]"}], "]"}]}], "}"}]}], "}"}], ",", "0.1"}], "]"}]}], "}"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "3", ",", "4"}], "}"}]}], "]"}], " ", "/.", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Plot91", "=", " ", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"Graphics3D", "[", RowBox[{"{", RowBox[{ RowBox[{"Lighter", "[", RowBox[{"Orange", ",", "0.5"}], "]"}], ",", RowBox[{"Tube", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"b1", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"b2", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"b3", "[", RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"b1", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"b2", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"b3", "[", RowBox[{"[", "i", "]"}], "]"}]}], "}"}]}], "}"}], ",", "0.1"}], "]"}]}], "}"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "5", ",", "6"}], "}"}]}], "]"}], " ", "/.", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Plot92", "=", " ", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"Graphics3D", "[", RowBox[{"{", RowBox[{ RowBox[{"Lighter", "[", RowBox[{"Orange", ",", "0.5"}], "]"}], ",", RowBox[{"Tube", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"b1", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"b2", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"b3", "[", RowBox[{"[", "1", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"b1", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"b2", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"b3", "[", RowBox[{"[", "i", "]"}], "]"}]}], "}"}]}], "}"}], ",", "0.1"}], "]"}]}], "}"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "7", ",", "8"}], "}"}]}], "]"}], " ", "/.", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Plot93", "=", " ", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"Graphics3D", "[", RowBox[{"{", RowBox[{ RowBox[{"Lighter", "[", RowBox[{"Orange", ",", "0.5"}], "]"}], ",", RowBox[{"Tube", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"b1", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"b2", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"b3", "[", RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"b1", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"b2", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"b3", "[", RowBox[{"[", "i", "]"}], "]"}]}], "}"}]}], "}"}], ",", "0.1"}], "]"}]}], "}"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "7", ",", "8"}], "}"}]}], "]"}], " ", "/.", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Plot94", "=", " ", RowBox[{ RowBox[{"Graphics3D", "[", RowBox[{"{", RowBox[{ RowBox[{"Lighter", "[", RowBox[{"Green", ",", "0.2"}], "]"}], ",", RowBox[{"Tube", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"b1", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"b2", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"b3", "[", RowBox[{"[", "2", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"b1", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"b2", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"b3", "[", RowBox[{"[", "1", "]"}], "]"}]}], "}"}]}], "}"}], ",", "0.1"}], "]"}]}], "}"}], "]"}], " ", "/.", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{"Show", "[", RowBox[{ "plot3", ",", "plot3aa", ",", "plot3b", ",", "plot3c", ",", "Plot9", ",", "Plot91", ",", "Plot92", ",", "Plot93", ",", "Plot94", ",", "plot3a", ",", RowBox[{"{", RowBox[{ RowBox[{"AspectRatio", " ", "\[Rule]", " ", "Automatic"}], ",", RowBox[{"(*", RowBox[{ RowBox[{"Axes", " ", "->", " ", "True"}], ","}], "*)"}], RowBox[{"Boxed", "\[Rule]", "False"}], ",", RowBox[{"SphericalRegion", "\[Rule]", "True"}]}], "}"}]}], "]"}], "\[IndentingNewLine]"}], "Input", CellChangeTimes->{{3.579683275360855*^9, 3.5796833411461706`*^9}, { 3.5796833720342245`*^9, 3.5796834557283716`*^9}, {3.579683488488429*^9, 3.5796835178788805`*^9}, {3.5796835721981764`*^9, 3.5796835803413906`*^9}, {3.579683651695916*^9, 3.579683739929671*^9}, { 3.5796837704589243`*^9, 3.5796838321258326`*^9}, {3.5796838663210926`*^9, 3.5796839172551823`*^9}, 3.5796839627448626`*^9, {3.585907377399228*^9, 3.585907378732304*^9}, {3.5859074197196484`*^9, 3.585907443942034*^9}, { 3.5859074991431913`*^9, 3.5859075881742835`*^9}, {3.585907731278469*^9, 3.5859077984813128`*^9}, {3.5859078312301855`*^9, 3.5859078362284718`*^9}, {3.5859078949488297`*^9, 3.5859079070525227`*^9}, {3.5859079760454683`*^9, 3.585908009973409*^9}, { 3.5859081325404196`*^9, 3.5859081526435695`*^9}, {3.585909245704187*^9, 3.585909252103553*^9}, {3.585909286629528*^9, 3.585909287092554*^9}, 3.5859094203001733`*^9, {3.585909506540106*^9, 3.5859095194268427`*^9}, 3.5859095497555776`*^9, {3.585909629659148*^9, 3.5859096302821836`*^9}, 3.585909866579699*^9, {3.5859100302510605`*^9, 3.5859101007060905`*^9}, { 3.585910144808613*^9, 3.585910160991538*^9}, 3.585910196896592*^9, { 3.5859102510476894`*^9, 3.5859103174724884`*^9}, {3.585910361746021*^9, 3.5859103727486506`*^9}, {3.5859104120368977`*^9, 3.5859104396994796`*^9}, {3.5859104768236027`*^9, 3.5859105651306543`*^9}, {3.5859106239850197`*^9, 3.585910647152345*^9}, 3.5859107391156054`*^9, {3.585910799548062*^9, 3.5859108112507315`*^9}, { 3.5859109377849684`*^9, 3.585910978616304*^9}, {3.5859110087210255`*^9, 3.5859110272970886`*^9}}], Cell[BoxData[ Graphics3DBox[{ {RGBColor[0.7, 0.7, 0.7], Opacity[0.4], SphereBox[{-1.65, 0, 0}, 0.3183387910180002]}, {RGBColor[0.7, 0.7, 0.7], Opacity[0.4], SphereBox[{1.65, 0, 0}, 0.3183387910180002]}, {RGBColor[0.7, 0.7, 0.7], Opacity[0.4], SphereBox[{-2.5059972462251365`, 1.5346516871748261`, 0}, 1.438899692083698]}, {RGBColor[0.7, 0.7, 0.7], Opacity[0.4], SphereBox[{-2.5059972462251365`, -1.5346516871748261`, 0}, 1.438899692083698]}, {RGBColor[0.7, 0.7, 0.7], Opacity[0.4], SphereBox[{2.5059972462251365`, -1.5346516871748261`, 0}, 1.438899692083698]}, {RGBColor[0.7, 0.7, 0.7], Opacity[0.4], SphereBox[{2.5059972462251365`, 1.5346516871748261`, 0}, 1.438899692083698]}, {RGBColor[1., 1., 1.], Opacity[0.8], SphereBox[{-2.7455472608147504`, 1.9641225010984673`, 0}, 0.3]}, {RGBColor[1., 1., 1.], Opacity[0.8], SphereBox[{-2.7455472608147504`, -1.9641225010984673`, 0}, 0.3]}, {RGBColor[1., 1., 1.], Opacity[0.8], SphereBox[{2.7455472608147504`, -1.9641225010984673`, 0}, 0.3]}, {RGBColor[1., 1., 1.], Opacity[0.8], SphereBox[{2.7455472608147504`, 1.9641225010984673`, 0}, 0.3]}, {RGBColor[1., 1., 1.], Opacity[0.8], SphereBox[{0, 0, 1.8344339983370048`}, 0.3]}, {RGBColor[1., 1., 1.], Opacity[0.8], SphereBox[{0, 0, -1.8344339983370048`}, 0.3]}, {RGBColor[0.7, 0.7, 1.], Opacity[0.6], SphereBox[{0, 0, 1.227140802818195}, 1.6444267225363753`]}, {RGBColor[0.7, 0.7, 1.], Opacity[0.6], SphereBox[{0, 0, -1.227140802818195}, 1.6444267225363753`]}, {RGBColor[0.2, 1., 0.2], Opacity[0.6], SphereBox[{-1.65, 0, 0}, 0.3183387910180002]}, {RGBColor[0.2, 1., 0.2], Opacity[0.6], SphereBox[{1.65, 0, 0}, 0.3183387910180002]}, {RGBColor[1., 0.75, 0.5], TubeBox[{{-1.65, 0, 0}, {-2.7455472608147504`, 1.9641225010984673`, 0}}, 0.1]}, {RGBColor[1., 0.75, 0.5], TubeBox[{{-1.65, 0, 0}, {-2.7455472608147504`, -1.9641225010984673`, 0}}, 0.1]}, {RGBColor[1., 0.75, 0.5], TubeBox[{{1.65, 0, 0}, {2.7455472608147504`, -1.9641225010984673`, 0}}, 0.1]}, {RGBColor[1., 0.75, 0.5], TubeBox[{{1.65, 0, 0}, {2.7455472608147504`, 1.9641225010984673`, 0}}, 0.1]}, {RGBColor[1., 0.75, 0.5], TubeBox[{{-1.65, 0, 0}, {0, 0, 1.8344339983370048`}}, 0.1]}, {RGBColor[1., 0.75, 0.5], TubeBox[{{-1.65, 0, 0}, {0, 0, -1.8344339983370048`}}, 0.1]}, {RGBColor[1., 0.75, 0.5], TubeBox[{{1.65, 0, 0}, {0, 0, 1.8344339983370048`}}, 0.1]}, {RGBColor[1., 0.75, 0.5], TubeBox[{{1.65, 0, 0}, {0, 0, -1.8344339983370048`}}, 0.1]}, {RGBColor[0.2, 1., 0.2], TubeBox[{{1.65, 0, 0}, {-1.65, 0, 0}}, 0.1]}, {RGBColor[0., 0., 0.], SphereBox[{-1.65, 0, 0}, 0.08]}, {RGBColor[0., 0., 0.], SphereBox[{1.65, 0, 0}, 0.08]}, {RGBColor[0., 0., 0.], SphereBox[{-2.7455472608147504`, 1.9641225010984673`, 0}, 0.08]}, {RGBColor[0., 0., 0.], SphereBox[{-2.7455472608147504`, -1.9641225010984673`, 0}, 0.08]}, {RGBColor[0., 0., 0.], SphereBox[{2.7455472608147504`, -1.9641225010984673`, 0}, 0.08]}, {RGBColor[0., 0., 0.], SphereBox[{2.7455472608147504`, 1.9641225010984673`, 0}, 0.08]}, {RGBColor[0., 0., 0.], SphereBox[{0, 0, 1.8344339983370048`}, 0.08]}, {RGBColor[0., 0., 0.], SphereBox[{0, 0, -1.8344339983370048`}, 0.08]}}, AspectRatio->Automatic, AutomaticImageSize->True, Boxed->False, ImageSize->{360., 360.}, SphericalRegion->True, ViewPoint->{-0.19104866644798024`, -3.136985822287172, 1.254121348920321}, ViewVertical->{-0.0019408960092615933`, -0.33217829931778015`, 1.3300151391459756`}]], "Output", CellChangeTimes->{ 3.585910573159113*^9, {3.5859106321914897`*^9, 3.5859106525156517`*^9}, 3.5859108175570917`*^9, {3.5859109616853356`*^9, 3.5859109852256823`*^9}, { 3.585911016404465*^9, 3.585911032554389*^9}}] }, Open ]] }, WindowToolbars->"EditBar", WindowSize->{1042, 1002}, WindowMargins->{{Automatic, 175}, {Automatic, 53}}, PrintingCopies->1, PrintingPageRange->{Automatic, Automatic}, PrivateNotebookOptions->{"VersionedStylesheet"->{"Default.nb"[8.] -> False}}, CellLabelAutoDelete->True, FrontEndVersion->"9.0 for Microsoft Windows (64-bit) (January 25, 2013)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[557, 20, 1673, 31, 92, "Input"], Cell[2233, 53, 2378, 43, 92, "Input"], Cell[4614, 98, 652, 12, 31, "Input"], Cell[5269, 112, 1255, 30, 92, "Input"], Cell[6527, 144, 5469, 162, 392, "Input"], Cell[11999, 308, 4017, 100, 192, "Input"], Cell[16019, 410, 3787, 113, 212, "Input"], Cell[CellGroupData[{ Cell[19831, 527, 2665, 54, 72, "Input"], Cell[22499, 583, 2530, 40, 31, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[25066, 628, 1494, 26, 72, "Input"], Cell[26563, 656, 867, 12, 31, "Output"], Cell[27433, 670, 882, 13, 31, "Output"], Cell[28318, 685, 866, 12, 31, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[29221, 702, 1298, 20, 31, "Input"], Cell[30522, 724, 878, 12, 31, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[31437, 741, 1296, 20, 31, "Input"], Cell[32736, 763, 2011, 28, 31, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[34784, 796, 1506, 25, 31, "Input"], Cell[36293, 823, 2002, 28, 31, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[38332, 856, 1588, 27, 31, "Input"], Cell[39923, 885, 1989, 28, 31, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[41949, 918, 578, 14, 31, "Input"], Cell[42530, 934, 1992, 28, 31, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[44559, 967, 1563, 25, 31, "Input"], Cell[46125, 994, 2014, 28, 31, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[48176, 1027, 1683, 27, 31, "Input"], Cell[49862, 1056, 2033, 28, 31, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[51932, 1089, 3139, 87, 132, "Input"], Cell[55074, 1178, 3322, 52, 379, "Output"], Cell[58399, 1232, 2045, 29, 31, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[60481, 1266, 3214, 88, 132, "Input"], Cell[63698, 1356, 3272, 55, 383, "Output"], Cell[66973, 1413, 2028, 29, 31, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[69038, 1447, 13035, 345, 552, "Input"], Cell[82076, 1794, 3918, 82, 407, "Output"] }, Open ]] } ] *) (* End of internal cache information *)