(* 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[ 164048, 4217] NotebookOptionsPosition[ 160233, 4103] NotebookOutlinePosition[ 160579, 4118] CellTagsIndexPosition[ 160536, 4115] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[BoxData[ StyleBox[ RowBox[{"(*", " ", RowBox[{ RowBox[{"NH3", " ", "with", " ", RowBox[{"N", "[", RowBox[{ RowBox[{"1", "s2", " ", "3", "sp"}], ",", " ", RowBox[{"1", " ", "sp", " ", "LP"}]}], "]"}]}], ",", RowBox[{ RowBox[{"~", "no"}], " ", "overlap"}], ",", " ", RowBox[{ "parametrized", " ", "ES", " ", "02.01", ".2013", "\[IndentingNewLine]", RowBox[{"Variables", ":", "\[IndentingNewLine]", " ", RowBox[{ RowBox[{"P", " ", "N"}], "-", RowBox[{"core", " ", "radius", "\[IndentingNewLine]", RowBox[{"Q", " ", "N"}]}], "-", RowBox[{ "H", " ", "cloud", " ", "radius", "\[IndentingNewLine]", "R", " ", "lone", " ", "pair", " ", "radius", "\[IndentingNewLine]", "S", " ", "distance", " ", "proton", " ", "to", " ", "center", " ", "of", " ", "NH", " ", "cloud", "\[IndentingNewLine]", RowBox[{"\[Alpha]", " ", "angle", " ", "lonePair"}]}], "-", "N", "-", RowBox[{ "H", "\[IndentingNewLine]", "\[Beta]", " ", "angle", " ", "H"}], "-", "N", "-", "H"}]}]}], ",", " ", RowBox[{ RowBox[{"not", " ", "a", " ", RowBox[{"variable", ":", " ", "cos\[Beta]"}]}], " ", "=", " ", RowBox[{"1.0", "-", RowBox[{"1.5", "*", RowBox[{ RowBox[{"Sin", "[", "\[Alpha]", "]"}], "^", "2"}]}]}]}]}], " ", "*)"}], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]]], "Input", CellChangeTimes->{{3.5661428484442425`*^9, 3.5661428783026953`*^9}, { 3.5661430328053665`*^9, 3.5661430369081736`*^9}, {3.5661430998854847`*^9, 3.5661431383537183`*^9}, {3.566143180801393*^9, 3.5661432248645825`*^9}, { 3.566143288683233*^9, 3.5661432893462706`*^9}, 3.566143344240477*^9, 3.566143438043442*^9, 3.566143491551536*^9}], Cell[BoxData[{ RowBox[{ RowBox[{"Clear", "[", "wc", "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Z", "=", "7.0"}], ";", " ", RowBox[{"w3", "=", RowBox[{"Sqrt", "[", "3", "]"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", "core", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Tc", "=", RowBox[{"2.25", "*", RowBox[{"k1", "/", RowBox[{"P", "^", "2"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vnec", "=", RowBox[{ RowBox[{"-", "3.0"}], "*", RowBox[{"Z", "/", "P"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Veec", "=", RowBox[{"3.0", "*", RowBox[{"s1", "/", "P"}]}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"XH", " ", "bond"}], "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Tb", "=", RowBox[{"2.25", "*", RowBox[{"k2", "/", RowBox[{"Q", "^", "2"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vneb", "=", RowBox[{ RowBox[{"-", "2.0"}], "*", RowBox[{ RowBox[{"(", RowBox[{"1.5", "-", RowBox[{"0.5", "*", RowBox[{ RowBox[{"(", RowBox[{"S", "/", "Q"}], ")"}], "^", "2"}]}]}], ")"}], "/", "Q"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Veeb", "=", RowBox[{"3.0", "*", RowBox[{"s2", "/", "Q"}]}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"lone", " ", "pair"}], "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Tp", "=", RowBox[{"2.25", "*", RowBox[{"k3", "/", RowBox[{"R", "^", "2"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Veep", "=", RowBox[{"3.0", "*", RowBox[{"s3", "/", "R"}]}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"core", "-", "bond"}], "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vnecb", "=", RowBox[{ RowBox[{ RowBox[{"-", "2.0"}], "*", RowBox[{"Z", "/", RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}]}]}], "-", RowBox[{"2", "/", RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Veecb", "=", RowBox[{"4.0", "/", RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Vnncb", "=", RowBox[{"Z", "/", RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}]}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"core", "-", RowBox[{"lone", " ", "pair"}]}], "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vnecp", "=", RowBox[{ RowBox[{"-", "2.0"}], "*", RowBox[{"Z", "/", RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Veecp", "=", RowBox[{"4.0", "/", RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}]}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"bond", "-", "bond"}], "*)"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ RowBox[{"d1", " ", "w", RowBox[{"(", "NH", ")"}]}], "-", RowBox[{"w", RowBox[{"(", "NH", ")"}]}]}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"d1", "=", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}], "*", RowBox[{"Sin", "[", "\[Alpha]", "]"}], "*", "w3"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ RowBox[{"d3", " ", "H"}], "-", "H"}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"d3", "=", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}], "*", RowBox[{"Sin", "[", "\[Alpha]", "]"}], "*", "w3"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ RowBox[{ RowBox[{"d5", " ", "w", RowBox[{"(", "NH", ")"}]}], "-", "H"}], ",", " ", RowBox[{"\[Beta]", " ", "is", " ", "angle", " ", "HNH"}]}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"cos\[Beta]", "=", RowBox[{"1.0", "-", RowBox[{"1.5", "*", RowBox[{ RowBox[{"Sin", "[", "\[Alpha]", "]"}], "^", "2"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"d5", "=", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}], "^", "2"}], "+", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}], "^", "2"}], "-", RowBox[{"2", "*", RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}], "*", RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}], "*", "cos\[Beta]"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Veebb", "=", RowBox[{"4.0", "/", "d1"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vnebb", "=", RowBox[{ RowBox[{"-", "4.0"}], "/", "d5"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Vnnbb", "=", RowBox[{"1.0", "/", "d3"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"bond", "-", RowBox[{"lone", " ", "pair"}]}], "*)"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ RowBox[{"d2", " ", "w", RowBox[{"(", "NH", ")"}]}], "-", RowBox[{"w", RowBox[{"(", "LP", ")"}]}]}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"d2", "=", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}], "^", "2"}], "+", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}], "^", "2"}], "-", RowBox[{"2", "*", RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}], "*", RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}], "*", RowBox[{"Cos", "[", "\[Alpha]", "]"}]}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Veebp", " ", "=", " ", RowBox[{"4", "/", "d2"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ RowBox[{"d4", " ", "w", RowBox[{"(", "LP", ")"}]}], "-", "H"}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"d4", "=", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}], "^", "2"}], "+", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}], "^", "2"}], "-", RowBox[{"2", "*", RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}], "*", RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}], "*", RowBox[{"Cos", "[", "\[Alpha]", "]"}]}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Vnebp", " ", "=", RowBox[{ RowBox[{"-", "2.0"}], "/", "d4"}]}], " ", ";"}], StyleBox["\[IndentingNewLine]", FontFamily->"Arial Narrow"], StyleBox["\[IndentingNewLine]", FontFamily->"Arial Narrow"], StyleBox[ RowBox[{"(*", RowBox[{ RowBox[{"lone", " ", "pair"}], "-", RowBox[{"lone", " ", "pair"}]}], "*)"}], "Text", FontFamily->"Arial Narrow"], "\[IndentingNewLine]", StyleBox[ RowBox[{"(*", " ", RowBox[{ RowBox[{"none", " ", "for", " ", "NH3"}], ",", " ", RowBox[{"but", " ", "for", " ", "H2O", " ", "and", " ", "HF"}]}], " ", "*)"}], "Text", FontFamily->"Arial Narrow"], "\[IndentingNewLine]", StyleBox[ RowBox[{"(*", " ", RowBox[{ RowBox[{ "H2O", " ", "has", " ", "two", " ", "independent", " ", "angles", " ", "as", " ", RowBox[{"variables", ":", " ", RowBox[{"H", "-", "O", "-", RowBox[{"H", " ", "and", " ", "LP"}], "-", "O", "-", "LP"}]}]}], ",", " ", RowBox[{ RowBox[{"or", " ", "LP"}], "-", "O", "-", "H"}]}], " ", "*)"}], "Text", FontFamily->"Arial Narrow"], "\[IndentingNewLine]", StyleBox[ RowBox[{"(*", " ", RowBox[{ RowBox[{ "HF", " ", "is", " ", "identical", " ", "to", " ", "NH3", " ", "with", " ", RowBox[{"LP", "/", "bond"}], " ", "reversed"}], ",", " ", RowBox[{"one", " ", "angle"}]}], " ", "*)"}], "Text", FontFamily->"Arial Narrow"], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"(*", RowBox[{"sum", " ", "terms"}], "*)"}]}]}], "Input", CellChangeTimes->{{3.5661428484442425`*^9, 3.5661428783026953`*^9}, { 3.5661430328053665`*^9, 3.5661430369081736`*^9}, {3.5661430998854847`*^9, 3.5661431383537183`*^9}, {3.566143180801393*^9, 3.5661432248645825`*^9}, { 3.566143288683233*^9, 3.5661432893462706`*^9}, 3.566143344240477*^9, 3.566143438043442*^9}], Cell[BoxData[{ RowBox[{ RowBox[{"T", "=", RowBox[{"Tc", "+", " ", RowBox[{"3", "*", "Tb"}], "+", " ", "Tp"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vne", "=", RowBox[{"Vnec", "+", RowBox[{"3", "*", "Vneb"}], "+", " ", RowBox[{"3", "*", "Vnecb"}], "+", "Vnecp", "+", RowBox[{"3", "*", "Vnebb"}], "+", RowBox[{"3", "*", "Vnebp"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vee", "=", RowBox[{"Veec", "+", RowBox[{"3", "*", "Veeb"}], "+", " ", "Veep", "+", RowBox[{"3", "*", "Veecb"}], "+", "Veecp", "+", RowBox[{"3", "*", "Veebb"}], "+", RowBox[{"3", "*", "Veebp"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vnn", "=", " ", RowBox[{ RowBox[{"3", "*", "Vnncb"}], "+", " ", RowBox[{"3", "*", "Vnnbb"}]}]}], ";"}], "\[IndentingNewLine]"}], "Input", CellChangeTimes->{{3.5661428484442425`*^9, 3.5661428783026953`*^9}, { 3.5661430328053665`*^9, 3.5661430369081736`*^9}, {3.5661430998854847`*^9, 3.5661431383537183`*^9}, {3.566143180801393*^9, 3.5661432248645825`*^9}, { 3.566143288683233*^9, 3.5661432893462706`*^9}, 3.566143344240477*^9, { 3.5661435118159714`*^9, 3.5661435187891836`*^9}, {3.5661470921646585`*^9, 3.5661471085134873`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ RowBox[{"Parameters", ":", " ", "G3"}], "//", RowBox[{ RowBox[{ RowBox[{ RowBox[{"6", "-", RowBox[{"311", "G", RowBox[{"(", RowBox[{"d", ",", "p"}], ")"}], " ", "can", " ", "be", " ", "read", " ", "for", " ", "LiH"}]}], " ", "..."}], " ", "HF", "\[IndentingNewLine]", "in", " ", "table", " ", "from", " ", RowBox[{"www", ".", "kimball"}]}], "-", RowBox[{ RowBox[{ RowBox[{"model", ".", "org"}], "/", "Kimball"}], "/", RowBox[{"params", ".", "xls"}]}]}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"par", "=", RowBox[{"{", RowBox[{ RowBox[{"k1", "\[Rule]", "0.95894892"}], ",", RowBox[{"k2", "\[Rule]", "1.34927913"}], ",", RowBox[{"k3", "\[Rule]", "1.23272312"}], ",", RowBox[{"s1", "\[Rule]", "0.42320792"}], ",", RowBox[{"s2", "\[Rule]", "0.29390809"}], ",", RowBox[{"s3", "\[Rule]", "0.39204779"}]}], "}"}]}], ";"}]}]], "Input", CellChangeTimes->{{3.5661428484442425`*^9, 3.5661428783026953`*^9}, { 3.5661430328053665`*^9, 3.5661430369081736`*^9}, {3.5661430998854847`*^9, 3.5661431383537183`*^9}, {3.566143180801393*^9, 3.5661432248645825`*^9}, { 3.566143288683233*^9, 3.5661432893462706`*^9}, {3.566143344240477*^9, 3.566143371150524*^9}}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"func", "=", RowBox[{ RowBox[{"T", "+", "Vne", "+", "Vee", "+", "Vnn"}], "/.", "par"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"t", "=", RowBox[{"FindMinimum", "[", RowBox[{"func", ",", RowBox[{"{", RowBox[{"P", ",", "0.2125"}], "}"}], ",", RowBox[{"{", RowBox[{"Q", ",", "1.12"}], "}"}], ",", RowBox[{"{", RowBox[{"R", ",", " ", "1.355"}], "}"}], ",", RowBox[{"{", RowBox[{"S", ",", "0.58"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"\[Alpha]", ",", "1.97"}], "}"}], ",", RowBox[{"MaxIterations", "\[Rule]", "500"}]}], "]"}]}], "\[IndentingNewLine]"}], "Input", CellChangeTimes->{{3.5661428484442425`*^9, 3.5661428783026953`*^9}, { 3.5661430328053665`*^9, 3.5661430369081736`*^9}, {3.5661430998854847`*^9, 3.5661431383537183`*^9}, {3.566143180801393*^9, 3.5661432248645825`*^9}, { 3.566143288683233*^9, 3.5661432893462706`*^9}, {3.566143344240477*^9, 3.5661433819769435`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"-", "56.50693567896597`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.21255884296750083`"}], ",", RowBox[{"Q", "\[Rule]", "1.1215796462307126`"}], ",", RowBox[{"R", "\[Rule]", "1.3545219025935082`"}], ",", RowBox[{"S", "\[Rule]", "0.5794206077585493`"}], ",", RowBox[{"\[Alpha]", "\[Rule]", "1.9706412301331597`"}]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{{3.5661466976555657`*^9, 3.5661467183100023`*^9}, 3.566146844592224*^9, 3.566146877617482*^9, 3.566146938785189*^9, 3.566146993166885*^9, 3.5661477546897445`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"u", "=", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"Vne", "/.", "u"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vee", "/.", "u"}], "/.", "par"}], "\[IndentingNewLine]", RowBox[{"Vnn", "/.", "u"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"Vne", "+", "Vee", "+", "Vnn"}], ")"}], "/", "T"}], "/.", "u"}], "/.", "par"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}], "*", "0.529177"}], "/.", "u"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{"angle", " ", "HNH"}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"180", "/", "\[Pi]"}], "*", RowBox[{"ArcCos", "[", "cos\[Beta]", "]"}]}], "/.", "u"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ RowBox[{"angle", " ", "LP"}], "-", "N", "-", "H"}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"180", "/", "\[Pi]"}], "*", "\[Alpha]"}], " ", "/.", " ", "u"}], "\[IndentingNewLine]", RowBox[{"d1", " ", "/.", "u"}], "\[IndentingNewLine]", RowBox[{"d3", " ", "/.", "u"}], "\[IndentingNewLine]", RowBox[{"d5", " ", "/.", "u"}], "\[IndentingNewLine]", RowBox[{"d2", " ", "/.", "u"}], "\[IndentingNewLine]", RowBox[{"d4", " ", "/.", "u"}], "\[IndentingNewLine]"}], "Input", CellChangeTimes->{{3.5661428484442425`*^9, 3.5661428783026953`*^9}, { 3.5661430328053665`*^9, 3.5661430369081736`*^9}, {3.5661430998854847`*^9, 3.5661431383537183`*^9}, {3.566143180801393*^9, 3.5661432248645825`*^9}, { 3.566143288683233*^9, 3.5661432893462706`*^9}, {3.566143344240477*^9, 3.5661433961729684`*^9}}], Cell[BoxData[ RowBox[{"-", "156.31115252225578`"}]], "Output", CellChangeTimes->{{3.5661466976555657`*^9, 3.566146718325602*^9}, 3.5661468446078243`*^9, 3.566146877617482*^9, 3.5661469388007894`*^9, 3.5661469931824846`*^9, 3.566147754700745*^9}], Cell[BoxData["31.340309455040384`"], "Output", CellChangeTimes->{{3.5661466976555657`*^9, 3.566146718325602*^9}, 3.5661468446078243`*^9, 3.566146877617482*^9, 3.5661469388007894`*^9, 3.5661469931824846`*^9, 3.566147754701745*^9}], Cell[BoxData["11.956971691023773`"], "Output", CellChangeTimes->{{3.5661466976555657`*^9, 3.566146718325602*^9}, 3.5661468446078243`*^9, 3.566146877617482*^9, 3.5661469388007894`*^9, 3.5661469931824846`*^9, 3.5661477547027454`*^9}], Cell[BoxData[ RowBox[{"-", "1.999999999676859`"}]], "Output", CellChangeTimes->{{3.5661466976555657`*^9, 3.566146718325602*^9}, 3.5661468446078243`*^9, 3.566146877617482*^9, 3.5661469388007894`*^9, 3.5661469931824846`*^9, 3.5661477547027454`*^9}], Cell[BoxData["1.0126114622502889`"], "Output", CellChangeTimes->{{3.5661466976555657`*^9, 3.566146718325602*^9}, 3.5661468446078243`*^9, 3.566146877617482*^9, 3.5661469388007894`*^9, 3.5661469931824846`*^9, 3.566147754703745*^9}], Cell[BoxData["105.8248111138388`"], "Output", CellChangeTimes->{{3.5661466976555657`*^9, 3.566146718325602*^9}, 3.5661468446078243`*^9, 3.566146877617482*^9, 3.5661469388007894`*^9, 3.5661469931824846`*^9, 3.5661477547047453`*^9}], Cell[BoxData["112.90942542109883`"], "Output", CellChangeTimes->{{3.5661466976555657`*^9, 3.566146718325602*^9}, 3.5661468446078243`*^9, 3.566146877617482*^9, 3.5661469388007894`*^9, 3.5661469931824846`*^9, 3.5661477547057457`*^9}], Cell[BoxData["2.128523276784143`"], "Output", CellChangeTimes->{{3.5661466976555657`*^9, 3.566146718325602*^9}, 3.5661468446078243`*^9, 3.566146877617482*^9, 3.5661469388007894`*^9, 3.5661469931824846`*^9, 3.5661477547067456`*^9}], Cell[BoxData["3.0529477354501084`"], "Output", CellChangeTimes->{{3.5661466976555657`*^9, 3.566146718325602*^9}, 3.5661468446078243`*^9, 3.566146877617482*^9, 3.5661469388007894`*^9, 3.5661469931824846`*^9, 3.5661477547067456`*^9}], Cell[BoxData["2.6141917600677806`"], "Output", CellChangeTimes->{{3.5661466976555657`*^9, 3.566146718325602*^9}, 3.5661468446078243`*^9, 3.566146877617482*^9, 3.5661469388007894`*^9, 3.5661469931824846`*^9, 3.5661477547077456`*^9}], Cell[BoxData["2.421443002175422`"], "Output", CellChangeTimes->{{3.5661466976555657`*^9, 3.566146718325602*^9}, 3.5661468446078243`*^9, 3.566146877617482*^9, 3.5661469388007894`*^9, 3.5661469931824846`*^9, 3.566147754708746*^9}], Cell[BoxData["2.9072482741647896`"], "Output", CellChangeTimes->{{3.5661466976555657`*^9, 3.566146718325602*^9}, 3.5661468446078243`*^9, 3.566146877617482*^9, 3.5661469388007894`*^9, 3.5661469931824846`*^9, 3.566147754708746*^9}] }, Open ]], Cell[BoxData[{ RowBox[{ RowBox[{"wd", " ", "=", " ", RowBox[{ RowBox[{"Sin", "[", "\[Alpha]", "]"}], " ", "/.", " ", "u"}]}], ";", " ", RowBox[{"wc", "=", RowBox[{ RowBox[{"Cos", "[", "\[Alpha]", "]"}], " ", "/.", " ", "u"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"(*", " ", RowBox[{ RowBox[{"Projection", " ", "onto", " ", "HHH"}], "-", "plane"}], " ", "*)"}]}]}], "Input", CellChangeTimes->{{3.5661428484442425`*^9, 3.5661428783026953`*^9}, { 3.5661430328053665`*^9, 3.5661430369081736`*^9}, {3.5661430998854847`*^9, 3.5661431383537183`*^9}, {3.566143180801393*^9, 3.5661432248645825`*^9}, { 3.566143288683233*^9, 3.5661432893462706`*^9}, {3.566143344240477*^9, 3.566143410119393*^9}}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"plot1", "=", RowBox[{ RowBox[{"Graphics", "[", RowBox[{"{", RowBox[{ RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", "P"}], "]"}], ",", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", "0.1"}], "]"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", "R"}], "]"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "wd"}], "*", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}], "/", "2"}]}], ",", RowBox[{"wd", "*", "w3", "*", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}], "/", "2"}]}]}], "}"}], ",", "Q"}], "]"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "wd"}], "*", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}], "/", "2"}]}], ",", RowBox[{ RowBox[{"-", "wd"}], "*", "w3", "*", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}], "/", "2"}]}]}], "}"}], ",", "Q"}], "]"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"wd", "*", RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}]}], ",", "0"}], "}"}], ",", "Q"}], "]"}], ",", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "wd"}], "*", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}], "/", "2"}]}], ",", RowBox[{"wd", "*", "w3", "*", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}], "/", "2"}]}]}], "}"}], ",", "0.1"}], "]"}], ",", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "wd"}], "*", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}], "/", "2"}]}], ",", RowBox[{ RowBox[{"-", "wd"}], "*", "w3", "*", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}], "/", "2"}]}]}], "}"}], ",", "0.1"}], "]"}], ",", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"wd", "*", RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}]}], ",", "0"}], "}"}], ",", "0.1"}], "]"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Thickness", "[", "0.012", "]"}], ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"wd", "*", RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}]}], ",", "0"}], "}"}]}], "}"}], "]"}], ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "wd"}], "*", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}], "/", "2"}]}], ",", RowBox[{"wd", "*", "w3", "*", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}], "/", "2"}]}]}], "}"}]}], "}"}], "]"}], ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "wd"}], "*", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}], "/", "2"}]}], ",", RowBox[{ RowBox[{"-", "wd"}], "*", "w3", "*", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}], "/", "2"}]}]}], "}"}]}], "}"}], "]"}]}], "}"}]}], "}"}], "]"}], "/.", "u"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"Show", "[", RowBox[{"plot1", ",", RowBox[{"{", RowBox[{ RowBox[{"AspectRatio", "\[Rule]", "Automatic"}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "3"}], ",", "3"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "3"}], ",", "3"}], "}"}]}], "}"}]}], ",", RowBox[{"Axes", "\[Rule]", "True"}], ",", RowBox[{"GridLines", "\[Rule]", "Automatic"}], ",", RowBox[{"Frame", "\[Rule]", "True"}]}], "}"}]}], "]"}], "\[IndentingNewLine]"}], "Input", CellChangeTimes->{{3.5661428484442425`*^9, 3.5661428783026953`*^9}, { 3.5661430328053665`*^9, 3.5661430369081736`*^9}, {3.5661430998854847`*^9, 3.5661431383537183`*^9}, {3.566143180801393*^9, 3.5661432248645825`*^9}, { 3.566143288683233*^9, 3.5661432893462706`*^9}, {3.566143344240477*^9, 3.566143410119393*^9}}], Cell[BoxData[ GraphicsBox[{CircleBox[{0, 0}, 0.21255884296750083`], DiskBox[{0, 0}, 0.1], CircleBox[{0, 0}, 1.3545219025935082`], CircleBox[{-0.6144517434138547, 1.0642616383920716`}, 1.1215796462307126`], CircleBox[{-0.6144517434138547, -1.0642616383920716`}, 1.1215796462307126`], CircleBox[{1.2289034868277093`, 0}, 1.1215796462307126`], DiskBox[{-0.8813100984419893, 1.5264738677250542`}, 0.1], DiskBox[{-0.8813100984419893, -1.5264738677250542`}, 0.1], DiskBox[{1.7626201968839785`, 0}, 0.1], {Thickness[0.012], LineBox[{{0, 0}, {1.7626201968839785`, 0}}], LineBox[{{0, 0}, {-0.8813100984419893, 1.5264738677250542`}}], LineBox[{{0, 0}, {-0.8813100984419893, -1.5264738677250542`}}]}}, AspectRatio->Automatic, Axes->True, Frame->True, GridLines->Automatic, PlotRange->{{-3, 3}, {-3, 3}}]], "Output", CellChangeTimes->{{3.566146697702366*^9, 3.566146718356802*^9}, 3.566146844639024*^9, 3.566146877648682*^9, 3.5661469388319893`*^9, 3.566146993213685*^9, 3.566147754734747*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"plot3", " ", "=", " ", RowBox[{ RowBox[{"Graphics3D", "[", RowBox[{"{", RowBox[{ RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", " ", "P"}], "]"}], ",", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", " ", "0.1"}], "]"}], ",", " ", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0", ",", RowBox[{"P", "+", "R"}]}], "}"}], ",", " ", "R"}], "]"}], ",", " ", "\[IndentingNewLine]", RowBox[{"{", RowBox[{ RowBox[{"Opacity", "[", "0.5", "]"}], ",", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "wd"}], "*", RowBox[{ RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q"}], ")"}], "/", "2"}]}], ",", " ", RowBox[{"wd", "*", "w3", "*", RowBox[{ RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q"}], ")"}], "/", "2"}]}], ",", RowBox[{"wc", "*", RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}]}]}], "}"}], ",", " ", "Q"}], "]"}], ",", " ", "\[IndentingNewLine]", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "wd"}], "*", RowBox[{ RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q"}], ")"}], "/", "2"}]}], ",", " ", RowBox[{ RowBox[{"-", "wd"}], "*", "w3", "*", RowBox[{ RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q"}], ")"}], "/", "2"}]}], ",", RowBox[{"wc", "*", RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}]}]}], "}"}], ",", " ", "Q"}], "]"}], ",", " ", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"wd", "*", RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q"}], ")"}]}], ",", " ", "0", ",", RowBox[{"wc", "*", RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}]}]}], "}"}], ",", " ", "Q"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Darker", "[", "Red", "]"}], ",", "\[IndentingNewLine]", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "wd"}], "*", RowBox[{ RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q", " ", "+", " ", "S"}], ")"}], "/", "2"}]}], ",", " ", RowBox[{"wd", "*", "w3", "*", RowBox[{ RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q", " ", "+", " ", "S"}], ")"}], "/", "2"}]}], ",", RowBox[{"wc", "*", RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}]}]}], "}"}], ",", " ", "0.06"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "wd"}], "*", RowBox[{ RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q", " ", "+", " ", "S"}], ")"}], "/", "2"}]}], ",", " ", RowBox[{ RowBox[{"-", "wd"}], "*", "w3", "*", RowBox[{ RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q", " ", "+", " ", "S"}], ")"}], "/", "2"}]}], ",", RowBox[{"wc", "*", RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}]}]}], "}"}], ",", " ", "0.06"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"wd", "*", RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q", " ", "+", " ", "S"}], ")"}]}], ",", " ", "0", ",", RowBox[{"wc", "*", RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}]}]}], "}"}], ",", " ", "0.06"}], "]"}]}], "}"}]}], "}"}], "]"}], " ", "/.", " ", "u"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"Show", "[", RowBox[{"plot3", ",", RowBox[{"{", RowBox[{ RowBox[{"AspectRatio", " ", "->", " ", "Automatic"}], ",", " ", RowBox[{"Boxed", " ", "\[Rule]", " ", "False"}]}], "}"}]}], "]"}]}], "Input", CellChangeTimes->{{3.5661464147646685`*^9, 3.566146439287912*^9}, { 3.566146505400828*^9, 3.566146609468611*^9}, 3.5661466450210733`*^9, { 3.5661467117267904`*^9, 3.5661467122883916`*^9}, {3.566146817869377*^9, 3.5661468378998117`*^9}, {3.5661468705194693`*^9, 3.5661468715646715`*^9}, {3.5661469228107615`*^9, 3.5661469312971764`*^9}, {3.5661469641976337`*^9, 3.5661469774888573`*^9}, {3.5661477373577533`*^9, 3.5661477492774353`*^9}}], Cell[BoxData[ Graphics3DBox[{SphereBox[{0, 0, 0}, 0.21255884296750083`], SphereBox[{0, 0, 0}, 0.1], SphereBox[{0, 0, 1.5670807455610092`}, 1.3545219025935082`], {Opacity[0.5], SphereBox[{-0.6144517434138547, 1.0642616383920716`, -0.5193474062925403}, 1.1215796462307126`], SphereBox[{-0.6144517434138547, -1.0642616383920716`, \ -0.5193474062925403}, 1.1215796462307126`], SphereBox[{1.2289034868277093`, 0, -0.5193474062925403}, 1.1215796462307126`]}, {RGBColor[ NCache[ Rational[2, 3], 0.6666666666666666], 0, 0], SphereBox[{-0.8813100984419893, 1.5264738677250542`, -0.744901643898485}, 0.06], SphereBox[{-0.8813100984419893, -1.5264738677250542`, \ -0.744901643898485}, 0.06], SphereBox[{1.7626201968839785`, 0, -0.744901643898485}, 0.06]}}, AspectRatio->Automatic, AutomaticImageSize->True, Boxed->False, ImageSize->{374.28863510256, 341.3514456692233}, ViewPoint->{-2.522546445024874, 2.255104243364661, 0.0355567750375356}, ViewVertical->{-0.5133091384817696, 0.3755035440019775, 0.811894035861656}]], "Output", CellChangeTimes->{3.566146877679882*^9, 3.5661469388475895`*^9, 3.566146993244885*^9, 3.566147755012763*^9}] }, Open ]], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{ RowBox[{"Projection", " ", "onto", " ", "LP"}], "-", "HH", "-", "plane"}], " ", "*)"}]], "Input", CellChangeTimes->{{3.5661428484442425`*^9, 3.5661428783026953`*^9}, { 3.5661430328053665`*^9, 3.5661430369081736`*^9}, {3.5661430998854847`*^9, 3.5661431383537183`*^9}, {3.566143180801393*^9, 3.5661432248645825`*^9}, { 3.566143288683233*^9, 3.5661432893462706`*^9}, {3.566143344240477*^9, 3.5661434217882133`*^9}}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"plot2", "=", RowBox[{ RowBox[{"Graphics", "[", RowBox[{"{", RowBox[{ RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", "P"}], "]"}], ",", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", "0.1"}], "]"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", RowBox[{"P", "+", "R"}]}], "}"}], ",", "R"}], "]"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", RowBox[{"P", "+", "R"}]}], "}"}], ",", "0.08"}], "]"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"wd", "*", RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}]}], ",", RowBox[{"wc", "*", RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}]}]}], "}"}], ",", "Q"}], "]"}], ",", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"wd", "*", RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}]}], ",", RowBox[{"wc", "*", RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}]}]}], "}"}], ",", "0.1"}], "]"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Dashing", "[", RowBox[{"{", "0.02", "}"}], "]"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "wd"}], "*", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}], "/", "2"}]}], ",", RowBox[{"wc", "*", RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}]}]}], "}"}], ",", "Q"}], "]"}]}], "}"}], ",", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "wd"}], "*", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}], "/", "2"}]}], ",", RowBox[{"wc", "*", RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}]}]}], "}"}], ",", "0.1"}], "]"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Thickness", "[", "0.012", "]"}], ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"wd", "*", RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}]}], ",", RowBox[{"wc", "*", RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}]}]}], "}"}]}], "}"}], "]"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Dashing", "[", RowBox[{"{", "0.02", "}"}], "]"}], ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "wd"}], "*", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}], "/", "2"}]}], ",", RowBox[{"wc", "*", RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}]}]}], "}"}]}], "}"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Dashing", "[", RowBox[{"{", RowBox[{"0.02", ",", "0.04"}], "}"}], "]"}], ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", RowBox[{"P", "+", "R"}]}], "}"}]}], "}"}], "]"}]}], "}"}]}], " ", "}"}]}], "}"}], "]"}], "/.", "u"}]}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{"Show", "[", RowBox[{"plot2", ",", RowBox[{"{", RowBox[{ RowBox[{"AspectRatio", "\[Rule]", "Automatic"}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "3"}], ",", "3"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "3"}], ",", "3"}], "}"}]}], "}"}]}], ",", RowBox[{"Axes", "\[Rule]", "True"}], ",", RowBox[{"GridLines", "\[Rule]", "Automatic"}], ",", RowBox[{"Frame", "\[Rule]", "True"}]}], "}"}]}], "]"}]}], "Input", CellChangeTimes->{{3.5661428484442425`*^9, 3.5661428783026953`*^9}, { 3.5661430328053665`*^9, 3.5661430369081736`*^9}, {3.5661430998854847`*^9, 3.5661431383537183`*^9}, {3.566143180801393*^9, 3.5661432248645825`*^9}, { 3.566143288683233*^9, 3.5661432893462706`*^9}, {3.566143344240477*^9, 3.5661434217882133`*^9}}], Cell[BoxData[ GraphicsBox[{CircleBox[{0, 0}, 0.21255884296750083`], DiskBox[{0, 0}, 0.1], CircleBox[{0, 1.5670807455610092`}, 1.3545219025935082`], CircleBox[{0, 1.5670807455610092`}, 0.08], CircleBox[{1.2289034868277093`, -0.5193474062925403}, 1.1215796462307126`], DiskBox[{1.7626201968839785`, -0.744901643898485}, 0.1], {Dashing[{0.02}], CircleBox[{-0.6144517434138547, -0.5193474062925403}, 1.1215796462307126`]}, DiskBox[{-0.8813100984419893, -0.744901643898485}, 0.1], {Thickness[0.012], LineBox[{{0, 0}, {1.7626201968839785`, -0.744901643898485}}], {Dashing[{0.02}], LineBox[{{0, 0}, {-0.8813100984419893, -0.744901643898485}}]}, {Dashing[{0.02, 0.04}], LineBox[{{0, 0}, {0, 1.5670807455610092`}}]}}}, AspectRatio->Automatic, Axes->True, Frame->True, GridLines->Automatic, PlotRange->{{-3, 3}, {-3, 3}}]], "Output", CellChangeTimes->{{3.5661466978895664`*^9, 3.5661467184192023`*^9}, 3.5661468448574247`*^9, 3.5661468777266817`*^9, 3.5661469388943896`*^9, 3.566146993276085*^9, 3.566147755049765*^9}] }, Open ]], Cell[BoxData[ StyleBox[ RowBox[{"(*", " ", RowBox[{ RowBox[{ StyleBox[ RowBox[{ StyleBox["H", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "2O"}]], " ", "with", " ", RowBox[{ StyleBox[ StyleBox[ StyleBox[ StyleBox["O", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "[", RowBox[{ RowBox[{"1", "s2", " ", StyleBox[ StyleBox[ StyleBox[ StyleBox[ StyleBox["2", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "sp"}], ",", " ", RowBox[{ StyleBox[ StyleBox[ StyleBox[ StyleBox["2", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "sp", " ", "LP"}]}], "]"}]}], ",", RowBox[{ RowBox[{"~", "no"}], " ", "overlap"}], ",", " ", RowBox[{"parametrized", " ", "ES", " ", StyleBox[ RowBox[{"0", StyleBox[ StyleBox[ StyleBox["5", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], ".01"}]], ".2013", "\[IndentingNewLine]", RowBox[{"Variables", ":", "\[IndentingNewLine]", StyleBox[ RowBox[{" ", StyleBox[" ", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]]}]], StyleBox[" ", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], RowBox[{ RowBox[{ StyleBox[ RowBox[{" ", " "}], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], StyleBox[ RowBox[{ StyleBox[" ", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], " "}]], RowBox[{"P", StyleBox[ RowBox[{" ", StyleBox[" ", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]]}]], StyleBox[ StyleBox[ StyleBox[ StyleBox["O", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]]}]}], "-", RowBox[{"core", " ", "radius", "\[IndentingNewLine]", "Q", StyleBox[ RowBox[{" ", StyleBox[" ", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]]}]], StyleBox[ StyleBox[ StyleBox[ StyleBox["O", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]]}], "-", RowBox[{"H", " ", "cloud", " ", "radius", "\[IndentingNewLine]", "R", StyleBox[ RowBox[{" ", StyleBox[" ", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]]}]], "lone", " ", "pair", " ", "radius", "\[IndentingNewLine]", "S", StyleBox[ RowBox[{" ", StyleBox[" ", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]]}]], "distance", " ", "proton", " ", "to", " ", "center", " ", "of", " ", StyleBox[ RowBox[{ StyleBox[ StyleBox[ StyleBox[ StyleBox[ StyleBox["O", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "H"}]], " ", "cloud", "\[IndentingNewLine]", "\[Alpha]", StyleBox[ RowBox[{" ", StyleBox[" ", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]]}]], "angle", " ", "lonePair"}], "-", StyleBox[ StyleBox[ StyleBox[ StyleBox[ StyleBox[ StyleBox[ StyleBox["O", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "-", RowBox[{ StyleBox[ RowBox[{ StyleBox[ StyleBox[ StyleBox[ StyleBox["l", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], StyleBox[ StyleBox[ StyleBox[ StyleBox[ StyleBox[ StyleBox[ StyleBox["o", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "n", StyleBox[ StyleBox[ StyleBox[ StyleBox[ StyleBox["e", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Pa", StyleBox[ StyleBox[ StyleBox[ StyleBox[ StyleBox["i", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], StyleBox[ StyleBox[ StyleBox[ StyleBox[ StyleBox[ StyleBox["r", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]]}]], "\[IndentingNewLine]", "\[Beta]", StyleBox[ RowBox[{" ", StyleBox[" ", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]]}]], "angle", " ", StyleBox[ StyleBox[ StyleBox[ StyleBox["H", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]]}], "-", StyleBox[ StyleBox[ StyleBox[ StyleBox[ StyleBox[ StyleBox[ StyleBox["O", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "-", "H"}]}]}]}], " ", "*)"}], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]]], "Input", CellChangeTimes->{{3.5661428484442425`*^9, 3.5661428783026953`*^9}, { 3.5661430328053665`*^9, 3.5661430369081736`*^9}, {3.5661430998854847`*^9, 3.5661431383537183`*^9}, {3.566143180801393*^9, 3.5661432248645825`*^9}, { 3.566143288683233*^9, 3.5661432893462706`*^9}, 3.566143344240477*^9, 3.566143438043442*^9, 3.566143491551536*^9, {3.5661922351554785`*^9, 3.5661923427020674`*^9}, {3.566201924306572*^9, 3.566201971387455*^9}, 3.566202322310071*^9, {3.5663083624775496`*^9, 3.5663083893563967`*^9}, { 3.5663654716175947`*^9, 3.566365471929596*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"Clear", "[", RowBox[{ "wc", ",", "k1", ",", "k2", ",", "k3", ",", "s1", ",", "s2", ",", "s3", ",", "par"}], "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Z", "=", "8.0"}], ";", " ", RowBox[{"w3", "=", RowBox[{"Sqrt", "[", "3", "]"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", "core", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Tc", "=", RowBox[{"2.25", "*", RowBox[{"k1", "/", RowBox[{"P", "^", "2"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vnec", "=", RowBox[{ RowBox[{"-", "3.0"}], "*", RowBox[{"Z", "/", "P"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Veec", "=", RowBox[{"3.0", "*", RowBox[{"s1", "/", "P"}]}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"XH", " ", "bond"}], "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Tb", "=", RowBox[{"2.25", "*", RowBox[{"k2", "/", RowBox[{"Q", "^", "2"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vneb", "=", RowBox[{ RowBox[{"-", "2.0"}], "*", RowBox[{ RowBox[{"(", RowBox[{"1.5", "-", RowBox[{"0.5", "*", RowBox[{ RowBox[{"(", RowBox[{"S", "/", "Q"}], ")"}], "^", "2"}]}]}], ")"}], "/", "Q"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Veeb", "=", RowBox[{"3.0", "*", RowBox[{"s2", "/", "Q"}]}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"lone", " ", "pair"}], "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Tp", "=", RowBox[{"2.25", "*", RowBox[{"k3", "/", RowBox[{"R", "^", "2"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Veep", "=", RowBox[{"3.0", "*", RowBox[{"s3", "/", "R"}]}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"core", "-", "bond"}], "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vnecb", "=", RowBox[{ RowBox[{ RowBox[{"-", "2.0"}], "*", RowBox[{"Z", "/", RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}]}]}], "-", RowBox[{"2", "/", RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Veecb", "=", RowBox[{"4.0", "/", RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Vnncb", "=", RowBox[{"Z", "/", RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}]}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"core", "-", RowBox[{"lone", " ", "pair"}]}], "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vnecp", "=", RowBox[{ RowBox[{"-", "2.0"}], "*", RowBox[{"Z", "/", RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Veecp", "=", RowBox[{"4.0", "/", RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}]}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"bond", "-", "bond"}], "*)"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ RowBox[{"d1", " ", "w", RowBox[{"(", "NH", ")"}]}], "-", RowBox[{"w", RowBox[{"(", "NH", ")"}]}]}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"d1", "=", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}], "*", RowBox[{"Sqrt", "[", RowBox[{"2", "*", RowBox[{"(", RowBox[{"1", "-", RowBox[{"Cos", "[", "\[Beta]", "]"}]}], ")"}]}], "]"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Veebb", "=", RowBox[{"4.0", "/", "d1"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ RowBox[{"d3", " ", "H"}], "-", "H"}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"d3", "=", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}], "*", RowBox[{"Sqrt", "[", RowBox[{"2", "*", RowBox[{"(", RowBox[{"1", "-", RowBox[{"Cos", "[", "\[Beta]", "]"}]}], ")"}]}], "]"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Vnnbb", "=", RowBox[{"1.0", "/", "d3"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ RowBox[{ RowBox[{"d5", " ", "w", RowBox[{"(", "OH", ")"}]}], "-", "H"}], ",", " ", RowBox[{"\[Beta]", " ", "is", " ", "angle", " ", "HOH"}]}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"d5", "=", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}], "^", "2"}], "+", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}], "^", "2"}], "-", RowBox[{"2", "*", RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}], "*", RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}], "*", RowBox[{"Cos", "[", "\[Beta]", "]"}]}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Vnebb", "=", RowBox[{ RowBox[{"-", "2.0"}], "/", "d5"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"bond", "-", RowBox[{"lone", " ", "pair"}]}], "*)"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ RowBox[{"d2", " ", "w", RowBox[{"(", "OH", ")"}]}], "-", RowBox[{"w", RowBox[{"(", "LP", ")"}]}]}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"d2", "=", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}], "^", "2"}], "+", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}], "^", "2"}], "+", RowBox[{"2", "*", RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}], "*", RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}], "*", RowBox[{"Cos", "[", RowBox[{"\[Alpha]", "/", "2"}], "]"}], "*", RowBox[{"Cos", "[", RowBox[{"\[Beta]", "/", "2"}], "]"}]}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Veebp", " ", "=", " ", RowBox[{"4", "/", "d2"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ RowBox[{"d4", " ", "w", RowBox[{"(", "LP", ")"}]}], "-", "H"}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"d4", "=", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}], "^", "2"}], "+", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}], "^", "2"}], "+", RowBox[{"2", "*", RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}], "*", RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}], "*", RowBox[{"Cos", "[", RowBox[{"\[Alpha]", "/", "2"}], "]"}], "*", RowBox[{"Cos", "[", RowBox[{"\[Beta]", "/", "2"}], "]"}]}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Vnebp", " ", "=", RowBox[{ RowBox[{"-", "2.0"}], "/", "d4"}]}], " ", ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", StyleBox[ RowBox[{"(*", RowBox[{ RowBox[{"lone", " ", "pair"}], "-", RowBox[{"lone", " ", "pair"}]}], "*)"}], "Text", FontFamily->"Arial Narrow"], "\[IndentingNewLine]", StyleBox[ RowBox[{"(*", " ", RowBox[{ RowBox[{"none", " ", "for", " ", "NH3"}], ",", " ", RowBox[{"but", " ", "for", " ", "H2O", " ", "and", " ", "HF"}]}], " ", "*)"}], "Text", FontFamily->"Arial Narrow"], "\[IndentingNewLine]", StyleBox[ RowBox[{"(*", " ", RowBox[{ RowBox[{ "H2O", " ", "has", " ", "two", " ", "independent", " ", "angles", " ", "as", " ", RowBox[{"variables", ":", " ", RowBox[{"H", "-", "O", "-", RowBox[{"H", " ", "and", " ", "LP"}], "-", "O", "-", StyleBox[ RowBox[{"L", StyleBox["P", "Text", FontFamily->"Arial Narrow"]}], "Text", FontFamily->"Arial Narrow"]}]}]}], ","}], " ", "*)"}], "Text", FontFamily->"Arial Narrow"]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ StyleBox["d6", "Input"], StyleBox[" ", "Input"], StyleBox["=", "Input"], StyleBox[" ", "Input"], RowBox[{ StyleBox[ RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}], "Input"], "*", RowBox[{"Sqrt", "[", RowBox[{"2", "*", RowBox[{"(", RowBox[{"1", "-", RowBox[{"Cos", "[", "\[Alpha]", "]"}]}], ")"}]}], "]"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Veepp", "=", RowBox[{"4", "/", "d6"}]}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"(*", RowBox[{"sum", " ", "terms"}], "*)"}]}]}], "Input", CellChangeTimes->{{3.5661428484442425`*^9, 3.5661428783026953`*^9}, { 3.5661430328053665`*^9, 3.5661430369081736`*^9}, {3.5661430998854847`*^9, 3.5661431383537183`*^9}, {3.566143180801393*^9, 3.5661432248645825`*^9}, { 3.566143288683233*^9, 3.5661432893462706`*^9}, 3.566143344240477*^9, 3.566143438043442*^9, {3.5661923476784763`*^9, 3.566192348645678*^9}, { 3.566192420296604*^9, 3.566192420858205*^9}, {3.5661924684694886`*^9, 3.566192482369113*^9}, {3.5661925519608355`*^9, 3.5661925980901165`*^9}, { 3.566192662409029*^9, 3.5661926728298473`*^9}, {3.566194733031866*^9, 3.5661947591463118`*^9}, {3.566201172322851*^9, 3.5662012049425087`*^9}, { 3.566201275782233*^9, 3.5662012764842343`*^9}, {3.566201607563616*^9, 3.5662016091548185`*^9}, {3.5662016860473537`*^9, 3.5662017455770583`*^9}, {3.566201840690425*^9, 3.5662018510176435`*^9}, { 3.5662021072324934`*^9, 3.56620213362774*^9}, {3.566202199116655*^9, 3.566202199381855*^9}, {3.5663049674119854`*^9, 3.5663050328385005`*^9}, { 3.566305081167385*^9, 3.566305085270192*^9}, {3.56630512948067*^9, 3.566305132819076*^9}, {3.5663056216303344`*^9, 3.56630564170757*^9}, 3.566305710550491*^9, {3.5663057513601627`*^9, 3.566305753419366*^9}, { 3.566307371173808*^9, 3.5663074141830835`*^9}, 3.566307457207959*^9, 3.5663085639675035`*^9, {3.5663156573399625`*^9, 3.5663156653739767`*^9}, { 3.5663157659941535`*^9, 3.5663157671173553`*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"T", "=", RowBox[{"Tc", "+", " ", RowBox[{"2", "*", "Tb"}], "+", " ", RowBox[{"2", "*", "Tp"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vne", "=", RowBox[{"Vnec", "+", RowBox[{"2", "*", "Vneb"}], "+", " ", RowBox[{"2", "*", "Vnecb"}], "+", RowBox[{"2", "*", "Vnecp"}], "+", " ", RowBox[{"2", "*", "Vnebb"}], "+", RowBox[{"4", "*", "Vnebp"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vee", "=", RowBox[{"Veec", "+", RowBox[{"2", "*", "Veeb"}], "+", " ", RowBox[{"2", "*", "Veep"}], "+", RowBox[{"2", "*", "Veecb"}], "+", RowBox[{"2", "*", "Veecp"}], "+", " ", "Veebb", "+", RowBox[{"4", "*", "Veebp"}], " ", "+", " ", "Veepp"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ "Vnn", "=", " ", RowBox[{ RowBox[{"2", "*", "Vnncb"}], "+", " ", "Vnnbb"}]}], ";"}], "\[IndentingNewLine]"}], "Input", CellChangeTimes->{{3.5661428484442425`*^9, 3.5661428783026953`*^9}, { 3.5661430328053665`*^9, 3.5661430369081736`*^9}, {3.5661430998854847`*^9, 3.5661431383537183`*^9}, {3.566143180801393*^9, 3.5661432248645825`*^9}, { 3.566143288683233*^9, 3.5661432893462706`*^9}, 3.566143344240477*^9, { 3.5661435118159714`*^9, 3.5661435187891836`*^9}, {3.5661470921646585`*^9, 3.5661471085134873`*^9}, {3.566194770659132*^9, 3.566194777023943*^9}, { 3.566200993811738*^9, 3.566201133478783*^9}, {3.5662012401361704`*^9, 3.5662012446913786`*^9}, {3.5662014429677267`*^9, 3.566201443654128*^9}, { 3.566201792735941*^9, 3.5662017938591433`*^9}, 3.5662022754943886`*^9, { 3.5662025519112744`*^9, 3.566202565436498*^9}, {3.566305456238844*^9, 3.566305473289674*^9}, {3.5663055060497313`*^9, 3.5663055070949335`*^9}, { 3.566308222654504*^9, 3.5663082634017754`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ RowBox[{"Parameters", ":", " ", "G3"}], "//", RowBox[{ RowBox[{ RowBox[{ RowBox[{"6", "-", RowBox[{"311", "G", RowBox[{"(", RowBox[{"d", ",", "p"}], ")"}], " ", "can", " ", "be", " ", "read", " ", "for", " ", "LiH"}]}], " ", "..."}], " ", "HF", "\[IndentingNewLine]", "in", " ", "table", " ", "from", " ", RowBox[{"www", ".", "kimball"}]}], "-", RowBox[{ RowBox[{ RowBox[{"model", ".", "org"}], "/", "Kimball"}], "/", RowBox[{"params", ".", "xls"}]}]}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"par", "=", RowBox[{"{", RowBox[{ RowBox[{"k1", "\[Rule]", "0.92438558"}], ",", RowBox[{"k2", "\[Rule]", "1.42313501"}], ",", RowBox[{"k3", "\[Rule]", "1.32139847"}], ",", RowBox[{"s1", "\[Rule]", "0.43732345"}], ",", RowBox[{"s2", "\[Rule]", "0.31779336"}], ",", RowBox[{"s3", "\[Rule]", "0.44741336"}]}], "}"}]}], ";"}]}]], "Input", CellChangeTimes->{{3.5661428484442425`*^9, 3.5661428783026953`*^9}, { 3.5661430328053665`*^9, 3.5661430369081736`*^9}, {3.5661430998854847`*^9, 3.5661431383537183`*^9}, {3.566143180801393*^9, 3.5661432248645825`*^9}, { 3.566143288683233*^9, 3.5661432893462706`*^9}, {3.566143344240477*^9, 3.566143371150524*^9}, {3.5663152612084665`*^9, 3.5663152707400837`*^9}, { 3.566315729474489*^9, 3.566315733312096*^9}, {3.566315805368623*^9, 3.5663158109066324`*^9}, {3.566316505903453*^9, 3.566316526277089*^9}}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"func", "=", RowBox[{ RowBox[{"T", "+", "Vne", "+", "Vee", "+", "Vnn"}], "/.", "par"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"t", "=", RowBox[{"FindMinimum", "[", RowBox[{"func", ",", RowBox[{"{", RowBox[{"P", ",", "0.18"}], "}"}], ",", RowBox[{"{", RowBox[{"Q", ",", "1.0"}], "}"}], ",", RowBox[{"{", RowBox[{"R", ",", " ", "1.1"}], "}"}], ",", RowBox[{"{", RowBox[{"S", ",", "0.6"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"\[Alpha]", ",", "1.97"}], "}"}], ",", RowBox[{"{", RowBox[{"\[Beta]", ",", "1.97"}], "}"}], ",", RowBox[{"MaxIterations", "\[Rule]", "500"}], ",", " ", RowBox[{"Method", " ", "\[Rule]", " ", "\"\\""}]}], "]"}]}]}], "Input", CellChangeTimes->{{3.5661428484442425`*^9, 3.5661428783026953`*^9}, { 3.5661430328053665`*^9, 3.5661430369081736`*^9}, {3.5661430998854847`*^9, 3.5661431383537183`*^9}, {3.566143180801393*^9, 3.5661432248645825`*^9}, { 3.566143288683233*^9, 3.5661432893462706`*^9}, {3.566143344240477*^9, 3.5661433819769435`*^9}, {3.5662024426798825`*^9, 3.5662024734587364`*^9}, {3.566305156858718*^9, 3.566305166265535*^9}, { 3.5663079842392845`*^9, 3.566307984972486*^9}, 3.5663654313383245`*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"-", "76.3817728692884`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.17819997281079103`"}], ",", RowBox[{"Q", "\[Rule]", "1.0173065142262625`"}], ",", RowBox[{"R", "\[Rule]", "1.12517354005205`"}], ",", RowBox[{"S", "\[Rule]", "0.5897462840993366`"}], ",", RowBox[{"\[Alpha]", "\[Rule]", "2.1112470512532235`"}], ",", RowBox[{"\[Beta]", "\[Rule]", "1.8036924790864421`"}]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{ 3.5663165372283077`*^9, 3.5663175724229746`*^9, 3.5663176162590513`*^9, 3.566317659330727*^9, 3.566320373812928*^9, 3.566320436821439*^9, 3.566320474401905*^9, 3.566321519556941*^9, 3.5663217066792693`*^9, 3.5663218529451265`*^9, 3.5663626992899256`*^9, 3.566362965987994*^9, 3.566363116185058*^9, 3.5663631564331284`*^9, 3.566363452786849*^9, 3.5663635926410947`*^9, 3.5663642788550997`*^9, 3.566364436430977*^9, { 3.5663645631343994`*^9, 3.5663645911676483`*^9}, 3.566365060759673*^9, 3.5663651610678496`*^9, 3.566365211206338*^9, 3.566365290938078*^9, { 3.5663653439157705`*^9, 3.5663653720114202`*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"u", "=", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"Vne", "/.", "u"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vee", "/.", "u"}], "/.", "par"}], "\[IndentingNewLine]", RowBox[{"Vnn", "/.", "u"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"Vne", "+", "Vee", "+", "Vnn"}], ")"}], "/", "T"}], "/.", "u"}], "/.", "par"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}], "*", "0.529177"}], "/.", "u"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{"angle", " ", "HNH"}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"180", "/", "\[Pi]"}], "*", "\[Beta]"}], "/.", "u"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ RowBox[{"angle", " ", "LP"}], "-", "N", "-", "LP"}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"180", "/", "\[Pi]"}], "*", "\[Alpha]"}], " ", "/.", " ", "u"}], "\[IndentingNewLine]", RowBox[{"d1", " ", "/.", "u"}], "\[IndentingNewLine]", RowBox[{"d3", " ", "/.", "u"}], "\[IndentingNewLine]", RowBox[{"d5", " ", "/.", "u"}], "\[IndentingNewLine]", RowBox[{"d2", " ", "/.", "u"}], "\[IndentingNewLine]", RowBox[{"d4", " ", "/.", "u"}]}], "Input", CellChangeTimes->{{3.5661428484442425`*^9, 3.5661428783026953`*^9}, { 3.5661430328053665`*^9, 3.5661430369081736`*^9}, {3.5661430998854847`*^9, 3.5661431383537183`*^9}, {3.566143180801393*^9, 3.5661432248645825`*^9}, { 3.566143288683233*^9, 3.5661432893462706`*^9}, {3.566143344240477*^9, 3.5661433961729684`*^9}, {3.5663078436050377`*^9, 3.5663078735570903`*^9}, 3.5663654431007447`*^9}], Cell[BoxData[ RowBox[{"-", "198.35190730103852`"}]], "Output", CellChangeTimes->{ 3.566316537243908*^9, 3.566317572438575*^9, 3.5663176162590513`*^9, 3.5663176593463273`*^9, 3.5663203738285284`*^9, 3.566320436837039*^9, 3.5663204744175053`*^9, 3.5663215195725408`*^9, 3.5663217066948695`*^9, 3.5663218529607267`*^9, 3.566362699305526*^9, 3.566362966003594*^9, 3.566363116185058*^9, 3.5663631564487286`*^9, 3.566363452802449*^9, 3.566363592656695*^9, 3.5663642788707*^9, 3.566364436446577*^9, { 3.5663645631343994`*^9, 3.5663645911676483`*^9}, 3.566365060759673*^9, 3.5663651610678496`*^9, 3.5663652112219377`*^9, 3.5663652909536777`*^9, { 3.5663653439157705`*^9, 3.5663653720114202`*^9}}], Cell[BoxData["36.269025578002974`"], "Output", CellChangeTimes->{ 3.566316537243908*^9, 3.566317572438575*^9, 3.5663176162590513`*^9, 3.5663176593463273`*^9, 3.5663203738285284`*^9, 3.566320436837039*^9, 3.5663204744175053`*^9, 3.5663215195725408`*^9, 3.5663217066948695`*^9, 3.5663218529607267`*^9, 3.566362699305526*^9, 3.566362966003594*^9, 3.566363116185058*^9, 3.5663631564487286`*^9, 3.566363452802449*^9, 3.566363592656695*^9, 3.5663642788707*^9, 3.566364436446577*^9, { 3.5663645631343994`*^9, 3.5663645911676483`*^9}, 3.566365060759673*^9, 3.5663651610678496`*^9, 3.5663652112219377`*^9, 3.5663652909536777`*^9, { 3.5663653439157705`*^9, 3.5663653720114202`*^9}}], Cell[BoxData["9.319335984458716`"], "Output", CellChangeTimes->{ 3.566316537243908*^9, 3.566317572438575*^9, 3.5663176162590513`*^9, 3.5663176593463273`*^9, 3.5663203738285284`*^9, 3.566320436837039*^9, 3.5663204744175053`*^9, 3.5663215195725408`*^9, 3.5663217066948695`*^9, 3.5663218529607267`*^9, 3.566362699305526*^9, 3.566362966003594*^9, 3.566363116185058*^9, 3.5663631564487286`*^9, 3.566363452802449*^9, 3.566363592656695*^9, 3.5663642788707*^9, 3.566364436446577*^9, { 3.5663645631343994`*^9, 3.5663645911676483`*^9}, 3.566365060759673*^9, 3.5663651610678496`*^9, 3.5663652112219377`*^9, 3.5663652909536777`*^9, { 3.5663653439157705`*^9, 3.56636537202702*^9}}], Cell[BoxData[ RowBox[{"-", "2.`"}]], "Output", CellChangeTimes->{ 3.566316537243908*^9, 3.566317572438575*^9, 3.5663176162590513`*^9, 3.5663176593463273`*^9, 3.5663203738285284`*^9, 3.566320436837039*^9, 3.5663204744175053`*^9, 3.5663215195725408`*^9, 3.5663217066948695`*^9, 3.5663218529607267`*^9, 3.566362699305526*^9, 3.566362966003594*^9, 3.566363116185058*^9, 3.5663631564487286`*^9, 3.566363452802449*^9, 3.566363592656695*^9, 3.5663642788707*^9, 3.566364436446577*^9, { 3.5663645631343994`*^9, 3.5663645911676483`*^9}, 3.566365060759673*^9, 3.5663651610678496`*^9, 3.5663652112219377`*^9, 3.5663652909536777`*^9, { 3.5663653439157705`*^9, 3.56636537202702*^9}}], Cell[BoxData["0.9447147056716415`"], "Output", CellChangeTimes->{ 3.566316537243908*^9, 3.566317572438575*^9, 3.5663176162590513`*^9, 3.5663176593463273`*^9, 3.5663203738285284`*^9, 3.566320436837039*^9, 3.5663204744175053`*^9, 3.5663215195725408`*^9, 3.5663217066948695`*^9, 3.5663218529607267`*^9, 3.566362699305526*^9, 3.566362966003594*^9, 3.566363116185058*^9, 3.5663631564487286`*^9, 3.566363452802449*^9, 3.566363592656695*^9, 3.5663642788707*^9, 3.566364436446577*^9, { 3.5663645631343994`*^9, 3.5663645911676483`*^9}, 3.566365060759673*^9, 3.5663651610678496`*^9, 3.5663652112219377`*^9, 3.5663652909536777`*^9, { 3.5663653439157705`*^9, 3.56636537202702*^9}}], Cell[BoxData["103.34396659114164`"], "Output", CellChangeTimes->{ 3.566316537243908*^9, 3.566317572438575*^9, 3.5663176162590513`*^9, 3.5663176593463273`*^9, 3.5663203738285284`*^9, 3.566320436837039*^9, 3.5663204744175053`*^9, 3.5663215195725408`*^9, 3.5663217066948695`*^9, 3.5663218529607267`*^9, 3.566362699305526*^9, 3.566362966003594*^9, 3.566363116185058*^9, 3.5663631564487286`*^9, 3.566363452802449*^9, 3.566363592656695*^9, 3.5663642788707*^9, 3.566364436446577*^9, { 3.5663645631343994`*^9, 3.5663645911676483`*^9}, 3.566365060759673*^9, 3.5663651610678496`*^9, 3.5663652112219377`*^9, 3.5663652909536777`*^9, { 3.5663653439157705`*^9, 3.56636537202702*^9}}], Cell[BoxData["120.96554554624991`"], "Output", CellChangeTimes->{ 3.566316537243908*^9, 3.566317572438575*^9, 3.5663176162590513`*^9, 3.5663176593463273`*^9, 3.5663203738285284`*^9, 3.566320436837039*^9, 3.5663204744175053`*^9, 3.5663215195725408`*^9, 3.5663217066948695`*^9, 3.5663218529607267`*^9, 3.566362699305526*^9, 3.566362966003594*^9, 3.566363116185058*^9, 3.5663631564487286`*^9, 3.566363452802449*^9, 3.566363592656695*^9, 3.5663642788707*^9, 3.566364436446577*^9, { 3.5663645631343994`*^9, 3.5663645911676483`*^9}, 3.566365060759673*^9, 3.5663651610678496`*^9, 3.5663652112219377`*^9, 3.5663652909536777`*^9, { 3.5663653439157705`*^9, 3.56636537202702*^9}}], Cell[BoxData["1.8756856345332227`"], "Output", CellChangeTimes->{ 3.566316537243908*^9, 3.566317572438575*^9, 3.5663176162590513`*^9, 3.5663176593463273`*^9, 3.5663203738285284`*^9, 3.566320436837039*^9, 3.5663204744175053`*^9, 3.5663215195725408`*^9, 3.5663217066948695`*^9, 3.5663218529607267`*^9, 3.566362699305526*^9, 3.566362966003594*^9, 3.566363116185058*^9, 3.5663631564487286`*^9, 3.566363452802449*^9, 3.566363592656695*^9, 3.5663642788707*^9, 3.566364436446577*^9, { 3.5663645631343994`*^9, 3.5663645911676483`*^9}, 3.566365060759673*^9, 3.5663651610678496`*^9, 3.5663652112219377`*^9, 3.5663652909536777`*^9, { 3.5663653439157705`*^9, 3.56636537202702*^9}}], Cell[BoxData["2.8009659614062543`"], "Output", CellChangeTimes->{ 3.566316537243908*^9, 3.566317572438575*^9, 3.5663176162590513`*^9, 3.5663176593463273`*^9, 3.5663203738285284`*^9, 3.566320436837039*^9, 3.5663204744175053`*^9, 3.5663215195725408`*^9, 3.5663217066948695`*^9, 3.5663218529607267`*^9, 3.566362699305526*^9, 3.566362966003594*^9, 3.566363116185058*^9, 3.5663631564487286`*^9, 3.566363452802449*^9, 3.566363592656695*^9, 3.5663642788707*^9, 3.566364436446577*^9, { 3.5663645631343994`*^9, 3.5663645911676483`*^9}, 3.566365060759673*^9, 3.5663651610678496`*^9, 3.5663652112219377`*^9, 3.5663652909536777`*^9, { 3.5663653439157705`*^9, 3.56636537202702*^9}}], Cell[BoxData["2.366755647766626`"], "Output", CellChangeTimes->{ 3.566316537243908*^9, 3.566317572438575*^9, 3.5663176162590513`*^9, 3.5663176593463273`*^9, 3.5663203738285284`*^9, 3.566320436837039*^9, 3.5663204744175053`*^9, 3.5663215195725408`*^9, 3.5663217066948695`*^9, 3.5663218529607267`*^9, 3.566362699305526*^9, 3.566362966003594*^9, 3.566363116185058*^9, 3.5663631564487286`*^9, 3.566363452802449*^9, 3.566363592656695*^9, 3.5663642788707*^9, 3.566364436446577*^9, { 3.5663645631343994`*^9, 3.5663645911676483`*^9}, 3.566365060759673*^9, 3.5663651610678496`*^9, 3.5663652112219377`*^9, 3.5663652909536777`*^9, { 3.5663653439157705`*^9, 3.56636537202702*^9}}], Cell[BoxData["2.019953840865406`"], "Output", CellChangeTimes->{ 3.566316537243908*^9, 3.566317572438575*^9, 3.5663176162590513`*^9, 3.5663176593463273`*^9, 3.5663203738285284`*^9, 3.566320436837039*^9, 3.5663204744175053`*^9, 3.5663215195725408`*^9, 3.5663217066948695`*^9, 3.5663218529607267`*^9, 3.566362699305526*^9, 3.566362966003594*^9, 3.566363116185058*^9, 3.5663631564487286`*^9, 3.566363452802449*^9, 3.566363592656695*^9, 3.5663642788707*^9, 3.566364436446577*^9, { 3.5663645631343994`*^9, 3.5663645911676483`*^9}, 3.566365060759673*^9, 3.5663651610678496`*^9, 3.5663652112219377`*^9, 3.5663652909536777`*^9, { 3.5663653439157705`*^9, 3.56636537202702*^9}}], Cell[BoxData["2.5115384660834215`"], "Output", CellChangeTimes->{ 3.566316537243908*^9, 3.566317572438575*^9, 3.5663176162590513`*^9, 3.5663176593463273`*^9, 3.5663203738285284`*^9, 3.566320436837039*^9, 3.5663204744175053`*^9, 3.5663215195725408`*^9, 3.5663217066948695`*^9, 3.5663218529607267`*^9, 3.566362699305526*^9, 3.566362966003594*^9, 3.566363116185058*^9, 3.5663631564487286`*^9, 3.566363452802449*^9, 3.566363592656695*^9, 3.5663642788707*^9, 3.566364436446577*^9, { 3.5663645631343994`*^9, 3.5663645911676483`*^9}, 3.566365060759673*^9, 3.5663651610678496`*^9, 3.5663652112219377`*^9, 3.5663652909536777`*^9, { 3.5663653439157705`*^9, 3.56636537202702*^9}}] }, Open ]], Cell[BoxData[ RowBox[{"(*", " ", RowBox[{ RowBox[{ RowBox[{"Projection", " ", "onto", " ", "HOH"}], "-", "plane"}], ",", " ", RowBox[{"LonePairs", " ", "dashed"}]}], " ", "*)"}]], "Input", CellChangeTimes->{{3.5661428484442425`*^9, 3.5661428783026953`*^9}, { 3.5661430328053665`*^9, 3.5661430369081736`*^9}, {3.5661430998854847`*^9, 3.5661431383537183`*^9}, {3.566143180801393*^9, 3.5661432248645825`*^9}, { 3.566143288683233*^9, 3.5661432893462706`*^9}, {3.566143344240477*^9, 3.566143410119393*^9}, {3.5663626853591013`*^9, 3.5663626922387133`*^9}, { 3.5663637103901014`*^9, 3.566363711466503*^9}, {3.566365546279326*^9, 3.5663655869329977`*^9}, 3.566369639264944*^9}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"plot1", "=", RowBox[{ RowBox[{"Graphics", "[", RowBox[{"{", RowBox[{ RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", "P"}], "]"}], ",", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", "0.1"}], "]"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Dashing", "[", RowBox[{"{", RowBox[{"0.02", ",", "0.02"}], "}"}], "]"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}]}], "*", RowBox[{"Cos", "[", RowBox[{"\[Alpha]", "/", "2"}], "]"}], "*", RowBox[{"Cos", "[", RowBox[{"\[Beta]", "/", "2"}], "]"}]}], ",", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}], "*", RowBox[{"Cos", "[", RowBox[{"\[Alpha]", "/", "2"}], "]"}], "*", RowBox[{"Sin", "[", RowBox[{"\[Beta]", "/", "2"}], "]"}]}]}], "}"}], ",", "R"}], "]"}]}], "}"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}]}], "*", RowBox[{"Sin", "[", RowBox[{"\[Beta]", "-", RowBox[{"\[Pi]", "/", "2"}]}], "]"}]}], ",", RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}]}], "*", RowBox[{"Cos", "[", RowBox[{"\[Beta]", "-", RowBox[{"\[Pi]", "/", "2"}]}], "]"}]}]}], "}"}], ",", "Q"}], "]"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"P", "+", "Q"}], ",", "0"}], "}"}], ",", "Q"}], "]"}], ",", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}]}], "*", RowBox[{"Sin", "[", RowBox[{"\[Beta]", "-", RowBox[{"\[Pi]", "/", "2"}]}], "]"}]}], ",", RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}]}], "*", RowBox[{"Cos", "[", RowBox[{"\[Beta]", "-", RowBox[{"\[Pi]", "/", "2"}]}], "]"}]}]}], "}"}], ",", "0.1"}], "]"}], ",", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"P", "+", "Q", "+", "S"}], ",", "0"}], "}"}], ",", "0.1"}], "]"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Thickness", "[", "0.012", "]"}], ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "+", "Q", "+", "S"}], ",", "0"}], "}"}]}], "}"}], "]"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Dashing", "[", RowBox[{"{", RowBox[{"0.01", ",", "0.01"}], "}"}], "]"}], ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}]}], "*", RowBox[{"Cos", "[", RowBox[{"\[Alpha]", "/", "2"}], "]"}], "*", RowBox[{"Cos", "[", RowBox[{"\[Beta]", "/", "2"}], "]"}]}], ",", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}], "*", RowBox[{"Cos", "[", RowBox[{"\[Alpha]", "/", "2"}], "]"}], "*", RowBox[{"Sin", "[", RowBox[{"\[Beta]", "/", "2"}], "]"}]}]}], "}"}]}], "}"}], "]"}]}], "}"}], ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}]}], "*", RowBox[{"Sin", "[", RowBox[{"\[Beta]", "-", RowBox[{"\[Pi]", "/", "2"}]}], "]"}]}], ",", RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}]}], "*", RowBox[{"Cos", "[", RowBox[{"\[Beta]", "-", RowBox[{"\[Pi]", "/", "2"}]}], "]"}]}]}], "}"}]}], "}"}], "]"}]}], "}"}]}], "}"}], "]"}], "/.", "u"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"Show", "[", RowBox[{"plot1", ",", RowBox[{"{", RowBox[{ RowBox[{"AspectRatio", "\[Rule]", "Automatic"}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "2"}], ",", "2.5"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "2.5"}], ",", "2"}], "}"}]}], "}"}]}], ",", RowBox[{"Axes", "\[Rule]", "True"}], ",", RowBox[{"GridLines", "\[Rule]", "Automatic"}], ",", RowBox[{"Frame", "\[Rule]", "True"}]}], "}"}]}], "]"}], "\[IndentingNewLine]"}], "Input", CellChangeTimes->{{3.5661428484442425`*^9, 3.5661428783026953`*^9}, { 3.5661430328053665`*^9, 3.5661430369081736`*^9}, {3.5661430998854847`*^9, 3.5661431383537183`*^9}, {3.566143180801393*^9, 3.5661432248645825`*^9}, { 3.566143288683233*^9, 3.5661432893462706`*^9}, {3.566143344240477*^9, 3.566143410119393*^9}, 3.566308742431816*^9, {3.5663176430442986`*^9, 3.566317682730768*^9}, {3.5663204234522157`*^9, 3.566320431267829*^9}, { 3.5663204642618876`*^9, 3.5663204702366977`*^9}, {3.566362812374524*^9, 3.5663628459301834`*^9}, {3.5663629092506943`*^9, 3.566362925334323*^9}, { 3.5663630995242286`*^9, 3.5663631074646425`*^9}, {3.566363142408704*^9, 3.5663631503023176`*^9}, {3.566363335505843*^9, 3.5663634250188007`*^9}, { 3.5663635259509773`*^9, 3.566363546730214*^9}, 3.5663635877426863`*^9, { 3.5663641657237015`*^9, 3.5663642640662737`*^9}, 3.5663644285061626`*^9, 3.5663644816710563`*^9, {3.5663645292199397`*^9, 3.566364585130438*^9}, { 3.5663650531156597`*^9, 3.566365054753663*^9}, {3.566365147714226*^9, 3.5663651764494762`*^9}, {3.566365323167734*^9, 3.566365365833809*^9}}], Cell[BoxData[ GraphicsBox[{CircleBox[{0, 0}, 0.17819997281079103`], DiskBox[{0, 0}, 0.1], {Dashing[{0.02, 0.02}], CircleBox[{-0.39823926301175055`, 0.5037517326518361}, 1.12517354005205]}, CircleBox[{-0.2759186527435787, -1.1632302685263332`}, 1.0173065142262625`], CircleBox[{1.1955064870370535`, 0}, 1.0173065142262625`], DiskBox[{-0.41203000130891465`, -1.7370546148211736`}, 0.1], DiskBox[{1.7852527711363901`, 0}, 0.1], {Thickness[0.012], LineBox[{{0, 0}, {1.7852527711363901`, 0}}], {Dashing[{0.01, 0.01}], LineBox[{{0, 0}, {-0.39823926301175055`, 0.5037517326518361}}]}, LineBox[{{0, 0}, {-0.41203000130891465`, -1.7370546148211736`}}]}}, AspectRatio->Automatic, Axes->True, Frame->True, GridLines->Automatic, PlotRange->{{-2, 2.5}, {-2.5, 2}}]], "Output", CellChangeTimes->{ 3.5663165372751083`*^9, 3.5663175724697742`*^9, 3.566317616305852*^9, 3.566317659377527*^9, 3.566320373875328*^9, 3.566320436883839*^9, 3.5663204744487047`*^9, 3.5663215196037407`*^9, 3.5663217067260695`*^9, 3.566321852991926*^9, 3.5663626993367257`*^9, 3.566362966034794*^9, 3.566363116216258*^9, 3.566363156464329*^9, 3.566363452833649*^9, 3.5663635926722946`*^9, 3.5663642789019003`*^9, 3.5663644364777765`*^9, { 3.566364563181199*^9, 3.566364591214449*^9}, 3.5663650608064737`*^9, 3.56636516111465*^9, 3.5663652112531376`*^9, 3.5663652909848776`*^9, { 3.566365343946971*^9, 3.56636537205822*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"plot3", " ", "=", " ", RowBox[{ RowBox[{"Graphics3D", "[", RowBox[{"{", RowBox[{ RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", " ", "P"}], "]"}], ",", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", " ", "0.1"}], "]"}], ",", " ", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}], "*", RowBox[{"Sin", "[", RowBox[{"\[Alpha]", "/", "2"}], "]"}]}], ",", "0", ",", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}], "*", RowBox[{"Cos", "[", RowBox[{"\[Alpha]", "/", "2"}], "]"}]}]}], "}"}], ",", " ", "R"}], "]"}], ",", " ", "\[IndentingNewLine]", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}]}], "*", RowBox[{"Sin", "[", RowBox[{"\[Alpha]", "/", "2"}], "]"}]}], ",", "0", ",", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}], "*", RowBox[{"Cos", "[", RowBox[{"\[Alpha]", "/", "2"}], "]"}]}]}], "}"}], ",", " ", "R"}], "]"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Opacity", "[", "0.5", "]"}], ",", " ", "\[IndentingNewLine]", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", RowBox[{ RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q"}], ")"}], "*", RowBox[{"Sin", "[", RowBox[{"\[Beta]", "/", "2"}], "]"}]}], ",", " ", RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q"}], ")"}]}], "*", RowBox[{"Cos", "[", RowBox[{"\[Beta]", "/", "2"}], "]"}]}]}], "}"}], ",", " ", "Q"}], "]"}], ",", " ", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q"}], ")"}]}], "*", RowBox[{"Sin", "[", RowBox[{"\[Beta]", "/", "2"}], "]"}]}], ",", " ", RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q"}], ")"}]}], "*", RowBox[{"Cos", "[", RowBox[{"\[Beta]", "/", "2"}], "]"}]}]}], "}"}], ",", " ", "Q"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Darker", "[", "Red", "]"}], ",", "\[IndentingNewLine]", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", RowBox[{ RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q", "+", "S"}], ")"}], "*", RowBox[{"Sin", "[", RowBox[{"\[Beta]", "/", "2"}], "]"}]}], ",", " ", RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q", "+", "S"}], ")"}]}], "*", RowBox[{"Cos", "[", RowBox[{"\[Beta]", "/", "2"}], "]"}]}]}], "}"}], ",", " ", "0.06"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q", "+", "S"}], ")"}]}], "*", RowBox[{"Sin", "[", RowBox[{"\[Beta]", "/", "2"}], "]"}]}], ",", " ", RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q", "+", "S"}], ")"}]}], "*", RowBox[{"Cos", "[", RowBox[{"\[Beta]", "/", "2"}], "]"}]}]}], "}"}], ",", " ", "0.06"}], "]"}]}], "}"}]}], "}"}], "]"}], " ", "/.", " ", "u"}]}], ";"}], " "}], "\[IndentingNewLine]", RowBox[{"Show", "[", RowBox[{"plot3", ",", RowBox[{"{", RowBox[{ RowBox[{"AspectRatio", " ", "->", " ", "Automatic"}], ",", " ", RowBox[{"Boxed", " ", "\[Rule]", " ", "False"}]}], "}"}]}], "]"}]}], "Input", CellChangeTimes->{{3.5661464147646685`*^9, 3.566146439287912*^9}, { 3.566146505400828*^9, 3.566146609468611*^9}, 3.5661466450210733`*^9, { 3.5661467117267904`*^9, 3.5661467122883916`*^9}, {3.566146817869377*^9, 3.5661468378998117`*^9}, {3.5661468705194693`*^9, 3.5661468715646715`*^9}, {3.5661469228107615`*^9, 3.5661469312971764`*^9}, {3.5661469641976337`*^9, 3.5661469774888573`*^9}, {3.5661477373577533`*^9, 3.5661477492774353`*^9}, {3.5663088744236484`*^9, 3.566308893580482*^9}, { 3.566317471849598*^9, 3.5663175599429526`*^9}, {3.5663175959010158`*^9, 3.56631760969144*^9}, {3.566321253732474*^9, 3.5663213099081726`*^9}, { 3.5663213595474596`*^9, 3.5663214564548297`*^9}, {3.5663214993705053`*^9, 3.566321505719716*^9}, {3.5663216132507052`*^9, 3.5663216271815295`*^9}, { 3.56632166145479*^9, 3.566321670892806*^9}, 3.5663217012504597`*^9, { 3.566321797674229*^9, 3.566321848187118*^9}, 3.566370359268608*^9}], Cell[BoxData[ Graphics3DBox[{SphereBox[{0, 0, 0}, 0.17819997281079103`], SphereBox[{0, 0, 0}, 0.1], SphereBox[{1.1342055348474336`, 0, 0.6421528780235041}, 1.12517354005205], SphereBox[{-1.1342055348474336`, 0, 0.6421528780235041}, 1.12517354005205], {Opacity[0.5], SphereBox[{0, 0.9378428172666112, -0.7414085315459367}, 1.0173065142262625`], SphereBox[{0, -0.9378428172666112, -0.7414085315459367}, 1.0173065142262625`]}, {RGBColor[ NCache[ Rational[2, 3], 0.6666666666666666], 0, 0], SphereBox[{0, 1.400482980703127, -1.107147179742172}, 0.06], SphereBox[{0, -1.400482980703127, -1.107147179742172}, 0.06]}}, AspectRatio->Automatic, Boxed->False, ImageSize->{451.47486168087, 326.}, ViewPoint->{2.33056220775643, 2.452889092660757, 0.04260158310888008}, ViewVertical->{0.19640236579812537`, 0.5622913296981552, 1.0909410783465248`}]], "Output", CellChangeTimes->{ 3.5663165372907085`*^9, 3.5663175725009747`*^9, 3.5663176163214517`*^9, 3.566317659408727*^9, 3.5663203738909283`*^9, 3.566320436899439*^9, 3.566320474479905*^9, 3.566321519634941*^9, 3.56632170675727*^9, 3.5663218530231266`*^9, 3.5663626996643267`*^9, 3.566362966050394*^9, 3.566363116247458*^9, 3.5663631564955287`*^9, 3.566363452849249*^9, 3.566363592703495*^9, 3.5663642789175*^9, 3.566364436508977*^9, { 3.5663645631967993`*^9, 3.566364591245649*^9}, 3.5663650608220735`*^9, 3.5663651611302495`*^9, 3.566365211268738*^9, 3.566365291000478*^9, { 3.566365343978171*^9, 3.56636537207382*^9}}] }, Open ]], Cell[BoxData[ StyleBox[ RowBox[{"(*", " ", RowBox[{ RowBox[{ StyleBox[ RowBox[{"H", StyleBox[ StyleBox["F", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]]}], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], " ", "with", " ", RowBox[{ StyleBox[ StyleBox[ StyleBox["F", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "[", RowBox[{ RowBox[{"1", "s2", " ", StyleBox[ StyleBox[ StyleBox[ StyleBox[ StyleBox["1", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "sp"}], ",", " ", RowBox[{ StyleBox[ StyleBox[ StyleBox[ StyleBox[ StyleBox[ StyleBox["3", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "p", " ", "LP"}]}], "]"}]}], ",", RowBox[{ RowBox[{"~", "no"}], " ", "overlap"}], ",", " ", RowBox[{ RowBox[{"parametrized", " ", "ES", " ", StyleBox[ RowBox[{"0", StyleBox[ StyleBox[ StyleBox["5", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], ".01"}]], ".2013", "\[IndentingNewLine]", RowBox[{"Variables", ":", "\[IndentingNewLine]", StyleBox[ RowBox[{" ", StyleBox[" ", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]]}]], StyleBox[" ", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], RowBox[{ RowBox[{ StyleBox[ RowBox[{" ", " "}], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], StyleBox[ RowBox[{ StyleBox[" ", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], " "}]], "P", StyleBox[ RowBox[{" ", StyleBox[" ", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]]}]], StyleBox[ StyleBox[ StyleBox[ StyleBox[ StyleBox["F", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]]}], "-", RowBox[{"core", " ", "radius", "\[IndentingNewLine]", "Q", StyleBox[ RowBox[{" ", StyleBox[" ", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]]}]], StyleBox[ StyleBox[ StyleBox[ StyleBox[ StyleBox["F", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]]}], "-", RowBox[{"H", " ", "cloud", " ", "radius", "\[IndentingNewLine]", "R", StyleBox[ RowBox[{" ", StyleBox[" ", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]]}]], "lone", " ", "pair", " ", "radius", "\[IndentingNewLine]", "S", StyleBox[ RowBox[{" ", StyleBox[" ", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]]}]], "distance", " ", "proton", " ", "to", " ", "center", " ", "of", " ", StyleBox[ RowBox[{ StyleBox[ StyleBox[ StyleBox[ StyleBox[ StyleBox[ StyleBox[ StyleBox["F", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "H"}]], " ", "cloud", "\[IndentingNewLine]", "\[Alpha]", " ", "angle", " ", "lonePair"}], "-", StyleBox[ StyleBox[ StyleBox[ StyleBox[ StyleBox[ StyleBox[ StyleBox["F", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "-", RowBox[{"H", StyleBox["\[IndentingNewLine]", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "\[Beta]", " ", "angle", " ", StyleBox[ StyleBox[ RowBox[{ StyleBox[ StyleBox["L", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], StyleBox[ StyleBox[ StyleBox["P", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]]}], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]]}], "-", StyleBox[ StyleBox[ StyleBox[ StyleBox[ StyleBox[ StyleBox[ StyleBox["F", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "-", RowBox[{ StyleBox[ RowBox[{"L", StyleBox["P", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]]}], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], StyleBox[" ", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "not", " ", "a", " ", RowBox[{"variable", ":", " ", "cos\[Beta]"}]}]}]}]}], " ", "=", " ", RowBox[{"1.0", "-", RowBox[{"1.5", "*", RowBox[{ RowBox[{"Sin", "[", "\[Alpha]", "]"}], "^", "2"}]}]}]}]}], " ", StyleBox[ StyleBox[ StyleBox[ StyleBox[ StyleBox["*)", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]]}], "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]]], "Input", CellChangeTimes->{{3.5661428484442425`*^9, 3.5661428783026953`*^9}, { 3.5661430328053665`*^9, 3.5661430369081736`*^9}, {3.5661430998854847`*^9, 3.5661431383537183`*^9}, {3.566143180801393*^9, 3.5661432248645825`*^9}, { 3.566143288683233*^9, 3.5661432893462706`*^9}, 3.566143344240477*^9, 3.566143438043442*^9, 3.566143491551536*^9, {3.5661922351554785`*^9, 3.5661923427020674`*^9}, {3.566201924306572*^9, 3.566201971387455*^9}, 3.566202322310071*^9, {3.5663083624775496`*^9, 3.5663083893563967`*^9}, { 3.5663654716175947`*^9, 3.566365471929596*^9}, {3.566381164724308*^9, 3.5663812477008533`*^9}, {3.566381498674094*^9, 3.5663815066145086`*^9}, { 3.5663815503101854`*^9, 3.5663816389963408`*^9}, {3.566381697512044*^9, 3.566381781892592*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"Clear", "[", RowBox[{ "wc", ",", "k1", ",", "k2", ",", "k3", ",", "s1", ",", "s2", ",", "s3", ",", "par"}], "]"}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Z", "=", "9.0"}], ";", " ", RowBox[{"w3", "=", RowBox[{"Sqrt", "[", "3", "]"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", "core", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Tc", "=", RowBox[{"2.25", "*", RowBox[{"k1", "/", RowBox[{"P", "^", "2"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vnec", "=", RowBox[{ RowBox[{"-", "3.0"}], "*", RowBox[{"Z", "/", "P"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Veec", "=", RowBox[{"3.0", "*", RowBox[{"s1", "/", "P"}]}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"XH", " ", "bond"}], "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Tb", "=", RowBox[{"2.25", "*", RowBox[{"k2", "/", RowBox[{"Q", "^", "2"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vneb", "=", RowBox[{ RowBox[{"-", "2.0"}], "*", RowBox[{ RowBox[{"(", RowBox[{"1.5", "-", RowBox[{"0.5", "*", RowBox[{ RowBox[{"(", RowBox[{"S", "/", "Q"}], ")"}], "^", "2"}]}]}], ")"}], "/", "Q"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Veeb", "=", RowBox[{"3.0", "*", RowBox[{"s2", "/", "Q"}]}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"lone", " ", "pair"}], "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Tp", "=", RowBox[{"2.25", "*", RowBox[{"k3", "/", RowBox[{"R", "^", "2"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Veep", "=", RowBox[{"3.0", "*", RowBox[{"s3", "/", "R"}]}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"core", "-", "bond"}], "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vnecb", "=", RowBox[{ RowBox[{ RowBox[{"-", "2.0"}], "*", RowBox[{"Z", "/", RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}]}]}], "-", RowBox[{"2", "/", RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Veecb", "=", RowBox[{"4.0", "/", RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Vnncb", "=", RowBox[{"Z", "/", RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}]}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"core", "-", RowBox[{"lone", " ", "pair"}]}], "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vnecp", "=", RowBox[{ RowBox[{"-", "2.0"}], "*", RowBox[{"Z", "/", RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}]}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Veecp", "=", RowBox[{"4.0", "/", RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}]}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"bond", "-", RowBox[{"lone", " ", "pair"}]}], "*)"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ RowBox[{"d2", " ", "w", RowBox[{"(", "FH", ")"}]}], "-", RowBox[{"w", RowBox[{"(", "LP", ")"}]}]}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"d2", "=", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}], "^", "2"}], "+", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}], "^", "2"}], "-", RowBox[{"2", "*", RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}], "*", RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}], "*", RowBox[{"Cos", "[", "\[Alpha]", "]"}]}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Veebp", " ", "=", " ", RowBox[{"4", "/", "d2"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ RowBox[{"d4", " ", "w", RowBox[{"(", "LP", ")"}]}], "-", "H"}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"d4", "=", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}], "^", "2"}], "+", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}], "^", "2"}], "-", RowBox[{"2", "*", RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}], "*", RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}], "*", RowBox[{"Cos", "[", "\[Alpha]", "]"}]}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Vnebp", " ", "=", RowBox[{ RowBox[{"-", "2.0"}], "/", "d4"}]}], " ", ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", StyleBox[ RowBox[{"(*", RowBox[{ RowBox[{"lone", " ", "pair"}], "-", RowBox[{"lone", " ", "pair"}]}], "*)"}], "Text", FontFamily->"Arial Narrow"], StyleBox["\[IndentingNewLine]", FontFamily->"Arial Narrow"], StyleBox[ RowBox[{"(*", " ", RowBox[{ RowBox[{"none", " ", "for", " ", "NH3"}], ",", " ", RowBox[{"but", " ", "for", " ", "H2O", " ", "and", " ", "HF"}]}], " ", "*)"}], "Text", FontFamily->"Arial Narrow"], "\[IndentingNewLine]", StyleBox[ RowBox[{"(*", " ", RowBox[{ StyleBox[ RowBox[{"H", StyleBox[ StyleBox["F", "Text", FontFamily->"Arial Narrow"], "Text", FontFamily->"Arial Narrow"]}]], " ", "has", " ", StyleBox[ RowBox[{ StyleBox["o", "Text", FontFamily->"Arial Narrow"], "ne"}]], StyleBox[" ", "Text", FontFamily->"Arial Narrow"], "angle", StyleBox[" ", "Text", FontFamily->"Arial Narrow"], "as", " ", RowBox[{"variable", ":", " ", RowBox[{"H", "-", StyleBox[ StyleBox[ StyleBox[ StyleBox["F", "Text", FontFamily->"Arial Narrow"], "Text", FontFamily->"Arial Narrow"], "Text", FontFamily->"Arial Narrow"], "Text", FontFamily->"Arial Narrow"], StyleBox["-", "Text", FontFamily->"Arial Narrow"], StyleBox[ StyleBox[ RowBox[{ StyleBox[ StyleBox[ StyleBox["L", "Text", FontFamily->"Arial Narrow"], "Text", FontFamily->"Arial Narrow"], "Text", FontFamily->"Arial Narrow"], StyleBox[ StyleBox[ StyleBox[ StyleBox["P", "Text", FontFamily->"Arial Narrow"], "Text", FontFamily->"Arial Narrow"], "Text", FontFamily->"Arial Narrow"], "Text", FontFamily->"Arial Narrow"]}], "Text", FontFamily->"Arial Narrow"], "Text", FontFamily->"Arial Narrow"]}]}]}], " ", "*)"}], "Text", FontFamily->"Arial Narrow"], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ RowBox[{"d1", " ", "w", RowBox[{"(", "LP", ")"}]}], "-", RowBox[{"w", RowBox[{"(", "LP", ")"}]}]}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"d1", "=", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}], "*", RowBox[{"Sin", "[", "\[Alpha]", "]"}], "*", "w3"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Veepp", " ", "=", " ", RowBox[{"4.0", "/", "d1"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"(*", RowBox[{"sum", " ", "terms"}], "*)"}]}]}], "Input", CellChangeTimes->{{3.5661428484442425`*^9, 3.5661428783026953`*^9}, { 3.5661430328053665`*^9, 3.5661430369081736`*^9}, {3.5661430998854847`*^9, 3.5661431383537183`*^9}, {3.566143180801393*^9, 3.5661432248645825`*^9}, { 3.566143288683233*^9, 3.5661432893462706`*^9}, 3.566143344240477*^9, 3.566143438043442*^9, {3.5661923476784763`*^9, 3.566192348645678*^9}, { 3.566192420296604*^9, 3.566192420858205*^9}, {3.5661924684694886`*^9, 3.566192482369113*^9}, {3.5661925519608355`*^9, 3.5661925980901165`*^9}, { 3.566192662409029*^9, 3.5661926728298473`*^9}, {3.566194733031866*^9, 3.5661947591463118`*^9}, {3.566201172322851*^9, 3.5662012049425087`*^9}, { 3.566201275782233*^9, 3.5662012764842343`*^9}, {3.566201607563616*^9, 3.5662016091548185`*^9}, {3.5662016860473537`*^9, 3.5662017455770583`*^9}, {3.566201840690425*^9, 3.5662018510176435`*^9}, { 3.5662021072324934`*^9, 3.56620213362774*^9}, {3.566202199116655*^9, 3.566202199381855*^9}, {3.5663049674119854`*^9, 3.5663050328385005`*^9}, { 3.566305081167385*^9, 3.566305085270192*^9}, {3.56630512948067*^9, 3.566305132819076*^9}, {3.5663056216303344`*^9, 3.56630564170757*^9}, 3.566305710550491*^9, {3.5663057513601627`*^9, 3.566305753419366*^9}, { 3.566307371173808*^9, 3.5663074141830835`*^9}, 3.566307457207959*^9, 3.5663085639675035`*^9, {3.5663156573399625`*^9, 3.5663156653739767`*^9}, { 3.5663157659941535`*^9, 3.5663157671173553`*^9}, {3.5663812557036676`*^9, 3.5663812564368687`*^9}, {3.5663813096797624`*^9, 3.566381328181395*^9}, { 3.5663813738582754`*^9, 3.566381461171629*^9}, {3.5663818620767326`*^9, 3.5663818862099752`*^9}, {3.566381924929243*^9, 3.566381974958531*^9}, { 3.5663820659690905`*^9, 3.56638212237879*^9}, {3.566382155014047*^9, 3.5663821831876965`*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"T", "=", RowBox[{"Tc", "+", " ", "Tb", "+", " ", RowBox[{"3", "*", "Tp"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vne", "=", RowBox[{"Vnec", "+", "Vneb", "+", " ", "Vnecb", "+", RowBox[{"3", "*", "Vnecp"}], "+", " ", RowBox[{"3", "*", "Vnebp"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vee", "=", RowBox[{"Veec", "+", "Veeb", "+", " ", RowBox[{"3", "*", "Veep"}], "+", "Veecb", "+", RowBox[{"3", "*", "Veecp"}], "+", " ", RowBox[{"3", "*", "Veebp"}], " ", "+", RowBox[{"3", "*", "Veepp"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ "Vnn", "=", " ", "Vnncb"}], ";"}], "\[IndentingNewLine]"}], "Input", CellChangeTimes->{{3.5661428484442425`*^9, 3.5661428783026953`*^9}, { 3.5661430328053665`*^9, 3.5661430369081736`*^9}, {3.5661430998854847`*^9, 3.5661431383537183`*^9}, {3.566143180801393*^9, 3.5661432248645825`*^9}, { 3.566143288683233*^9, 3.5661432893462706`*^9}, 3.566143344240477*^9, { 3.5661435118159714`*^9, 3.5661435187891836`*^9}, {3.5661470921646585`*^9, 3.5661471085134873`*^9}, {3.566194770659132*^9, 3.566194777023943*^9}, { 3.566200993811738*^9, 3.566201133478783*^9}, {3.5662012401361704`*^9, 3.5662012446913786`*^9}, {3.5662014429677267`*^9, 3.566201443654128*^9}, { 3.566201792735941*^9, 3.5662017938591433`*^9}, 3.5662022754943886`*^9, { 3.5662025519112744`*^9, 3.566202565436498*^9}, {3.566305456238844*^9, 3.566305473289674*^9}, {3.5663055060497313`*^9, 3.5663055070949335`*^9}, { 3.566308222654504*^9, 3.5663082634017754`*^9}, {3.566382191440111*^9, 3.5663823356935644`*^9}, {3.5663823743660326`*^9, 3.566382374693633*^9}, { 3.5663824375617437`*^9, 3.5663824722562046`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ RowBox[{"Parameters", ":", " ", "G3"}], "//", RowBox[{ RowBox[{ RowBox[{ RowBox[{"6", "-", RowBox[{"311", "G", RowBox[{"(", RowBox[{"d", ",", "p"}], ")"}], " ", "can", " ", "be", " ", "read", " ", "for", " ", "LiH"}]}], " ", "..."}], " ", "HF", "\[IndentingNewLine]", "in", " ", "table", " ", "from", " ", RowBox[{"www", ".", "kimball"}]}], "-", RowBox[{ RowBox[{ RowBox[{"model", ".", "org"}], "/", "Kimball"}], "/", RowBox[{"params", ".", "xls"}]}]}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"par", "=", RowBox[{"{", RowBox[{ RowBox[{"k1", "\[Rule]", "0.91058613"}], ",", RowBox[{"k2", "\[Rule]", "1.41656890"}], ",", RowBox[{"k3", "\[Rule]", "1.21497466"}], ",", RowBox[{"s1", "\[Rule]", "0.51017912"}], ",", RowBox[{"s2", "\[Rule]", "0.44166435"}], ",", RowBox[{"s3", "\[Rule]", "0.36888374"}]}], "}"}]}], ";"}]}]], "Input", CellChangeTimes->{{3.5661428484442425`*^9, 3.5661428783026953`*^9}, { 3.5661430328053665`*^9, 3.5661430369081736`*^9}, {3.5661430998854847`*^9, 3.5661431383537183`*^9}, {3.566143180801393*^9, 3.5661432248645825`*^9}, { 3.566143288683233*^9, 3.5661432893462706`*^9}, {3.566143344240477*^9, 3.566143371150524*^9}, {3.5663152612084665`*^9, 3.5663152707400837`*^9}, { 3.566315729474489*^9, 3.566315733312096*^9}, {3.566315805368623*^9, 3.5663158109066324`*^9}, {3.566316505903453*^9, 3.566316526277089*^9}, { 3.566382621319267*^9, 3.5663826949981966`*^9}}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"func", "=", RowBox[{ RowBox[{"T", "+", "Vne", "+", "Vee", "+", "Vnn"}], "/.", "par"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"t", "=", RowBox[{"FindMinimum", "[", RowBox[{"func", ",", RowBox[{"{", RowBox[{"P", ",", "0.16"}], "}"}], ",", RowBox[{"{", RowBox[{"Q", ",", "1.0"}], "}"}], ",", RowBox[{"{", RowBox[{"R", ",", " ", "0.83"}], "}"}], ",", RowBox[{"{", RowBox[{"S", ",", "0.6"}], "}"}], ",", " ", RowBox[{"{", RowBox[{"\[Alpha]", ",", "1.82"}], "}"}], ",", RowBox[{"MaxIterations", "\[Rule]", "500"}], ",", " ", RowBox[{"Method", " ", "\[Rule]", " ", "\"\\""}]}], "]"}]}]}], "Input", CellChangeTimes->{{3.5661428484442425`*^9, 3.5661428783026953`*^9}, { 3.5661430328053665`*^9, 3.5661430369081736`*^9}, {3.5661430998854847`*^9, 3.5661431383537183`*^9}, {3.566143180801393*^9, 3.5661432248645825`*^9}, { 3.566143288683233*^9, 3.5661432893462706`*^9}, {3.566143344240477*^9, 3.5661433819769435`*^9}, {3.5662024426798825`*^9, 3.5662024734587364`*^9}, {3.566305156858718*^9, 3.566305166265535*^9}, { 3.5663079842392845`*^9, 3.566307984972486*^9}, 3.5663654313383245`*^9, 3.5663824984954505`*^9, {3.5663827052318144`*^9, 3.56638270819582*^9}, { 3.566385995585103*^9, 3.5663860074879236`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"-", "100.39980413479705`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.15508308756218367`"}], ",", RowBox[{"Q", "\[Rule]", "0.9829704623880443`"}], ",", RowBox[{"R", "\[Rule]", "0.829677130509701`"}], ",", RowBox[{"S", "\[Rule]", "0.5833733146453088`"}], ",", RowBox[{"\[Alpha]", "\[Rule]", "1.8164601591901128`"}]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{ 3.5663165372283077`*^9, 3.5663175724229746`*^9, 3.5663176162590513`*^9, 3.566317659330727*^9, 3.566320373812928*^9, 3.566320436821439*^9, 3.566320474401905*^9, 3.566321519556941*^9, 3.5663217066792693`*^9, 3.5663218529451265`*^9, 3.5663626992899256`*^9, 3.566362965987994*^9, 3.566363116185058*^9, 3.5663631564331284`*^9, 3.566363452786849*^9, 3.5663635926410947`*^9, 3.5663642788550997`*^9, 3.566364436430977*^9, { 3.5663645631343994`*^9, 3.5663645911676483`*^9}, 3.566365060759673*^9, 3.5663651610678496`*^9, 3.566365211206338*^9, 3.566365290938078*^9, { 3.5663653439157705`*^9, 3.5663653720114202`*^9}, 3.5663827425158796`*^9, 3.566383042399929*^9, 3.566383724812604*^9, 3.5663841404922447`*^9, 3.566384221284786*^9, 3.5663848937707677`*^9, 3.566385052329446*^9, 3.5663851086767454`*^9, 3.5663851390655985`*^9, 3.5663854423925314`*^9, 3.566385496836627*^9, 3.5663856065204196`*^9, 3.5663856438824854`*^9, 3.5663857054869933`*^9, 3.566385797480355*^9, 3.566385928317785*^9, 3.56638601658274*^9, 3.5663861224757257`*^9, 3.5663862556219597`*^9, 3.5663870661444235`*^9, 3.5663871844238315`*^9, 3.5663872179950905`*^9, 3.5663872767915936`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"u", "=", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"Vne", "/.", "u"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vee", "/.", "u"}], "/.", "par"}], "\[IndentingNewLine]", RowBox[{"Vnn", "/.", "u"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"Vne", "+", "Vee", "+", "Vnn"}], ")"}], "/", "T"}], "/.", "u"}], "/.", "par"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"P", "+", "Q", "+", "S"}], ")"}], "*", "0.529177"}], "/.", "u"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ RowBox[{"angle", " ", "LP"}], "-", "FH"}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"180", "/", "\[Pi]"}], "*", "\[Alpha]"}], " ", "/.", " ", "u"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{ RowBox[{"angle", " ", "LP"}], "-", "F", "-", "LP"}], " ", "*)"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"180", "/", "\[Pi]"}], "*", RowBox[{"ArcCos", "[", RowBox[{"1.0", "-", RowBox[{"1.5", "*", RowBox[{ RowBox[{"Sin", "[", "\[Alpha]", "]"}], "^", "2"}]}]}], "]"}]}], " ", "/.", " ", "u"}], "\[IndentingNewLine]", RowBox[{"d1", " ", "/.", "u"}], "\[IndentingNewLine]", RowBox[{"d2", " ", "/.", "u"}], "\[IndentingNewLine]", RowBox[{"d4", " ", "/.", "u"}]}], "Input", CellChangeTimes->{{3.5661428484442425`*^9, 3.5661428783026953`*^9}, { 3.5661430328053665`*^9, 3.5661430369081736`*^9}, {3.5661430998854847`*^9, 3.5661431383537183`*^9}, {3.566143180801393*^9, 3.5661432248645825`*^9}, { 3.566143288683233*^9, 3.5661432893462706`*^9}, {3.566143344240477*^9, 3.5661433961729684`*^9}, {3.5663078436050377`*^9, 3.5663078735570903`*^9}, 3.5663654431007447`*^9, {3.5663830104354725`*^9, 3.566383031199109*^9}, { 3.5663836080676546`*^9, 3.566383713471384*^9}}], Cell[BoxData[ RowBox[{"-", "251.3586523772973`"}]], "Output", CellChangeTimes->{ 3.566316537243908*^9, 3.566317572438575*^9, 3.5663176162590513`*^9, 3.5663176593463273`*^9, 3.5663203738285284`*^9, 3.566320436837039*^9, 3.5663204744175053`*^9, 3.5663215195725408`*^9, 3.5663217066948695`*^9, 3.5663218529607267`*^9, 3.566362699305526*^9, 3.566362966003594*^9, 3.566363116185058*^9, 3.5663631564487286`*^9, 3.566363452802449*^9, 3.566363592656695*^9, 3.5663642788707*^9, 3.566364436446577*^9, { 3.5663645631343994`*^9, 3.5663645911676483`*^9}, 3.566365060759673*^9, 3.5663651610678496`*^9, 3.5663652112219377`*^9, 3.5663652909536777`*^9, { 3.5663653439157705`*^9, 3.5663653720114202`*^9}, 3.56638274253148*^9, 3.5663830424155293`*^9, 3.5663837248282037`*^9, 3.566384140507845*^9, 3.5663842213003864`*^9, 3.5663848937863674`*^9, 3.566385052345046*^9, 3.566385108692345*^9, 3.5663851390811987`*^9, 3.566385442408131*^9, 3.5663854968522267`*^9, 3.5663856065360193`*^9, 3.566385643898085*^9, 3.5663857055025935`*^9, 3.566385797480355*^9, 3.5663859283333845`*^9, 3.56638601659834*^9, 3.566386122491326*^9, 3.56638625563756*^9, 3.5663870661444235`*^9, 3.566387184439431*^9, 3.5663872180106907`*^9, 3.5663872768071938`*^9}], Cell[BoxData["45.33082269224218`"], "Output", CellChangeTimes->{ 3.566316537243908*^9, 3.566317572438575*^9, 3.5663176162590513`*^9, 3.5663176593463273`*^9, 3.5663203738285284`*^9, 3.566320436837039*^9, 3.5663204744175053`*^9, 3.5663215195725408`*^9, 3.5663217066948695`*^9, 3.5663218529607267`*^9, 3.566362699305526*^9, 3.566362966003594*^9, 3.566363116185058*^9, 3.5663631564487286`*^9, 3.566363452802449*^9, 3.566363592656695*^9, 3.5663642788707*^9, 3.566364436446577*^9, { 3.5663645631343994`*^9, 3.5663645911676483`*^9}, 3.566365060759673*^9, 3.5663651610678496`*^9, 3.5663652112219377`*^9, 3.5663652909536777`*^9, { 3.5663653439157705`*^9, 3.5663653720114202`*^9}, 3.56638274253148*^9, 3.5663830424155293`*^9, 3.5663837248282037`*^9, 3.566384140507845*^9, 3.5663842213003864`*^9, 3.5663848937863674`*^9, 3.566385052345046*^9, 3.566385108692345*^9, 3.5663851390811987`*^9, 3.566385442408131*^9, 3.5663854968522267`*^9, 3.5663856065360193`*^9, 3.566385643898085*^9, 3.5663857055025935`*^9, 3.566385797480355*^9, 3.5663859283333845`*^9, 3.56638601659834*^9, 3.566386122491326*^9, 3.56638625563756*^9, 3.5663870661444235`*^9, 3.566387184439431*^9, 3.5663872180106907`*^9, 3.5663872768071938`*^9}], Cell[BoxData["5.228220951527105`"], "Output", CellChangeTimes->{ 3.566316537243908*^9, 3.566317572438575*^9, 3.5663176162590513`*^9, 3.5663176593463273`*^9, 3.5663203738285284`*^9, 3.566320436837039*^9, 3.5663204744175053`*^9, 3.5663215195725408`*^9, 3.5663217066948695`*^9, 3.5663218529607267`*^9, 3.566362699305526*^9, 3.566362966003594*^9, 3.566363116185058*^9, 3.5663631564487286`*^9, 3.566363452802449*^9, 3.566363592656695*^9, 3.5663642788707*^9, 3.566364436446577*^9, { 3.5663645631343994`*^9, 3.5663645911676483`*^9}, 3.566365060759673*^9, 3.5663651610678496`*^9, 3.5663652112219377`*^9, 3.5663652909536777`*^9, { 3.5663653439157705`*^9, 3.5663653720114202`*^9}, 3.56638274253148*^9, 3.5663830424155293`*^9, 3.5663837248282037`*^9, 3.566384140507845*^9, 3.5663842213003864`*^9, 3.5663848937863674`*^9, 3.566385052345046*^9, 3.566385108692345*^9, 3.5663851390811987`*^9, 3.566385442408131*^9, 3.5663854968522267`*^9, 3.5663856065360193`*^9, 3.566385643898085*^9, 3.5663857055025935`*^9, 3.566385797480355*^9, 3.5663859283333845`*^9, 3.56638601659834*^9, 3.566386122491326*^9, 3.56638625563756*^9, 3.5663870661444235`*^9, 3.566387184439431*^9, 3.5663872180106907`*^9, 3.5663872768071938`*^9}], Cell[BoxData[ RowBox[{"-", "1.999999995379135`"}]], "Output", CellChangeTimes->{ 3.566316537243908*^9, 3.566317572438575*^9, 3.5663176162590513`*^9, 3.5663176593463273`*^9, 3.5663203738285284`*^9, 3.566320436837039*^9, 3.5663204744175053`*^9, 3.5663215195725408`*^9, 3.5663217066948695`*^9, 3.5663218529607267`*^9, 3.566362699305526*^9, 3.566362966003594*^9, 3.566363116185058*^9, 3.5663631564487286`*^9, 3.566363452802449*^9, 3.566363592656695*^9, 3.5663642788707*^9, 3.566364436446577*^9, { 3.5663645631343994`*^9, 3.5663645911676483`*^9}, 3.566365060759673*^9, 3.5663651610678496`*^9, 3.5663652112219377`*^9, 3.5663652909536777`*^9, { 3.5663653439157705`*^9, 3.5663653720114202`*^9}, 3.56638274253148*^9, 3.5663830424155293`*^9, 3.5663837248282037`*^9, 3.566384140507845*^9, 3.5663842213003864`*^9, 3.5663848937863674`*^9, 3.566385052345046*^9, 3.566385108692345*^9, 3.5663851390811987`*^9, 3.566385442408131*^9, 3.5663854968522267`*^9, 3.5663856065360193`*^9, 3.566385643898085*^9, 3.5663857055025935`*^9, 3.566385797480355*^9, 3.5663859283333845`*^9, 3.56638601659834*^9, 3.566386122491326*^9, 3.56638625563756*^9, 3.5663870661444235`*^9, 3.566387184439431*^9, 3.5663872180106907`*^9, 3.5663872768071938`*^9}], Cell[BoxData["0.9109395039260724`"], "Output", CellChangeTimes->{ 3.566316537243908*^9, 3.566317572438575*^9, 3.5663176162590513`*^9, 3.5663176593463273`*^9, 3.5663203738285284`*^9, 3.566320436837039*^9, 3.5663204744175053`*^9, 3.5663215195725408`*^9, 3.5663217066948695`*^9, 3.5663218529607267`*^9, 3.566362699305526*^9, 3.566362966003594*^9, 3.566363116185058*^9, 3.5663631564487286`*^9, 3.566363452802449*^9, 3.566363592656695*^9, 3.5663642788707*^9, 3.566364436446577*^9, { 3.5663645631343994`*^9, 3.5663645911676483`*^9}, 3.566365060759673*^9, 3.5663651610678496`*^9, 3.5663652112219377`*^9, 3.5663652909536777`*^9, { 3.5663653439157705`*^9, 3.5663653720114202`*^9}, 3.56638274253148*^9, 3.5663830424155293`*^9, 3.5663837248282037`*^9, 3.566384140507845*^9, 3.5663842213003864`*^9, 3.5663848937863674`*^9, 3.566385052345046*^9, 3.566385108692345*^9, 3.5663851390811987`*^9, 3.566385442408131*^9, 3.5663854968522267`*^9, 3.5663856065360193`*^9, 3.566385643898085*^9, 3.5663857055025935`*^9, 3.566385797480355*^9, 3.5663859283333845`*^9, 3.56638601659834*^9, 3.566386122491326*^9, 3.56638625563756*^9, 3.5663870661444235`*^9, 3.566387184439431*^9, 3.5663872180106907`*^9, 3.5663872768071938`*^9}], Cell[BoxData["104.07550077525512`"], "Output", CellChangeTimes->{ 3.566316537243908*^9, 3.566317572438575*^9, 3.5663176162590513`*^9, 3.5663176593463273`*^9, 3.5663203738285284`*^9, 3.566320436837039*^9, 3.5663204744175053`*^9, 3.5663215195725408`*^9, 3.5663217066948695`*^9, 3.5663218529607267`*^9, 3.566362699305526*^9, 3.566362966003594*^9, 3.566363116185058*^9, 3.5663631564487286`*^9, 3.566363452802449*^9, 3.566363592656695*^9, 3.5663642788707*^9, 3.566364436446577*^9, { 3.5663645631343994`*^9, 3.5663645911676483`*^9}, 3.566365060759673*^9, 3.5663651610678496`*^9, 3.5663652112219377`*^9, 3.5663652909536777`*^9, { 3.5663653439157705`*^9, 3.5663653720114202`*^9}, 3.56638274253148*^9, 3.5663830424155293`*^9, 3.5663837248282037`*^9, 3.566384140507845*^9, 3.5663842213003864`*^9, 3.5663848937863674`*^9, 3.566385052345046*^9, 3.566385108692345*^9, 3.5663851390811987`*^9, 3.566385442408131*^9, 3.5663854968522267`*^9, 3.5663856065360193`*^9, 3.566385643898085*^9, 3.5663857055025935`*^9, 3.566385797480355*^9, 3.5663859283333845`*^9, 3.56638601659834*^9, 3.566386122491326*^9, 3.56638625563756*^9, 3.5663870661444235`*^9, 3.566387184439431*^9, 3.5663872180106907`*^9, 3.5663872768071938`*^9}], Cell[BoxData["114.28529518872062`"], "Output", CellChangeTimes->{ 3.566316537243908*^9, 3.566317572438575*^9, 3.5663176162590513`*^9, 3.5663176593463273`*^9, 3.5663203738285284`*^9, 3.566320436837039*^9, 3.5663204744175053`*^9, 3.5663215195725408`*^9, 3.5663217066948695`*^9, 3.5663218529607267`*^9, 3.566362699305526*^9, 3.566362966003594*^9, 3.566363116185058*^9, 3.5663631564487286`*^9, 3.566363452802449*^9, 3.566363592656695*^9, 3.5663642788707*^9, 3.566364436446577*^9, { 3.5663645631343994`*^9, 3.5663645911676483`*^9}, 3.566365060759673*^9, 3.5663651610678496`*^9, 3.5663652112219377`*^9, 3.5663652909536777`*^9, { 3.5663653439157705`*^9, 3.5663653720114202`*^9}, 3.56638274253148*^9, 3.5663830424155293`*^9, 3.5663837248282037`*^9, 3.566384140507845*^9, 3.5663842213003864`*^9, 3.5663848937863674`*^9, 3.566385052345046*^9, 3.566385108692345*^9, 3.5663851390811987`*^9, 3.566385442408131*^9, 3.5663854968522267`*^9, 3.5663856065360193`*^9, 3.566385643898085*^9, 3.5663857055025935`*^9, 3.566385797480355*^9, 3.5663859283333845`*^9, 3.56638601659834*^9, 3.566386122491326*^9, 3.56638625563756*^9, 3.5663870661444235`*^9, 3.566387184439431*^9, 3.5663872180106907`*^9, 3.5663872768071938`*^9}], Cell[BoxData["1.6544443144832426`"], "Output", CellChangeTimes->{ 3.566316537243908*^9, 3.566317572438575*^9, 3.5663176162590513`*^9, 3.5663176593463273`*^9, 3.5663203738285284`*^9, 3.566320436837039*^9, 3.5663204744175053`*^9, 3.5663215195725408`*^9, 3.5663217066948695`*^9, 3.5663218529607267`*^9, 3.566362699305526*^9, 3.566362966003594*^9, 3.566363116185058*^9, 3.5663631564487286`*^9, 3.566363452802449*^9, 3.566363592656695*^9, 3.5663642788707*^9, 3.566364436446577*^9, { 3.5663645631343994`*^9, 3.5663645911676483`*^9}, 3.566365060759673*^9, 3.5663651610678496`*^9, 3.5663652112219377`*^9, 3.5663652909536777`*^9, { 3.5663653439157705`*^9, 3.5663653720114202`*^9}, 3.56638274253148*^9, 3.5663830424155293`*^9, 3.5663837248282037`*^9, 3.566384140507845*^9, 3.5663842213003864`*^9, 3.5663848937863674`*^9, 3.566385052345046*^9, 3.566385108692345*^9, 3.5663851390811987`*^9, 3.566385442408131*^9, 3.5663854968522267`*^9, 3.5663856065360193`*^9, 3.566385643898085*^9, 3.5663857055025935`*^9, 3.566385797480355*^9, 3.5663859283333845`*^9, 3.56638601659834*^9, 3.566386122491326*^9, 3.56638625563756*^9, 3.5663870661444235`*^9, 3.566387184439431*^9, 3.5663872180106907`*^9, 3.5663872768071938`*^9}], Cell[BoxData["1.676315147540139`"], "Output", CellChangeTimes->{ 3.566316537243908*^9, 3.566317572438575*^9, 3.5663176162590513`*^9, 3.5663176593463273`*^9, 3.5663203738285284`*^9, 3.566320436837039*^9, 3.5663204744175053`*^9, 3.5663215195725408`*^9, 3.5663217066948695`*^9, 3.5663218529607267`*^9, 3.566362699305526*^9, 3.566362966003594*^9, 3.566363116185058*^9, 3.5663631564487286`*^9, 3.566363452802449*^9, 3.566363592656695*^9, 3.5663642788707*^9, 3.566364436446577*^9, { 3.5663645631343994`*^9, 3.5663645911676483`*^9}, 3.566365060759673*^9, 3.5663651610678496`*^9, 3.5663652112219377`*^9, 3.5663652909536777`*^9, { 3.5663653439157705`*^9, 3.5663653720114202`*^9}, 3.56638274253148*^9, 3.5663830424155293`*^9, 3.5663837248282037`*^9, 3.566384140507845*^9, 3.5663842213003864`*^9, 3.5663848937863674`*^9, 3.566385052345046*^9, 3.566385108692345*^9, 3.5663851390811987`*^9, 3.566385442408131*^9, 3.5663854968522267`*^9, 3.5663856065360193`*^9, 3.566385643898085*^9, 3.5663857055025935`*^9, 3.566385797480355*^9, 3.5663859283333845`*^9, 3.56638601659834*^9, 3.566386122491326*^9, 3.56638625563756*^9, 3.5663870661444235`*^9, 3.566387184439431*^9, 3.5663872180106907`*^9, 3.5663872768071938`*^9}], Cell[BoxData["2.1811936669618626`"], "Output", CellChangeTimes->{ 3.566316537243908*^9, 3.566317572438575*^9, 3.5663176162590513`*^9, 3.5663176593463273`*^9, 3.5663203738285284`*^9, 3.566320436837039*^9, 3.5663204744175053`*^9, 3.5663215195725408`*^9, 3.5663217066948695`*^9, 3.5663218529607267`*^9, 3.566362699305526*^9, 3.566362966003594*^9, 3.566363116185058*^9, 3.5663631564487286`*^9, 3.566363452802449*^9, 3.566363592656695*^9, 3.5663642788707*^9, 3.566364436446577*^9, { 3.5663645631343994`*^9, 3.5663645911676483`*^9}, 3.566365060759673*^9, 3.5663651610678496`*^9, 3.5663652112219377`*^9, 3.5663652909536777`*^9, { 3.5663653439157705`*^9, 3.5663653720114202`*^9}, 3.56638274253148*^9, 3.5663830424155293`*^9, 3.5663837248282037`*^9, 3.566384140507845*^9, 3.5663842213003864`*^9, 3.5663848937863674`*^9, 3.566385052345046*^9, 3.566385108692345*^9, 3.5663851390811987`*^9, 3.566385442408131*^9, 3.5663854968522267`*^9, 3.5663856065360193`*^9, 3.566385643898085*^9, 3.5663857055025935`*^9, 3.566385797480355*^9, 3.5663859283333845`*^9, 3.56638601659834*^9, 3.566386122491326*^9, 3.56638625563756*^9, 3.5663870661444235`*^9, 3.566387184439431*^9, 3.5663872180106907`*^9, 3.5663872768071938`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"wd", " ", "=", " ", RowBox[{ RowBox[{"Sin", "[", "\[Alpha]", "]"}], " ", "/.", " ", "u"}]}], ";", " ", RowBox[{"wc", "=", RowBox[{ RowBox[{"Cos", "[", "\[Alpha]", "]"}], " ", "/.", " ", "u"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"plot3", " ", "=", " ", RowBox[{ RowBox[{"Graphics3D", "[", RowBox[{"{", RowBox[{ RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", " ", "P"}], "]"}], ",", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", " ", "0.1"}], "]"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Opacity", "[", "0.5", "]"}], ",", " ", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0", ",", RowBox[{"P", "+", "Q"}]}], "}"}], ",", " ", "Q"}], "]"}]}], "}"}], ",", " ", "\[IndentingNewLine]", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "wd"}], "*", RowBox[{ RowBox[{"(", RowBox[{"P", " ", "+", " ", "R"}], ")"}], "/", "2"}]}], ",", " ", RowBox[{"wd", "*", "w3", "*", RowBox[{ RowBox[{"(", RowBox[{"P", " ", "+", " ", "R"}], ")"}], "/", "2"}]}], ",", RowBox[{"wc", "*", RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}]}]}], "}"}], ",", " ", "R"}], "]"}], ",", " ", "\[IndentingNewLine]", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "wd"}], "*", RowBox[{ RowBox[{"(", RowBox[{"P", " ", "+", " ", "R"}], ")"}], "/", "2"}]}], ",", " ", RowBox[{ RowBox[{"-", "wd"}], "*", "w3", "*", RowBox[{ RowBox[{"(", RowBox[{"P", " ", "+", " ", "R"}], ")"}], "/", "2"}]}], ",", RowBox[{"wc", "*", RowBox[{"(", RowBox[{"P", "+", "Q"}], ")"}]}]}], "}"}], ",", " ", "R"}], "]"}], ",", " ", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"wd", "*", RowBox[{"(", RowBox[{"P", " ", "+", " ", "R"}], ")"}]}], ",", " ", "0", ",", RowBox[{"wc", "*", RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}]}]}], "}"}], ",", " ", "R"}], "]"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Darker", "[", "Red", "]"}], ",", "\[IndentingNewLine]", RowBox[{"Sphere", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0", ",", RowBox[{"P", "+", "Q", "+", "S"}]}], "}"}], ",", " ", "0.06"}], "]"}]}], "}"}]}], "}"}], "]"}], " ", "/.", " ", "u"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"Show", "[", RowBox[{"plot3", ",", RowBox[{"{", RowBox[{ RowBox[{"AspectRatio", " ", "->", " ", "Automatic"}], ",", " ", RowBox[{"Boxed", " ", "\[Rule]", " ", "False"}]}], "}"}]}], "]"}], "\[IndentingNewLine]"}], "Input", CellChangeTimes->{{3.5661428484442425`*^9, 3.5661428783026953`*^9}, { 3.5661430328053665`*^9, 3.5661430369081736`*^9}, {3.5661430998854847`*^9, 3.5661431383537183`*^9}, {3.566143180801393*^9, 3.5661432248645825`*^9}, { 3.566143288683233*^9, 3.5661432893462706`*^9}, {3.566143344240477*^9, 3.566143410119393*^9}, 3.566308742431816*^9, {3.5663176430442986`*^9, 3.566317682730768*^9}, {3.5663204234522157`*^9, 3.566320431267829*^9}, { 3.5663204642618876`*^9, 3.5663204702366977`*^9}, {3.566362812374524*^9, 3.5663628459301834`*^9}, {3.5663629092506943`*^9, 3.566362925334323*^9}, { 3.5663630995242286`*^9, 3.5663631074646425`*^9}, {3.566363142408704*^9, 3.5663631503023176`*^9}, {3.566363335505843*^9, 3.5663634250188007`*^9}, { 3.5663635259509773`*^9, 3.566363546730214*^9}, 3.5663635877426863`*^9, { 3.5663641657237015`*^9, 3.5663642640662737`*^9}, 3.5663644285061626`*^9, 3.5663644816710563`*^9, {3.5663645292199397`*^9, 3.566364585130438*^9}, { 3.5663650531156597`*^9, 3.566365054753663*^9}, {3.566365147714226*^9, 3.5663651764494762`*^9}, {3.566365323167734*^9, 3.566365365833809*^9}, { 3.566383849395733*^9, 3.5663838785677843`*^9}, {3.566383973837152*^9, 3.5663841070301857`*^9}, 3.566384191894335*^9, {3.5663843074125376`*^9, 3.566384326803372*^9}, {3.566384842134677*^9, 3.5663848872655563`*^9}, { 3.566386201848665*^9, 3.5663862026910667`*^9}}], Cell[BoxData[ Graphics3DBox[{SphereBox[{0, 0, 0}, 0.15508308756218367`], SphereBox[{0, 0, 0}, 0.1], {Opacity[0.5], SphereBox[{0, 0, 1.138053549950228}, 0.9829704623880443]}, SphereBox[{-0.477596935163073, 0.8272221572416213, -0.23949396065108544`}, 0.829677130509701], SphereBox[{-0.477596935163073, -0.8272221572416213, -0.2767749418678407}, 0.829677130509701], SphereBox[{0.955193870326146, 0, -0.23949396065108544`}, 0.829677130509701], {RGBColor[ NCache[ Rational[2, 3], 0.6666666666666666], 0, 0], SphereBox[{0, 0, 1.7214268645955368`}, 0.06]}}, AspectRatio->Automatic, Boxed->False, ImageSize->{368.6432719409886, 407.}, ViewPoint->{1.5246376839886455`, -2.699389920119084, 1.3560140086744012`}, ViewVertical->{-0.3374982216564296, -0.9474048594238383, \ -0.05850779766207861}]], "Output", CellChangeTimes->{3.566384221331587*^9, 3.566384893817568*^9, 3.5663850523762465`*^9, 3.566385108723545*^9, 3.5663851391123986`*^9, 3.566385442439331*^9, 3.5663854968834267`*^9, 3.5663856065672197`*^9, 3.566385643929285*^9, 3.5663857055337934`*^9, 3.5663857975115547`*^9, 3.566385928364585*^9, 3.5663860166295395`*^9, 3.566386122522526*^9, 3.5663862556687593`*^9, 3.5663870664252243`*^9, 3.5663871844706316`*^9, 3.5663872180262904`*^9, 3.5663872768383937`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"plot2", "=", RowBox[{ RowBox[{"Graphics", "[", RowBox[{"{", RowBox[{ RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", "P"}], "]"}], ",", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", "0.1"}], "]"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", RowBox[{"P", "+", "Q"}]}], "}"}], ",", "Q"}], "]"}], ",", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", RowBox[{"P", "+", "Q", "+", "S"}]}], "}"}], ",", "0.1"}], "]"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"wd", "*", RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}]}], ",", RowBox[{"wc", "*", RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}]}]}], "}"}], ",", "R"}], "]"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Dashing", "[", RowBox[{"{", "0.02", "}"}], "]"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "wd"}], "*", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}], "/", "2"}]}], ",", RowBox[{"wc", "*", RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}]}]}], "}"}], ",", "R"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Thickness", "[", "0.01", "]"}], ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", RowBox[{"P", "+", "Q", "+", "S"}]}], "}"}]}], "}"}], "]"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Dashing", "[", RowBox[{"{", "0.005", "}"}], "]"}], ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "wd"}], "*", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}], "/", "2"}]}], ",", RowBox[{"wc", "*", RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}]}]}], "}"}]}], "}"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Dashing", "[", RowBox[{"{", "0.005", "}"}], "]"}], ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"wd", "*", RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}]}], ",", RowBox[{"wc", "*", RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}]}]}], "}"}]}], "}"}], "]"}]}], "}"}]}], " ", "}"}]}], "}"}], "]"}], "/.", "u"}]}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{"Show", "[", RowBox[{"plot2", ",", RowBox[{"{", RowBox[{ RowBox[{"AspectRatio", "\[Rule]", "Automatic"}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "2.5"}], ",", "2"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "1.5"}], ",", "2.5"}], "}"}]}], "}"}]}], ",", RowBox[{"Axes", "\[Rule]", "True"}], ",", RowBox[{"GridLines", "\[Rule]", "Automatic"}], ",", RowBox[{"Frame", "\[Rule]", "True"}]}], "}"}]}], "]"}]}], "Input", CellChangeTimes->{{3.5661428484442425`*^9, 3.5661428783026953`*^9}, { 3.5661430328053665`*^9, 3.5661430369081736`*^9}, {3.5661430998854847`*^9, 3.5661431383537183`*^9}, {3.566143180801393*^9, 3.5661432248645825`*^9}, { 3.566143288683233*^9, 3.5661432893462706`*^9}, {3.566143344240477*^9, 3.566143410119393*^9}, 3.566308742431816*^9, {3.5663176430442986`*^9, 3.566317682730768*^9}, {3.5663204234522157`*^9, 3.566320431267829*^9}, { 3.5663204642618876`*^9, 3.5663204702366977`*^9}, {3.566362812374524*^9, 3.5663628459301834`*^9}, {3.5663629092506943`*^9, 3.566362925334323*^9}, { 3.5663630995242286`*^9, 3.5663631074646425`*^9}, {3.566363142408704*^9, 3.5663631503023176`*^9}, {3.566363335505843*^9, 3.5663634250188007`*^9}, { 3.5663635259509773`*^9, 3.566363546730214*^9}, 3.5663635877426863`*^9, { 3.5663641657237015`*^9, 3.5663642640662737`*^9}, 3.5663644285061626`*^9, 3.5663644816710563`*^9, {3.5663645292199397`*^9, 3.566364585130438*^9}, { 3.5663650531156597`*^9, 3.566365054753663*^9}, {3.566365147714226*^9, 3.5663651764494762`*^9}, {3.566365323167734*^9, 3.566365365833809*^9}, { 3.566383849395733*^9, 3.5663838785677843`*^9}, {3.566383973837152*^9, 3.5663841070301857`*^9}, 3.566384191894335*^9, {3.5663851922772923`*^9, 3.5663854368389215`*^9}, {3.5663854838262043`*^9, 3.5663855197842674`*^9}, {3.56638554987672*^9, 3.566385586770785*^9}, { 3.5663856269408555`*^9, 3.566385636316472*^9}, {3.566385700650985*^9, 3.5663857193086176`*^9}, {3.5663857518814745`*^9, 3.5663857852811337`*^9}, {3.5663859175693655`*^9, 3.566385920658171*^9}, { 3.5663860574704113`*^9, 3.5663860852384605`*^9}, {3.566386249023148*^9, 3.566386251846753*^9}, {3.566387051105997*^9, 3.566387057907609*^9}, { 3.5663871787454214`*^9, 3.566387213112282*^9}, {3.566387262330368*^9, 3.566387271799585*^9}}], Cell[BoxData[ GraphicsBox[{CircleBox[{0, 0}, 0.15508308756218367`], DiskBox[{0, 0}, 0.1], CircleBox[{0, 1.138053549950228}, 0.9829704623880443], DiskBox[{0, 1.7214268645955368`}, 0.1], CircleBox[{0.955193870326146, -0.23949396065108544`}, 0.829677130509701], {Dashing[{0.02}], CircleBox[{-0.477596935163073, -0.23949396065108544`}, 0.829677130509701]}, {Thickness[0.01], LineBox[{{0, 0}, {0, 1.7214268645955368`}}], {Dashing[{0.005}], LineBox[{{0, 0}, {-0.477596935163073, -0.23949396065108544`}}]}, {Dashing[{0.005}], LineBox[{{0, 0}, {0.955193870326146, -0.23949396065108544`}}]}}}, AspectRatio->Automatic, Axes->True, Frame->True, GridLines->Automatic, PlotRange->{{-2.5, 2}, {-1.5, 2.5}}]], "Output", CellChangeTimes->{3.566384140539045*^9, 3.5663842213939867`*^9, 3.5663848938799677`*^9, 3.5663850524386463`*^9, 3.566385108785945*^9, 3.5663851391591988`*^9, 3.5663854425017314`*^9, 3.566385496914627*^9, 3.5663856065984197`*^9, 3.566385643976085*^9, 3.5663857055805936`*^9, 3.5663857975739555`*^9, 3.566385928426985*^9, 3.5663860166919403`*^9, 3.566386122569326*^9, 3.56638625571556*^9, 3.566387066472024*^9, 3.5663871845174313`*^9, 3.5663872180730906`*^9, 3.566387276885194*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"plot1", "=", RowBox[{ RowBox[{"Graphics", "[", RowBox[{"{", RowBox[{ RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", "P"}], "]"}], ",", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", "0.1"}], "]"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", "Q"}], "]"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "wd"}], "*", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}], "/", "2"}]}], ",", RowBox[{"wd", "*", "w3", "*", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}], "/", "2"}]}]}], "}"}], ",", "R"}], "]"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", "wd"}], "*", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}], "/", "2"}]}], ",", RowBox[{ RowBox[{"-", "wd"}], "*", "w3", "*", RowBox[{ RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}], "/", "2"}]}]}], "}"}], ",", "R"}], "]"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"wd", "*", RowBox[{"(", RowBox[{"P", "+", "R"}], ")"}]}], ",", "0"}], "}"}], ",", "R"}], "]"}]}], "}"}], "]"}], "/.", "u"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"Show", "[", RowBox[{"plot1", ",", RowBox[{"{", RowBox[{ RowBox[{"AspectRatio", "\[Rule]", "Automatic"}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", "1.5"}], ",", "2"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "2"}], ",", "2"}], "}"}]}], "}"}]}], ",", RowBox[{"Axes", "\[Rule]", "True"}], ",", RowBox[{"GridLines", "\[Rule]", "Automatic"}], ",", RowBox[{"Frame", "\[Rule]", "True"}]}], "}"}]}], "]"}], "\[IndentingNewLine]"}], "Input", CellChangeTimes->{{3.5661464147646685`*^9, 3.566146439287912*^9}, { 3.566146505400828*^9, 3.566146609468611*^9}, 3.5661466450210733`*^9, { 3.5661467117267904`*^9, 3.5661467122883916`*^9}, {3.566146817869377*^9, 3.5661468378998117`*^9}, {3.5661468705194693`*^9, 3.5661468715646715`*^9}, {3.5661469228107615`*^9, 3.5661469312971764`*^9}, {3.5661469641976337`*^9, 3.5661469774888573`*^9}, {3.5661477373577533`*^9, 3.5661477492774353`*^9}, {3.5663088744236484`*^9, 3.566308893580482*^9}, { 3.566317471849598*^9, 3.5663175599429526`*^9}, {3.5663175959010158`*^9, 3.56631760969144*^9}, {3.566321253732474*^9, 3.5663213099081726`*^9}, { 3.5663213595474596`*^9, 3.5663214564548297`*^9}, {3.5663214993705053`*^9, 3.566321505719716*^9}, {3.5663216132507052`*^9, 3.5663216271815295`*^9}, { 3.56632166145479*^9, 3.566321670892806*^9}, 3.5663217012504597`*^9, { 3.566321797674229*^9, 3.566321848187118*^9}, 3.5663838550117435`*^9, 3.566383890408205*^9, {3.5663849894145355`*^9, 3.566385046136235*^9}, { 3.5663850934823184`*^9, 3.5663850968363247`*^9}, {3.5663851324199867`*^9, 3.56638513435439*^9}, {3.566386098108483*^9, 3.5663861172809167`*^9}}], Cell[BoxData[ GraphicsBox[{CircleBox[{0, 0}, 0.15508308756218367`], DiskBox[{0, 0}, 0.1], CircleBox[{0, 0}, 0.9829704623880443], CircleBox[{-0.477596935163073, 0.8272221572416213}, 0.829677130509701], CircleBox[{-0.477596935163073, -0.8272221572416213}, 0.829677130509701], CircleBox[{0.955193870326146, 0}, 0.829677130509701]}, AspectRatio->Automatic, Axes->True, Frame->True, GridLines->Automatic, PlotRange->{{-1.5, 2}, {-2, 2}}]], "Output", CellChangeTimes->{3.5663841406170444`*^9, 3.566384221425187*^9, 3.5663848939267683`*^9, 3.5663850524542465`*^9, 3.5663851088171453`*^9, 3.566385139174799*^9, 3.5663854425173316`*^9, 3.5663854969458275`*^9, 3.5663856066296196`*^9, 3.566385643991685*^9, 3.5663857055961933`*^9, 3.566385797605155*^9, 3.566385928458185*^9, 3.56638601670754*^9, 3.566386122600526*^9, 3.5663862557311597`*^9, 3.566387066487624*^9, 3.5663871845330315`*^9, 3.5663872181042905`*^9, 3.566387276900794*^9}] }, Open ]] }, WindowSize->{707, 817}, WindowMargins->{{392, Automatic}, {Automatic, 15}}, FrontEndVersion->"9.0 for Microsoft Windows (64-bit) (November 20, 2012)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[557, 20, 1850, 42, 173, "Input"], Cell[2410, 64, 8795, 281, 1052, "Input"], Cell[11208, 347, 1365, 30, 112, "Input"], Cell[12576, 379, 1388, 34, 92, "Input"], Cell[CellGroupData[{ Cell[13989, 417, 1006, 25, 92, "Input"], Cell[14998, 444, 639, 14, 31, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[15674, 463, 1786, 46, 332, "Input"], Cell[17463, 511, 255, 4, 31, "Output"], Cell[17721, 517, 238, 3, 31, "Output"], Cell[17962, 522, 240, 3, 31, "Output"], Cell[18205, 527, 256, 4, 31, "Output"], Cell[18464, 533, 238, 3, 31, "Output"], Cell[18705, 538, 239, 3, 31, "Output"], Cell[18947, 543, 240, 3, 31, "Output"], Cell[19190, 548, 239, 3, 31, "Output"], Cell[19432, 553, 240, 3, 31, "Output"], Cell[19675, 558, 240, 3, 31, "Output"], Cell[19918, 563, 237, 3, 31, "Output"], Cell[20158, 568, 238, 3, 31, "Output"] }, Open ]], Cell[20411, 574, 751, 18, 52, "Input"], Cell[CellGroupData[{ Cell[21187, 596, 5711, 165, 292, "Input"], Cell[26901, 763, 1048, 20, 379, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[27986, 788, 5377, 142, 272, "Input"], Cell[33366, 932, 1218, 26, 356, "Output"] }, Open ]], Cell[34599, 961, 475, 9, 31, "Input"], Cell[CellGroupData[{ Cell[35099, 974, 5434, 152, 252, "Input"], Cell[40536, 1128, 1088, 22, 379, "Output"] }, Open ]], Cell[41639, 1153, 13470, 397, 173, "Input"], Cell[55112, 1552, 10484, 317, 1072, "Input"], Cell[65599, 1871, 1966, 40, 112, "Input"], Cell[67568, 1913, 1584, 36, 92, "Input"], Cell[CellGroupData[{ Cell[69177, 1953, 1301, 30, 72, "Input"], Cell[70481, 1985, 1171, 22, 52, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[71689, 2012, 1808, 46, 312, "Input"], Cell[73500, 2060, 722, 11, 31, "Output"], Cell[74225, 2073, 705, 10, 31, "Output"], Cell[74933, 2085, 701, 10, 31, "Output"], Cell[75637, 2097, 703, 11, 31, "Output"], Cell[76343, 2110, 702, 10, 31, "Output"], Cell[77048, 2122, 702, 10, 31, "Output"], Cell[77753, 2134, 702, 10, 31, "Output"], Cell[78458, 2146, 702, 10, 31, "Output"], Cell[79163, 2158, 702, 10, 31, "Output"], Cell[79868, 2170, 701, 10, 31, "Output"], Cell[80572, 2182, 701, 10, 31, "Output"], Cell[81276, 2194, 702, 10, 31, "Output"] }, Open ]], Cell[81993, 2207, 707, 12, 31, "Input"], Cell[CellGroupData[{ Cell[82725, 2223, 7121, 184, 352, "Input"], Cell[89849, 2409, 1477, 27, 381, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[91363, 2441, 5805, 143, 232, "Input"], Cell[97171, 2586, 1572, 30, 341, "Output"] }, Open ]], Cell[98758, 2619, 12449, 362, 173, "Input"], Cell[111210, 2983, 9516, 275, 832, "Input"], Cell[120729, 3260, 1858, 34, 112, "Input"], Cell[122590, 3296, 1635, 37, 92, "Input"], Cell[CellGroupData[{ Cell[124250, 3337, 1359, 30, 72, "Input"], Cell[125612, 3369, 1675, 29, 31, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[127324, 3403, 1985, 52, 272, "Input"], Cell[129312, 3457, 1280, 19, 31, "Output"], Cell[130595, 3478, 1263, 18, 31, "Output"], Cell[131861, 3498, 1263, 18, 31, "Output"], Cell[133127, 3518, 1280, 19, 31, "Output"], Cell[134410, 3539, 1264, 18, 31, "Output"], Cell[135677, 3559, 1264, 18, 31, "Output"], Cell[136944, 3579, 1264, 18, 31, "Output"], Cell[138211, 3599, 1264, 18, 31, "Output"], Cell[139478, 3619, 1263, 18, 31, "Output"], Cell[140744, 3639, 1264, 18, 31, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[142045, 3662, 4805, 117, 212, "Input"], Cell[146853, 3781, 1341, 26, 422, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[148231, 3812, 6110, 148, 212, "Input"], Cell[154344, 3962, 1267, 24, 335, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[155648, 3991, 3596, 90, 172, "Input"], Cell[159247, 4083, 970, 17, 423, "Output"] }, Open ]] } ] *) (* End of internal cache information *)