(* 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[ 107721, 2537] NotebookOptionsPosition[ 104125, 2439] NotebookOutlinePosition[ 104637, 2458] CellTagsIndexPosition[ 104594, 2455] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{"F2", " ", "with", " ", "tetrahedric", " ", RowBox[{"clouds", ":", " ", RowBox[{"Linnet", " ", "cannot", " ", "represent", " ", RowBox[{"F2", " ", "!"}], " ", "18.05", ".2013"}]}]}], " ", "*)"}], "\n", RowBox[{ RowBox[{ RowBox[{"Clear", "[", RowBox[{ "a", ",", "z", ",", "sig1", ",", "sig2", ",", "sig3", ",", "sig4", ",", "k1", ",", "k2", ",", "k3", ",", "k4", ",", "k5", ",", "nc", ",", "vee", ",", "vne", ",", "vnn", ",", "R1", ",", "R2", ",", "R3", ",", "R4", ",", "R5", ",", "\n", "xw", ",", "yw", ",", "zw", ",", "xn", ",", "yn", ",", "zn", ",", "ow", ",", "cn", ",", "rr", ",", "d1", ",", "d2", ",", "d3", ",", "d4", ",", "d5", ",", "i", ",", "j", ",", "d", ",", "t", ",", "sq", ",", "Ekin", ",", "Epot", ",", "pi", ",", "p", ",", "cc", ",", "dv", ",", "dw", ",", "dcw", ",", "dcc", ",", "tw", ",", "tv", ",", "w4", ",", "w5", ",", "q"}], "]"}], ";"}], "\n", RowBox[{ RowBox[{"bohr", "=", "0.529177"}], ";"}], " ", "\n", RowBox[{"(*", RowBox[{ RowBox[{"R1", "=", "0.1486"}], ";", RowBox[{"\[Alpha]", "=", "1.23"}], ";", "\[IndentingNewLine]", RowBox[{"R3", "=", "0.96"}], ";", " ", RowBox[{"R2", "=", "1.166"}], ";"}], "*)"}], "\[IndentingNewLine]"}]}]], "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.5773415047046976`*^9, 3.5773415137059135`*^9}, { 3.5773416513449545`*^9, 3.577341669752987*^9}, {3.5773417239318824`*^9, 3.577341893956581*^9}, 3.5773419407566633`*^9, {3.5773420175555983`*^9, 3.577342056165666*^9}, {3.5773420929973307`*^9, 3.5773421056489525`*^9}, { 3.5773421396414127`*^9, 3.5773423470593767`*^9}, {3.577342382081438*^9, 3.577342441158742*^9}, {3.5773425028880506`*^9, 3.577342523682887*^9}, { 3.5773425980794177`*^9, 3.5773426626791315`*^9}, {3.5773426933487854`*^9, 3.5773426941131864`*^9}, {3.577342738354864*^9, 3.577342889035529*^9}, { 3.5773429412644205`*^9, 3.57734295234044*^9}, {3.577342991824109*^9, 3.5773430231333647`*^9}, {3.5773430561274223`*^9, 3.577343143737176*^9}, { 3.5773431993668737`*^9, 3.577343258334977*^9}, 3.5773434426249013`*^9, { 3.577343480042567*^9, 3.5773435125686235`*^9}, {3.5773436063903885`*^9, 3.5773436272008247`*^9}, {3.5773436617392855`*^9, 3.5773437808610945`*^9}, {3.577343868193248*^9, 3.577343879659268*^9}, { 3.5773439214829416`*^9, 3.577344118838888*^9}, {3.5773441503509436`*^9, 3.577344170740179*^9}, {3.5773446920694942`*^9, 3.577344693286297*^9}, 3.5773526894956493`*^9, {3.5773528388503113`*^9, 3.577352838897112*^9}, { 3.577353457328998*^9, 3.577353503364679*^9}, {3.5773556679320807`*^9, 3.5773557165573664`*^9}, 3.5773607112650003`*^9, {3.577363145384075*^9, 3.5773631641821084`*^9}, 3.5773713130222435`*^9, {3.5773735851350346`*^9, 3.577373695474028*^9}, {3.577373964730501*^9, 3.5773739677725067`*^9}, 3.5773740596254673`*^9, {3.577374092245125*^9, 3.5773741045067463`*^9}, { 3.5773741853616886`*^9, 3.577374230414568*^9}, 3.5773742637986264`*^9, { 3.5773746139412413`*^9, 3.5773746339092765`*^9}, {3.5773750495252066`*^9, 3.5773750633000307`*^9}, 3.5773754561555204`*^9, {3.5773757271123962`*^9, 3.577375728063998*^9}, 3.577375963078411*^9, {3.5773761205918875`*^9, 3.5773762048008356`*^9}, 3.5774291695117903`*^9, {3.5774424964722877`*^9, 3.5774424977046895`*^9}, {3.5774428197828555`*^9, 3.577442825336465*^9}, { 3.577442927267044*^9, 3.5774429281406455`*^9}, {3.5774429766411304`*^9, 3.577442992833959*^9}, {3.5774430583072743`*^9, 3.577443066715689*^9}, { 3.5774436977055974`*^9, 3.5774437122604227`*^9}, 3.5774438419902506`*^9, { 3.577446912512444*^9, 3.577446913651246*^9}, {3.577447169460495*^9, 3.577447179070112*^9}, {3.577449454744052*^9, 3.57744946514927*^9}, { 3.577449615377534*^9, 3.5774496616940155`*^9}, {3.577450578242425*^9, 3.5774506245901065`*^9}, {3.577451618688603*^9, 3.5774516195310044`*^9}, { 3.5774516668614874`*^9, 3.5774516695602922`*^9}, 3.577455410676982*^9, { 3.5774557532847834`*^9, 3.5774558029396706`*^9}, {3.5774565586985984`*^9, 3.5774565589793987`*^9}, 3.577457075496306*^9, {3.577457117273179*^9, 3.577457148941235*^9}, {3.577678508293971*^9, 3.5776785428012314`*^9}, { 3.5776816489512606`*^9, 3.5776816557372727`*^9}, 3.5776817445014286`*^9, { 3.5776818073071394`*^9, 3.5776818117219467`*^9}, 3.577681885338476*^9, { 3.577682125485298*^9, 3.577682126124899*^9}, 3.5776821818169966`*^9, { 3.577682250082717*^9, 3.5776822522355204`*^9}, {3.5776823382072716`*^9, 3.5776823388624725`*^9}, 3.577682369656927*^9, {3.577682527513604*^9, 3.5776825496656427`*^9}, {3.577682803010088*^9, 3.5776828162233114`*^9}, { 3.577684829687648*^9, 3.5776848354440575`*^9}, 3.5776850034095526`*^9, { 3.5776850335332055`*^9, 3.577685052362439*^9}, {3.5776850854032965`*^9, 3.5776852462551794`*^9}, {3.5776852913080587`*^9, 3.577685293601262*^9}, { 3.577685324068116*^9, 3.577685477260385*^9}, 3.5776855432329006`*^9, 3.5776855922013865`*^9, {3.577685645615881*^9, 3.577685645803081*^9}, { 3.5776856793431396`*^9, 3.5776856800139413`*^9}, {3.5776857127739983`*^9, 3.577685783754123*^9}, {3.5776858189009852`*^9, 3.5776858890075083`*^9}, { 3.5776859215491652`*^9, 3.577685994728894*^9}, 3.577686107969493*^9, { 3.5776872286366615`*^9, 3.5776873647301006`*^9}, {3.5776874640554748`*^9, 3.5776874877519164`*^9}, {3.5776877853380394`*^9, 3.577687807490078*^9}, { 3.577697174278837*^9, 3.57769723305974*^9}, {3.577698938949665*^9, 3.5776989710545216`*^9}, 3.5776992675018425`*^9, {3.5776996873921795`*^9, 3.577699689232983*^9}, {3.5777703272886305`*^9, 3.577770335821845*^9}, { 3.577773226787723*^9, 3.577773228893727*^9}, {3.577775199551588*^9, 3.5777752073828015`*^9}, {3.577788184156394*^9, 3.5777881865899982`*^9}, { 3.5778008491634393`*^9, 3.5778008555282507`*^9}, {3.577801589103939*^9, 3.577801632799616*^9}, {3.5778052419439545`*^9, 3.5778052431295567`*^9}, { 3.577809188251686*^9, 3.5778092039453135`*^9}, {3.577809267531025*^9, 3.577809279043845*^9}, {3.5778094313625126`*^9, 3.577809438866126*^9}, { 3.577809500751435*^9, 3.577809520345069*^9}, 3.577854175078371*^9, { 3.57785477798823*^9, 3.577854780874235*^9}, {3.5778554270897703`*^9, 3.5778554580714245`*^9}, {3.577855503155504*^9, 3.577855507398711*^9}, { 3.5778846533117895`*^9, 3.577884661548604*^9}, {3.5778883011732435`*^9, 3.577888342294916*^9}, {3.5778884221358557`*^9, 3.5778884303570704`*^9}, 3.5778924889237995`*^9, 3.577892597967991*^9, {3.5778933300304766`*^9, 3.577893333743283*^9}, 3.5778934445034776`*^9, {3.5778936859919014`*^9, 3.57789369068751*^9}, 3.5778968416902447`*^9, {3.5778968733115*^9, 3.5778968737015004`*^9}, {3.5778969790484858`*^9, 3.5778969839936943`*^9}, {3.5778970399041924`*^9, 3.5778970599970274`*^9}, 3.577897095767891*^9, 3.577897163082009*^9, {3.5778974648177385`*^9, 3.5778974651297393`*^9}, 3.577897572161527*^9}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"nc", "=", "9"}], ";", " ", RowBox[{"(*", " ", RowBox[{"number", " ", "of", " ", "clouds"}], " ", "*)"}], " ", RowBox[{"z", "=", "9.0"}], ";", " ", RowBox[{"pi", "=", "0.0"}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"sig1", "=", "0.325"}], ";", " ", RowBox[{"sig2", "=", "0.325"}], ";", " ", RowBox[{"sig3", "=", "sig2"}], ";", " ", RowBox[{"(*", " ", RowBox[{ RowBox[{"screening", " ", RowBox[{"const", ".", " ", "from"}], " ", "e"}], "-", RowBox[{"e", " ", "interaction", " ", "in", " ", "doubly", " ", RowBox[{"occ", ".", " ", "clouds"}]}]}], " ", "*)"}], "\n", RowBox[{"c", "=", RowBox[{"{", RowBox[{ RowBox[{"k1", "->", "0.8828"}], ",", " ", RowBox[{"k2", "->", "1.33"}], ",", " ", RowBox[{"k3", "->", "1.2"}]}], "}"}]}], ";", RowBox[{"(*", " ", RowBox[{ RowBox[{ "parameters", " ", "for", " ", "kinetic", " ", "energy", " ", "of", " ", "clouds"}], ";"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{"cloud", " ", "charges"}], " ", "*)"}], "\n", RowBox[{"ow", "=", RowBox[{"{", RowBox[{ RowBox[{"-", "2"}], ",", RowBox[{"-", "2"}], ",", RowBox[{"-", "2"}], ",", RowBox[{"-", "2"}], ",", RowBox[{"-", "2"}], ",", RowBox[{"-", "2"}], ",", RowBox[{"-", "2"}], ",", RowBox[{"-", "2"}], ",", RowBox[{"-", "2"}]}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{"nucleus", " ", "charges"}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"cn", "=", RowBox[{"{", " ", RowBox[{ "z", ",", "z", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0"}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{"core", " ", "contributions"}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Ekin", " ", "=", " ", RowBox[{ RowBox[{"2", "*", RowBox[{"(", RowBox[{"2.25", "*", RowBox[{"k1", "/", RowBox[{"R1", "^", "2"}]}]}], ")"}]}], " ", "/.", " ", "c"}]}], ";"}], "\n", RowBox[{ RowBox[{"vee", "=", RowBox[{"2", "*", RowBox[{"(", RowBox[{"3.0", "*", RowBox[{"sig1", "/", "R1"}]}], ")"}]}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"vne", "=", RowBox[{ RowBox[{"-", "2"}], "*", RowBox[{"(", RowBox[{"3.0", "*", RowBox[{"z", "/", "R1"}]}], ")"}]}]}], ";"}], "\n", RowBox[{"(*", " ", RowBox[{ RowBox[{"bonding", " ", "clouds"}], ",", " ", RowBox[{"Ekin", " ", "and", " ", "internal", " ", "vee"}]}], " ", "*)"}]}], "\n", RowBox[{ RowBox[{"Ekin", " ", "=", " ", RowBox[{ RowBox[{"Ekin", " ", "+", " ", RowBox[{"2.25", "*", RowBox[{"(", RowBox[{ RowBox[{"k2", "/", RowBox[{"R2", "^", "2"}]}], "+", RowBox[{"6", "*", RowBox[{"k3", "/", RowBox[{"R3", "^", "2"}]}]}]}], ")"}]}]}], "/.", " ", "c"}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"vee", "=", RowBox[{ RowBox[{"vee", "+", RowBox[{"3.0", "*", RowBox[{"(", RowBox[{ RowBox[{"sig2", "/", "R2"}], "+", RowBox[{"6", "*", RowBox[{"sig3", "/", "R3"}]}]}], ")"}]}]}], " ", "/.", " ", "c"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ "nuclear", " ", "coordinates", " ", "in", " ", "terms", " ", "of", " ", "radii"}], " ", "*)"}]}], "\n", RowBox[{ RowBox[{ RowBox[{"d1", "=", RowBox[{"R1", "+", "R2"}]}], ";", " ", RowBox[{"d2", "=", RowBox[{"R1", "+", "R3"}]}], ";", " ", RowBox[{"d3", "=", RowBox[{"d1", "+", RowBox[{"d2", "*", RowBox[{"Cos", "[", "\[Alpha]", "]"}]}]}]}], ";", " ", RowBox[{"d4", "=", RowBox[{"d2", "*", RowBox[{"Sin", "[", "\[Alpha]", "]"}]}]}], ";", " ", RowBox[{"d5", "=", RowBox[{"d4", "*", RowBox[{ RowBox[{"Sqrt", "[", "3", "]"}], "/", "2"}]}]}], ";"}], "\n", RowBox[{"(*", RowBox[{ RowBox[{"p", "=", RowBox[{"d1", "+", "d3"}]}], ";", " ", RowBox[{"q", "=", "d1"}], ";", " ", RowBox[{"w4", "=", RowBox[{"pi", "+", "R4"}]}], ";", " ", RowBox[{"w5", "=", RowBox[{"pi", "+", "R5"}]}], ";"}], "*)"}]}], "\n", RowBox[{ RowBox[{"xn", "=", RowBox[{"{", RowBox[{ RowBox[{"-", "d1"}], ",", "d1", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0"}], "}"}]}], ";"}], "\n", RowBox[{ RowBox[{"yn", "=", RowBox[{"{", RowBox[{ "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0"}], "}"}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"zn", "=", RowBox[{"{", RowBox[{ "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0", ",", "0"}], "}"}]}], ";"}], "\n", RowBox[{"(*", " ", RowBox[{ "cloud", " ", "coordinates", " ", "in", " ", "terms", " ", "of", " ", "radii"}], " ", "*)"}], " "}], "\n", RowBox[{ RowBox[{"xw", "=", RowBox[{"{", RowBox[{"d1", ",", RowBox[{"-", "d1"}], ",", "0", ",", "d3", ",", " ", "d3", ",", " ", "d3", ",", " ", RowBox[{"-", "d3"}], ",", RowBox[{"-", "d3"}], ",", RowBox[{"-", "d3"}]}], "}"}]}], ";"}], "\n", RowBox[{ RowBox[{"yw", "=", RowBox[{"{", RowBox[{"0", ",", "0", ",", "0", ",", "d4", ",", " ", RowBox[{ RowBox[{"-", "d4"}], "/", "2"}], ",", RowBox[{ RowBox[{"-", "d4"}], "/", "2"}], ",", RowBox[{"-", "d4"}], ",", " ", RowBox[{"d4", "/", "2"}], ",", RowBox[{"d4", "/", "2"}]}], "}"}]}], ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"zw", "=", RowBox[{"{", RowBox[{"0", ",", "0", ",", "0", ",", "0", ",", " ", "d5", ",", RowBox[{"-", "d5"}], ",", " ", "0", ",", " ", "d5", ",", " ", RowBox[{"-", "d5"}]}], "}"}]}], ";"}], "\n", RowBox[{"(*", " ", RowBox[{"cloud", " ", "distances", " ", "and", " ", "reciprocal"}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"tw", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xw", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"yw", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"zw", "[", RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "nc"}], "}"}]}], "]"}]}], ";"}], " ", "\n", RowBox[{"(*", RowBox[{"MatrixForm", "[", "tw", "]"}], "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"dw", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{ RowBox[{"sq", "=", RowBox[{"Norm", "[", RowBox[{ RowBox[{"tw", "[", RowBox[{"[", "i", "]"}], "]"}], "-", RowBox[{"tw", "[", RowBox[{"[", "j", "]"}], "]"}]}], "]"}]}], ";", RowBox[{"If", "[", RowBox[{"(*", RowBox[{"(", RowBox[{ RowBox[{"sq", ">", "0"}], " ", "&&"}]}], "*)"}], " ", RowBox[{ RowBox[{"i", "\[NotEqual]", "j"}], ",", RowBox[{"1", "/", "sq"}], ",", "0"}], "]"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "nc"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", "nc"}], "}"}]}], "]"}]}], ";"}], "\n", RowBox[{"(*", RowBox[{"MatrixForm", "[", "dw", "]"}], "*)"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ RowBox[{"nuclear", " ", "distances", " ", "and", " ", "cloud"}], "-", RowBox[{"nucleus", " ", "reciprocal", " ", "distances"}]}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"tv", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xn", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"yn", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"zn", "[", RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "nc"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"MatrixForm", "[", "tv", "]"}], "*)"}]}], "\n", RowBox[{ RowBox[{ RowBox[{"dcw", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{ RowBox[{"sq", "=", RowBox[{"Norm", "[", RowBox[{ RowBox[{"tw", "[", RowBox[{"[", "i", "]"}], "]"}], "-", RowBox[{"tv", "[", RowBox[{"[", "j", "]"}], "]"}]}], "]"}]}], ";", RowBox[{"If", "[", RowBox[{"(*", RowBox[{"(", RowBox[{ RowBox[{"sq", ">", " ", "0.0"}], " ", "&&"}]}], "*)"}], " ", RowBox[{ RowBox[{"i", "\[NotEqual]", "j"}], ",", RowBox[{"1", "/", "sq"}], ",", "0"}], "]"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "nc"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", "nc"}], "}"}]}], "]"}]}], ";"}], "\n", RowBox[{"(*", RowBox[{"MatrixForm", "[", "dcw", "]"}], "*)"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{"reciprocal", " ", "nuclear", " ", "distances"}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"dcc", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{ RowBox[{"sq", "=", RowBox[{"Norm", "[", RowBox[{ RowBox[{"tv", "[", RowBox[{"[", "i", "]"}], "]"}], "-", RowBox[{"tv", "[", RowBox[{"[", "j", "]"}], "]"}]}], "]"}]}], ";", RowBox[{"If", "[", RowBox[{"(*", RowBox[{"(", RowBox[{ RowBox[{"sq", ">", " ", "0.0"}], " ", "&&"}]}], "*)"}], " ", RowBox[{ RowBox[{"i", "\[NotEqual]", "j"}], ",", RowBox[{"1", "/", "sq"}], ",", "0"}], "]"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "nc"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", "nc"}], "}"}]}], "]"}]}], ";"}], "\n", RowBox[{"(*", RowBox[{"MatrixForm", "[", "dcc", "]"}], "*)"}], "\n", RowBox[{"(*", " ", RowBox[{ RowBox[{"sum", " ", "of", " ", "cloud"}], "-", RowBox[{"nucleus", " ", "potential", " ", "energies"}]}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"vne", "=", RowBox[{"vne", "+", RowBox[{"cn", ".", RowBox[{"(", RowBox[{"ow", ".", "dcw"}], ")"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ RowBox[{ RowBox[{"sum", " ", "of", " ", "cloud"}], "-", RowBox[{"cloud", " ", "potential", " ", "energies"}]}], ",", " ", RowBox[{"symm", ".", " ", "matrix"}], ",", " ", RowBox[{"hence", " ", "0.5"}]}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"vee", "=", RowBox[{"vee", "+", RowBox[{"0.5", "*", RowBox[{"ow", ".", RowBox[{"(", RowBox[{"ow", ".", "dw"}], ")"}]}]}]}]}], ";"}], "\n", RowBox[{"(*", " ", RowBox[{ RowBox[{ RowBox[{"sum", " ", "of", " ", "nucleus"}], "-", RowBox[{"nucleus", " ", "potential", " ", "energies"}]}], ",", " ", RowBox[{"symm", " ", "matrix", " ", "hence", " ", "0.5"}]}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"vnn", " ", "=", RowBox[{"0.5", "*", RowBox[{"cn", ".", RowBox[{"(", RowBox[{"cn", ".", "dcc"}], ")"}]}]}]}], ";"}], "\[IndentingNewLine]"}], "\n", RowBox[{ RowBox[{"Epot", "=", RowBox[{ RowBox[{"vne", "+", "vee", "+", "vnn"}], " ", "/.", " ", "c"}]}], ";"}], "\n", RowBox[{"func", "=", RowBox[{ RowBox[{"Ekin", "+", "Epot"}], " ", "/.", " ", "c"}]}], "\n", RowBox[{"t", " ", "=", " ", RowBox[{"FindMinimum", "[", RowBox[{"func", ",", RowBox[{"{", RowBox[{"R1", ",", "0.148"}], "}"}], ",", RowBox[{"{", RowBox[{"R2", ",", "1.166"}], "}"}], ",", RowBox[{"{", RowBox[{"R3", ",", "0.96"}], "}"}], ",", RowBox[{"{", RowBox[{"\[Alpha]", ",", "1.23"}], "}"}], ",", RowBox[{"Method", " ", "\[Rule]", "\"\\""}], ",", " ", RowBox[{"MaxIterations", " ", "->", " ", "500"}]}], "]"}]}]}], "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.5773415047046976`*^9, 3.5773415137059135`*^9}, { 3.5773416513449545`*^9, 3.577341669752987*^9}, {3.5773417239318824`*^9, 3.577341893956581*^9}, 3.5773419407566633`*^9, {3.5773420175555983`*^9, 3.577342056165666*^9}, {3.5773420929973307`*^9, 3.5773421056489525`*^9}, { 3.5773421396414127`*^9, 3.5773423470593767`*^9}, {3.577342382081438*^9, 3.577342441158742*^9}, {3.5773425028880506`*^9, 3.577342523682887*^9}, { 3.5773425980794177`*^9, 3.5773426626791315`*^9}, {3.5773426933487854`*^9, 3.5773426941131864`*^9}, {3.577342738354864*^9, 3.577342889035529*^9}, { 3.5773429412644205`*^9, 3.57734295234044*^9}, {3.577342991824109*^9, 3.5773430231333647`*^9}, {3.5773430561274223`*^9, 3.577343143737176*^9}, { 3.5773431993668737`*^9, 3.577343258334977*^9}, 3.5773434426249013`*^9, { 3.577343480042567*^9, 3.5773435125686235`*^9}, {3.5773436063903885`*^9, 3.5773436272008247`*^9}, {3.5773436617392855`*^9, 3.5773437808610945`*^9}, {3.577343868193248*^9, 3.577343879659268*^9}, { 3.5773439214829416`*^9, 3.577344118838888*^9}, {3.5773441503509436`*^9, 3.577344170740179*^9}, {3.5773446920694942`*^9, 3.577344693286297*^9}, 3.5773526894956493`*^9, {3.5773528388503113`*^9, 3.577352838897112*^9}, { 3.577353457328998*^9, 3.577353503364679*^9}, {3.5773556679320807`*^9, 3.5773557165573664`*^9}, 3.5773607112650003`*^9, {3.577363145384075*^9, 3.5773631641821084`*^9}, 3.5773713130222435`*^9, {3.5773735851350346`*^9, 3.577373695474028*^9}, {3.577373964730501*^9, 3.5773739677725067`*^9}, 3.5773740596254673`*^9, {3.577374092245125*^9, 3.5773741045067463`*^9}, { 3.5773741853616886`*^9, 3.577374230414568*^9}, 3.5773742637986264`*^9, { 3.5773746139412413`*^9, 3.5773746339092765`*^9}, {3.5773750495252066`*^9, 3.5773750633000307`*^9}, 3.5773754561555204`*^9, {3.5773757271123962`*^9, 3.577375728063998*^9}, 3.577375963078411*^9, {3.5773761205918875`*^9, 3.5773762048008356`*^9}, 3.5774291695117903`*^9, {3.5774424964722877`*^9, 3.5774424977046895`*^9}, {3.5774428197828555`*^9, 3.577442825336465*^9}, { 3.577442927267044*^9, 3.5774429281406455`*^9}, {3.5774429766411304`*^9, 3.577442992833959*^9}, {3.5774430583072743`*^9, 3.577443066715689*^9}, { 3.5774436977055974`*^9, 3.5774437122604227`*^9}, 3.5774438419902506`*^9, { 3.577446912512444*^9, 3.577446913651246*^9}, {3.577447169460495*^9, 3.577447179070112*^9}, {3.577449454744052*^9, 3.57744946514927*^9}, { 3.577449615377534*^9, 3.5774496616940155`*^9}, {3.577450578242425*^9, 3.5774506245901065`*^9}, {3.577451618688603*^9, 3.5774516195310044`*^9}, { 3.5774516668614874`*^9, 3.5774516695602922`*^9}, 3.577455410676982*^9, { 3.5774557532847834`*^9, 3.5774558029396706`*^9}, {3.5774565586985984`*^9, 3.5774565589793987`*^9}, 3.577457075496306*^9, {3.577457117273179*^9, 3.577457148941235*^9}, {3.577678508293971*^9, 3.5776785428012314`*^9}, { 3.5776816489512606`*^9, 3.5776816557372727`*^9}, 3.5776817445014286`*^9, { 3.5776818073071394`*^9, 3.5776818117219467`*^9}, 3.577681885338476*^9, { 3.577682125485298*^9, 3.577682126124899*^9}, 3.5776821818169966`*^9, { 3.577682250082717*^9, 3.5776822522355204`*^9}, {3.5776823382072716`*^9, 3.5776823388624725`*^9}, 3.577682369656927*^9, {3.577682527513604*^9, 3.5776825496656427`*^9}, {3.577682803010088*^9, 3.5776828162233114`*^9}, { 3.577684829687648*^9, 3.5776848354440575`*^9}, 3.5776850034095526`*^9, { 3.5776850335332055`*^9, 3.577685052362439*^9}, {3.5776850854032965`*^9, 3.5776852462551794`*^9}, {3.5776852913080587`*^9, 3.577685293601262*^9}, { 3.577685324068116*^9, 3.577685477260385*^9}, 3.5776855432329006`*^9, 3.5776855922013865`*^9, {3.577685645615881*^9, 3.577685645803081*^9}, { 3.5776856793431396`*^9, 3.5776856800139413`*^9}, {3.5776857127739983`*^9, 3.577685783754123*^9}, {3.5776858189009852`*^9, 3.5776858890075083`*^9}, { 3.5776859215491652`*^9, 3.577685994728894*^9}, 3.577686107969493*^9, { 3.5776872286366615`*^9, 3.5776873647301006`*^9}, {3.5776874640554748`*^9, 3.5776874877519164`*^9}, {3.5776877853380394`*^9, 3.577687807490078*^9}, { 3.577697174278837*^9, 3.57769723305974*^9}, {3.577698938949665*^9, 3.5776989710545216`*^9}, 3.5776992675018425`*^9, {3.5776993025083036`*^9, 3.577699341367972*^9}, {3.5776994640153875`*^9, 3.577699496120244*^9}, { 3.577699693788191*^9, 3.5776997265794487`*^9}, {3.5776997671863203`*^9, 3.5776997677167206`*^9}, {3.577699803175583*^9, 3.577699803783984*^9}, { 3.5776998675880966`*^9, 3.5776998682744975`*^9}, {3.577699911829774*^9, 3.5776999680366726`*^9}, {3.577770340813854*^9, 3.5777703575682836`*^9}, { 3.57777040109236*^9, 3.577770466144474*^9}, {3.577770587325487*^9, 3.5777706167159386`*^9}, {3.577770670614033*^9, 3.5777706894432664`*^9}, { 3.5777707698926077`*^9, 3.577770801389063*^9}, {3.577770851387151*^9, 3.577770892009622*^9}, {3.5777713137875633`*^9, 3.5777713421328125`*^9}, 3.5777713898844967`*^9, {3.577771482189859*^9, 3.5777715441375675`*^9}, { 3.5777716319345217`*^9, 3.5777717046306496`*^9}, {3.57777208670652*^9, 3.5777720884069233`*^9}, {3.5777722492900057`*^9, 3.5777722501324077`*^9}, {3.577772298820093*^9, 3.577772347928979*^9}, { 3.577772387069448*^9, 3.5777724011094728`*^9}, {3.5777724628699813`*^9, 3.5777724640867834`*^9}, {3.5777725348173075`*^9, 3.577772535472509*^9}, { 3.5777726008678236`*^9, 3.577772767304516*^9}, {3.577772817692604*^9, 3.577772840390644*^9}, {3.577772881746317*^9, 3.5777728836963205`*^9}, 3.5777729249271927`*^9, {3.5777729796208887`*^9, 3.5777730045497327`*^9}, 3.5777730453750043`*^9, 3.5777730778542614`*^9, {3.5777731197559347`*^9, 3.577773144934379*^9}, {3.577773188598856*^9, 3.577773204058483*^9}, { 3.5777735446850815`*^9, 3.577773565573518*^9}, {3.5777736641032915`*^9, 3.5777736802649193`*^9}, {3.5777737110749736`*^9, 3.577773757360255*^9}, { 3.5777741482345414`*^9, 3.5777741523061485`*^9}, {3.577774405822194*^9, 3.577774457239884*^9}, 3.577774494149549*^9, 3.577774526738006*^9, { 3.5777745718376856`*^9, 3.577774610338553*^9}, {3.5777746431142106`*^9, 3.5777747969148808`*^9}, 3.5777749559103603`*^9, {3.5777749904956207`*^9, 3.577774991213222*^9}, {3.577775052318529*^9, 3.5777751512539034`*^9}, 3.577775256444888*^9, {3.577775310498983*^9, 3.577775359670269*^9}, { 3.5777754112595596`*^9, 3.5777754119771614`*^9}, 3.5777754682308598`*^9, { 3.577775501942519*^9, 3.577775520007351*^9}, {3.5777812436886044`*^9, 3.5777813396911726`*^9}, {3.577781428736129*^9, 3.5777814336813374`*^9}, { 3.5777878209721565`*^9, 3.5777878541378145`*^9}, {3.5777879075211086`*^9, 3.577787952636388*^9}, {3.5777879932276587`*^9, 3.5777880083440857`*^9}, { 3.577788064004983*^9, 3.5777880923190327`*^9}, {3.5777881342675066`*^9, 3.5777881746091776`*^9}, {3.57778828427737*^9, 3.577788288692178*^9}, { 3.5777883311242523`*^9, 3.5777883470518804`*^9}, {3.5777884558932714`*^9, 3.5777885753894815`*^9}, {3.577788608305539*^9, 3.5777886580852265`*^9}, { 3.577788700283301*^9, 3.57778870563411*^9}, {3.577788765912616*^9, 3.57778885355357*^9}, {3.5777954490267544`*^9, 3.577795455719166*^9}, { 3.5777954887756243`*^9, 3.577795537166909*^9}, {3.5777955732653723`*^9, 3.577795590207002*^9}, {3.5777956205022554`*^9, 3.577795640610691*^9}, { 3.5777957491556816`*^9, 3.577795809886588*^9}, 3.57779591919598*^9, { 3.577795962376856*^9, 3.57779603864539*^9}, {3.577796068847043*^9, 3.5777961433527737`*^9}, {3.5777962295117254`*^9, 3.577796231243328*^9}, 3.5777962883862286`*^9, {3.577796331676305*^9, 3.577796332549906*^9}, 3.57779642014406*^9, {3.5777964865689764`*^9, 3.57779652826785*^9}, { 3.5777965709027247`*^9, 3.5777968014399295`*^9}, {3.577800521219663*^9, 3.577800573729355*^9}, {3.577800657719903*^9, 3.5778006618695097`*^9}, { 3.57780097488406*^9, 3.577800986833681*^9}, {3.5778011009166813`*^9, 3.577801178792018*^9}, {3.5778012158264832`*^9, 3.57780124821214*^9}, { 3.577801344339509*^9, 3.5778013721387577`*^9}, {3.5778014595613112`*^9, 3.5778014797321463`*^9}, 3.5778052599775867`*^9, 3.5778062355408998`*^9, 3.577809012361377*^9, {3.5778093581983843`*^9, 3.577809377916819*^9}, 3.5778509988815928`*^9, {3.5778539201895237`*^9, 3.5778539359923515`*^9}, { 3.5778539765368223`*^9, 3.5778540720401907`*^9}, 3.577854102943845*^9, { 3.577854197823211*^9, 3.5778542710029397`*^9}, {3.5778543785807285`*^9, 3.577854500011342*^9}, {3.577854532428199*^9, 3.5778546388515863`*^9}, { 3.577854680051258*^9, 3.577854727428541*^9}, {3.5778554835306697`*^9, 3.577855532202755*^9}, {3.57787732657312*^9, 3.5778774264444957`*^9}, { 3.5778774894686065`*^9, 3.5778776046436086`*^9}, {3.57788370263052*^9, 3.5778837650618296`*^9}, {3.577883801253893*^9, 3.577883802439495*^9}, { 3.5778838363695545`*^9, 3.577883925133711*^9}, {3.5778839773002024`*^9, 3.577884171801344*^9}, {3.577884202517798*^9, 3.5778842301922464`*^9}, { 3.577884359110873*^9, 3.577884440589816*^9}, {3.577884512412342*^9, 3.5778845286051707`*^9}, {3.577884580568862*^9, 3.5778845858104715`*^9}, { 3.5778883619509506`*^9, 3.5778883699849644`*^9}, {3.577890350798444*^9, 3.5778903689424753`*^9}, {3.5778927844039183`*^9, 3.577892803935152*^9}, { 3.5778929458486013`*^9, 3.5778929587342243`*^9}, 3.577893153843767*^9, { 3.5778932178818793`*^9, 3.5778932411571207`*^9}, 3.5778933148672495`*^9, { 3.5778934240206413`*^9, 3.5778934291686506`*^9}, {3.5778936211113873`*^9, 3.5778936614530582`*^9}, {3.5778937553340235`*^9, 3.5778937598892317`*^9}, {3.5778973974724207`*^9, 3.577897401107227*^9}, 3.577897435505287*^9, {3.577897505767811*^9, 3.577897507624214*^9}, { 3.577897614609202*^9, 3.577897657275277*^9}, {3.577897705572962*^9, 3.5778977077101655`*^9}, {3.577897928902954*^9, 3.5778979754378357`*^9}}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"Power", "::", "infy"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Infinite expression \[NoBreak]\\!\\(1\\/0\\)\[NoBreak] \ encountered. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/message/General/infy\\\", ButtonNote -> \ \\\"Power::infy\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.5778979939706683`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"Power", "::", "infy"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Infinite expression \[NoBreak]\\!\\(1\\/0\\)\[NoBreak] \ encountered. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/message/General/infy\\\", ButtonNote -> \ \\\"Power::infy\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.5778979940174685`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"Power", "::", "infy"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Infinite expression \[NoBreak]\\!\\(1\\/0\\)\[NoBreak] \ encountered. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/message/General/infy\\\", ButtonNote -> \ \\\"Power::infy\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.577897994064268*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"General", "::", "stop"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Further output of \[NoBreak]\\!\\(\\*StyleBox[\\(Power :: \ infy\\), \\\"MessageName\\\"]\\)\[NoBreak] will be suppressed during this \ calculation. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/message/General/stop\\\", ButtonNote -> \ \\\"General::stop\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.5778979940954685`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"Power", "::", "infy"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Infinite expression \[NoBreak]\\!\\(1\\/0\\)\[NoBreak] \ encountered. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/message/General/infy\\\", ButtonNote -> \ \\\"Power::infy\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.577897994142268*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"Power", "::", "infy"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Infinite expression \[NoBreak]\\!\\(1\\/0\\)\[NoBreak] \ encountered. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/message/General/infy\\\", ButtonNote -> \ \\\"Power::infy\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.5778979942046685`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"Power", "::", "infy"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Infinite expression \[NoBreak]\\!\\(1\\/0\\)\[NoBreak] \ encountered. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/message/General/infy\\\", ButtonNote -> \ \\\"Power::infy\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.5778979942514687`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"General", "::", "stop"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Further output of \[NoBreak]\\!\\(\\*StyleBox[\\(Power :: \ infy\\), \\\"MessageName\\\"]\\)\[NoBreak] will be suppressed during this \ calculation. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/message/General/stop\\\", ButtonNote -> \ \\\"General::stop\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.577897994282669*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"Infinity", "::", "indet"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Indeterminate expression \[NoBreak]\\!\\(0\\\\ \ ComplexInfinity\\)\[NoBreak] encountered. \\!\\(\\*ButtonBox[\\\"\ \[RightSkeleton]\\\", ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/message/General/indet\\\", ButtonNote -> \ \\\"Infinity::indet\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.5778979943294687`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"Infinity", "::", "indet"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Indeterminate expression \[NoBreak]\\!\\(0\\\\ \ ComplexInfinity\\)\[NoBreak] encountered. \\!\\(\\*ButtonBox[\\\"\ \[RightSkeleton]\\\", ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/message/General/indet\\\", ButtonNote -> \ \\\"Infinity::indet\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.577897994376269*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"Infinity", "::", "indet"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Indeterminate expression \[NoBreak]\\!\\(0\\\\ \ ComplexInfinity\\)\[NoBreak] encountered. \\!\\(\\*ButtonBox[\\\"\ \[RightSkeleton]\\\", ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/message/General/indet\\\", ButtonNote -> \ \\\"Infinity::indet\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.577897994423069*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"General", "::", "stop"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Further output of \[NoBreak]\\!\\(\\*StyleBox[\\(Infinity \ :: indet\\), \\\"MessageName\\\"]\\)\[NoBreak] will be suppressed during this \ calculation. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/message/General/stop\\\", ButtonNote -> \ \\\"General::stop\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.577897994454269*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"Infinity", "::", "indet"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Indeterminate expression \[NoBreak]\\!\\(0\\\\ \ ComplexInfinity\\)\[NoBreak] encountered. \\!\\(\\*ButtonBox[\\\"\ \[RightSkeleton]\\\", ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/message/General/indet\\\", ButtonNote -> \ \\\"Infinity::indet\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.5778979945166693`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"Infinity", "::", "indet"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Indeterminate expression \[NoBreak]\\!\\(0\\\\ \ ComplexInfinity\\)\[NoBreak] encountered. \\!\\(\\*ButtonBox[\\\"\ \[RightSkeleton]\\\", ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/message/General/indet\\\", ButtonNote -> \ \\\"Infinity::indet\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.577897994563469*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"Infinity", "::", "indet"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Indeterminate expression \[NoBreak]\\!\\(0\\\\ \ ComplexInfinity\\)\[NoBreak] encountered. \\!\\(\\*ButtonBox[\\\"\ \[RightSkeleton]\\\", ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/message/General/indet\\\", ButtonNote -> \ \\\"Infinity::indet\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.5778979946102695`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"General", "::", "stop"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Further output of \[NoBreak]\\!\\(\\*StyleBox[\\(Infinity \ :: indet\\), \\\"MessageName\\\"]\\)\[NoBreak] will be suppressed during this \ calculation. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/message/General/stop\\\", ButtonNote -> \ \\\"General::stop\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.5778979946414695`*^9}], Cell[BoxData["Indeterminate"], "Output", CellChangeTimes->{3.5778979946414695`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"FindMinimum", "::", "nrnum"}], "MessageName"], RowBox[{ ":", " "}], "\<\"The function value \[NoBreak]\\!\\(Indeterminate\\)\ \[NoBreak] is not a real number at \[NoBreak]\\!\\({R1, R2, R3, \[Alpha]}\\)\ \[NoBreak] = \[NoBreak]\\!\\({0.148`, 1.166`, 0.96`, 1.23`}\\)\[NoBreak]. \\!\ \\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", ButtonStyle->\\\"Link\\\", \ ButtonFrame->None, ButtonData:>\\\"paclet:ref/FindMinimum\\\", ButtonNote -> \ \\\"FindMinimum::nrnum\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.5778979947194695`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"FindMinimum", "::", "nrnum"}], "MessageName"], RowBox[{ ":", " "}], "\<\"The function value \[NoBreak]\\!\\(Indeterminate\\)\ \[NoBreak] is not a real number at \[NoBreak]\\!\\({R1, R2, R3, \[Alpha]}\\)\ \[NoBreak] = \[NoBreak]\\!\\({0.148`, 1.166`, 0.96`, 1.23`}\\)\[NoBreak]. \\!\ \\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", ButtonStyle->\\\"Link\\\", \ ButtonFrame->None, ButtonData:>\\\"paclet:ref/FindMinimum\\\", ButtonNote -> \ \\\"FindMinimum::nrnum\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.5778979947974696`*^9}], Cell[BoxData[ RowBox[{"FindMinimum", "[", RowBox[{"func", ",", RowBox[{"{", RowBox[{"R1", ",", "0.148`"}], "}"}], ",", RowBox[{"{", RowBox[{"R2", ",", "1.166`"}], "}"}], ",", RowBox[{"{", RowBox[{"R3", ",", "0.96`"}], "}"}], ",", RowBox[{"{", RowBox[{"\[Alpha]", ",", "1.23`"}], "}"}], ",", RowBox[{"Method", "\[Rule]", "\<\"Newton\"\>"}], ",", RowBox[{"MaxIterations", "\[Rule]", "500"}]}], "]"}]], "Output", CellChangeTimes->{3.5778979947974696`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", "func", " ", "*)"}], "\[IndentingNewLine]", RowBox[{ 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[{"Vir", " ", "=", " ", RowBox[{ RowBox[{ RowBox[{"-", "Epot"}], "/", "Ekin"}], " ", "/.", " ", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}]}], " ", " ", " ", "\n", RowBox[{ RowBox[{"2", "*", "d1", "*", "bohr"}], " ", "/.", " ", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}], " ", " ", " "}]}]}]], "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.5773415047046976`*^9, 3.5773415137059135`*^9}, { 3.5773416513449545`*^9, 3.577341669752987*^9}, {3.5773417239318824`*^9, 3.577341893956581*^9}, 3.5773419407566633`*^9, {3.5773420175555983`*^9, 3.577342056165666*^9}, {3.5773420929973307`*^9, 3.5773421056489525`*^9}, { 3.5773421396414127`*^9, 3.5773423470593767`*^9}, {3.577342382081438*^9, 3.577342441158742*^9}, {3.5773425028880506`*^9, 3.577342523682887*^9}, { 3.5773425980794177`*^9, 3.5773426626791315`*^9}, {3.5773426933487854`*^9, 3.5773426941131864`*^9}, {3.577342738354864*^9, 3.577342889035529*^9}, { 3.5773429412644205`*^9, 3.57734295234044*^9}, {3.577342991824109*^9, 3.5773430231333647`*^9}, {3.5773430561274223`*^9, 3.577343143737176*^9}, { 3.5773431993668737`*^9, 3.577343258334977*^9}, 3.5773434426249013`*^9, { 3.577343480042567*^9, 3.5773435125686235`*^9}, {3.5773436063903885`*^9, 3.5773436272008247`*^9}, {3.5773436617392855`*^9, 3.5773437808610945`*^9}, {3.577343868193248*^9, 3.577343879659268*^9}, { 3.5773439214829416`*^9, 3.577344086609232*^9}, {3.5773444013153844`*^9, 3.5773444098953996`*^9}, {3.5773527214601054`*^9, 3.577352754235763*^9}, { 3.577352839084312*^9, 3.577352839115512*^9}, {3.5773562861299667`*^9, 3.577356291496376*^9}, {3.5773717712106485`*^9, 3.5773718737028284`*^9}, { 3.5774560203260527`*^9, 3.577456044802496*^9}, {3.577678590303315*^9, 3.577678629615384*^9}, {3.577687576656472*^9, 3.577687617044944*^9}, { 3.577703686849204*^9, 3.5777036868648043`*^9}, {3.5777751762451468`*^9, 3.577775176260747*^9}, {3.577795838871439*^9, 3.577795838918239*^9}, { 3.577806042521761*^9, 3.577806042552961*^9}, {3.577888393541006*^9, 3.577888393556606*^9}, {3.577893410713818*^9, 3.577893410760618*^9}, { 3.5778966896835775`*^9, 3.5778966896991777`*^9}, {3.5778975884947557`*^9, 3.5778975955303683`*^9}}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"FindMinimum", "::", "nrnum"}], "MessageName"], RowBox[{ ":", " "}], "\<\"The function value \[NoBreak]\\!\\(Indeterminate\\)\ \[NoBreak] is not a real number at \[NoBreak]\\!\\({R1, R2, R3, \[Alpha]}\\)\ \[NoBreak] = \[NoBreak]\\!\\({0.148`, 1.166`, 0.96`, 1.23`}\\)\[NoBreak]. \\!\ \\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", ButtonStyle->\\\"Link\\\", \ ButtonFrame->None, ButtonData:>\\\"paclet:ref/FindMinimum\\\", ButtonNote -> \ \\\"FindMinimum::nrnum\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.5778979948754697`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.148`}\\)\[NoBreak] is neither a \ list of replacement rules nor a valid dispatch table, and so cannot be used \ for replacing. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/message/ReplaceAll/reps\\\", ButtonNote -> \ \\\"ReplaceAll::reps\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.5778979949066696`*^9}], Cell[BoxData[ RowBox[{ RowBox[{ FractionBox["3.9726`", SuperscriptBox["R1", "2"]], "+", RowBox[{"2.25`", " ", RowBox[{"(", RowBox[{ FractionBox["1.33`", SuperscriptBox["R2", "2"]], "+", FractionBox["7.199999999999999`", SuperscriptBox["R3", "2"]]}], ")"}]}]}], "/.", "\[VeryThinSpace]", RowBox[{"{", RowBox[{"R1", ",", "0.148`"}], "}"}]}]], "Output", CellChangeTimes->{3.5778979949066696`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"FindMinimum", "::", "nrnum"}], "MessageName"], RowBox[{ ":", " "}], "\<\"The function value \[NoBreak]\\!\\(Indeterminate\\)\ \[NoBreak] is not a real number at \[NoBreak]\\!\\({R1, R2, R3, \[Alpha]}\\)\ \[NoBreak] = \[NoBreak]\\!\\({0.148`, 1.166`, 0.96`, 1.23`}\\)\[NoBreak]. \\!\ \\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", ButtonStyle->\\\"Link\\\", \ ButtonFrame->None, ButtonData:>\\\"paclet:ref/FindMinimum\\\", ButtonNote -> \ \\\"FindMinimum::nrnum\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.5778979949846697`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.148`}\\)\[NoBreak] is neither a \ list of replacement rules nor a valid dispatch table, and so cannot be used \ for replacing. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/message/ReplaceAll/reps\\\", ButtonNote -> \ \\\"ReplaceAll::reps\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.57789799501587*^9}], Cell[BoxData[ RowBox[{"Indeterminate", "/.", "\[VeryThinSpace]", RowBox[{"{", RowBox[{"R1", ",", "0.148`"}], "}"}]}]], "Output", CellChangeTimes->{3.57789799501587*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"FindMinimum", "::", "nrnum"}], "MessageName"], RowBox[{ ":", " "}], "\<\"The function value \[NoBreak]\\!\\(Indeterminate\\)\ \[NoBreak] is not a real number at \[NoBreak]\\!\\({R1, R2, R3, \[Alpha]}\\)\ \[NoBreak] = \[NoBreak]\\!\\({0.148`, 1.166`, 0.96`, 1.23`}\\)\[NoBreak]. \\!\ \\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", ButtonStyle->\\\"Link\\\", \ ButtonFrame->None, ButtonData:>\\\"paclet:ref/FindMinimum\\\", ButtonNote -> \ \\\"FindMinimum::nrnum\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.57789799507827*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.148`}\\)\[NoBreak] is neither a \ list of replacement rules nor a valid dispatch table, and so cannot be used \ for replacing. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/message/ReplaceAll/reps\\\", ButtonNote -> \ \\\"ReplaceAll::reps\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.5778979951094704`*^9}], Cell[BoxData[ RowBox[{ RowBox[{ FractionBox["1.9500000000000002`", "R1"], "+", RowBox[{"3.`", " ", RowBox[{"(", RowBox[{ FractionBox["0.325`", "R2"], "+", FractionBox["1.9500000000000002`", "R3"]}], ")"}]}], "+", RowBox[{"0.5`", " ", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", RowBox[{"(", RowBox[{ RowBox[{"-", FractionBox["2", RowBox[{"Abs", "[", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", "R1"}], "-", RowBox[{"2", " ", "R2"}]}], "]"}]]}], "-", FractionBox["2", RowBox[{"Abs", "[", RowBox[{ RowBox[{"-", "R1"}], "-", "R2"}], "]"}]], "-", FractionBox["6.000000000000001`", RowBox[{"Abs", "[", RowBox[{"R1", "+", "R3"}], "]"}]], "-", FractionBox["6", SqrtBox[ RowBox[{ RowBox[{"0.8882851417666464`", " ", SuperscriptBox[ RowBox[{"Abs", "[", RowBox[{"R1", "+", "R3"}], "]"}], "2"]}], "+", SuperscriptBox[ RowBox[{"Abs", "[", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", "R1"}], "-", RowBox[{"2", " ", "R2"}], "-", RowBox[{"0.3342377271245026`", " ", RowBox[{"(", RowBox[{"R1", "+", "R3"}], ")"}]}]}], "]"}], "2"]}]]]}], ")"}]}], "-", RowBox[{"4", " ", RowBox[{"(", RowBox[{ RowBox[{"-", FractionBox["2.0000000000000004`", RowBox[{"Abs", "[", RowBox[{"R1", "+", "R3"}], "]"}]]}], "-", FractionBox["4", SqrtBox[ RowBox[{"0.`", "\[VeryThinSpace]", "+", RowBox[{"2.664855425299939`", " ", SuperscriptBox[ RowBox[{"Abs", "[", RowBox[{"R1", "+", "R3"}], "]"}], "2"]}]}]]], "-", FractionBox["2", SqrtBox[ RowBox[{ RowBox[{"0.8882851417666464`", " ", SuperscriptBox[ RowBox[{"Abs", "[", RowBox[{"R1", "+", "R3"}], "]"}], "2"]}], "+", SuperscriptBox[ RowBox[{"Abs", "[", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", "R1"}], "-", RowBox[{"2", " ", "R2"}], "-", RowBox[{"0.6684754542490052`", " ", RowBox[{"(", RowBox[{"R1", "+", "R3"}], ")"}]}]}], "]"}], "2"]}]]], "-", FractionBox["2", SqrtBox[ RowBox[{"0.`", "\[VeryThinSpace]", "+", RowBox[{"0.8882851417666465`", " ", SuperscriptBox[ RowBox[{"Abs", "[", RowBox[{"R1", "+", "R3"}], "]"}], "2"]}], "+", SuperscriptBox[ RowBox[{"Abs", "[", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", "R1"}], "-", RowBox[{"2", " ", "R2"}], "-", RowBox[{"0.6684754542490052`", " ", RowBox[{"(", RowBox[{"R1", "+", "R3"}], ")"}]}]}], "]"}], "2"]}]]], "-", FractionBox["2", SqrtBox[ RowBox[{ RowBox[{"3.5531405670665857`", " ", SuperscriptBox[ RowBox[{"Abs", "[", RowBox[{"R1", "+", "R3"}], "]"}], "2"]}], "+", SuperscriptBox[ RowBox[{"Abs", "[", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", "R1"}], "-", RowBox[{"2", " ", "R2"}], "-", RowBox[{"0.6684754542490052`", " ", RowBox[{"(", RowBox[{"R1", "+", "R3"}], ")"}]}]}], "]"}], "2"]}]]], "-", FractionBox["2", SqrtBox[ RowBox[{ RowBox[{"0.8882851417666464`", " ", SuperscriptBox[ RowBox[{"Abs", "[", RowBox[{"R1", "+", "R3"}], "]"}], "2"]}], "+", SuperscriptBox[ RowBox[{"Abs", "[", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", "R1"}], "-", RowBox[{"2", " ", "R2"}], "-", RowBox[{"0.3342377271245026`", " ", RowBox[{"(", RowBox[{"R1", "+", "R3"}], ")"}]}]}], "]"}], "2"]}]]], "-", FractionBox["2", SqrtBox[ RowBox[{ RowBox[{"0.8882851417666464`", " ", SuperscriptBox[ RowBox[{"Abs", "[", RowBox[{"R1", "+", "R3"}], "]"}], "2"]}], "+", SuperscriptBox[ RowBox[{"Abs", "[", RowBox[{ RowBox[{"-", "R1"}], "-", "R2", "-", RowBox[{"0.3342377271245026`", " ", RowBox[{"(", RowBox[{"R1", "+", "R3"}], ")"}]}]}], "]"}], "2"]}]]]}], ")"}]}], "-", RowBox[{"2", " ", RowBox[{"(", RowBox[{ RowBox[{"-", FractionBox["2.`", RowBox[{"Abs", "[", RowBox[{"R1", "+", "R3"}], "]"}]]}], "-", FractionBox["4", SqrtBox[ RowBox[{"0.`", "\[VeryThinSpace]", "+", RowBox[{"2.6648554252999395`", " ", SuperscriptBox[ RowBox[{"Abs", "[", RowBox[{"R1", "+", "R3"}], "]"}], "2"]}]}]]], "-", FractionBox["4", SqrtBox[ RowBox[{ RowBox[{"0.8882851417666464`", " ", SuperscriptBox[ RowBox[{"Abs", "[", RowBox[{"R1", "+", "R3"}], "]"}], "2"]}], "+", SuperscriptBox[ RowBox[{"Abs", "[", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", "R1"}], "-", RowBox[{"2", " ", "R2"}], "-", RowBox[{"0.6684754542490052`", " ", RowBox[{"(", RowBox[{"R1", "+", "R3"}], ")"}]}]}], "]"}], "2"]}]]], "-", FractionBox["2", SqrtBox[ RowBox[{ RowBox[{"3.553140567066586`", " ", SuperscriptBox[ RowBox[{"Abs", "[", RowBox[{"R1", "+", "R3"}], "]"}], "2"]}], "+", SuperscriptBox[ RowBox[{"Abs", "[", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", "R1"}], "-", RowBox[{"2", " ", "R2"}], "-", RowBox[{"0.6684754542490052`", " ", RowBox[{"(", RowBox[{"R1", "+", "R3"}], ")"}]}]}], "]"}], "2"]}]]], "-", FractionBox["2", SqrtBox[ RowBox[{ RowBox[{"0.8882851417666465`", " ", SuperscriptBox[ RowBox[{"Abs", "[", RowBox[{"R1", "+", "R3"}], "]"}], "2"]}], "+", SuperscriptBox[ RowBox[{"Abs", "[", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", "R1"}], "-", RowBox[{"2", " ", "R2"}], "-", RowBox[{"0.3342377271245026`", " ", RowBox[{"(", RowBox[{"R1", "+", "R3"}], ")"}]}]}], "]"}], "2"]}]]], "-", FractionBox["2", SqrtBox[ RowBox[{ RowBox[{"0.8882851417666465`", " ", SuperscriptBox[ RowBox[{"Abs", "[", RowBox[{"R1", "+", "R3"}], "]"}], "2"]}], "+", SuperscriptBox[ RowBox[{"Abs", "[", RowBox[{ RowBox[{"-", "R1"}], "-", "R2", "-", RowBox[{"0.3342377271245026`", " ", RowBox[{"(", RowBox[{"R1", "+", "R3"}], ")"}]}]}], "]"}], "2"]}]]]}], ")"}]}], "-", RowBox[{"2", " ", RowBox[{"(", RowBox[{ RowBox[{"-", FractionBox["2", RowBox[{"Abs", "[", RowBox[{ RowBox[{"-", "R1"}], "-", "R2"}], "]"}]]}], "-", FractionBox["2", RowBox[{"Abs", "[", RowBox[{"R1", "+", "R2"}], "]"}]], "-", FractionBox["6", SqrtBox[ RowBox[{ RowBox[{"0.8882851417666464`", " ", SuperscriptBox[ RowBox[{"Abs", "[", RowBox[{"R1", "+", "R3"}], "]"}], "2"]}], "+", SuperscriptBox[ RowBox[{"Abs", "[", RowBox[{ RowBox[{"-", "R1"}], "-", "R2", "-", RowBox[{"0.3342377271245026`", " ", RowBox[{"(", RowBox[{"R1", "+", "R3"}], ")"}]}]}], "]"}], "2"]}]]], "-", FractionBox["6", SqrtBox[ RowBox[{ RowBox[{"0.8882851417666464`", " ", SuperscriptBox[ RowBox[{"Abs", "[", RowBox[{"R1", "+", "R3"}], "]"}], "2"]}], "+", SuperscriptBox[ RowBox[{"Abs", "[", RowBox[{"R1", "+", "R2", "+", RowBox[{"0.3342377271245026`", " ", RowBox[{"(", RowBox[{"R1", "+", "R3"}], ")"}]}]}], "]"}], "2"]}]]]}], ")"}]}], "-", RowBox[{"2", " ", RowBox[{"(", RowBox[{ RowBox[{"-", FractionBox["2", RowBox[{"Abs", "[", RowBox[{"R1", "+", "R2"}], "]"}]]}], "-", FractionBox["2", RowBox[{"Abs", "[", RowBox[{ RowBox[{"2", " ", "R1"}], "+", RowBox[{"2", " ", "R2"}]}], "]"}]], "-", FractionBox["6.000000000000001`", RowBox[{"Abs", "[", RowBox[{"R1", "+", "R3"}], "]"}]], "-", FractionBox["6", SqrtBox[ RowBox[{ RowBox[{"0.8882851417666464`", " ", SuperscriptBox[ RowBox[{"Abs", "[", RowBox[{"R1", "+", "R3"}], "]"}], "2"]}], "+", SuperscriptBox[ RowBox[{"Abs", "[", RowBox[{ RowBox[{"2", " ", "R1"}], "+", RowBox[{"2", " ", "R2"}], "+", RowBox[{"0.3342377271245026`", " ", RowBox[{"(", RowBox[{"R1", "+", "R3"}], ")"}]}]}], "]"}], "2"]}]]]}], ")"}]}], "-", RowBox[{"4", " ", RowBox[{"(", RowBox[{ RowBox[{"-", FractionBox["2.0000000000000004`", RowBox[{"Abs", "[", RowBox[{"R1", "+", "R3"}], "]"}]]}], "-", FractionBox["4", SqrtBox[ RowBox[{"0.`", "\[VeryThinSpace]", "+", RowBox[{"2.664855425299939`", " ", SuperscriptBox[ RowBox[{"Abs", "[", RowBox[{"R1", "+", "R3"}], "]"}], "2"]}]}]]], "-", FractionBox["2", SqrtBox[ RowBox[{ RowBox[{"0.8882851417666464`", " ", SuperscriptBox[ RowBox[{"Abs", "[", RowBox[{"R1", "+", "R3"}], "]"}], "2"]}], "+", SuperscriptBox[ RowBox[{"Abs", "[", RowBox[{"R1", "+", "R2", "+", RowBox[{"0.3342377271245026`", " ", RowBox[{"(", RowBox[{"R1", "+", "R3"}], ")"}]}]}], "]"}], "2"]}]]], "-", FractionBox["2", SqrtBox[ RowBox[{ RowBox[{"0.8882851417666464`", " ", SuperscriptBox[ RowBox[{"Abs", "[", RowBox[{"R1", "+", "R3"}], "]"}], "2"]}], "+", SuperscriptBox[ RowBox[{"Abs", "[", RowBox[{ RowBox[{"2", " ", "R1"}], "+", RowBox[{"2", " ", "R2"}], "+", RowBox[{"0.3342377271245026`", " ", RowBox[{"(", RowBox[{"R1", "+", "R3"}], ")"}]}]}], "]"}], "2"]}]]], "-", FractionBox["2", SqrtBox[ RowBox[{ RowBox[{"0.8882851417666464`", " ", SuperscriptBox[ RowBox[{"Abs", "[", RowBox[{"R1", "+", "R3"}], "]"}], "2"]}], "+", SuperscriptBox[ RowBox[{"Abs", "[", RowBox[{ RowBox[{"2", " ", "R1"}], "+", RowBox[{"2", " ", "R2"}], "+", RowBox[{"0.6684754542490052`", " ", RowBox[{"(", RowBox[{"R1", "+", "R3"}], ")"}]}]}], "]"}], "2"]}]]], "-", FractionBox["2", SqrtBox[ RowBox[{"0.`", "\[VeryThinSpace]", "+", RowBox[{"0.8882851417666465`", " ", SuperscriptBox[ RowBox[{"Abs", "[", RowBox[{"R1", "+", "R3"}], "]"}], "2"]}], "+", SuperscriptBox[ RowBox[{"Abs", "[", RowBox[{ RowBox[{"2", " ", "R1"}], "+", RowBox[{"2", " ", "R2"}], "+", RowBox[{"0.6684754542490052`", " ", RowBox[{"(", RowBox[{"R1", "+", "R3"}], ")"}]}]}], "]"}], "2"]}]]], "-", FractionBox["2", SqrtBox[ RowBox[{ RowBox[{"3.5531405670665857`", " ", SuperscriptBox[ RowBox[{"Abs", "[", RowBox[{"R1", "+", "R3"}], "]"}], "2"]}], "+", SuperscriptBox[ RowBox[{"Abs", "[", RowBox[{ RowBox[{"2", " ", "R1"}], "+", RowBox[{"2", " ", "R2"}], "+", RowBox[{"0.6684754542490052`", " ", RowBox[{"(", RowBox[{"R1", "+", "R3"}], ")"}]}]}], "]"}], "2"]}]]]}], ")"}]}], "-", RowBox[{"2", " ", RowBox[{"(", RowBox[{ RowBox[{"-", FractionBox["2.`", RowBox[{"Abs", "[", RowBox[{"R1", "+", "R3"}], "]"}]]}], "-", FractionBox["4", SqrtBox[ RowBox[{"0.`", "\[VeryThinSpace]", "+", RowBox[{"2.6648554252999395`", " ", SuperscriptBox[ RowBox[{"Abs", "[", RowBox[{"R1", "+", "R3"}], "]"}], "2"]}]}]]], "-", FractionBox["2", SqrtBox[ RowBox[{ RowBox[{"0.8882851417666465`", " ", SuperscriptBox[ RowBox[{"Abs", "[", RowBox[{"R1", "+", "R3"}], "]"}], "2"]}], "+", SuperscriptBox[ RowBox[{"Abs", "[", RowBox[{"R1", "+", "R2", "+", RowBox[{"0.3342377271245026`", " ", RowBox[{"(", RowBox[{"R1", "+", "R3"}], ")"}]}]}], "]"}], "2"]}]]], "-", FractionBox["2", SqrtBox[ RowBox[{ RowBox[{"0.8882851417666465`", " ", SuperscriptBox[ RowBox[{"Abs", "[", RowBox[{"R1", "+", "R3"}], "]"}], "2"]}], "+", SuperscriptBox[ RowBox[{"Abs", "[", RowBox[{ RowBox[{"2", " ", "R1"}], "+", RowBox[{"2", " ", "R2"}], "+", RowBox[{"0.3342377271245026`", " ", RowBox[{"(", RowBox[{"R1", "+", "R3"}], ")"}]}]}], "]"}], "2"]}]]], "-", FractionBox["4", SqrtBox[ RowBox[{ RowBox[{"0.8882851417666464`", " ", SuperscriptBox[ RowBox[{"Abs", "[", RowBox[{"R1", "+", "R3"}], "]"}], "2"]}], "+", SuperscriptBox[ RowBox[{"Abs", "[", RowBox[{ RowBox[{"2", " ", "R1"}], "+", RowBox[{"2", " ", "R2"}], "+", RowBox[{"0.6684754542490052`", " ", RowBox[{"(", RowBox[{"R1", "+", "R3"}], ")"}]}]}], "]"}], "2"]}]]], "-", FractionBox["2", SqrtBox[ RowBox[{ RowBox[{"3.553140567066586`", " ", SuperscriptBox[ RowBox[{"Abs", "[", RowBox[{"R1", "+", "R3"}], "]"}], "2"]}], "+", SuperscriptBox[ RowBox[{"Abs", "[", RowBox[{ RowBox[{"2", " ", "R1"}], "+", RowBox[{"2", " ", "R2"}], "+", RowBox[{"0.6684754542490052`", " ", RowBox[{"(", RowBox[{"R1", "+", "R3"}], ")"}]}]}], "]"}], "2"]}]]]}], ")"}]}]}], ")"}]}]}], "/.", "\[VeryThinSpace]", RowBox[{"{", RowBox[{"R1", ",", "0.148`"}], "}"}]}]], "Output", CellChangeTimes->{3.57789799512507*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"FindMinimum", "::", "nrnum"}], "MessageName"], RowBox[{ ":", " "}], "\<\"The function value \[NoBreak]\\!\\(Indeterminate\\)\ \[NoBreak] is not a real number at \[NoBreak]\\!\\({R1, R2, R3, \[Alpha]}\\)\ \[NoBreak] = \[NoBreak]\\!\\({0.148`, 1.166`, 0.96`, 1.23`}\\)\[NoBreak]. \\!\ \\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", ButtonStyle->\\\"Link\\\", \ ButtonFrame->None, ButtonData:>\\\"paclet:ref/FindMinimum\\\", ButtonNote -> \ \\\"FindMinimum::nrnum\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.5778979951874704`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.148`}\\)\[NoBreak] is neither a \ list of replacement rules nor a valid dispatch table, and so cannot be used \ for replacing. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/message/ReplaceAll/reps\\\", ButtonNote -> \ \\\"ReplaceAll::reps\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.5778979952186704`*^9}], Cell[BoxData[ RowBox[{"Indeterminate", "/.", "\[VeryThinSpace]", RowBox[{"{", RowBox[{"R1", ",", "0.148`"}], "}"}]}]], "Output", CellChangeTimes->{3.5778979952186704`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"FindMinimum", "::", "nrnum"}], "MessageName"], RowBox[{ ":", " "}], "\<\"The function value \[NoBreak]\\!\\(Indeterminate\\)\ \[NoBreak] is not a real number at \[NoBreak]\\!\\({R1, R2, R3, \[Alpha]}\\)\ \[NoBreak] = \[NoBreak]\\!\\({0.148`, 1.166`, 0.96`, 1.23`}\\)\[NoBreak]. \\!\ \\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", ButtonStyle->\\\"Link\\\", \ ButtonFrame->None, ButtonData:>\\\"paclet:ref/FindMinimum\\\", ButtonNote -> \ \\\"FindMinimum::nrnum\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.5778979952966704`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.148`}\\)\[NoBreak] is neither a \ list of replacement rules nor a valid dispatch table, and so cannot be used \ for replacing. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/message/ReplaceAll/reps\\\", ButtonNote -> \ \\\"ReplaceAll::reps\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.5778979953278704`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.148`}\\)\[NoBreak] is neither a \ list of replacement rules nor a valid dispatch table, and so cannot be used \ for replacing. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/message/ReplaceAll/reps\\\", ButtonNote -> \ \\\"ReplaceAll::reps\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.5778979953434706`*^9}], Cell[BoxData[ RowBox[{"Indeterminate", "/.", "\[VeryThinSpace]", RowBox[{"{", RowBox[{"R1", ",", "0.148`"}], "}"}]}]], "Output", CellChangeTimes->{3.5778979953434706`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"FindMinimum", "::", "nrnum"}], "MessageName"], RowBox[{ ":", " "}], "\<\"The function value \[NoBreak]\\!\\(Indeterminate\\)\ \[NoBreak] is not a real number at \[NoBreak]\\!\\({R1, R2, R3, \[Alpha]}\\)\ \[NoBreak] = \[NoBreak]\\!\\({0.148`, 1.166`, 0.96`, 1.23`}\\)\[NoBreak]. \\!\ \\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", ButtonStyle->\\\"Link\\\", \ ButtonFrame->None, ButtonData:>\\\"paclet:ref/FindMinimum\\\", ButtonNote -> \ \\\"FindMinimum::nrnum\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.5778979954214706`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.148`}\\)\[NoBreak] is neither a \ list of replacement rules nor a valid dispatch table, and so cannot be used \ for replacing. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/message/ReplaceAll/reps\\\", ButtonNote -> \ \\\"ReplaceAll::reps\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.577897995452671*^9}], Cell[BoxData[ RowBox[{ RowBox[{"1.058354`", " ", RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}]}], "/.", "\[VeryThinSpace]", RowBox[{"{", RowBox[{"R1", ",", "0.148`"}], "}"}]}]], "Output", CellChangeTimes->{3.577897995452671*^9}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{"rr", "=", RowBox[{"{", RowBox[{ "R1", ",", "R1", ",", "R2", ",", "R3", ",", "R3", ",", "R3", ",", "R3", ",", "R3", ",", "R3"}], "}"}]}], " ", RowBox[{"(*", " ", RowBox[{"/.", " ", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}], " ", " ", " ", "*)"}], ";"}]], "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.5773415047046976`*^9, 3.5773415137059135`*^9}, { 3.5773416513449545`*^9, 3.577341669752987*^9}, {3.5773417239318824`*^9, 3.577341893956581*^9}, 3.5773419407566633`*^9, {3.5773420175555983`*^9, 3.577342056165666*^9}, {3.5773420929973307`*^9, 3.5773421056489525`*^9}, { 3.5773421396414127`*^9, 3.5773423470593767`*^9}, {3.577342382081438*^9, 3.577342441158742*^9}, {3.5773425028880506`*^9, 3.577342523682887*^9}, { 3.5773425980794177`*^9, 3.5773426626791315`*^9}, {3.5773426933487854`*^9, 3.5773426941131864`*^9}, {3.577342738354864*^9, 3.577342889035529*^9}, { 3.5773429412644205`*^9, 3.57734295234044*^9}, {3.577342991824109*^9, 3.5773430231333647`*^9}, {3.5773430561274223`*^9, 3.577343143737176*^9}, { 3.5773431993668737`*^9, 3.577343258334977*^9}, 3.5773434426249013`*^9, { 3.577343480042567*^9, 3.5773435125686235`*^9}, {3.5773436063903885`*^9, 3.5773436272008247`*^9}, {3.5773436617392855`*^9, 3.5773437808610945`*^9}, {3.577343868193248*^9, 3.577343879659268*^9}, { 3.5773439214829416`*^9, 3.577344086609232*^9}, {3.5773444013153844`*^9, 3.5773444098953996`*^9}, {3.5773527214601054`*^9, 3.577352754235763*^9}, { 3.577352839084312*^9, 3.577352839115512*^9}, {3.5773562861299667`*^9, 3.577356291496376*^9}, {3.5773717712106485`*^9, 3.577371904793683*^9}, { 3.577372601444106*^9, 3.577372628525754*^9}, 3.577375878432662*^9, 3.577442909545413*^9, {3.5774471930165367`*^9, 3.5774472038429556`*^9}, { 3.5774495179865627`*^9, 3.577449527440179*^9}, {3.577455105602846*^9, 3.577455114822462*^9}, {3.577455475619896*^9, 3.577455493403927*^9}, 3.577678642797407*^9, 3.5777036868804045`*^9, 3.577775176260747*^9, 3.577795838918239*^9, 3.5778060425685606`*^9, {3.5778846007864976`*^9, 3.5778846245921392`*^9}, 3.577888393556606*^9, 3.577893410760618*^9, 3.5778966896991777`*^9, 3.5778975955303683`*^9}], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{ RowBox[{"projection", " ", "on", " ", "xy"}], "-", RowBox[{"plane", " ", "of", " ", "molecule"}]}], " ", "*)"}]], "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.5773413317003937`*^9, 3.5773413362088013`*^9}, { 3.577352839115512*^9, 3.577352839131112*^9}, {3.5773715990175457`*^9, 3.5773716034635534`*^9}, {3.5773716711988726`*^9, 3.5773716785932856`*^9}, {3.5773719148557005`*^9, 3.5773720320119066`*^9}, {3.5773720625879602`*^9, 3.577372089326407*^9}, { 3.5773721384508934`*^9, 3.577372247229884*^9}, {3.5773723618276854`*^9, 3.577372407098965*^9}, {3.577372640428575*^9, 3.57737267169103*^9}, { 3.5773727054338894`*^9, 3.5773727803452206`*^9}, {3.577372855490553*^9, 3.5773728739765854`*^9}, {3.5773729151762576`*^9, 3.577372935331493*^9}, 3.577373068493327*^9}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"plot0", "=", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"Graphics3D", "[", RowBox[{"{", RowBox[{ RowBox[{"Darker", "[", RowBox[{"Green", ",", "0.3"}], "]"}], ",", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xw", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"yw", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"zw", "[", RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", RowBox[{"rr", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}], "}"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "2"}], "}"}]}], "]"}], " ", "/.", " ", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}]}], " ", " ", " ", ";"}], "\n", RowBox[{ RowBox[{"plot1", "=", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"Graphics3D", "[", RowBox[{"{", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xw", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"yw", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"zw", "[", RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", RowBox[{"rr", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}], "}"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "3", ",", "3"}], "}"}]}], "]"}], " ", "/.", " ", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}]}], " ", " ", " ", ";"}], "\n", RowBox[{ RowBox[{ RowBox[{"plot2", "=", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"Graphics3D", "[", RowBox[{"{", RowBox[{ RowBox[{"Opacity", "[", "0.5", "]"}], ",", RowBox[{"Lighter", "[", RowBox[{"Green", ",", "0.4"}], "]"}], ",", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xw", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"yw", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"zw", "[", RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", RowBox[{"rr", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}], "}"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "4", ",", "9"}], "}"}]}], "]"}], " ", "/.", " ", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}]}], " ", " ", " ", ";"}], "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{"plot2a", "=", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"Graphics3D", "[", RowBox[{"{", RowBox[{ RowBox[{"Opacity", "[", "0.5", "]"}], ",", RowBox[{"Lighter", "[", RowBox[{"Cyan", ",", "0.3"}], "]"}], ",", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xw", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"yw", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"zw", "[", RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", RowBox[{"rr", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}], "}"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "10", ",", "nc"}], "}"}]}], "]"}], " ", "/.", " ", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}]}], " ", " ", " ", ";"}], "*)"}], "\n"}], "\[IndentingNewLine]", RowBox[{"Show", "[", RowBox[{"plot0", ",", "plot1", ",", "plot2", ",", RowBox[{"(*", RowBox[{"plot2a", ","}], "*)"}], RowBox[{"{", RowBox[{ RowBox[{"AspectRatio", " ", "\[Rule]", " ", "Automatic"}], ",", RowBox[{"Axes", " ", "->", " ", "True"}], ",", RowBox[{"PlotRange", " ", "\[Rule]", " ", "All"}], ",", " ", RowBox[{"AxesLabel", " ", "->", " ", RowBox[{"{", RowBox[{"xl", ",", " ", "yl", ",", " ", "zl"}], "}"}]}], ",", RowBox[{"SphericalRegion", "\[Rule]", "True"}]}], "}"}]}], "]"}], "\n"}], "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.5773413317003937`*^9, 3.5773413362088013`*^9}, { 3.577352839115512*^9, 3.577352839131112*^9}, {3.5773715990175457`*^9, 3.5773716034635534`*^9}, {3.5773716711988726`*^9, 3.5773716785932856`*^9}, {3.5773719148557005`*^9, 3.5773720320119066`*^9}, {3.5773720625879602`*^9, 3.577372089326407*^9}, { 3.5773721384508934`*^9, 3.577372247229884*^9}, {3.5773723618276854`*^9, 3.577372407098965*^9}, {3.577372640428575*^9, 3.57737267169103*^9}, { 3.5773727054338894`*^9, 3.5773727803452206`*^9}, {3.577372855490553*^9, 3.5773728739765854`*^9}, {3.5773729151762576`*^9, 3.577372935331493*^9}, 3.577373068493327*^9, 3.577375863940237*^9, {3.577442873244149*^9, 3.5774428788913593`*^9}, {3.5774433556593966`*^9, 3.5774433846442475`*^9}, {3.5774434296971264`*^9, 3.5774434305863285`*^9}, {3.5774434802880154`*^9, 3.5774435257464952`*^9}, 3.577443563108561*^9, 3.5774436149318523`*^9, {3.577443975932086*^9, 3.5774439769772882`*^9}, {3.5774443030334606`*^9, 3.5774444499077187`*^9}, {3.577444584567155*^9, 3.5774445963139753`*^9}, { 3.577444660835689*^9, 3.577444676014516*^9}, {3.577444738118225*^9, 3.577444777929495*^9}, 3.5774449044145164`*^9, {3.5774460837297883`*^9, 3.577446116255845*^9}, {3.577446165957533*^9, 3.5774461885775723`*^9}, { 3.577446237358858*^9, 3.5774463455450478`*^9}, {3.577446401939147*^9, 3.5774464300347967`*^9}, 3.5774465521986113`*^9, {3.577446588484275*^9, 3.5774466383887625`*^9}, 3.5774507241494813`*^9, {3.577451828867772*^9, 3.5774518404585924`*^9}, {3.5774520011856747`*^9, 3.577452007347686*^9}, { 3.5774555692980604`*^9, 3.5774555760684724`*^9}, {3.5774556601526203`*^9, 3.5774556771410503`*^9}, {3.5774566838108177`*^9, 3.577456688662427*^9}, { 3.577678653530226*^9, 3.5776786628122425`*^9}, {3.57767878666086*^9, 3.5776787972532787`*^9}, {3.577687666855831*^9, 3.577687709365906*^9}, { 3.5777036868804045`*^9, 3.5777036870988045`*^9}, {3.577771101829991*^9, 3.5777711034523935`*^9}, {3.577771732835499*^9, 3.5777717348791027`*^9}, { 3.577775176276347*^9, 3.577775176494747*^9}, {3.577795838918239*^9, 3.5777958391678395`*^9}, {3.5778019220553236`*^9, 3.5778019563129835`*^9}, {3.577806042584161*^9, 3.5778060428025613`*^9}, { 3.5778842585374966`*^9, 3.5778843062891803`*^9}, {3.5778883935722055`*^9, 3.5778883937906065`*^9}, {3.577893410760618*^9, 3.5778934107918186`*^9}, { 3.5778966897147775`*^9, 3.577896689948778*^9}, {3.577897595545968*^9, 3.5778975955615683`*^9}}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"FindMinimum", "::", "nrnum"}], "MessageName"], RowBox[{ ":", " "}], "\<\"The function value \[NoBreak]\\!\\(Indeterminate\\)\ \[NoBreak] is not a real number at \[NoBreak]\\!\\({R1, R2, R3, \[Alpha]}\\)\ \[NoBreak] = \[NoBreak]\\!\\({0.148`, 1.166`, 0.96`, 1.23`}\\)\[NoBreak]. \\!\ \\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", ButtonStyle->\\\"Link\\\", \ ButtonFrame->None, ButtonData:>\\\"paclet:ref/FindMinimum\\\", ButtonNote -> \ \\\"FindMinimum::nrnum\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.577897995561871*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.148`}\\)\[NoBreak] is neither a \ list of replacement rules nor a valid dispatch table, and so cannot be used \ for replacing. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/message/ReplaceAll/reps\\\", ButtonNote -> \ \\\"ReplaceAll::reps\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.577897995593071*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.148`}\\)\[NoBreak] is neither a \ list of replacement rules nor a valid dispatch table, and so cannot be used \ for replacing. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/message/ReplaceAll/reps\\\", ButtonNote -> \ \\\"ReplaceAll::reps\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.5778979956242714`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"FindMinimum", "::", "nrnum"}], "MessageName"], RowBox[{ ":", " "}], "\<\"The function value \[NoBreak]\\!\\(Indeterminate\\)\ \[NoBreak] is not a real number at \[NoBreak]\\!\\({R1, R2, R3, \[Alpha]}\\)\ \[NoBreak] = \[NoBreak]\\!\\({0.148`, 1.166`, 0.96`, 1.23`}\\)\[NoBreak]. \\!\ \\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", ButtonStyle->\\\"Link\\\", \ ButtonFrame->None, ButtonData:>\\\"paclet:ref/FindMinimum\\\", ButtonNote -> \ \\\"FindMinimum::nrnum\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.5778979957022715`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.148`}\\)\[NoBreak] is neither a \ list of replacement rules nor a valid dispatch table, and so cannot be used \ for replacing. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/message/ReplaceAll/reps\\\", ButtonNote -> \ \\\"ReplaceAll::reps\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.577897995717871*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.148`}\\)\[NoBreak] is neither a \ list of replacement rules nor a valid dispatch table, and so cannot be used \ for replacing. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/message/ReplaceAll/reps\\\", ButtonNote -> \ \\\"ReplaceAll::reps\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.577897995749071*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"FindMinimum", "::", "nrnum"}], "MessageName"], RowBox[{ ":", " "}], "\<\"The function value \[NoBreak]\\!\\(Indeterminate\\)\ \[NoBreak] is not a real number at \[NoBreak]\\!\\({R1, R2, R3, \[Alpha]}\\)\ \[NoBreak] = \[NoBreak]\\!\\({0.148`, 1.166`, 0.96`, 1.23`}\\)\[NoBreak]. \\!\ \\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", ButtonStyle->\\\"Link\\\", \ ButtonFrame->None, ButtonData:>\\\"paclet:ref/FindMinimum\\\", ButtonNote -> \ \\\"FindMinimum::nrnum\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.577897995827071*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.148`}\\)\[NoBreak] is neither a \ list of replacement rules nor a valid dispatch table, and so cannot be used \ for replacing. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/message/ReplaceAll/reps\\\", ButtonNote -> \ \\\"ReplaceAll::reps\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.5778979958582716`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.148`}\\)\[NoBreak] is neither a \ list of replacement rules nor a valid dispatch table, and so cannot be used \ for replacing. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/message/ReplaceAll/reps\\\", ButtonNote -> \ \\\"ReplaceAll::reps\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.577897995873872*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.148`}\\)\[NoBreak] is neither a \ list of replacement rules nor a valid dispatch table, and so cannot be used \ for replacing. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/message/ReplaceAll/reps\\\", ButtonNote -> \ \\\"ReplaceAll::reps\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.5778979959050713`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.148`}\\)\[NoBreak] is neither a \ list of replacement rules nor a valid dispatch table, and so cannot be used \ for replacing. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/message/ReplaceAll/reps\\\", ButtonNote -> \ \\\"ReplaceAll::reps\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.5778979959362717`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"Show", "::", "gcomb"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Could not combine the graphics objects in \ \[NoBreak]\\!\\(Show[\\(\\(\\(\\(\\(\\({\\*Graphics3DBox[List[RGBColor[0.`, \ 0.7`, 0.`], SphereBox[List[Plus[Skeleton[2]], 0, 0], R1]]], \ \\*Graphics3DBox[List[RGBColor[0.`, 0.7`, 0.`], \ SphereBox[List[Plus[Skeleton[2]], 0, 0], R1]]]}\\)\\)/.\[VeryThinSpace] \ \\(\\({R1, 0.148`}\\)\\)\\)\\), \ \\(\\(\\(\\({\\*Graphics3DBox[List[SphereBox[List[0, 0, 0], R2]]]}\\)\\)/.\ \[VeryThinSpace] \\(\\({R1, 0.148`}\\)\\)\\)\\), \ \\(\\(\\(\\({\\*Graphics3DBox[List[Opacity[0.5`], RGBColor[0.4`, 1.`, 0.4`], \ SphereBox[List[Plus[Skeleton[3]], Times[Skeleton[2]], 0], R3]]], \ \\*Graphics3DBox[List[Opacity[0.5`], RGBColor[0.4`, 1.`, 0.4`], \ SphereBox[List[Plus[Skeleton[3]], Times[Skeleton[2]], Times[Skeleton[2]]], \ R3]]], \\*Graphics3DBox[List[Opacity[0.5`], RGBColor[0.4`, 1.`, 0.4`], \ SphereBox[List[Plus[Skeleton[3]], Times[Skeleton[2]], Times[Skeleton[2]]], \ R3]]], \\*Graphics3DBox[List[Opacity[0.5`], RGBColor[0.4`, 1.`, 0.4`], \ SphereBox[List[Plus[Skeleton[3]], Times[Skeleton[2]], 0], R3]]], \ \\*Graphics3DBox[List[Opacity[0.5`], RGBColor[0.4`, 1.`, 0.4`], \ SphereBox[List[Plus[Skeleton[3]], Times[Skeleton[2]], Times[Skeleton[2]]], \ R3]]], \\*Graphics3DBox[List[Opacity[0.5`], RGBColor[0.4`, 1.`, 0.4`], \ SphereBox[List[Plus[Skeleton[3]], Times[Skeleton[2]], Times[Skeleton[2]]], \ R3]]]}\\)\\)/.\[VeryThinSpace] \\(\\({R1, 0.148`}\\)\\)\\)\\), \ \\(\\({\\(\\(AspectRatio \[Rule] Automatic\\)\\), \\(\\(Axes \[Rule] \ True\\)\\), \\(\\(PlotRange \[Rule] All\\)\\), \\(\\(AxesLabel \[Rule] \ \\(\\({xl, yl, zl}\\)\\)\\)\\), \\(\\(SphericalRegion \[Rule] \ True\\)\\)}\\)\\)\\)\\)]\\)\[NoBreak]. \ \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", ButtonStyle->\\\"Link\\\", \ ButtonFrame->None, ButtonData:>\\\"paclet:ref/message/Show/gcomb\\\", \ ButtonNote -> \\\"Show::gcomb\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.5778979959674716`*^9}], Cell[BoxData[ RowBox[{"Show", "[", RowBox[{ RowBox[{ RowBox[{"{", RowBox[{ Graphics3DBox[ {RGBColor[0., 0.7, 0.], SphereBox[{$CellContext`R1 + $CellContext`R2, 0, 0}, $CellContext`R1]}], ",", Graphics3DBox[ {RGBColor[0., 0.7, 0.], SphereBox[{-$CellContext`R1 - $CellContext`R2, 0, 0}, $CellContext`R1]}]}], "}"}], "/.", "\[VeryThinSpace]", RowBox[{"{", RowBox[{"R1", ",", "0.148`"}], "}"}]}], ",", RowBox[{ RowBox[{"{", Graphics3DBox[SphereBox[{0, 0, 0}, $CellContext`R2]], "}"}], "/.", "\[VeryThinSpace]", RowBox[{"{", RowBox[{"R1", ",", "0.148`"}], "}"}]}], ",", RowBox[{ RowBox[{"{", RowBox[{ Graphics3DBox[ {RGBColor[0.4, 1., 0.4], Opacity[0.5], SphereBox[{$CellContext`R1 + $CellContext`R2 + 0.3342377271245026 ($CellContext`R1 + $CellContext`R3), 0.9424888019316975 ($CellContext`R1 + $CellContext`R3), 0}, $CellContext`R3]}], ",", Graphics3DBox[ {RGBColor[0.4, 1., 0.4], Opacity[0.5], SphereBox[{$CellContext`R1 + $CellContext`R2 + 0.3342377271245026 ($CellContext`R1 + $CellContext`R3), \ (-0.47124440096584874`) ($CellContext`R1 + $CellContext`R3), 0.8162192452552101 ($CellContext`R1 + $CellContext`R3)}, \ $CellContext`R3]}], ",", Graphics3DBox[ {RGBColor[0.4, 1., 0.4], Opacity[0.5], SphereBox[{$CellContext`R1 + $CellContext`R2 + 0.3342377271245026 ($CellContext`R1 + $CellContext`R3), \ (-0.47124440096584874`) ($CellContext`R1 + $CellContext`R3), \ (-0.8162192452552101) ($CellContext`R1 + $CellContext`R3)}, \ $CellContext`R3]}], ",", Graphics3DBox[ {RGBColor[0.4, 1., 0.4], Opacity[0.5], SphereBox[{-$CellContext`R1 - $CellContext`R2 - 0.3342377271245026 ($CellContext`R1 + $CellContext`R3), \ (-0.9424888019316975) ($CellContext`R1 + $CellContext`R3), 0}, $CellContext`R3]}], ",", Graphics3DBox[ {RGBColor[0.4, 1., 0.4], Opacity[0.5], SphereBox[{-$CellContext`R1 - $CellContext`R2 - 0.3342377271245026 ($CellContext`R1 + $CellContext`R3), 0.47124440096584874` ($CellContext`R1 + $CellContext`R3), 0.8162192452552101 ($CellContext`R1 + $CellContext`R3)}, \ $CellContext`R3]}], ",", Graphics3DBox[ {RGBColor[0.4, 1., 0.4], Opacity[0.5], SphereBox[{-$CellContext`R1 - $CellContext`R2 - 0.3342377271245026 ($CellContext`R1 + $CellContext`R3), 0.47124440096584874` ($CellContext`R1 + $CellContext`R3), \ (-0.8162192452552101) ($CellContext`R1 + $CellContext`R3)}, \ $CellContext`R3]}]}], "}"}], "/.", "\[VeryThinSpace]", RowBox[{"{", RowBox[{"R1", ",", "0.148`"}], "}"}]}], ",", RowBox[{"{", RowBox[{ RowBox[{"AspectRatio", "\[Rule]", "Automatic"}], ",", RowBox[{"Axes", "\[Rule]", "True"}], ",", RowBox[{"PlotRange", "\[Rule]", "All"}], ",", RowBox[{"AxesLabel", "\[Rule]", RowBox[{"{", RowBox[{"xl", ",", "yl", ",", "zl"}], "}"}]}], ",", RowBox[{"SphericalRegion", "\[Rule]", "True"}]}], "}"}]}], "]"}]], "Output", CellChangeTimes->{3.5778979959830713`*^9}] }, Open ]], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{ RowBox[{"projection", " ", "on", " ", "xy"}], "-", "plane"}], " ", "*)"}]], "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.5773413317003937`*^9, 3.5773413362088013`*^9}, { 3.577352839115512*^9, 3.577352839131112*^9}, {3.5773715990175457`*^9, 3.5773716034635534`*^9}, {3.5773716711988726`*^9, 3.5773716785932856`*^9}, { 3.5773719148557005`*^9, 3.5773720320119066`*^9}, {3.5773720625879602`*^9, 3.577372089326407*^9}, {3.5773721384508934`*^9, 3.577372247229884*^9}, { 3.5773723618276854`*^9, 3.577372407098965*^9}, {3.577372640428575*^9, 3.57737267169103*^9}, {3.5773727054338894`*^9, 3.5773727803452206`*^9}, { 3.577372855490553*^9, 3.5773728739765854`*^9}, {3.5773729151762576`*^9, 3.577372935331493*^9}, {3.577373068493327*^9, 3.577373077494543*^9}, { 3.5774504640814247`*^9, 3.5774504745490427`*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"plot3", "=", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"Graphics", "[", RowBox[{"{", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"xw", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"yw", "[", RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", RowBox[{"rr", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}], "}"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "nc"}], "}"}]}], "]"}], " ", "/.", " ", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}]}], " ", " ", " ", ";", RowBox[{"Show", "[", RowBox[{"plot3", ",", RowBox[{"{", RowBox[{ RowBox[{"AspectRatio", " ", "\[Rule]", " ", "Automatic"}], ",", RowBox[{"Axes", " ", "->", " ", "True"}], ",", RowBox[{"GridLines", " ", "->", " ", "Automatic"}], ",", " ", RowBox[{"PlotRange", " ", "\[Rule]", " ", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "4"}], ",", "4"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "3"}], ",", "3"}], "}"}]}], "}"}]}], ",", " ", RowBox[{"Frame", " ", "->", " ", "True"}]}], "}"}]}], "]"}]}], "\[IndentingNewLine]"}]], "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.5773413317003937`*^9, 3.5773413362088013`*^9}, { 3.577352839115512*^9, 3.577352839131112*^9}, {3.5773715990175457`*^9, 3.5773716034635534`*^9}, {3.5773716711988726`*^9, 3.5773716785932856`*^9}, {3.5773719148557005`*^9, 3.5773720320119066`*^9}, {3.5773720625879602`*^9, 3.577372089326407*^9}, { 3.5773721384508934`*^9, 3.577372247229884*^9}, {3.5773723618276854`*^9, 3.577372407098965*^9}, {3.577372640428575*^9, 3.57737267169103*^9}, { 3.5773727054338894`*^9, 3.5773727803452206`*^9}, {3.577372855490553*^9, 3.5773728739765854`*^9}, {3.5773729151762576`*^9, 3.577372935331493*^9}, { 3.577373068493327*^9, 3.577373077494543*^9}, 3.577375870741849*^9, { 3.5774428830877666`*^9, 3.577442884476169*^9}, {3.5774502155573883`*^9, 3.5774503400456066`*^9}, {3.5774503829300823`*^9, 3.5774503991853104`*^9}, 3.577678826113329*^9, {3.5776877166823187`*^9, 3.577687719771124*^9}, 3.5777036871144047`*^9, 3.5777751765103474`*^9, 3.5777958391678395`*^9, { 3.5778019656730003`*^9, 3.577801968293805*^9}, 3.577806042818161*^9, 3.577888393806206*^9, 3.5778934107918186`*^9, 3.577896689948778*^9, 3.5778975955615683`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"FindMinimum", "::", "nrnum"}], "MessageName"], RowBox[{ ":", " "}], "\<\"The function value \[NoBreak]\\!\\(Indeterminate\\)\ \[NoBreak] is not a real number at \[NoBreak]\\!\\({R1, R2, R3, \[Alpha]}\\)\ \[NoBreak] = \[NoBreak]\\!\\({0.148`, 1.166`, 0.96`, 1.23`}\\)\[NoBreak]. \\!\ \\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", ButtonStyle->\\\"Link\\\", \ ButtonFrame->None, ButtonData:>\\\"paclet:ref/FindMinimum\\\", ButtonNote -> \ \\\"FindMinimum::nrnum\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.577897996076672*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.148`}\\)\[NoBreak] is neither a \ list of replacement rules nor a valid dispatch table, and so cannot be used \ for replacing. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/message/ReplaceAll/reps\\\", ButtonNote -> \ \\\"ReplaceAll::reps\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.577897996107872*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.148`}\\)\[NoBreak] is neither a \ list of replacement rules nor a valid dispatch table, and so cannot be used \ for replacing. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/message/ReplaceAll/reps\\\", ButtonNote -> \ \\\"ReplaceAll::reps\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.577897996139072*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"Show", "::", "gtype"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\(ReplaceAll\\)\[NoBreak] is not a type of \ graphics. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/message/Show/gtype\\\", ButtonNote -> \ \\\"Show::gtype\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.5778979961546717`*^9}], Cell[BoxData[ RowBox[{"Show", "[", RowBox[{ RowBox[{ RowBox[{"{", RowBox[{ GraphicsBox[ CircleBox[{$CellContext`R1 + $CellContext`R2, 0}, $CellContext`R1]], ",", GraphicsBox[ CircleBox[{-$CellContext`R1 - $CellContext`R2, 0}, $CellContext`R1]], ",", GraphicsBox[CircleBox[{0, 0}, $CellContext`R2]], ",", GraphicsBox[ CircleBox[{$CellContext`R1 + $CellContext`R2 + 0.3342377271245026 ($CellContext`R1 + $CellContext`R3), 0.9424888019316975 ($CellContext`R1 + $CellContext`R3)}, \ $CellContext`R3]], ",", GraphicsBox[ CircleBox[{$CellContext`R1 + $CellContext`R2 + 0.3342377271245026 ($CellContext`R1 + $CellContext`R3), \ (-0.47124440096584874`) ($CellContext`R1 + $CellContext`R3)}, \ $CellContext`R3]], ",", GraphicsBox[ CircleBox[{$CellContext`R1 + $CellContext`R2 + 0.3342377271245026 ($CellContext`R1 + $CellContext`R3), \ (-0.47124440096584874`) ($CellContext`R1 + $CellContext`R3)}, \ $CellContext`R3]], ",", GraphicsBox[ CircleBox[{-$CellContext`R1 - $CellContext`R2 - 0.3342377271245026 ($CellContext`R1 + $CellContext`R3), \ (-0.9424888019316975) ($CellContext`R1 + $CellContext`R3)}, $CellContext`R3]], ",", GraphicsBox[ CircleBox[{-$CellContext`R1 - $CellContext`R2 - 0.3342377271245026 ($CellContext`R1 + $CellContext`R3), 0.47124440096584874` ($CellContext`R1 + $CellContext`R3)}, \ $CellContext`R3]], ",", GraphicsBox[ CircleBox[{-$CellContext`R1 - $CellContext`R2 - 0.3342377271245026 ($CellContext`R1 + $CellContext`R3), 0.47124440096584874` ($CellContext`R1 + $CellContext`R3)}, \ $CellContext`R3]]}], "}"}], "/.", "\[VeryThinSpace]", RowBox[{"{", RowBox[{"R1", ",", "0.148`"}], "}"}]}], ",", RowBox[{"{", RowBox[{ RowBox[{"AspectRatio", "\[Rule]", "Automatic"}], ",", RowBox[{"Axes", "\[Rule]", "True"}], ",", RowBox[{"GridLines", "\[Rule]", "Automatic"}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "4"}], ",", "4"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "3"}], ",", "3"}], "}"}]}], "}"}]}], ",", RowBox[{"Frame", "\[Rule]", "True"}]}], "}"}]}], "]"}]], "Output", CellChangeTimes->{3.577897996170272*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ RowBox[{"projection", " ", "on", " ", "xz"}], "-", "plane"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"plot4", "=", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"Graphics", "[", RowBox[{"{", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"yw", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"zw", "[", RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", RowBox[{"rr", "[", RowBox[{"[", "i", "]"}], "]"}]}], "]"}], "}"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "6"}], "}"}]}], "]"}], " ", "/.", " ", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}]}], " ", " ", " ", ";"}], "\n", RowBox[{"Show", "[", RowBox[{"plot4", ",", RowBox[{"{", RowBox[{ RowBox[{"AspectRatio", " ", "\[Rule]", " ", "Automatic"}], ",", RowBox[{"Axes", " ", "->", " ", "True"}], ",", RowBox[{"GridLines", " ", "->", " ", "Automatic"}], ",", " ", RowBox[{"PlotRange", " ", "\[Rule]", " ", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "4"}], ",", "4"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "3"}], ",", "3"}], "}"}]}], "}"}]}], ",", " ", RowBox[{"Frame", " ", "->", " ", "True"}]}], "}"}]}], "]"}]}]}]], "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.5773413317003937`*^9, 3.5773413362088013`*^9}, { 3.577352839115512*^9, 3.577352839131112*^9}, {3.5773715990175457`*^9, 3.5773716034635534`*^9}, {3.5773716711988726`*^9, 3.5773716785932856`*^9}, {3.5773719148557005`*^9, 3.5773720320119066`*^9}, {3.5773720625879602`*^9, 3.577372089326407*^9}, { 3.5773721384508934`*^9, 3.577372247229884*^9}, {3.5773723618276854`*^9, 3.577372407098965*^9}, {3.577372640428575*^9, 3.57737267169103*^9}, { 3.5773727054338894`*^9, 3.5773727803452206`*^9}, {3.577372855490553*^9, 3.5773728739765854`*^9}, {3.5773729151762576`*^9, 3.577372935331493*^9}, { 3.577373068493327*^9, 3.577373077494543*^9}, 3.577375870741849*^9, { 3.5774428830877666`*^9, 3.577442884476169*^9}, {3.5774502155573883`*^9, 3.5774503400456066`*^9}, {3.5774503829300823`*^9, 3.5774503991853104`*^9}, 3.577450442132186*^9, {3.577450485905863*^9, 3.577450491911873*^9}, 3.5776788392953525`*^9, {3.577687724622732*^9, 3.57768772900634*^9}, 3.5777036871144047`*^9, 3.5777751765103474`*^9, 3.5777958391678395`*^9, { 3.577801978543023*^9, 3.5778019826146297`*^9}, 3.577806042818161*^9, 3.577888393806206*^9, 3.5778934107918186`*^9, 3.577896689964378*^9, 3.5778971958420663`*^9, 3.5778972409885454`*^9, 3.5778975955771685`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"FindMinimum", "::", "nrnum"}], "MessageName"], RowBox[{ ":", " "}], "\<\"The function value \[NoBreak]\\!\\(Indeterminate\\)\ \[NoBreak] is not a real number at \[NoBreak]\\!\\({R1, R2, R3, \[Alpha]}\\)\ \[NoBreak] = \[NoBreak]\\!\\({0.148`, 1.166`, 0.96`, 1.23`}\\)\[NoBreak]. \\!\ \\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", ButtonStyle->\\\"Link\\\", \ ButtonFrame->None, ButtonData:>\\\"paclet:ref/FindMinimum\\\", ButtonNote -> \ \\\"FindMinimum::nrnum\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.577897996263872*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.148`}\\)\[NoBreak] is neither a \ list of replacement rules nor a valid dispatch table, and so cannot be used \ for replacing. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/message/ReplaceAll/reps\\\", ButtonNote -> \ \\\"ReplaceAll::reps\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.5778979962794724`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.148`}\\)\[NoBreak] is neither a \ list of replacement rules nor a valid dispatch table, and so cannot be used \ for replacing. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/message/ReplaceAll/reps\\\", ButtonNote -> \ \\\"ReplaceAll::reps\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.5778979963106723`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"Show", "::", "gtype"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\(ReplaceAll\\)\[NoBreak] is not a type of \ graphics. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/message/Show/gtype\\\", ButtonNote -> \ \\\"Show::gtype\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.577897996341872*^9}], Cell[BoxData[ RowBox[{"Show", "[", RowBox[{ RowBox[{ RowBox[{"{", RowBox[{ GraphicsBox[CircleBox[{0, 0}, $CellContext`R1]], ",", GraphicsBox[CircleBox[{0, 0}, $CellContext`R1]], ",", GraphicsBox[CircleBox[{0, 0}, $CellContext`R2]], ",", GraphicsBox[ CircleBox[{ 0.9424888019316975 ($CellContext`R1 + $CellContext`R3), 0}, $CellContext`R3]], ",", GraphicsBox[ CircleBox[{(-0.47124440096584874`) ($CellContext`R1 + $CellContext`R3), 0.8162192452552101 ($CellContext`R1 + $CellContext`R3)}, \ $CellContext`R3]], ",", GraphicsBox[ CircleBox[{(-0.47124440096584874`) ($CellContext`R1 + \ $CellContext`R3), (-0.8162192452552101) ($CellContext`R1 + $CellContext`R3)}, \ $CellContext`R3]]}], "}"}], "/.", "\[VeryThinSpace]", RowBox[{"{", RowBox[{"R1", ",", "0.148`"}], "}"}]}], ",", RowBox[{"{", RowBox[{ RowBox[{"AspectRatio", "\[Rule]", "Automatic"}], ",", RowBox[{"Axes", "\[Rule]", "True"}], ",", RowBox[{"GridLines", "\[Rule]", "Automatic"}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "4"}], ",", "4"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "3"}], ",", "3"}], "}"}]}], "}"}]}], ",", RowBox[{"Frame", "\[Rule]", "True"}]}], "}"}]}], "]"}]], "Output", CellChangeTimes->{3.5778979963574724`*^9}] }, Open ]], Cell[BoxData[" "], "Input", CellChangeTimes->{3.577701861864399*^9}] }, WindowToolbars->"EditBar", WindowSize->{1023, 861}, WindowMargins->{{Automatic, 374}, {Automatic, 5}}, PrintingCopies->1, PrintingPageRange->{Automatic, Automatic}, PrivateNotebookOptions->{"VersionedStylesheet"->{"Default.nb"[8.] -> False}}, 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, 7088, 106, 172, "Input"], Cell[CellGroupData[{ Cell[7670, 130, 22389, 510, 952, "Input"], Cell[30062, 642, 447, 10, 40, "Message"], Cell[30512, 654, 447, 10, 40, "Message"], Cell[30962, 666, 445, 10, 40, "Message"], Cell[31410, 678, 530, 11, 21, "Message"], Cell[31943, 691, 445, 10, 40, "Message"], Cell[32391, 703, 447, 10, 40, "Message"], Cell[32841, 715, 447, 10, 40, "Message"], Cell[33291, 727, 528, 11, 21, "Message"], Cell[33822, 740, 477, 10, 21, "Message"], Cell[34302, 752, 475, 10, 21, "Message"], Cell[34780, 764, 475, 10, 21, "Message"], Cell[35258, 776, 532, 11, 21, "Message"], Cell[35793, 789, 477, 10, 21, "Message"], Cell[36273, 801, 475, 10, 21, "Message"], Cell[36751, 813, 477, 10, 21, "Message"], Cell[37231, 825, 534, 11, 21, "Message"], Cell[37768, 838, 84, 1, 31, "Output"], Cell[37855, 841, 587, 11, 21, "Message"], Cell[38445, 854, 587, 11, 21, "Message"], Cell[39035, 867, 501, 13, 31, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[39573, 885, 2879, 52, 152, "Input"], Cell[42455, 939, 587, 11, 21, "Message"], Cell[43045, 952, 540, 11, 21, "Message"], Cell[43588, 965, 456, 14, 46, "Output"], Cell[44047, 981, 587, 11, 21, "Message"], Cell[44637, 994, 537, 11, 21, "Message"], Cell[45177, 1007, 177, 4, 31, "Output"], Cell[45357, 1013, 584, 11, 21, "Message"], Cell[45944, 1026, 540, 11, 21, "Message"], Cell[46487, 1039, 16735, 450, 1287, "Output"], Cell[63225, 1491, 587, 11, 21, "Message"], Cell[63815, 1504, 540, 11, 21, "Message"], Cell[64358, 1517, 180, 4, 31, "Output"], Cell[64541, 1523, 587, 11, 21, "Message"], Cell[65131, 1536, 540, 11, 21, "Message"], Cell[65674, 1549, 540, 11, 21, "Message"], Cell[66217, 1562, 180, 4, 31, "Output"], Cell[66400, 1568, 587, 11, 21, "Message"], Cell[66990, 1581, 538, 11, 21, "Message"], Cell[67531, 1594, 248, 7, 31, "Output"] }, Open ]], Cell[67794, 1604, 2282, 37, 31, "Input"], Cell[70079, 1643, 859, 15, 31, "Input"], Cell[CellGroupData[{ Cell[70963, 1662, 6953, 159, 232, "Input"], Cell[77919, 1823, 585, 11, 21, "Message"], Cell[78507, 1836, 538, 11, 21, "Message"], Cell[79048, 1849, 540, 11, 21, "Message"], Cell[79591, 1862, 587, 11, 21, "Message"], Cell[80181, 1875, 538, 11, 21, "Message"], Cell[80722, 1888, 538, 11, 21, "Message"], Cell[81263, 1901, 585, 11, 21, "Message"], Cell[81851, 1914, 540, 11, 21, "Message"], Cell[82394, 1927, 538, 11, 21, "Message"], Cell[82935, 1940, 540, 11, 21, "Message"], Cell[83478, 1953, 540, 11, 21, "Message"], Cell[84021, 1966, 2028, 33, 264, "Message"], Cell[86052, 2001, 3259, 77, 275, "Output"] }, Open ]], Cell[89326, 2081, 885, 15, 31, "Input"], Cell[CellGroupData[{ Cell[90236, 2100, 2656, 59, 92, "Input"], Cell[92895, 2161, 585, 11, 21, "Message"], Cell[93483, 2174, 538, 11, 21, "Message"], Cell[94024, 2187, 538, 11, 21, "Message"], Cell[94565, 2200, 444, 10, 21, "Message"], Cell[95012, 2212, 2452, 60, 256, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[97501, 2277, 2934, 65, 92, "Input"], Cell[100438, 2344, 585, 11, 21, "Message"], Cell[101026, 2357, 540, 11, 21, "Message"], Cell[101569, 2370, 540, 11, 21, "Message"], Cell[102112, 2383, 442, 10, 21, "Message"], Cell[102557, 2395, 1480, 38, 143, "Output"] }, Open ]], Cell[104052, 2436, 69, 1, 31, "Input"] } ] *) (* End of internal cache information *)