(* 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[ 63619, 1572] NotebookOptionsPosition[ 62058, 1523] NotebookOutlinePosition[ 62616, 1542] CellTagsIndexPosition[ 62573, 1539] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[BoxData[ StyleBox[ RowBox[{"(*", " ", RowBox[{ RowBox[{ StyleBox[ RowBox[{ StyleBox["H", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], 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]]}]], " ", "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[ StyleBox[ StyleBox[ StyleBox[ StyleBox[ StyleBox[ StyleBox["\[IndentingNewLine]", "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]], "\[Beta]", " ", "angle", " ", 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["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[{ StyleBox["L", "Text", FontFamily->"Arial Narrow", FontSize->14, FontColor->GrayLevel[0]], StyleBox[ 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"], "\[IndentingNewLine]", 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[ RowBox[{ StyleBox[ StyleBox[ StyleBox[ StyleBox[ StyleBox["L", "Text", FontFamily->"Arial Narrow"], "Text", FontFamily->"Arial Narrow"], "Text", FontFamily->"Arial Narrow"], "Text", FontFamily->"Arial Narrow"], "Text", FontFamily->"Arial Narrow"], StyleBox[ StyleBox[ 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"], StyleBox["\[IndentingNewLine]", FontFamily->"Arial Narrow"], 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}] }, 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}], 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}], 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.5663862556531596`*^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.5663862556531596`*^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.5663862556531596`*^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.5663862556531596`*^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.5663862556531596`*^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.5663862556531596`*^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.5663862556531596`*^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.5663862556531596`*^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}] }, 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.01", "}"}], "]"}], ",", 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.01", "}"}], "]"}], ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", RowBox[{"{", 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}}], 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.01}], LineBox[{{0, 0}, {-0.477596935163073, -0.23949396065108544`}}]}, {Dashing[{0.01}], LineBox[{{0, 0}, {0.9847602180718847, -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}] }, 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}] }, Open ]] }, WindowSize->{829, 1002}, WindowMargins->{{359, Automatic}, {Automatic, 16}}, PrintingCopies->1, PrintingPageRange->{Automatic, Automatic}, PrivateNotebookOptions->{"VersionedStylesheet"->{"Default.nb"[8.] -> False}}, FrontEndVersion->"9.0 for Microsoft Windows (64-bit) (November 20, 2012)", StyleDefinitions->FrontEnd`FileName[{"Utility"}, "Notepad.nb", CharacterEncoding -> "WindowsANSI"] ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[557, 20, 13922, 405, 194, "Input"], Cell[14482, 427, 9689, 280, 832, "Input"], Cell[24174, 709, 1858, 34, 112, "Input"], Cell[26035, 745, 1635, 37, 92, "Input"], Cell[CellGroupData[{ Cell[27695, 786, 1359, 30, 72, "Input"], Cell[29057, 818, 1571, 27, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[30665, 850, 1985, 52, 272, "Input"], Cell[32653, 904, 1178, 17, 30, "Output"], Cell[33834, 923, 1161, 16, 30, "Output"], Cell[34998, 941, 1164, 16, 30, "Output"], Cell[36165, 959, 1181, 17, 30, "Output"], Cell[37349, 978, 1165, 16, 30, "Output"], Cell[38517, 996, 1165, 16, 30, "Output"], Cell[39685, 1014, 1165, 16, 30, "Output"], Cell[40853, 1032, 1165, 16, 30, "Output"], Cell[42021, 1050, 1164, 16, 30, "Output"], Cell[43188, 1068, 1165, 16, 30, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[44390, 1089, 4805, 117, 192, "Input"], Cell[49198, 1208, 1242, 25, 454, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[50477, 1238, 5880, 144, 192, "Input"], Cell[56360, 1384, 1171, 23, 335, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[57568, 1412, 3596, 90, 152, "Input"], Cell[61167, 1504, 875, 16, 423, "Output"] }, Open ]] } ] *) (* End of internal cache information *)