(* 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[ 155694, 4381] NotebookOptionsPosition[ 149533, 4217] NotebookOutlinePosition[ 150093, 4238] CellTagsIndexPosition[ 150050, 4235] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell["\<\ (* H2C=CH2 with spherical bananas: Two fused methanes with CH4 data \ CH4Mt(min).nb 20.06.2012 TEST! *) Clear[z,sig1,sig2,sig3,sig4,k1,k2,k4,nc,R1,R2,R3,R4,w,p,vee,vne,vnn, xc,yc,zc,xn,yn,zn,oc,ch,rr,cs,ss,d1,d2,d3,d4,d5,pi,i,j,d,t,a]; z=6.0;\ \>", "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.566799684588899*^9, 3.56679969039211*^9}, 3.5668006247245502`*^9, {3.577339457029501*^9, 3.577339457575502*^9}, { 3.577339970067602*^9, 3.5773399713936043`*^9}}], Cell["\<\ pi=0.0; (* distance of bananas from C-C axis, pi=0 is regular tetrahedric *) nc=9; (* number of clouds *) sig1=0.3; sig2=0.3; sig4=sig2; sig3=sig2;(* screening const. from e-e \ interaction in doubly occ. clouds *) k1=1.027; k2=1.27; k4=1.656; k3=1.2; (* parameters for kinetic energy of \ clouds; k=1.0 Kimball's lowest value *) bohr=0.529177; rad=57.29578;\ \>", "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.566799684588899*^9, 3.56679969039211*^9}, { 3.5668006247245502`*^9, 3.5668006303405604`*^9}, {3.5668170547514086`*^9, 3.5668170590102158`*^9}, {3.5668172256185083`*^9, 3.566817256584563*^9}, { 3.566817851008007*^9, 3.566817851554008*^9}, {3.56681792693334*^9, 3.5668179459965734`*^9}, {3.566818031469124*^9, 3.566818055025165*^9}, { 3.5668182716003456`*^9, 3.5668183126440177`*^9}, {3.5668184559302692`*^9, 3.5668184565698705`*^9}}], Cell["\<\ (* C He-shells *) Ekin = 2*(2.25*k1/R1^2); vee=2*(3.0*sig1/R1); vne=-2*(3.0*z/R1); (* this is the common edge assumption *) (*R4=R2;*)\ \>", "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.566799684588899*^9, 3.56679969039211*^9}, { 3.5668006247245502`*^9, 3.56680063598777*^9}, {3.566802700978997*^9, 3.56680270817061*^9}}], Cell["\<\ (* bonding pairs *) Ekin = Ekin + 2.25*(4*k2/R2^2+k4/R4^2+k3/R3^2); vee=vee+3.0*(4*sig2/R2+sig3/R3); (* cloud occupation *) oc={-2,-2,-2,-2,-2,-2,-2,-1,-1}; (* nuclear charges for C1,C2,H3,H4,H5,H6, banana7, banana8 *) ch={6,6,1,1,1,1,0,0,0}; (* cloud radii in the same order *) rr={R1,R1,R2,R2,R2,R2,R3,R4,R4}; (* w is half angle between two C-H of CH4, i.e. 109.47\[Degree]/2 *) (*w=ArcCos[-1/3]/2;*) cs=Cos[w]; ss=Sin[w]; (* edge length of tetrahedron of 4 equal clouds *) a=R4; (*a=4*(R1+R2)/Sqrt[6];*) (* 4/Sqrt[6] is also Sqrt[8/3] *) (* x is C-C bond axis, xy plane of molecule *) (* nuclear coordinates in terms of radii; C nucleus assumed in center of \ C(1s) cloud *) d1=R1+R3; (* R1+R3 C-C sigma *) d2=d1+(R1+R2+p)*cs; d3=(R1+R2+p)*ss; xn={-d1,d1,-d2,-d2,d2,d2,0,0,0}; yn={0,0,d3,-d3,-d3,d3,0,0,0}; zn={0,0,0,0,0,0,0,0,0}; (* cloud coordinates in terms of radii *) d4=d1+(R1+R2)*cs; d5=(R1+R2)*ss; xc={-d1,d1,-d4,-d4,d4,d4,0,0,0}; yc={0,0,d5,-d5,-d5,d5,0,0,0}; zc={0,0,0,0,0,0,0,pi+a,-pi-a}; (* potential energy of protons in CH-clouds with eccentricity p *) vne=vne-4*(3-(p/R2)^2)/R2;\ \>", "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.566799684588899*^9, 3.56679969039211*^9}, { 3.5668006247245502`*^9, 3.566800645940588*^9}, {3.566803037549588*^9, 3.5668030837880697`*^9}, {3.5668036983355494`*^9, 3.5668037058079624`*^9}, {3.5668042399061003`*^9, 3.566804247784114*^9}, { 3.5668167100220027`*^9, 3.5668167195380197`*^9}, {3.5668183334856544`*^9, 3.5668184100505886`*^9}, {3.566818441531444*^9, 3.5668184420618453`*^9}, { 3.5668185025743513`*^9, 3.566818580777289*^9}, {3.5668186404941936`*^9, 3.56681867269265*^9}, {3.566819101974204*^9, 3.566819103939808*^9}, { 3.5668193038697586`*^9, 3.5668193198753867`*^9}, {3.566821030995592*^9, 3.566821037282403*^9}, 3.5773407322849407`*^9}], Cell["\<\ (* cc: sum of cloud-cloud potential energies *) For[i = 1, i < nc, i++, For[j = i+1, j < nc+1, j++, vee = vee + \ oc[[i]]*oc[[j]]/Sqrt[(xc[[i]]-xc[[j]])^2+(yc[[i]]-yc[[j]])^2+(zc[[i]]-zc[[j]])\ ^2]]]\ \>", "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.566799684588899*^9, 3.56679969039211*^9}, { 3.5668006247245502`*^9, 3.566800657812209*^9}}], Cell["\<\ (* nn: sum of nuclei-nuclei potential energies *) vnn = 0.0; For[i = 1, i < nc-3, i++, For[j = i+1, j < nc-2, j++, vnn = vnn + \ ch[[i]]*ch[[j]]/Sqrt[(xn[[i]]-xn[[j]])^2+(yn[[i]]-yn[[j]])^2]]]\ \>", "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.566799684588899*^9, 3.56679969039211*^9}, { 3.5668006247245502`*^9, 3.5668006625858173`*^9}, {3.5668191723771276`*^9, 3.566819178835539*^9}}], Cell["\<\ (* cn: sum of cloud-nuclei potential energies *) For[i = 1, i < nc+1, i++, For[j = 1, j < nc-2, \ j++,sq=Sqrt[(xc[[i]]-xn[[j]])^2+(yc[[i]]-yn[[j]])^2+(zc[[i]]-zn[[j]])^2]; If[(i \[NotEqual] j)&&(sq \[NotEqual] 0), vne = vne + oc[[i]]*ch[[j]]/sq]]]\ \>", "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.566799684588899*^9, 3.56679969039211*^9}, { 3.5668006247245502`*^9, 3.5668006669226246`*^9}, {3.5668194563132267`*^9, 3.5668194682316475`*^9}, {3.566819622905919*^9, 3.566819624512722*^9}, { 3.5668196983320513`*^9, 3.5668197141972795`*^9}, 3.5668208452148657`*^9, { 3.57734048694331*^9, 3.577340525880978*^9}, {3.577340561792241*^9, 3.5773405777666693`*^9}, {3.577340862217569*^9, 3.577340866445176*^9}}], Cell["\<\ Epot=vne+vee+vnn; func=Ekin+Epot; (* results of CH4 computation; if this is not available, decomment the \ minimize function *) (* R1=0.2623610; R2=1.2461360; p=0.53986226; *) (* minimization function for R1, R2,R3,R4,w, p *)\ \>", "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.566799684588899*^9, 3.56679969039211*^9}, { 3.5668006247245502`*^9, 3.5668006719614334`*^9}, {3.566802680886162*^9, 3.5668026869077725`*^9}, {3.5668027663743124`*^9, 3.5668027861551466`*^9}, {3.5668031435049744`*^9, 3.5668031532549915`*^9}, {3.566803746040433*^9, 3.5668037659304676`*^9}, { 3.5668040820650234`*^9, 3.5668040971502495`*^9}, {3.56680426798615*^9, 3.5668042822601748`*^9}, {3.566804358279108*^9, 3.5668043613679132`*^9}, { 3.5668048168263135`*^9, 3.566804817497115*^9}, {3.566817988179048*^9, 3.5668180095978856`*^9}, {3.5668186901646805`*^9, 3.566818702270302*^9}, { 3.566820740304682*^9, 3.5668207656859264`*^9}, {3.566820986363914*^9, 3.566820991605523*^9}, {3.5773389203105583`*^9, 3.5773389237581644`*^9}, { 3.5773391147180996`*^9, 3.577339118852107*^9}, {3.5773395247336197`*^9, 3.5773396163213806`*^9}, 3.57733970132593*^9, 3.5773398311181583`*^9}], Cell[CellGroupData[{ Cell["\<\ t = FindMinimum[func,{R1,0.2,0.28},{R2,1.1,1.2},{R3,1.1,1.3},{R4,1.5,1.9},{p,\ 0.5,0.6},{w,1.1,1.4},{Method -> Automatic, MaxIterations -> 800}]\ \>", "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.566799684588899*^9, 3.56679969039211*^9}, { 3.5668006247245502`*^9, 3.5668006719614334`*^9}, {3.566802680886162*^9, 3.5668026869077725`*^9}, {3.5668027663743124`*^9, 3.5668027861551466`*^9}, {3.5668031435049744`*^9, 3.5668031532549915`*^9}, {3.566803746040433*^9, 3.5668037659304676`*^9}, { 3.5668040820650234`*^9, 3.5668040971502495`*^9}, {3.56680426798615*^9, 3.5668042822601748`*^9}, {3.566804358279108*^9, 3.5668043613679132`*^9}, { 3.5668048168263135`*^9, 3.566804817497115*^9}, {3.566817988179048*^9, 3.5668180095978856`*^9}, {3.5668186901646805`*^9, 3.566818702270302*^9}, { 3.566820740304682*^9, 3.5668207656859264`*^9}, {3.566820986363914*^9, 3.566820991605523*^9}, {3.5773389203105583`*^9, 3.5773389237581644`*^9}, { 3.5773391147180996`*^9, 3.577339118852107*^9}, {3.5773395247336197`*^9, 3.5773396163213806`*^9}, 3.57733970132593*^9, {3.5773398311181583`*^9, 3.577339837451769*^9}, {3.5773411439072638`*^9, 3.5773411619876957`*^9}}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"FindMinimum", "::", "lsbrak"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Unable to bracket a minimum along the direction \[NoBreak]\ \\!\\({0, 0, 0.19999999999999996`, 0, 0, 0}\\)\[NoBreak] from the point \ \[NoBreak]\\!\\({0.31620422775675705`, 2.180312896149525`*^12, 1.1`, 1.5`, \ 0.5`, 1.1`}\\)\[NoBreak]. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/FindMinimum\\\", ButtonNote -> \ \\\"FindMinimum::lsbrak\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.577341184482935*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"FindMinimum", "::", "lsbrak"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Unable to bracket a minimum along the direction \[NoBreak]\ \\!\\({0, 0, 0.19999999999999996`, 0, 0, 0}\\)\[NoBreak] from the point \ \[NoBreak]\\!\\({0.31620422775675705`, 2.180312896149525`*^12, 1.1`, 1.5`, \ 0.5`, 1.1`}\\)\[NoBreak]. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/FindMinimum\\\", ButtonNote -> \ \\\"FindMinimum::lsbrak\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.577341184560935*^9}], Cell[BoxData[ RowBox[{"FindMinimum", "[", RowBox[{"func", ",", RowBox[{"{", RowBox[{"R1", ",", "0.2`", ",", "0.28`"}], "}"}], ",", RowBox[{"{", RowBox[{"R2", ",", "1.1`", ",", "1.2`"}], "}"}], ",", RowBox[{"{", RowBox[{"R3", ",", "1.1`", ",", "1.3`"}], "}"}], ",", RowBox[{"{", RowBox[{"R4", ",", "1.5`", ",", "1.9`"}], "}"}], ",", RowBox[{"{", RowBox[{"p", ",", "0.5`", ",", "0.6`"}], "}"}], ",", RowBox[{"{", RowBox[{"w", ",", "1.1`", ",", "1.4`"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Method", "\[Rule]", "Automatic"}], ",", RowBox[{"MaxIterations", "\[Rule]", "800"}]}], "}"}]}], "]"}]], "Output", CellChangeTimes->{3.577341184560935*^9}] }, Open ]], Cell["", "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.566799684588899*^9, 3.56679969039211*^9}, { 3.5668006247245502`*^9, 3.5668006719614334`*^9}, {3.566802680886162*^9, 3.5668026869077725`*^9}, {3.5668027663743124`*^9, 3.5668027861551466`*^9}, {3.5668031435049744`*^9, 3.5668031532549915`*^9}, {3.566803746040433*^9, 3.5668037659304676`*^9}, { 3.5668040820650234`*^9, 3.5668040971502495`*^9}, {3.56680426798615*^9, 3.5668042822601748`*^9}, {3.566804358279108*^9, 3.5668043613679132`*^9}, { 3.5668048168263135`*^9, 3.566804817497115*^9}, {3.566817988179048*^9, 3.5668180095978856`*^9}, {3.5668186901646805`*^9, 3.566818702270302*^9}, { 3.566820740304682*^9, 3.5668207656859264`*^9}, {3.566820986363914*^9, 3.566820991605523*^9}, {3.5773389203105583`*^9, 3.5773389237581644`*^9}, { 3.5773391147180996`*^9, 3.577339118852107*^9}, {3.5773395247336197`*^9, 3.5773396163213806`*^9}, 3.57733970132593*^9, 3.5773398311181583`*^9}], Cell[CellGroupData[{ Cell["\<\ vne /. t[[2]] vee /. t[[2]] vnn /. t[[2]] -Epot/Ekin /. t[[2]] 2*d1*bohr /. t[[2]] (R1+R2+p)*bohr /. t[[2]] (R1+R2)*bohr /. t[[2]] 2*w*rad /. t[[2]] 4*(R1*(R1+2*R4))^(3/2)/(R1+R4)^3 /. t[[2]]\ \>", "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.566799684588899*^9, 3.56679969039211*^9}, { 3.5668006247245502`*^9, 3.5668006905878663`*^9}}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"FindMinimum", "::", "lsbrak"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Unable to bracket a minimum along the direction \[NoBreak]\ \\!\\({0, 0, 0.19999999999999996`, 0, 0, 0}\\)\[NoBreak] from the point \ \[NoBreak]\\!\\({0.31620422775675705`, 2.180312896149525`*^12, 1.1`, 1.5`, \ 0.5`, 1.1`}\\)\[NoBreak]. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/FindMinimum\\\", ButtonNote -> \ \\\"FindMinimum::lsbrak\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.5773411846545353`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.2`, 0.28`}\\)\[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.577341184685735*^9}], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"-", FractionBox["36.`", "R1"]}], "-", FractionBox[ RowBox[{"4", " ", RowBox[{"(", RowBox[{"3", "-", FractionBox[ SuperscriptBox["p", "2"], SuperscriptBox["R2", "2"]]}], ")"}]}], "R2"]}], "/.", "\[VeryThinSpace]", RowBox[{"{", RowBox[{"R1", ",", "0.2`", ",", "0.28`"}], "}"}]}]], "Output", CellChangeTimes->{3.577341184685735*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"FindMinimum", "::", "lsbrak"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Unable to bracket a minimum along the direction \[NoBreak]\ \\!\\({0, 0, 0.19999999999999996`, 0, 0, 0}\\)\[NoBreak] from the point \ \[NoBreak]\\!\\({0.31620422775675705`, 2.180312896149525`*^12, 1.1`, 1.5`, \ 0.5`, 1.1`}\\)\[NoBreak]. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/FindMinimum\\\", ButtonNote -> \ \\\"FindMinimum::lsbrak\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.5773411847637353`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.2`, 0.28`}\\)\[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.5773411847793355`*^9}], Cell[BoxData[ RowBox[{ RowBox[{ FractionBox["1.7999999999999998`", "R1"], "+", RowBox[{"3.`", " ", RowBox[{"(", RowBox[{ FractionBox["1.2`", "R2"], "+", FractionBox["0.3`", "R3"]}], ")"}]}], "+", FractionBox["4", SqrtBox[ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", "R1"}], "-", RowBox[{"2", " ", "R3"}]}], ")"}], "2"]]], "+", FractionBox["4", SqrtBox[ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"-", "R1"}], "-", "R3"}], ")"}], "2"]]], "+", FractionBox["4", SqrtBox[ SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R3"}], ")"}], "2"]]], "+", FractionBox["2", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"-", "R1"}], "-", "R3"}], ")"}], "2"], "+", SuperscriptBox[ RowBox[{"(", RowBox[{"0.`", "\[VeryThinSpace]", "-", "R4"}], ")"}], "2"]}]]], "+", FractionBox["2", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R3"}], ")"}], "2"], "+", SuperscriptBox[ RowBox[{"(", RowBox[{"0.`", "\[VeryThinSpace]", "-", "R4"}], ")"}], "2"]}]]], "+", FractionBox["2", SqrtBox[ SuperscriptBox[ RowBox[{"(", RowBox[{"0.`", "\[VeryThinSpace]", "-", "R4"}], ")"}], "2"]]], "+", FractionBox["2", SqrtBox[ SuperscriptBox[ RowBox[{"(", RowBox[{"0.`", "\[VeryThinSpace]", "+", "R4"}], ")"}], "2"]]], "+", FractionBox["1", SqrtBox[ SuperscriptBox[ RowBox[{"(", RowBox[{"0.`", "\[VeryThinSpace]", "+", RowBox[{"2", " ", "R4"}]}], ")"}], "2"]]], "+", FractionBox["2", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"-", "R1"}], "-", "R3"}], ")"}], "2"], "+", SuperscriptBox[ RowBox[{"(", RowBox[{"0.`", "\[VeryThinSpace]", "+", "R4"}], ")"}], "2"]}]]], "+", FractionBox["2", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R3"}], ")"}], "2"], "+", SuperscriptBox[ RowBox[{"(", RowBox[{"0.`", "\[VeryThinSpace]", "+", "R4"}], ")"}], "2"]}]]], "+", FractionBox["8", SqrtBox[ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", "R1"}], "-", RowBox[{"2", " ", "R3"}], "-", RowBox[{"2", " ", RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"]]], "+", FractionBox["4", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]]], "+", FractionBox["16", SqrtBox[ RowBox[{ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Cos", "[", "w", "]"}], "2"]}], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}]]], "+", FractionBox["8", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", "R1"}], "-", RowBox[{"2", " ", "R3"}], "-", RowBox[{ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}]]], "+", FractionBox["8", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"-", "R1"}], "-", "R3", "-", RowBox[{ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}]]], "+", FractionBox["4", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"0.`", "\[VeryThinSpace]", "-", "R4"}], ")"}], "2"], "+", SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"-", "R1"}], "-", "R3", "-", RowBox[{ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}]]], "+", FractionBox["4", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"0.`", "\[VeryThinSpace]", "+", "R4"}], ")"}], "2"], "+", SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"-", "R1"}], "-", "R3", "-", RowBox[{ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}]]], "+", FractionBox["8", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R3", "+", RowBox[{ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}]]], "+", FractionBox["4", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"0.`", "\[VeryThinSpace]", "-", "R4"}], ")"}], "2"], "+", SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R3", "+", RowBox[{ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}]]], "+", FractionBox["4", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"0.`", "\[VeryThinSpace]", "+", "R4"}], ")"}], "2"], "+", SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R3", "+", RowBox[{ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}]]], "+", FractionBox["8", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"2", " ", "R1"}], "+", RowBox[{"2", " ", "R3"}], "+", RowBox[{ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}]]], "+", FractionBox["8", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", "R1"}], "-", RowBox[{"2", " ", "R3"}], "-", RowBox[{"2", " ", RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{"4", " ", SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}]]]}], "/.", "\[VeryThinSpace]", RowBox[{"{", RowBox[{"R1", ",", "0.2`", ",", "0.28`"}], "}"}]}]], "Output", CellChangeTimes->{3.5773411847949357`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"FindMinimum", "::", "lsbrak"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Unable to bracket a minimum along the direction \[NoBreak]\ \\!\\({0, 0, 0.19999999999999996`, 0, 0, 0}\\)\[NoBreak] from the point \ \[NoBreak]\\!\\({0.31620422775675705`, 2.180312896149525`*^12, 1.1`, 1.5`, \ 0.5`, 1.1`}\\)\[NoBreak]. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/FindMinimum\\\", ButtonNote -> \ \\\"FindMinimum::lsbrak\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.5773411848729353`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.2`, 0.28`}\\)\[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.5773411848885355`*^9}], Cell[BoxData[ RowBox[{ RowBox[{"0.`", "\[VeryThinSpace]", "+", FractionBox["36", SqrtBox[ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", "R1"}], "-", RowBox[{"2", " ", "R3"}]}], ")"}], "2"]]], "+", FractionBox["2", SqrtBox[ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", "R1"}], "-", RowBox[{"2", " ", "R3"}], "-", RowBox[{"2", " ", RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"]]], "+", FractionBox["1", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]]], "+", FractionBox["24", SqrtBox[ RowBox[{ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Cos", "[", "w", "]"}], "2"]}], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}]]], "+", FractionBox["12", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", "R1"}], "-", RowBox[{"2", " ", "R3"}], "-", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}]]], "+", FractionBox["12", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"2", " ", "R1"}], "+", RowBox[{"2", " ", "R3"}], "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}]]], "+", FractionBox["2", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", "R1"}], "-", RowBox[{"2", " ", "R3"}], "-", RowBox[{"2", " ", RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{"4", " ", SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}]]]}], "/.", "\[VeryThinSpace]", RowBox[{"{", RowBox[{"R1", ",", "0.2`", ",", "0.28`"}], "}"}]}]], "Output", CellChangeTimes->{3.5773411848885355`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"FindMinimum", "::", "lsbrak"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Unable to bracket a minimum along the direction \[NoBreak]\ \\!\\({0, 0, 0.19999999999999996`, 0, 0, 0}\\)\[NoBreak] from the point \ \[NoBreak]\\!\\({0.31620422775675705`, 2.180312896149525`*^12, 1.1`, 1.5`, \ 0.5`, 1.1`}\\)\[NoBreak]. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/FindMinimum\\\", ButtonNote -> \ \\\"FindMinimum::lsbrak\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.5773411849665356`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.2`, 0.28`}\\)\[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.577341184997736*^9}], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"0.`", "\[VeryThinSpace]", "-", FractionBox["34.2`", "R1"], "-", FractionBox[ RowBox[{"4", " ", RowBox[{"(", RowBox[{"3", "-", FractionBox[ SuperscriptBox["p", "2"], SuperscriptBox["R2", "2"]]}], ")"}]}], "R2"], "+", RowBox[{"3.`", " ", RowBox[{"(", RowBox[{ FractionBox["1.2`", "R2"], "+", FractionBox["0.3`", "R3"]}], ")"}]}], "+", FractionBox["40", SqrtBox[ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", "R1"}], "-", RowBox[{"2", " ", "R3"}]}], ")"}], "2"]]], "+", FractionBox["4", SqrtBox[ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"-", "R1"}], "-", "R3"}], ")"}], "2"]]], "+", FractionBox["4", SqrtBox[ SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R3"}], ")"}], "2"]]], "+", FractionBox["2", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"-", "R1"}], "-", "R3"}], ")"}], "2"], "+", SuperscriptBox[ RowBox[{"(", RowBox[{"0.`", "\[VeryThinSpace]", "-", "R4"}], ")"}], "2"]}]]], "+", FractionBox["2", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R3"}], ")"}], "2"], "+", SuperscriptBox[ RowBox[{"(", RowBox[{"0.`", "\[VeryThinSpace]", "-", "R4"}], ")"}], "2"]}]]], "+", FractionBox["2", SqrtBox[ SuperscriptBox[ RowBox[{"(", RowBox[{"0.`", "\[VeryThinSpace]", "-", "R4"}], ")"}], "2"]]], "+", FractionBox["2", SqrtBox[ SuperscriptBox[ RowBox[{"(", RowBox[{"0.`", "\[VeryThinSpace]", "+", "R4"}], ")"}], "2"]]], "+", FractionBox["1", SqrtBox[ SuperscriptBox[ RowBox[{"(", RowBox[{"0.`", "\[VeryThinSpace]", "+", RowBox[{"2", " ", "R4"}]}], ")"}], "2"]]], "+", FractionBox["2", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"-", "R1"}], "-", "R3"}], ")"}], "2"], "+", SuperscriptBox[ RowBox[{"(", RowBox[{"0.`", "\[VeryThinSpace]", "+", "R4"}], ")"}], "2"]}]]], "+", FractionBox["2", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R3"}], ")"}], "2"], "+", SuperscriptBox[ RowBox[{"(", RowBox[{"0.`", "\[VeryThinSpace]", "+", "R4"}], ")"}], "2"]}]]], "+", FractionBox["8", SqrtBox[ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", "R1"}], "-", RowBox[{"2", " ", "R3"}], "-", RowBox[{"2", " ", RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"]]], "+", FractionBox["2", SqrtBox[ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", "R1"}], "-", RowBox[{"2", " ", "R3"}], "-", RowBox[{"2", " ", RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"]]], "+", FractionBox["4", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]]], "+", FractionBox["1", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]]], "+", FractionBox["16", SqrtBox[ RowBox[{ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Cos", "[", "w", "]"}], "2"]}], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}]]], "+", FractionBox["8", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", "R1"}], "-", RowBox[{"2", " ", "R3"}], "-", RowBox[{ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}]]], "+", FractionBox["8", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"-", "R1"}], "-", "R3", "-", RowBox[{ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}]]], "+", FractionBox["4", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"0.`", "\[VeryThinSpace]", "-", "R4"}], ")"}], "2"], "+", SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"-", "R1"}], "-", "R3", "-", RowBox[{ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}]]], "+", FractionBox["4", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"0.`", "\[VeryThinSpace]", "+", "R4"}], ")"}], "2"], "+", SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"-", "R1"}], "-", "R3", "-", RowBox[{ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}]]], "+", FractionBox["8", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R3", "+", RowBox[{ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}]]], "+", FractionBox["4", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"0.`", "\[VeryThinSpace]", "-", "R4"}], ")"}], "2"], "+", SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R3", "+", RowBox[{ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}]]], "+", FractionBox["4", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"0.`", "\[VeryThinSpace]", "+", "R4"}], ")"}], "2"], "+", SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R3", "+", RowBox[{ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}]]], "+", FractionBox["8", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"2", " ", "R1"}], "+", RowBox[{"2", " ", "R3"}], "+", RowBox[{ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}]]], "+", FractionBox["8", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", "R1"}], "-", RowBox[{"2", " ", "R3"}], "-", RowBox[{"2", " ", RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{"4", " ", SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}]]], "+", FractionBox["24", SqrtBox[ RowBox[{ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Cos", "[", "w", "]"}], "2"]}], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}]]], "+", FractionBox["12", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", "R1"}], "-", RowBox[{"2", " ", "R3"}], "-", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}]]], "+", FractionBox["12", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"2", " ", "R1"}], "+", RowBox[{"2", " ", "R3"}], "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}]]], "+", FractionBox["2", SqrtBox[ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", "R1"}], "-", RowBox[{"2", " ", "R3"}], "-", RowBox[{"2", " ", RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{"4", " ", SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}]]]}], ")"}]}], "/", RowBox[{"(", RowBox[{ FractionBox["4.621499999999999`", SuperscriptBox["R1", "2"]], "+", RowBox[{"2.25`", " ", RowBox[{"(", RowBox[{ FractionBox["5.08`", SuperscriptBox["R2", "2"]], "+", FractionBox["1.2`", SuperscriptBox["R3", "2"]], "+", FractionBox["1.656`", SuperscriptBox["R4", "2"]]}], ")"}]}]}], ")"}]}], "/.", "\[VeryThinSpace]", RowBox[{"{", RowBox[{"R1", ",", "0.2`", ",", "0.28`"}], "}"}]}]], "Output", CellChangeTimes->{3.577341184997736*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"FindMinimum", "::", "lsbrak"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Unable to bracket a minimum along the direction \[NoBreak]\ \\!\\({0, 0, 0.19999999999999996`, 0, 0, 0}\\)\[NoBreak] from the point \ \[NoBreak]\\!\\({0.31620422775675705`, 2.180312896149525`*^12, 1.1`, 1.5`, \ 0.5`, 1.1`}\\)\[NoBreak]. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/FindMinimum\\\", ButtonNote -> \ \\\"FindMinimum::lsbrak\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.577341185075736*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.2`, 0.28`}\\)\[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.577341185106936*^9}], Cell[BoxData[ RowBox[{ RowBox[{"1.058354`", " ", RowBox[{"(", RowBox[{"R1", "+", "R3"}], ")"}]}], "/.", "\[VeryThinSpace]", RowBox[{"{", RowBox[{"R1", ",", "0.2`", ",", "0.28`"}], "}"}]}]], "Output", CellChangeTimes->{3.577341185106936*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"FindMinimum", "::", "lsbrak"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Unable to bracket a minimum along the direction \[NoBreak]\ \\!\\({0, 0, 0.19999999999999996`, 0, 0, 0}\\)\[NoBreak] from the point \ \[NoBreak]\\!\\({0.31620422775675705`, 2.180312896149525`*^12, 1.1`, 1.5`, \ 0.5`, 1.1`}\\)\[NoBreak]. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/FindMinimum\\\", ButtonNote -> \ \\\"FindMinimum::lsbrak\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.5773411851849365`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.2`, 0.28`}\\)\[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.577341185216136*^9}], Cell[BoxData[ RowBox[{ RowBox[{"0.529177`", " ", RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}]}], "/.", "\[VeryThinSpace]", RowBox[{"{", RowBox[{"R1", ",", "0.2`", ",", "0.28`"}], "}"}]}]], "Output", CellChangeTimes->{3.577341185216136*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"FindMinimum", "::", "lsbrak"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Unable to bracket a minimum along the direction \[NoBreak]\ \\!\\({0, 0, 0.19999999999999996`, 0, 0, 0}\\)\[NoBreak] from the point \ \[NoBreak]\\!\\({0.31620422775675705`, 2.180312896149525`*^12, 1.1`, 1.5`, \ 0.5`, 1.1`}\\)\[NoBreak]. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/FindMinimum\\\", ButtonNote -> \ \\\"FindMinimum::lsbrak\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.577341185294136*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.2`, 0.28`}\\)\[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.5773411853253365`*^9}], Cell[BoxData[ RowBox[{ RowBox[{"0.529177`", " ", RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}]}], "/.", "\[VeryThinSpace]", RowBox[{"{", RowBox[{"R1", ",", "0.2`", ",", "0.28`"}], "}"}]}]], "Output", CellChangeTimes->{3.5773411853253365`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"FindMinimum", "::", "lsbrak"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Unable to bracket a minimum along the direction \[NoBreak]\ \\!\\({0, 0, 0.19999999999999996`, 0, 0, 0}\\)\[NoBreak] from the point \ \[NoBreak]\\!\\({0.31620422775675705`, 2.180312896149525`*^12, 1.1`, 1.5`, \ 0.5`, 1.1`}\\)\[NoBreak]. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/FindMinimum\\\", ButtonNote -> \ \\\"FindMinimum::lsbrak\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.5773411854033365`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.2`, 0.28`}\\)\[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.577341185434537*^9}], Cell[BoxData[ RowBox[{ RowBox[{"114.59156`", " ", "w"}], "/.", "\[VeryThinSpace]", RowBox[{"{", RowBox[{"R1", ",", "0.2`", ",", "0.28`"}], "}"}]}]], "Output", CellChangeTimes->{3.577341185434537*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"FindMinimum", "::", "lsbrak"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Unable to bracket a minimum along the direction \[NoBreak]\ \\!\\({0, 0, 0.19999999999999996`, 0, 0, 0}\\)\[NoBreak] from the point \ \[NoBreak]\\!\\({0.31620422775675705`, 2.180312896149525`*^12, 1.1`, 1.5`, \ 0.5`, 1.1`}\\)\[NoBreak]. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/FindMinimum\\\", ButtonNote -> \ \\\"FindMinimum::lsbrak\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.577341185512537*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.2`, 0.28`}\\)\[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.5773411855437365`*^9}], Cell[BoxData[ RowBox[{ FractionBox[ RowBox[{"4", " ", SuperscriptBox[ RowBox[{"(", RowBox[{"R1", " ", RowBox[{"(", RowBox[{"R1", "+", RowBox[{"2", " ", "R4"}]}], ")"}]}], ")"}], RowBox[{"3", "/", "2"}]]}], SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R4"}], ")"}], "3"]], "/.", "\[VeryThinSpace]", RowBox[{"{", RowBox[{"R1", ",", "0.2`", ",", "0.28`"}], "}"}]}]], "Output", CellChangeTimes->{3.5773411855437365`*^9}] }, Open ]], Cell["\<\ (* Cvnf: sum of C1-nuclei forces *) Cvnf = 0.0; For[j = 2, j < 7, j++, Cvnf = Cvnf - \ (xn[[j]]-xn[[1]])*ch[[1]]*ch[[j]]/((xn[[1]]-xn[[j]])^2+(yn[[1]]-yn[[j]])^2+(\ zn[[1]]-zn[[j]])^2)^(3/2)]\ \>", "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.566799684588899*^9, 3.56679969039211*^9}, { 3.5668006247245502`*^9, 3.5668006962662764`*^9}}], Cell[CellGroupData[{ Cell["\<\ (* Cvef: sum of C1 - clouds forces *) Cvef = 0.0; For[j = 2, j < nc+1, j++, Cvef = Cvef - \ (xc[[j]]-xn[[1]])*oc[[j]]*ch[[1]]/((xc[[1]]-xn[[j]])^2+(yc[[1]]-yn[[j]])^2+(\ zc[[1]]-zn[[j]])^2)^(3/2)] Cforcetot = Cvnf+Cvef /. t[[2]] Cforcenn=Cvnf /. t[[2]] CForcene=Cvef /. t[[2]]\ \>", "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.566799684588899*^9, 3.56679969039211*^9}, { 3.5668006247245502`*^9, 3.566800704503091*^9}}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"FindMinimum", "::", "lsbrak"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Unable to bracket a minimum along the direction \[NoBreak]\ \\!\\({0, 0, 0.19999999999999996`, 0, 0, 0}\\)\[NoBreak] from the point \ \[NoBreak]\\!\\({0.31620422775675705`, 2.180312896149525`*^12, 1.1`, 1.5`, \ 0.5`, 1.1`}\\)\[NoBreak]. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/FindMinimum\\\", ButtonNote -> \ \\\"FindMinimum::lsbrak\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.577341185668537*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.2`, 0.28`}\\)\[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.5773411856997375`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.2`, 0.28`}\\)\[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.577341185730937*^9}], Cell[BoxData[ RowBox[{ RowBox[{"0.`", "\[VeryThinSpace]", "+", FractionBox[ RowBox[{"24", " ", RowBox[{"(", RowBox[{"R1", "+", "R3"}], ")"}]}], SuperscriptBox[ RowBox[{"(", SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"-", "R1"}], "-", "R3"}], ")"}], "2"], ")"}], RowBox[{"3", "/", "2"}]]], "-", FractionBox[ RowBox[{"24", " ", RowBox[{"(", RowBox[{ RowBox[{"2", " ", "R1"}], "+", RowBox[{"2", " ", "R3"}]}], ")"}]}], SuperscriptBox[ RowBox[{"(", SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", "R1"}], "-", RowBox[{"2", " ", "R3"}]}], ")"}], "2"], ")"}], RowBox[{"3", "/", "2"}]]], "-", FractionBox[ RowBox[{"24", " ", RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Cos", "[", "w", "]"}], "2"]}], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}], ")"}], RowBox[{"3", "/", "2"}]]], "+", FractionBox[ RowBox[{"12", " ", RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Cos", "[", "w", "]"}], "2"]}], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}], ")"}], RowBox[{"3", "/", "2"}]]], "+", FractionBox[ RowBox[{"24", " ", RowBox[{"(", RowBox[{ RowBox[{"2", " ", "R1"}], "+", RowBox[{"2", " ", "R3"}], "+", RowBox[{ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", "R1"}], "-", RowBox[{"2", " ", "R3"}], "-", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}], ")"}], RowBox[{"3", "/", "2"}]]], "-", FractionBox[ RowBox[{"12", " ", RowBox[{"(", RowBox[{ RowBox[{"2", " ", "R1"}], "+", RowBox[{"2", " ", "R3"}], "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", "R1"}], "-", RowBox[{"2", " ", "R3"}], "-", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}], ")"}], RowBox[{"3", "/", "2"}]]]}], "/.", "\[VeryThinSpace]", RowBox[{"{", RowBox[{"R1", ",", "0.2`", ",", "0.28`"}], "}"}]}]], "Output", CellChangeTimes->{3.577341185730937*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"FindMinimum", "::", "lsbrak"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Unable to bracket a minimum along the direction \[NoBreak]\ \\!\\({0, 0, 0.19999999999999996`, 0, 0, 0}\\)\[NoBreak] from the point \ \[NoBreak]\\!\\({0.31620422775675705`, 2.180312896149525`*^12, 1.1`, 1.5`, \ 0.5`, 1.1`}\\)\[NoBreak]. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/FindMinimum\\\", ButtonNote -> \ \\\"FindMinimum::lsbrak\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.577341185808937*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.2`, 0.28`}\\)\[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.5773411858401375`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.2`, 0.28`}\\)\[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.5773411858713374`*^9}], Cell[BoxData[ RowBox[{ RowBox[{"0.`", "\[VeryThinSpace]", "-", FractionBox[ RowBox[{"36", " ", RowBox[{"(", RowBox[{ RowBox[{"2", " ", "R1"}], "+", RowBox[{"2", " ", "R3"}]}], ")"}]}], SuperscriptBox[ RowBox[{"(", SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", "R1"}], "-", RowBox[{"2", " ", "R3"}]}], ")"}], "2"], ")"}], RowBox[{"3", "/", "2"}]]], "+", FractionBox[ RowBox[{"12", " ", RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Cos", "[", "w", "]"}], "2"]}], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}], ")"}], RowBox[{"3", "/", "2"}]]], "-", FractionBox[ RowBox[{"12", " ", RowBox[{"(", RowBox[{ RowBox[{"2", " ", "R1"}], "+", RowBox[{"2", " ", "R3"}], "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", "R1"}], "-", RowBox[{"2", " ", "R3"}], "-", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}], ")"}], RowBox[{"3", "/", "2"}]]]}], "/.", "\[VeryThinSpace]", RowBox[{"{", RowBox[{"R1", ",", "0.2`", ",", "0.28`"}], "}"}]}]], "Output", CellChangeTimes->{3.5773411858713374`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"FindMinimum", "::", "lsbrak"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Unable to bracket a minimum along the direction \[NoBreak]\ \\!\\({0, 0, 0.19999999999999996`, 0, 0, 0}\\)\[NoBreak] from the point \ \[NoBreak]\\!\\({0.31620422775675705`, 2.180312896149525`*^12, 1.1`, 1.5`, \ 0.5`, 1.1`}\\)\[NoBreak]. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/FindMinimum\\\", ButtonNote -> \ \\\"FindMinimum::lsbrak\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.5773411859493375`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.2`, 0.28`}\\)\[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.5773411859805374`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.2`, 0.28`}\\)\[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.577341186011738*^9}], Cell[BoxData[ RowBox[{ RowBox[{"0.`", "\[VeryThinSpace]", "+", FractionBox[ RowBox[{"24", " ", RowBox[{"(", RowBox[{"R1", "+", "R3"}], ")"}]}], SuperscriptBox[ RowBox[{"(", SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"-", "R1"}], "-", "R3"}], ")"}], "2"], ")"}], RowBox[{"3", "/", "2"}]]], "+", FractionBox[ RowBox[{"12", " ", RowBox[{"(", RowBox[{ RowBox[{"2", " ", "R1"}], "+", RowBox[{"2", " ", "R3"}]}], ")"}]}], SuperscriptBox[ RowBox[{"(", SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", "R1"}], "-", RowBox[{"2", " ", "R3"}]}], ")"}], "2"], ")"}], RowBox[{"3", "/", "2"}]]], "-", FractionBox[ RowBox[{"24", " ", RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Cos", "[", "w", "]"}], "2"]}], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}], ")"}], RowBox[{"3", "/", "2"}]]], "+", FractionBox[ RowBox[{"24", " ", RowBox[{"(", RowBox[{ RowBox[{"2", " ", "R1"}], "+", RowBox[{"2", " ", "R3"}], "+", RowBox[{ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", "R1"}], "-", RowBox[{"2", " ", "R3"}], "-", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}], ")"}], RowBox[{"3", "/", "2"}]]]}], "/.", "\[VeryThinSpace]", RowBox[{"{", RowBox[{"R1", ",", "0.2`", ",", "0.28`"}], "}"}]}]], "Output", CellChangeTimes->{3.577341186011738*^9}] }, Open ]], Cell["\<\ (* Hvnfx: sum of H3-nuclei forces *) Hvnfx = 0.0; For[j = 1, j < 8, j++, If[j != 3, Hvnfx = Hvnfx - \ (xn[[j]]-xn[[3]])*ch[[3]]*ch[[j]]/((xn[[3]]-xn[[j]])^2+(yn[[3]]-yn[[j]])^2+(\ zn[[3]]-zn[[j]])^2)^(3/2)]]\ \>", "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.566799684588899*^9, 3.56679969039211*^9}, { 3.5668006247245502`*^9, 3.5668007097759*^9}, {3.5668187369647627`*^9, 3.566818737651164*^9}}], Cell[CellGroupData[{ Cell["\<\ (* Hvefx: sum of H3 - clouds forces *) Hvefx = 0.0; For[j = 1, j < nc+1, j++, If[j != 3, Hvefx = Hvefx- \ (xc[[j]]-xn[[3]])*oc[[j]]*ch[[3]]/((xc[[j]]-xn[[3]])^2+(yc[[j]]-yn[[3]])^2+(\ zc[[j]]-zn[[3]])^2)^(3/2)]] Hforcetotx = Hvnfx+Hvefx/. t[[2]] Hforcennx=Hvnfx /. t[[2]] HForcenex=Hvefx /. t[[2]]\ \>", "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.566799684588899*^9, 3.56679969039211*^9}, { 3.5668006247245502`*^9, 3.566800722365122*^9}}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"FindMinimum", "::", "lsbrak"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Unable to bracket a minimum along the direction \[NoBreak]\ \\!\\({0, 0, 0.19999999999999996`, 0, 0, 0}\\)\[NoBreak] from the point \ \[NoBreak]\\!\\({0.31620422775675705`, 2.180312896149525`*^12, 1.1`, 1.5`, \ 0.5`, 1.1`}\\)\[NoBreak]. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/FindMinimum\\\", ButtonNote -> \ \\\"FindMinimum::lsbrak\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.577341186120938*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.2`, 0.28`}\\)\[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.5773411861365376`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.2`, 0.28`}\\)\[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.577341186167738*^9}], Cell[BoxData[ RowBox[{ RowBox[{"0.`", "\[VeryThinSpace]", "-", FractionBox[ RowBox[{ RowBox[{"2", " ", "R1"}], "+", RowBox[{"2", " ", "R3"}], "+", RowBox[{"2", " ", RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], SuperscriptBox[ RowBox[{"(", SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", "R1"}], "-", RowBox[{"2", " ", "R3"}], "-", RowBox[{"2", " ", RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], ")"}], RowBox[{"3", "/", "2"}]]], "-", FractionBox[ RowBox[{"4", " ", RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Cos", "[", "w", "]"}], "2"]}], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}], ")"}], RowBox[{"3", "/", "2"}]]], "-", FractionBox[ RowBox[{"6", " ", RowBox[{"(", RowBox[{ RowBox[{"2", " ", "R1"}], "+", RowBox[{"2", " ", "R3"}], "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", "R1"}], "-", RowBox[{"2", " ", "R3"}], "-", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}], ")"}], RowBox[{"3", "/", "2"}]]], "+", FractionBox[ RowBox[{"2", " ", RowBox[{"(", RowBox[{"R1", "+", "R3", "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R3", "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}], ")"}], RowBox[{"3", "/", "2"}]]], "+", FractionBox[ RowBox[{"R1", "+", "R3", "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"0.`", "\[VeryThinSpace]", "-", "R4"}], ")"}], "2"], "+", SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R3", "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}], ")"}], RowBox[{"3", "/", "2"}]]], "+", FractionBox[ RowBox[{"R1", "+", "R3", "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"0.`", "\[VeryThinSpace]", "+", "R4"}], ")"}], "2"], "+", SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R3", "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}], ")"}], RowBox[{"3", "/", "2"}]]], "+", FractionBox[ RowBox[{"2", " ", RowBox[{"(", RowBox[{ RowBox[{"2", " ", "R1"}], "+", RowBox[{"2", " ", "R3"}], "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"2", " ", "R1"}], "+", RowBox[{"2", " ", "R3"}], "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}], ")"}], RowBox[{"3", "/", "2"}]]], "-", FractionBox[ RowBox[{ RowBox[{"2", " ", "R1"}], "+", RowBox[{"2", " ", "R3"}], "+", RowBox[{"2", " ", RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", "R1"}], "-", RowBox[{"2", " ", "R3"}], "-", RowBox[{"2", " ", RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{"4", " ", SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}], ")"}], RowBox[{"3", "/", "2"}]]], "+", FractionBox[ RowBox[{"2", " ", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}]}], " ", RowBox[{"Cos", "[", "w", "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}]}], " ", RowBox[{"Cos", "[", "w", "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}]}], " ", RowBox[{"Sin", "[", "w", "]"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Sin", "[", "w", "]"}]}]}], ")"}], "2"]}], ")"}], RowBox[{"3", "/", "2"}]]], "+", FractionBox[ RowBox[{"2", " ", RowBox[{"(", RowBox[{ RowBox[{"2", " ", "R1"}], "+", RowBox[{"2", " ", "R3"}], "+", RowBox[{ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"2", " ", "R1"}], "+", RowBox[{"2", " ", "R3"}], "+", RowBox[{ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}]}], " ", RowBox[{"Sin", "[", "w", "]"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Sin", "[", "w", "]"}]}]}], ")"}], "2"]}], ")"}], RowBox[{"3", "/", "2"}]]], "+", FractionBox[ RowBox[{"2", " ", RowBox[{"(", RowBox[{ RowBox[{"2", " ", "R1"}], "+", RowBox[{"2", " ", "R3"}], "+", RowBox[{ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"2", " ", "R1"}], "+", RowBox[{"2", " ", "R3"}], "+", RowBox[{ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], " ", RowBox[{"Sin", "[", "w", "]"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Sin", "[", "w", "]"}]}]}], ")"}], "2"]}], ")"}], RowBox[{"3", "/", "2"}]]]}], "/.", "\[VeryThinSpace]", RowBox[{"{", RowBox[{"R1", ",", "0.2`", ",", "0.28`"}], "}"}]}]], "Output", CellChangeTimes->{3.577341186167738*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"FindMinimum", "::", "lsbrak"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Unable to bracket a minimum along the direction \[NoBreak]\ \\!\\({0, 0, 0.19999999999999996`, 0, 0, 0}\\)\[NoBreak] from the point \ \[NoBreak]\\!\\({0.31620422775675705`, 2.180312896149525`*^12, 1.1`, 1.5`, \ 0.5`, 1.1`}\\)\[NoBreak]. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/FindMinimum\\\", ButtonNote -> \ \\\"FindMinimum::lsbrak\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.577341186245738*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.2`, 0.28`}\\)\[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.577341186276938*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.2`, 0.28`}\\)\[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.577341186308138*^9}], Cell[BoxData[ RowBox[{ RowBox[{"0.`", "\[VeryThinSpace]", "-", FractionBox[ RowBox[{ RowBox[{"2", " ", "R1"}], "+", RowBox[{"2", " ", "R3"}], "+", RowBox[{"2", " ", RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], SuperscriptBox[ RowBox[{"(", SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", "R1"}], "-", RowBox[{"2", " ", "R3"}], "-", RowBox[{"2", " ", RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], ")"}], RowBox[{"3", "/", "2"}]]], "-", FractionBox[ RowBox[{"6", " ", RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Cos", "[", "w", "]"}], "2"]}], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}], ")"}], RowBox[{"3", "/", "2"}]]], "-", FractionBox[ RowBox[{"6", " ", RowBox[{"(", RowBox[{ RowBox[{"2", " ", "R1"}], "+", RowBox[{"2", " ", "R3"}], "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", "R1"}], "-", RowBox[{"2", " ", "R3"}], "-", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}], ")"}], RowBox[{"3", "/", "2"}]]], "-", FractionBox[ RowBox[{ RowBox[{"2", " ", "R1"}], "+", RowBox[{"2", " ", "R3"}], "+", RowBox[{"2", " ", RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", "R1"}], "-", RowBox[{"2", " ", "R3"}], "-", RowBox[{"2", " ", RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{"4", " ", SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}], ")"}], RowBox[{"3", "/", "2"}]]]}], "/.", "\[VeryThinSpace]", RowBox[{"{", RowBox[{"R1", ",", "0.2`", ",", "0.28`"}], "}"}]}]], "Output", CellChangeTimes->{3.577341186308138*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"FindMinimum", "::", "lsbrak"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Unable to bracket a minimum along the direction \[NoBreak]\ \\!\\({0, 0, 0.19999999999999996`, 0, 0, 0}\\)\[NoBreak] from the point \ \[NoBreak]\\!\\({0.31620422775675705`, 2.180312896149525`*^12, 1.1`, 1.5`, \ 0.5`, 1.1`}\\)\[NoBreak]. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/FindMinimum\\\", ButtonNote -> \ \\\"FindMinimum::lsbrak\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.5773411863861384`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.2`, 0.28`}\\)\[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.5773411864173384`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.2`, 0.28`}\\)\[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.5773411864485383`*^9}], Cell[BoxData[ RowBox[{ RowBox[{"0.`", "\[VeryThinSpace]", "+", FractionBox[ RowBox[{"2", " ", RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Cos", "[", "w", "]"}], "2"]}], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}], ")"}], RowBox[{"3", "/", "2"}]]], "+", FractionBox[ RowBox[{"2", " ", RowBox[{"(", RowBox[{"R1", "+", "R3", "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R3", "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}], ")"}], RowBox[{"3", "/", "2"}]]], "+", FractionBox[ RowBox[{"R1", "+", "R3", "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"0.`", "\[VeryThinSpace]", "-", "R4"}], ")"}], "2"], "+", SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R3", "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}], ")"}], RowBox[{"3", "/", "2"}]]], "+", FractionBox[ RowBox[{"R1", "+", "R3", "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"0.`", "\[VeryThinSpace]", "+", "R4"}], ")"}], "2"], "+", SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R3", "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}], ")"}], RowBox[{"3", "/", "2"}]]], "+", FractionBox[ RowBox[{"2", " ", RowBox[{"(", RowBox[{ RowBox[{"2", " ", "R1"}], "+", RowBox[{"2", " ", "R3"}], "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"2", " ", "R1"}], "+", RowBox[{"2", " ", "R3"}], "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}], ")"}], RowBox[{"3", "/", "2"}]]], "+", FractionBox[ RowBox[{"2", " ", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}]}], " ", RowBox[{"Cos", "[", "w", "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}]}], " ", RowBox[{"Cos", "[", "w", "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}]}], " ", RowBox[{"Sin", "[", "w", "]"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Sin", "[", "w", "]"}]}]}], ")"}], "2"]}], ")"}], RowBox[{"3", "/", "2"}]]], "+", FractionBox[ RowBox[{"2", " ", RowBox[{"(", RowBox[{ RowBox[{"2", " ", "R1"}], "+", RowBox[{"2", " ", "R3"}], "+", RowBox[{ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"2", " ", "R1"}], "+", RowBox[{"2", " ", "R3"}], "+", RowBox[{ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}]}], " ", RowBox[{"Sin", "[", "w", "]"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Sin", "[", "w", "]"}]}]}], ")"}], "2"]}], ")"}], RowBox[{"3", "/", "2"}]]], "+", FractionBox[ RowBox[{"2", " ", RowBox[{"(", RowBox[{ RowBox[{"2", " ", "R1"}], "+", RowBox[{"2", " ", "R3"}], "+", RowBox[{ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"2", " ", "R1"}], "+", RowBox[{"2", " ", "R3"}], "+", RowBox[{ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], " ", RowBox[{"Sin", "[", "w", "]"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Sin", "[", "w", "]"}]}]}], ")"}], "2"]}], ")"}], RowBox[{"3", "/", "2"}]]]}], "/.", "\[VeryThinSpace]", RowBox[{"{", RowBox[{"R1", ",", "0.2`", ",", "0.28`"}], "}"}]}]], "Output", CellChangeTimes->{3.5773411864485383`*^9}] }, Open ]], Cell["\<\ (* Hvnfy: sum of H3-nuclei forces *) Hvnfy = 0.0; For[j = 1, j < 8, j++, If[j != 3, Hvnfy = Hvnfy - \ (yn[[j]]-yn[[3]])*ch[[3]]*ch[[j]]/((xn[[3]]-xn[[j]])^2+(yn[[3]]-yn[[j]])^2+(\ zn[[3]]-zn[[j]])^2)^(3/2)]]\ \>", "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.566799684588899*^9, 3.56679969039211*^9}, { 3.5668006247245502`*^9, 3.566800727622331*^9}, {3.5668187604896045`*^9, 3.566818761051205*^9}}], Cell[CellGroupData[{ Cell["\<\ (* Hvefy: sum of H3 - clouds forces *) Hvefy = 0.0; For[j = 1, j < nc+1, j++, If[j != 3, Hvefy = Hvefy - \ (yc[[j]]-yn[[3]])*oc[[j]]*ch[[3]]/((xc[[j]]-xn[[3]])^2+(yc[[j]]-yn[[3]])^2+(\ zc[[j]]-zn[[3]])^2)^(3/2)]] Hforcetoty = Hvnfy+Hvefy /. t[[2]] Hforcenny=Hvnfy /. t[[2]] HForceney=Hvefy /. t[[2]]\ \>", "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.566799684588899*^9, 3.56679969039211*^9}, { 3.5668006247245502`*^9, 3.5668007324115396`*^9}}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"FindMinimum", "::", "lsbrak"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Unable to bracket a minimum along the direction \[NoBreak]\ \\!\\({0, 0, 0.19999999999999996`, 0, 0, 0}\\)\[NoBreak] from the point \ \[NoBreak]\\!\\({0.31620422775675705`, 2.180312896149525`*^12, 1.1`, 1.5`, \ 0.5`, 1.1`}\\)\[NoBreak]. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/FindMinimum\\\", ButtonNote -> \ \\\"FindMinimum::lsbrak\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.5773411865577383`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.2`, 0.28`}\\)\[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.5773411865889387`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.2`, 0.28`}\\)\[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.577341186604539*^9}], Cell[BoxData[ RowBox[{ RowBox[{"0.`", "\[VeryThinSpace]", "+", FractionBox[ RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Sin", "[", "w", "]"}]}], RowBox[{"4", " ", SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}], ")"}], RowBox[{"3", "/", "2"}]]}]], "+", FractionBox[ RowBox[{"4", " ", RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Sin", "[", "w", "]"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Cos", "[", "w", "]"}], "2"]}], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}], ")"}], RowBox[{"3", "/", "2"}]]], "+", FractionBox[ RowBox[{"6", " ", RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Sin", "[", "w", "]"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", "R1"}], "-", RowBox[{"2", " ", "R3"}], "-", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}], ")"}], RowBox[{"3", "/", "2"}]]], "-", FractionBox[ RowBox[{"2", " ", RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Sin", "[", "w", "]"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R3", "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}], ")"}], RowBox[{"3", "/", "2"}]]], "-", FractionBox[ RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Sin", "[", "w", "]"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"0.`", "\[VeryThinSpace]", "-", "R4"}], ")"}], "2"], "+", SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R3", "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}], ")"}], RowBox[{"3", "/", "2"}]]], "-", FractionBox[ RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Sin", "[", "w", "]"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"0.`", "\[VeryThinSpace]", "+", "R4"}], ")"}], "2"], "+", SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R3", "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}], ")"}], RowBox[{"3", "/", "2"}]]], "-", FractionBox[ RowBox[{"2", " ", RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Sin", "[", "w", "]"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"2", " ", "R1"}], "+", RowBox[{"2", " ", "R3"}], "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}], ")"}], RowBox[{"3", "/", "2"}]]], "+", FractionBox[ RowBox[{"2", " ", RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Sin", "[", "w", "]"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", "R1"}], "-", RowBox[{"2", " ", "R3"}], "-", RowBox[{"2", " ", RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{"4", " ", SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}], ")"}], RowBox[{"3", "/", "2"}]]], "+", FractionBox[ RowBox[{"2", " ", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}]}], " ", RowBox[{"Sin", "[", "w", "]"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Sin", "[", "w", "]"}]}]}], ")"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}]}], " ", RowBox[{"Cos", "[", "w", "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}]}], " ", RowBox[{"Sin", "[", "w", "]"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Sin", "[", "w", "]"}]}]}], ")"}], "2"]}], ")"}], RowBox[{"3", "/", "2"}]]], "+", FractionBox[ RowBox[{"2", " ", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}]}], " ", RowBox[{"Sin", "[", "w", "]"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Sin", "[", "w", "]"}]}]}], ")"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"2", " ", "R1"}], "+", RowBox[{"2", " ", "R3"}], "+", RowBox[{ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}]}], " ", RowBox[{"Sin", "[", "w", "]"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Sin", "[", "w", "]"}]}]}], ")"}], "2"]}], ")"}], RowBox[{"3", "/", "2"}]]], "+", FractionBox[ RowBox[{"2", " ", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], " ", RowBox[{"Sin", "[", "w", "]"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Sin", "[", "w", "]"}]}]}], ")"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"2", " ", "R1"}], "+", RowBox[{"2", " ", "R3"}], "+", RowBox[{ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], " ", RowBox[{"Sin", "[", "w", "]"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Sin", "[", "w", "]"}]}]}], ")"}], "2"]}], ")"}], RowBox[{"3", "/", "2"}]]]}], "/.", "\[VeryThinSpace]", RowBox[{"{", RowBox[{"R1", ",", "0.2`", ",", "0.28`"}], "}"}]}]], "Output", CellChangeTimes->{3.577341186620139*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"FindMinimum", "::", "lsbrak"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Unable to bracket a minimum along the direction \[NoBreak]\ \\!\\({0, 0, 0.19999999999999996`, 0, 0, 0}\\)\[NoBreak] from the point \ \[NoBreak]\\!\\({0.31620422775675705`, 2.180312896149525`*^12, 1.1`, 1.5`, \ 0.5`, 1.1`}\\)\[NoBreak]. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/FindMinimum\\\", ButtonNote -> \ \\\"FindMinimum::lsbrak\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.5773411866981387`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.2`, 0.28`}\\)\[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.577341186713739*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.2`, 0.28`}\\)\[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.577341186744939*^9}], Cell[BoxData[ RowBox[{ RowBox[{"0.`", "\[VeryThinSpace]", "+", FractionBox[ RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Sin", "[", "w", "]"}]}], RowBox[{"4", " ", SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}], ")"}], RowBox[{"3", "/", "2"}]]}]], "+", FractionBox[ RowBox[{"6", " ", RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Sin", "[", "w", "]"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Cos", "[", "w", "]"}], "2"]}], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}], ")"}], RowBox[{"3", "/", "2"}]]], "+", FractionBox[ RowBox[{"6", " ", RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Sin", "[", "w", "]"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", "R1"}], "-", RowBox[{"2", " ", "R3"}], "-", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}], ")"}], RowBox[{"3", "/", "2"}]]], "+", FractionBox[ RowBox[{"2", " ", RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Sin", "[", "w", "]"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", "2"}], " ", "R1"}], "-", RowBox[{"2", " ", "R3"}], "-", RowBox[{"2", " ", RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{"4", " ", SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}], ")"}], RowBox[{"3", "/", "2"}]]]}], "/.", "\[VeryThinSpace]", RowBox[{"{", RowBox[{"R1", ",", "0.2`", ",", "0.28`"}], "}"}]}]], "Output", CellChangeTimes->{3.577341186744939*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"FindMinimum", "::", "lsbrak"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Unable to bracket a minimum along the direction \[NoBreak]\ \\!\\({0, 0, 0.19999999999999996`, 0, 0, 0}\\)\[NoBreak] from the point \ \[NoBreak]\\!\\({0.31620422775675705`, 2.180312896149525`*^12, 1.1`, 1.5`, \ 0.5`, 1.1`}\\)\[NoBreak]. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/FindMinimum\\\", ButtonNote -> \ \\\"FindMinimum::lsbrak\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.577341186822939*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.2`, 0.28`}\\)\[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.5773411868541393`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.2`, 0.28`}\\)\[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.5773411868853393`*^9}], Cell[BoxData[ RowBox[{ RowBox[{"0.`", "\[VeryThinSpace]", "-", FractionBox[ RowBox[{"2", " ", RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Sin", "[", "w", "]"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Cos", "[", "w", "]"}], "2"]}], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}], ")"}], RowBox[{"3", "/", "2"}]]], "-", FractionBox[ RowBox[{"2", " ", RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Sin", "[", "w", "]"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R3", "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}], ")"}], RowBox[{"3", "/", "2"}]]], "-", FractionBox[ RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Sin", "[", "w", "]"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"0.`", "\[VeryThinSpace]", "-", "R4"}], ")"}], "2"], "+", SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R3", "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}], ")"}], RowBox[{"3", "/", "2"}]]], "-", FractionBox[ RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Sin", "[", "w", "]"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"0.`", "\[VeryThinSpace]", "+", "R4"}], ")"}], "2"], "+", SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R3", "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}], ")"}], RowBox[{"3", "/", "2"}]]], "-", FractionBox[ RowBox[{"2", " ", RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Sin", "[", "w", "]"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"2", " ", "R1"}], "+", RowBox[{"2", " ", "R3"}], "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}], ")"}], RowBox[{"3", "/", "2"}]]], "+", FractionBox[ RowBox[{"2", " ", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}]}], " ", RowBox[{"Sin", "[", "w", "]"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Sin", "[", "w", "]"}]}]}], ")"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}]}], " ", RowBox[{"Cos", "[", "w", "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}]}], " ", RowBox[{"Sin", "[", "w", "]"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Sin", "[", "w", "]"}]}]}], ")"}], "2"]}], ")"}], RowBox[{"3", "/", "2"}]]], "+", FractionBox[ RowBox[{"2", " ", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}]}], " ", RowBox[{"Sin", "[", "w", "]"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Sin", "[", "w", "]"}]}]}], ")"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"2", " ", "R1"}], "+", RowBox[{"2", " ", "R3"}], "+", RowBox[{ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}]}], " ", RowBox[{"Sin", "[", "w", "]"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Sin", "[", "w", "]"}]}]}], ")"}], "2"]}], ")"}], RowBox[{"3", "/", "2"}]]], "+", FractionBox[ RowBox[{"2", " ", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], " ", RowBox[{"Sin", "[", "w", "]"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Sin", "[", "w", "]"}]}]}], ")"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{"2", " ", "R1"}], "+", RowBox[{"2", " ", "R3"}], "+", RowBox[{ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", SuperscriptBox[ RowBox[{"(", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"R1", "+", "R2"}], ")"}], " ", RowBox[{"Sin", "[", "w", "]"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Sin", "[", "w", "]"}]}]}], ")"}], "2"]}], ")"}], RowBox[{"3", "/", "2"}]]]}], "/.", "\[VeryThinSpace]", RowBox[{"{", RowBox[{"R1", ",", "0.2`", ",", "0.28`"}], "}"}]}]], "Output", CellChangeTimes->{3.5773411868853393`*^9}] }, Open ]], Cell["\<\ (* Hvnfz: sum of H3-nuclei forces *) Hvnfz = 0.0; For[j = 7, j < 10, j++, If[j != 3, Hvnfz = Hvnfz - \ (zn[[j]]-zn[[3]])*ch[[3]]*ch[[j]]/((xn[[3]]-xn[[j]])^2+(yn[[3]]-yn[[j]])^2+(\ zn[[3]]-zn[[j]])^2)^(3/2)]]\ \>", "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.566799684588899*^9, 3.56679969039211*^9}, { 3.5668006247245502`*^9, 3.5668007404299536`*^9}, {3.5668020045937743`*^9, 3.566802043765443*^9}, {3.566802206068128*^9, 3.566802232900175*^9}, { 3.5668022898090754`*^9, 3.566802295066284*^9}, {3.56680237223962*^9, 3.56680239537446*^9}, {3.566818772064825*^9, 3.566818772720026*^9}}], Cell[CellGroupData[{ Cell["\<\ (* Hvefz: sum of H3 - clouds forces *) Hvefz = 0.0; For[j = 7, j < nc+1, j++, If[j != 3, Hvefz = Hvefz - \ (zc[[j]]-zn[[3]])*oc[[j]]*ch[[3]]/((xc[[j]]-xn[[3]])^2+(yc[[j]]-yn[[3]])^2+(\ zc[[j]]-zn[[3]])^2)^(3/2)]] Hforcetotz = Hvnfz+Hvefz /. t[[2]] Hforcennz=Hvnfz /. t[[2]] HForcenez=Hvefz /. t[[2]] \ \>", "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.566799684588899*^9, 3.56679969039211*^9}, { 3.5668006247245502`*^9, 3.5668007404299536`*^9}, {3.5668020045937743`*^9, 3.566802043765443*^9}, {3.566802206068128*^9, 3.566802232900175*^9}, { 3.5668022898090754`*^9, 3.566802295066284*^9}, {3.56680237223962*^9, 3.566802503638651*^9}}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"FindMinimum", "::", "lsbrak"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Unable to bracket a minimum along the direction \[NoBreak]\ \\!\\({0, 0, 0.19999999999999996`, 0, 0, 0}\\)\[NoBreak] from the point \ \[NoBreak]\\!\\({0.31620422775675705`, 2.180312896149525`*^12, 1.1`, 1.5`, \ 0.5`, 1.1`}\\)\[NoBreak]. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/FindMinimum\\\", ButtonNote -> \ \\\"FindMinimum::lsbrak\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.5773411869945393`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.2`, 0.28`}\\)\[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.5773411870257397`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.2`, 0.28`}\\)\[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.577341187056939*^9}], Cell[BoxData[ RowBox[{ RowBox[{"0.`", "\[VeryThinSpace]", "+", FractionBox[ RowBox[{"0.`", "\[VeryThinSpace]", "-", "R4"}], SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"0.`", "\[VeryThinSpace]", "-", "R4"}], ")"}], "2"], "+", SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R3", "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}], ")"}], RowBox[{"3", "/", "2"}]]], "+", FractionBox[ RowBox[{"0.`", "\[VeryThinSpace]", "+", "R4"}], SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"0.`", "\[VeryThinSpace]", "+", "R4"}], ")"}], "2"], "+", SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R3", "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}], ")"}], RowBox[{"3", "/", "2"}]]]}], "/.", "\[VeryThinSpace]", RowBox[{"{", RowBox[{"R1", ",", "0.2`", ",", "0.28`"}], "}"}]}]], "Output", CellChangeTimes->{3.577341187056939*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"FindMinimum", "::", "lsbrak"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Unable to bracket a minimum along the direction \[NoBreak]\ \\!\\({0, 0, 0.19999999999999996`, 0, 0, 0}\\)\[NoBreak] from the point \ \[NoBreak]\\!\\({0.31620422775675705`, 2.180312896149525`*^12, 1.1`, 1.5`, \ 0.5`, 1.1`}\\)\[NoBreak]. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/FindMinimum\\\", ButtonNote -> \ \\\"FindMinimum::lsbrak\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.5773411871349397`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.2`, 0.28`}\\)\[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.5773411871661396`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.2`, 0.28`}\\)\[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.5773411871973395`*^9}], Cell[BoxData[ RowBox[{"0.`", "/.", "\[VeryThinSpace]", RowBox[{"{", RowBox[{"R1", ",", "0.2`", ",", "0.28`"}], "}"}]}]], "Output", CellChangeTimes->{3.5773411871973395`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"FindMinimum", "::", "lsbrak"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Unable to bracket a minimum along the direction \[NoBreak]\ \\!\\({0, 0, 0.19999999999999996`, 0, 0, 0}\\)\[NoBreak] from the point \ \[NoBreak]\\!\\({0.31620422775675705`, 2.180312896149525`*^12, 1.1`, 1.5`, \ 0.5`, 1.1`}\\)\[NoBreak]. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/FindMinimum\\\", ButtonNote -> \ \\\"FindMinimum::lsbrak\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.57734118727534*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.2`, 0.28`}\\)\[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.57734118730654*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.2`, 0.28`}\\)\[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.5773411873221397`*^9}], Cell[BoxData[ RowBox[{ RowBox[{"0.`", "\[VeryThinSpace]", "+", FractionBox[ RowBox[{"0.`", "\[VeryThinSpace]", "-", "R4"}], SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"0.`", "\[VeryThinSpace]", "-", "R4"}], ")"}], "2"], "+", SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R3", "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}], ")"}], RowBox[{"3", "/", "2"}]]], "+", FractionBox[ RowBox[{"0.`", "\[VeryThinSpace]", "+", "R4"}], SuperscriptBox[ RowBox[{"(", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"0.`", "\[VeryThinSpace]", "+", "R4"}], ")"}], "2"], "+", SuperscriptBox[ RowBox[{"(", RowBox[{"R1", "+", "R3", "+", RowBox[{ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], " ", RowBox[{"Cos", "[", "w", "]"}]}]}], ")"}], "2"], "+", RowBox[{ SuperscriptBox[ RowBox[{"(", RowBox[{"p", "+", "R1", "+", "R2"}], ")"}], "2"], " ", SuperscriptBox[ RowBox[{"Sin", "[", "w", "]"}], "2"]}]}], ")"}], RowBox[{"3", "/", "2"}]]]}], "/.", "\[VeryThinSpace]", RowBox[{"{", RowBox[{"R1", ",", "0.2`", ",", "0.28`"}], "}"}]}]], "Output", CellChangeTimes->{3.57734118733774*^9}] }, Open ]], Cell[CellGroupData[{ Cell["\<\ (* Wvnfy: sum of w7-nuclei forces *) Wvnfy = 0.0; For[j = 1, j < 8, j++, Wvnfy = Wvnfy - \ (yn[[j]]-zc[[7]])*oc[[7]]*ch[[j]]/((xc[[7]]-xn[[j]])^2+(yc[[7]]-yn[[j]])^2+(\ zc[[7]]-zn[[j]])^2)^(3/2)]\ \>", "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.566799684588899*^9, 3.56679969039211*^9}, { 3.5668006247245502`*^9, 3.5668007404299536`*^9}, {3.5668020045937743`*^9, 3.566802043765443*^9}, {3.566802512827067*^9, 3.5668025183962765`*^9}, { 3.5668187803016396`*^9, 3.5668187808632402`*^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.57734118741574*^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.57734118746254*^9}] }, Open ]], Cell[CellGroupData[{ Cell["\<\ (* Wvefy: sum of W7 - clouds forces *) Wvefy = 0.0; For[j = 1, j < nc+1, j++, If[j != 7, Wvefy = Wvefy - \ (yc[[j]]-zc[[7]])*oc[[j]]*oc[[7]]/((xc[[j]]-xc[[7]])^2+(yc[[j]]-yc[[7]])^2+(\ zc[[j]]-zc[[7]])^2)^(3/2)]] Wforcetoty = Wvnfy+Wvefy /. t[[2]] Wforcenny=Wvnfy /. t[[2]] WForceney=Wvefy /. t[[2]]\ \>", "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.566799684588899*^9, 3.56679969039211*^9}, { 3.5668006247245502`*^9, 3.566800746217564*^9}}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"FindMinimum", "::", "lsbrak"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Unable to bracket a minimum along the direction \[NoBreak]\ \\!\\({0, 0, 0.19999999999999996`, 0, 0, 0}\\)\[NoBreak] from the point \ \[NoBreak]\\!\\({0.31620422775675705`, 2.180312896149525`*^12, 1.1`, 1.5`, \ 0.5`, 1.1`}\\)\[NoBreak]. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/FindMinimum\\\", ButtonNote -> \ \\\"FindMinimum::lsbrak\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.57734118757174*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.2`, 0.28`}\\)\[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.5773411875873404`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.2`, 0.28`}\\)\[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.5773411876185403`*^9}], Cell[BoxData[ RowBox[{"Indeterminate", "/.", "\[VeryThinSpace]", RowBox[{"{", RowBox[{"R1", ",", "0.2`", ",", "0.28`"}], "}"}]}]], "Output", CellChangeTimes->{3.5773411876185403`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"FindMinimum", "::", "lsbrak"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Unable to bracket a minimum along the direction \[NoBreak]\ \\!\\({0, 0, 0.19999999999999996`, 0, 0, 0}\\)\[NoBreak] from the point \ \[NoBreak]\\!\\({0.31620422775675705`, 2.180312896149525`*^12, 1.1`, 1.5`, \ 0.5`, 1.1`}\\)\[NoBreak]. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/FindMinimum\\\", ButtonNote -> \ \\\"FindMinimum::lsbrak\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.577341187696541*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.2`, 0.28`}\\)\[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.5773411877277403`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.2`, 0.28`}\\)\[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.5773411877589407`*^9}], Cell[BoxData[ RowBox[{"Indeterminate", "/.", "\[VeryThinSpace]", RowBox[{"{", RowBox[{"R1", ",", "0.2`", ",", "0.28`"}], "}"}]}]], "Output", CellChangeTimes->{3.5773411877589407`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"FindMinimum", "::", "lsbrak"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Unable to bracket a minimum along the direction \[NoBreak]\ \\!\\({0, 0, 0.19999999999999996`, 0, 0, 0}\\)\[NoBreak] from the point \ \[NoBreak]\\!\\({0.31620422775675705`, 2.180312896149525`*^12, 1.1`, 1.5`, \ 0.5`, 1.1`}\\)\[NoBreak]. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/FindMinimum\\\", ButtonNote -> \ \\\"FindMinimum::lsbrak\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.577341187836941*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.2`, 0.28`}\\)\[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.5773411878681407`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.2`, 0.28`}\\)\[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.5773411878993406`*^9}], Cell[BoxData[ RowBox[{"0.`", "/.", "\[VeryThinSpace]", RowBox[{"{", RowBox[{"R1", ",", "0.2`", ",", "0.28`"}], "}"}]}]], "Output", CellChangeTimes->{3.5773411878993406`*^9}] }, Open ]], Cell[CellGroupData[{ Cell["\<\ (* projection on xy-plane of molecule *) plot1=Graphics[{Circle[{xc[[1]],yc[[1]]},R1], \ Circle[{xc[[2]],yc[[2]]},R1],Circle[{xc[[3]],yc[[3]]},R2],Circle[{xc[[4]],yc[[\ 4]]},R2],Circle[{xc[[5]],yc[[5]]},R2],Circle[{xc[[6]],yc[[6]]},R2],Circle[{xc[\ [7]],yc[[7]]},R3],Circle[{xc[[8]],yc[[8]]},R4],Circle[{xc[[9]],yc[[9]]},R4],\ Disk[{xn[[1]],yn[[1]]},0.08], Disk[{xn[[2]],yn[[2]]},0.08], \ Disk[{xn[[3]],yn[[3]]},0.08], \ Disk[{xn[[4]],yn[[4]]},0.08],Disk[{xn[[5]],yn[[5]]},0.08],Disk[{xn[[6]],yn[[6]\ ]},0.08]} ] /. t[[2]]; Show[plot1,{AspectRatio \[Rule] Automatic,Axes -> True,GridLines -> \ Automatic, PlotRange \[Rule] {{-4,4},{-3,3}}, Frame -> True}]\ \>", "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.5668007580735846`*^9, 3.5668007676052017`*^9}, { 3.5668188223905134`*^9, 3.5668188450261526`*^9}}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"FindMinimum", "::", "lsbrak"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Unable to bracket a minimum along the direction \[NoBreak]\ \\!\\({0, 0, 0.19999999999999996`, 0, 0, 0}\\)\[NoBreak] from the point \ \[NoBreak]\\!\\({0.31620422775675705`, 2.180312896149525`*^12, 1.1`, 1.5`, \ 0.5`, 1.1`}\\)\[NoBreak]. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/FindMinimum\\\", ButtonNote -> \ \\\"FindMinimum::lsbrak\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.577341188008541*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.2`, 0.28`}\\)\[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.5773411880241413`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.2`, 0.28`}\\)\[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.5773411880553412`*^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.577341188086541*^9}], Cell[BoxData[ RowBox[{"Show", "[", RowBox[{ RowBox[{ GraphicsBox[{ CircleBox[{-$CellContext`R1 - $CellContext`R3, 0}, $CellContext`R1], CircleBox[{$CellContext`R1 + $CellContext`R3, 0}, $CellContext`R1], CircleBox[{-$CellContext`R1 - $CellContext`R3 - ($CellContext`R1 + \ $CellContext`R2) Cos[$CellContext`w], ($CellContext`R1 + $CellContext`R2) Sin[$CellContext`w]}, $CellContext`R2], CircleBox[{-$CellContext`R1 - $CellContext`R3 - ($CellContext`R1 + \ $CellContext`R2) Cos[$CellContext`w], -($CellContext`R1 + $CellContext`R2) Sin[$CellContext`w]}, $CellContext`R2], CircleBox[{$CellContext`R1 + $CellContext`R3 + ($CellContext`R1 + \ $CellContext`R2) Cos[$CellContext`w], -($CellContext`R1 + $CellContext`R2) Sin[$CellContext`w]}, $CellContext`R2], CircleBox[{$CellContext`R1 + $CellContext`R3 + ($CellContext`R1 + \ $CellContext`R2) Cos[$CellContext`w], ($CellContext`R1 + $CellContext`R2) Sin[$CellContext`w]}, $CellContext`R2], CircleBox[{0, 0}, $CellContext`R3], CircleBox[{0, 0}, $CellContext`R4], CircleBox[{0, 0}, $CellContext`R4], DiskBox[{-$CellContext`R1 - $CellContext`R3, 0}, 0.08], DiskBox[{$CellContext`R1 + $CellContext`R3, 0}, 0.08], DiskBox[{-$CellContext`R1 - $CellContext`R3 - ($CellContext`p + \ $CellContext`R1 + $CellContext`R2) Cos[$CellContext`w], ($CellContext`p + $CellContext`R1 + \ $CellContext`R2) Sin[$CellContext`w]}, 0.08], DiskBox[{-$CellContext`R1 - $CellContext`R3 - ($CellContext`p + \ $CellContext`R1 + $CellContext`R2) Cos[$CellContext`w], -($CellContext`p + $CellContext`R1 + \ $CellContext`R2) Sin[$CellContext`w]}, 0.08], DiskBox[{$CellContext`R1 + $CellContext`R3 + ($CellContext`p + \ $CellContext`R1 + $CellContext`R2) Cos[$CellContext`w], -($CellContext`p + $CellContext`R1 + \ $CellContext`R2) Sin[$CellContext`w]}, 0.08], DiskBox[{$CellContext`R1 + $CellContext`R3 + ($CellContext`p + \ $CellContext`R1 + $CellContext`R2) Cos[$CellContext`w], ($CellContext`p + $CellContext`R1 + \ $CellContext`R2) Sin[$CellContext`w]}, 0.08]}], "/.", "\[VeryThinSpace]", RowBox[{"{", RowBox[{"R1", ",", "0.2`", ",", "0.28`"}], "}"}]}], ",", 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.577341188086541*^9}] }, Open ]], Cell[CellGroupData[{ Cell["\<\ (* projection on xz-plane, perpendicular to molecular plane *) plot2=Graphics[{Circle[{xc[[1]],zc[[1]]},R1], \ Circle[{xc[[2]],zc[[2]]},R1],Circle[{xc[[3]],zc[[3]]},R2],Circle[{xc[[4]],zc[[\ 4]]},R2],Circle[{xc[[5]],zc[[5]]},R2],Circle[{xc[[6]],zc[[6]]},R2],Circle[{xc[\ [7]],zc[[7]]},R3],Circle[{xc[[8]],zc[[8]]},R4],Circle[{xc[[9]],zc[[9]]},R4],\ Disk[{xn[[1]],zn[[1]]},0.08], Disk[{xn[[2]],zn[[2]]},0.08], \ Disk[{xn[[3]],zn[[3]]},0.08], \ Disk[{xn[[4]],zn[[4]]},0.08],Disk[{xn[[5]],zn[[5]]},0.08],Disk[{xn[[6]],zn[[6]\ ]},0.08]} ] /. t[[2]]; Show[plot2,{AspectRatio \[Rule] Automatic,Axes -> True,GridLines -> \ Automatic, PlotRange \[Rule] {{-4,4},{-3,3}}, Frame -> True}]\ \>", "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.5668007580735846`*^9, 3.5668007725660105`*^9}, { 3.566818876475808*^9, 3.566818903120655*^9}}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"FindMinimum", "::", "lsbrak"}], "MessageName"], RowBox[{ ":", " "}], "\<\"Unable to bracket a minimum along the direction \[NoBreak]\ \\!\\({0, 0, 0.19999999999999996`, 0, 0, 0}\\)\[NoBreak] from the point \ \[NoBreak]\\!\\({0.31620422775675705`, 2.180312896149525`*^12, 1.1`, 1.5`, \ 0.5`, 1.1`}\\)\[NoBreak]. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \ ButtonStyle->\\\"Link\\\", ButtonFrame->None, \ ButtonData:>\\\"paclet:ref/FindMinimum\\\", ButtonNote -> \ \\\"FindMinimum::lsbrak\\\"]\\)\"\>"}]], "Message", "MSG", CellChangeTimes->{3.5773411881957417`*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.2`, 0.28`}\\)\[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.577341188211342*^9}], Cell[BoxData[ RowBox[{ StyleBox[ RowBox[{"ReplaceAll", "::", "reps"}], "MessageName"], RowBox[{ ":", " "}], "\<\"\[NoBreak]\\!\\({R1, 0.2`, 0.28`}\\)\[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.5773411882425413`*^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.5773411882737417`*^9}], Cell[BoxData[ RowBox[{"Show", "[", RowBox[{ RowBox[{ GraphicsBox[{ CircleBox[{-$CellContext`R1 - $CellContext`R3, 0}, $CellContext`R1], CircleBox[{$CellContext`R1 + $CellContext`R3, 0}, $CellContext`R1], CircleBox[{-$CellContext`R1 - $CellContext`R3 - ($CellContext`R1 + \ $CellContext`R2) Cos[$CellContext`w], 0}, $CellContext`R2], CircleBox[{-$CellContext`R1 - $CellContext`R3 - ($CellContext`R1 + \ $CellContext`R2) Cos[$CellContext`w], 0}, $CellContext`R2], CircleBox[{$CellContext`R1 + $CellContext`R3 + ($CellContext`R1 + \ $CellContext`R2) Cos[$CellContext`w], 0}, $CellContext`R2], CircleBox[{$CellContext`R1 + $CellContext`R3 + ($CellContext`R1 + \ $CellContext`R2) Cos[$CellContext`w], 0}, $CellContext`R2], CircleBox[{0, 0}, $CellContext`R3], CircleBox[{0, 0. + $CellContext`R4}, $CellContext`R4], CircleBox[{0, 0. - $CellContext`R4}, $CellContext`R4], DiskBox[{-$CellContext`R1 - $CellContext`R3, 0}, 0.08], DiskBox[{$CellContext`R1 + $CellContext`R3, 0}, 0.08], DiskBox[{-$CellContext`R1 - $CellContext`R3 - ($CellContext`p + \ $CellContext`R1 + $CellContext`R2) Cos[$CellContext`w], 0}, 0.08], DiskBox[{-$CellContext`R1 - $CellContext`R3 - ($CellContext`p + \ $CellContext`R1 + $CellContext`R2) Cos[$CellContext`w], 0}, 0.08], DiskBox[{$CellContext`R1 + $CellContext`R3 + ($CellContext`p + \ $CellContext`R1 + $CellContext`R2) Cos[$CellContext`w], 0}, 0.08], DiskBox[{$CellContext`R1 + $CellContext`R3 + ($CellContext`p + \ $CellContext`R1 + $CellContext`R2) Cos[$CellContext`w], 0}, 0.08]}], "/.", "\[VeryThinSpace]", RowBox[{"{", RowBox[{"R1", ",", "0.2`", ",", "0.28`"}], "}"}]}], ",", 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.5773411882737417`*^9}] }, Open ]], Cell["", "Input", PageWidth->WindowWidth, CellChangeTimes->{{3.5668007580735846`*^9, 3.5668007629251933`*^9}}] }, WindowToolbars->"EditBar", WindowSize->{1112, 964}, WindowMargins->{{Automatic, 238}, {Automatic, 0}}, PrintingCopies->1, PrintingPageRange->{Automatic, Automatic}, PrivateNotebookOptions->{"VersionedStylesheet"->{"Default.nb"[8.] -> False}}, CellLabelAutoDelete->True, Magnification->1.25, 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, 484, 10, 113, "Input"], Cell[1044, 32, 880, 16, 152, "Input"], Cell[1927, 50, 347, 11, 133, "Input"], Cell[2277, 63, 1862, 41, 556, "Input"], Cell[4142, 106, 371, 10, 113, "Input"], Cell[4516, 118, 417, 11, 113, "Input"], Cell[4936, 131, 753, 14, 113, "Input"], Cell[5692, 147, 1210, 21, 113, "Input"], Cell[CellGroupData[{ Cell[6927, 172, 1207, 18, 75, "Input"], Cell[8137, 192, 619, 12, 34, "Message"], Cell[8759, 206, 619, 12, 34, "Message"], Cell[9381, 220, 728, 19, 60, "Output"] }, Open ]], Cell[10124, 242, 977, 14, 36, "Input"], Cell[CellGroupData[{ Cell[11126, 260, 369, 13, 190, "Input"], Cell[11498, 275, 621, 12, 34, "Message"], Cell[12122, 289, 543, 11, 29, "Message"], Cell[12668, 302, 438, 15, 69, "Output"], Cell[13109, 319, 621, 12, 34, "Message"], Cell[13733, 333, 545, 11, 29, "Message"], Cell[14281, 346, 8655, 287, 728, "Output"], Cell[22939, 635, 621, 12, 34, "Message"], Cell[23563, 649, 545, 11, 29, "Message"], Cell[24111, 662, 3259, 105, 309, "Output"], Cell[27373, 769, 621, 12, 34, "Message"], Cell[27997, 783, 543, 11, 29, "Message"], Cell[28543, 796, 13335, 405, 1155, "Output"], Cell[41881, 1203, 619, 12, 34, "Message"], Cell[42503, 1217, 543, 11, 29, "Message"], Cell[43049, 1230, 260, 7, 36, "Output"], Cell[43312, 1239, 621, 12, 34, "Message"], Cell[43936, 1253, 543, 11, 29, "Message"], Cell[44482, 1266, 270, 7, 36, "Output"], Cell[44755, 1275, 619, 12, 34, "Message"], Cell[45377, 1289, 545, 11, 29, "Message"], Cell[45925, 1302, 262, 7, 36, "Output"], Cell[46190, 1311, 621, 12, 34, "Message"], Cell[46814, 1325, 543, 11, 29, "Message"], Cell[47360, 1338, 210, 5, 36, "Output"], Cell[47573, 1345, 619, 12, 34, "Message"], Cell[48195, 1359, 545, 11, 29, "Message"], Cell[48743, 1372, 498, 16, 63, "Output"] }, Open ]], Cell[49256, 1391, 365, 10, 113, "Input"], Cell[CellGroupData[{ Cell[49646, 1405, 448, 13, 171, "Input"], Cell[50097, 1420, 619, 12, 34, "Message"], Cell[50719, 1434, 545, 11, 29, "Message"], Cell[51267, 1447, 543, 11, 29, "Message"], Cell[51813, 1460, 4156, 133, 222, "Output"], Cell[55972, 1595, 619, 12, 34, "Message"], Cell[56594, 1609, 545, 11, 29, "Message"], Cell[57142, 1622, 545, 11, 29, "Message"], Cell[57690, 1635, 2255, 71, 116, "Output"], Cell[59948, 1708, 621, 12, 34, "Message"], Cell[60572, 1722, 545, 11, 29, "Message"], Cell[61120, 1735, 543, 11, 29, "Message"], Cell[61666, 1748, 2536, 82, 116, "Output"] }, Open ]], Cell[64217, 1833, 432, 12, 133, "Input"], Cell[CellGroupData[{ Cell[64674, 1849, 473, 14, 190, "Input"], Cell[65150, 1865, 619, 12, 34, "Message"], Cell[65772, 1879, 545, 11, 29, "Message"], Cell[66320, 1892, 543, 11, 29, "Message"], Cell[66866, 1905, 10980, 337, 566, "Output"], Cell[77849, 2244, 619, 12, 34, "Message"], Cell[78471, 2258, 543, 11, 29, "Message"], Cell[79017, 2271, 543, 11, 29, "Message"], Cell[79563, 2284, 3420, 105, 169, "Output"], Cell[82986, 2391, 621, 12, 34, "Message"], Cell[83610, 2405, 545, 11, 29, "Message"], Cell[84158, 2418, 545, 11, 29, "Message"], Cell[84706, 2431, 8441, 258, 460, "Output"] }, Open ]], Cell[93162, 2692, 434, 12, 133, "Input"], Cell[CellGroupData[{ Cell[93621, 2708, 477, 14, 190, "Input"], Cell[94101, 2724, 621, 12, 34, "Message"], Cell[94725, 2738, 545, 11, 29, "Message"], Cell[95273, 2751, 543, 11, 29, "Message"], Cell[95819, 2764, 10074, 310, 576, "Output"], Cell[105896, 3076, 621, 12, 34, "Message"], Cell[106520, 3090, 543, 11, 29, "Message"], Cell[107066, 3103, 543, 11, 29, "Message"], Cell[107612, 3116, 2966, 91, 172, "Output"], Cell[110581, 3209, 619, 12, 34, "Message"], Cell[111203, 3223, 545, 11, 29, "Message"], Cell[111751, 3236, 545, 11, 29, "Message"], Cell[112299, 3249, 7989, 245, 468, "Output"] }, Open ]], Cell[120303, 3497, 627, 14, 133, "Input"], Cell[CellGroupData[{ Cell[120955, 3515, 681, 17, 190, "Input"], Cell[121639, 3534, 621, 12, 34, "Message"], Cell[122263, 3548, 545, 11, 29, "Message"], Cell[122811, 3561, 543, 11, 29, "Message"], Cell[123357, 3574, 1715, 49, 116, "Output"], Cell[125075, 3625, 621, 12, 34, "Message"], Cell[125699, 3639, 545, 11, 29, "Message"], Cell[126247, 3652, 545, 11, 29, "Message"], Cell[126795, 3665, 182, 4, 36, "Output"], Cell[126980, 3671, 618, 12, 34, "Message"], Cell[127601, 3685, 542, 11, 29, "Message"], Cell[128146, 3698, 545, 11, 29, "Message"], Cell[128694, 3711, 1714, 49, 116, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[130445, 3765, 521, 12, 113, "Input"], Cell[130969, 3779, 444, 10, 47, "Message"], Cell[131416, 3791, 474, 10, 29, "Message"] }, Open ]], Cell[CellGroupData[{ Cell[131927, 3806, 475, 14, 190, "Input"], Cell[132405, 3822, 618, 12, 34, "Message"], Cell[133026, 3836, 545, 11, 29, "Message"], Cell[133574, 3849, 545, 11, 29, "Message"], Cell[134122, 3862, 192, 4, 36, "Output"], Cell[134317, 3868, 619, 12, 34, "Message"], Cell[134939, 3882, 545, 11, 29, "Message"], Cell[135487, 3895, 545, 11, 29, "Message"], Cell[136035, 3908, 192, 4, 36, "Output"], Cell[136230, 3914, 619, 12, 34, "Message"], Cell[136852, 3928, 545, 11, 29, "Message"], Cell[137400, 3941, 545, 11, 29, "Message"], Cell[137948, 3954, 182, 4, 36, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[138167, 3963, 835, 16, 229, "Input"], Cell[139005, 3981, 619, 12, 34, "Message"], Cell[139627, 3995, 545, 11, 29, "Message"], Cell[140175, 4008, 545, 11, 29, "Message"], Cell[140723, 4021, 442, 10, 29, "Message"], Cell[141168, 4033, 2855, 56, 178, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[144060, 4094, 853, 16, 229, "Input"], Cell[144916, 4112, 621, 12, 34, "Message"], Cell[145540, 4126, 543, 11, 29, "Message"], Cell[146086, 4139, 545, 11, 29, "Message"], Cell[146634, 4152, 444, 10, 29, "Message"], Cell[147081, 4164, 2321, 46, 178, "Output"] }, Open ]], Cell[149417, 4213, 112, 2, 36, "Input"] } ] *) (* End of internal cache information *)