(* 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[ 37520, 953] NotebookOptionsPosition[ 36308, 913] NotebookOutlinePosition[ 36792, 931] CellTagsIndexPosition[ 36749, 928] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[BoxData[ StyleBox[ RowBox[{"(*", " ", RowBox[{ RowBox[{"C2H6", " ", "Molecule", " ", "after", " ", "Kimball"}], ",", " ", RowBox[{"non", " ", "overlapping"}], ",", RowBox[{ RowBox[{"parameterized", " ", "with", " ", "G3MP2"}], "//", RowBox[{"6", "-", RowBox[{"31", "g", RowBox[{"(", "d", ")"}]}]}]}]}], " ", "*)"}], "Text"]], "Input", CellChangeTimes->{ 3.5788956762396965`*^9, {3.5788957137421627`*^9, 3.578895724880582*^9}, { 3.578895806827526*^9, 3.578895904780098*^9}, {3.5788960148694916`*^9, 3.5788960197835*^9}, {3.5788961272676888`*^9, 3.578896146970524*^9}, { 3.5788962074362297`*^9, 3.578896248542302*^9}, {3.5788963748741245`*^9, 3.578896407337781*^9}, 3.5789031386419754`*^9, {3.5789304190931296`*^9, 3.5789304212303333`*^9}, {3.5789317063449907`*^9, 3.5789317142698045`*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"Clear", "[", RowBox[{"k1", ",", "k2", ",", "k3", ",", "s1", ",", "s2", ",", "s3"}], "]"}], ";"}], "\[IndentingNewLine]", RowBox[{" ", RowBox[{ RowBox[{"n", "=", "2."}], ";"}], "\[IndentingNewLine]", StyleBox[ RowBox[{"(*", " ", RowBox[{"Parameterliste", " ", "c"}], " ", "*)"}], "Text"]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"c", " ", "=", " ", RowBox[{"{", RowBox[{ RowBox[{"k1", " ", "->", " ", "1.024905"}], ",", RowBox[{"k2", " ", "->", " ", "1.37426345"}], ",", RowBox[{"k3", " ", "->", " ", "1.20537762"}], ",", RowBox[{"s1", " ", "->", " ", "0.30582536"}], ",", RowBox[{"s2", " ", "->", " ", "0.30677632"}], ",", RowBox[{"s3", " ", "->", " ", "0.3418"}]}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{"c", " ", "=", " ", RowBox[{"{", RowBox[{ RowBox[{"k1", " ", "->", " ", "1.02246687"}], ",", RowBox[{"k2", " ", "->", " ", "1.37426345"}], ",", RowBox[{"k3", " ", "->", " ", "1.20537762"}], ",", RowBox[{"s1", " ", "->", " ", "0.30582536"}], ",", RowBox[{"s2", " ", "->", " ", "0.30677632"}], ",", RowBox[{"s3", " ", "->", " ", "0.35441063"}]}], "}"}]}], ";"}], "*)"}]}], "\[IndentingNewLine]", RowBox[{ StyleBox[ RowBox[{"(*", " ", RowBox[{"Terms", " ", "from", " ", RowBox[{ RowBox[{"ChemEdu", "/", "Kimball"}], "/", RowBox[{"Alkan", ".", "sys"}]}]}], " ", "*)"}], "Text"]}]}], "Input", CellChangeTimes->{ 3.5788956762396965`*^9, {3.5788957137421627`*^9, 3.578895724880582*^9}, { 3.578895806827526*^9, 3.578895904780098*^9}, {3.5788960148694916`*^9, 3.5788960197835*^9}, {3.5788961272676888`*^9, 3.578896146970524*^9}, { 3.5788962074362297`*^9, 3.578896248542302*^9}, {3.5788963748741245`*^9, 3.578896407337781*^9}, {3.578897895970396*^9, 3.5788979299784555`*^9}, { 3.578901683206219*^9, 3.5789017130022717`*^9}, 3.5789017480087333`*^9, 3.5789031559112062`*^9, {3.578930508668487*^9, 3.5789305512877617`*^9}, { 3.578930588213027*^9, 3.5789305892270284`*^9}, {3.578930623203888*^9, 3.5789307850073724`*^9}, {3.5789308197642336`*^9, 3.5789308212306356`*^9}, {3.5789308916959596`*^9, 3.5789309168276043`*^9}, {3.578930968635295*^9, 3.578931053967445*^9}, { 3.5789310902531085`*^9, 3.5789312758466344`*^9}}], Cell[BoxData[{ RowBox[{ RowBox[{"T", "=", RowBox[{ RowBox[{ RowBox[{"2.25", "*", "n", "*", RowBox[{"k1", "/", RowBox[{"P", "^", "2"}]}]}], " ", "+", RowBox[{ RowBox[{"(", RowBox[{"n", " ", "-", " ", "1.0"}], ")"}], "*", "2.25", "*", RowBox[{"k2", "/", RowBox[{"Q", "^", "2"}]}]}], " ", "+", " ", RowBox[{"4.5", "*", RowBox[{"(", RowBox[{"n", " ", "+", " ", "1.0"}], ")"}], "*", RowBox[{"k3", "/", RowBox[{"R", "^", "2"}]}]}]}], " ", "/.", "c"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"p", " ", "=", "X"}], ";", " ", RowBox[{"p2", "=", " ", RowBox[{"p", "*", "p"}]}], ";", " ", RowBox[{"z", "=", " ", "6."}], ";", " ", RowBox[{"ad", "=", " ", RowBox[{"8.", "/", "3."}]}], ";", " ", RowBox[{"zd", "=", " ", RowBox[{"2.", "/", " ", "3."}]}], ";", " ", RowBox[{"vd", " ", "=", RowBox[{"4.", "/", " ", "3."}]}], ";", " ", RowBox[{"q", "=", " ", RowBox[{"P", " ", "+", " ", "Q"}]}], ";", " ", RowBox[{"r", "=", " ", RowBox[{"P", " ", "+", " ", "R"}]}], ";", " ", RowBox[{"q2", "=", " ", RowBox[{"q", "*", "q"}]}], ";", " ", RowBox[{"r2", "=", " ", RowBox[{"r", "*", "r"}]}], ";", " ", RowBox[{"r2p2", "=", RowBox[{"r2", "*", "p2"}]}], ";", " ", RowBox[{"r2p", "=", " ", RowBox[{"r2", "*", "p"}]}], ";", " ", RowBox[{"qr", "=", " ", RowBox[{"q", "*", "r"}]}], ";", RowBox[{"qrp", "=", RowBox[{"qr", "*", "p"}]}], ";"}]}], "Input", CellChangeTimes->{ 3.5788956762396965`*^9, {3.5788957137421627`*^9, 3.578895724880582*^9}, { 3.578895806827526*^9, 3.578895904780098*^9}, {3.5788960148694916`*^9, 3.5788960197835*^9}, {3.5788961272676888`*^9, 3.578896146970524*^9}, { 3.5788962074362297`*^9, 3.578896248542302*^9}, {3.5788963748741245`*^9, 3.578896407337781*^9}, {3.578897895970396*^9, 3.5788979299784555`*^9}, { 3.578901683206219*^9, 3.578901728461899*^9}, {3.578902190753111*^9, 3.578902205167536*^9}, {3.578902454549574*^9, 3.578902456686778*^9}}], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"Vee", "=", RowBox[{ RowBox[{ RowBox[{"3.0", "*", "n", "*", RowBox[{"s1", "/", "P"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"n", " ", "-", " ", "1"}], ")"}], "*", "3.0", "*", RowBox[{"s2", "/", " ", "Q"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"n", " ", "+", " ", "1"}], ")"}], "*", "6.0", "*", RowBox[{"s3", "/", " ", "R"}]}], " ", "+", " ", RowBox[{"8.0", "/", " ", "q"}], "+", " ", RowBox[{"24.0", "/", " ", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{"zd", "*", "qr"}], "+", "q2", " ", "+", " ", "r2"}], "]"}]}], "+", " ", RowBox[{"24.0", "/", " ", "r"}], " ", "+", " ", RowBox[{"4.0", "/", RowBox[{"(", RowBox[{"2", "*", "q"}], ")"}]}], " ", "+", " ", RowBox[{"24.0", "/", " ", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{"vd", "*", "qr"}], " ", "+", " ", RowBox[{"4.0", "*", "q2"}], " ", "+", " ", "r2"}], "]"}]}], "+", " ", RowBox[{"24.0", "/", " ", RowBox[{"Sqrt", "[", RowBox[{"ad", "*", "r2"}], "]"}]}], " ", "+", " ", RowBox[{"12.0", "/", " ", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{"ad", "*", "qr"}], " ", "+", " ", RowBox[{"4.0", "*", "q2"}], " ", "+", RowBox[{"4.0", "*", "r2"}]}], "]"}]}], "+", " ", RowBox[{"24.0", "/", " ", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{"ad", "*", "qr"}], " ", "+", " ", RowBox[{"4.0", "*", "q2"}], " ", "+", " ", RowBox[{"vd", "*", "r2"}]}], "]"}]}]}], " ", "/.", "c"}]}], ";"}], "\[IndentingNewLine]"}]], "Input", CellChangeTimes->{ 3.5788956762396965`*^9, {3.5788957137421627`*^9, 3.578895724880582*^9}, { 3.578895806827526*^9, 3.578895904780098*^9}, {3.5788960148694916`*^9, 3.5788960197835*^9}, {3.5788961272676888`*^9, 3.578896146970524*^9}, { 3.5788962074362297`*^9, 3.578896248542302*^9}, {3.5788963748741245`*^9, 3.5788964162141967`*^9}, {3.5788976338275356`*^9, 3.578897654762772*^9}, { 3.57889796086651*^9, 3.578897963409314*^9}, {3.5789018226860642`*^9, 3.5789018266172714`*^9}, {3.578901856694124*^9, 3.5789019010138016`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"Vne", " ", "=", RowBox[{ RowBox[{ RowBox[{"-", "2.0"}], "*", RowBox[{"(", RowBox[{"n", " ", "+", " ", "1.0"}], ")"}], "*", RowBox[{ RowBox[{"(", RowBox[{"3.0", " ", "-", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"(", RowBox[{"p", " ", "-", " ", "1.0"}], ")"}], "*", RowBox[{"(", RowBox[{"1.0", " ", "+", " ", RowBox[{"P", "/", " ", "R"}]}], ")"}]}], ")"}], "^", "2"}]}], ")"}], "/", "R"}]}], " ", "-", " ", RowBox[{"3.0", "*", "n", "*", RowBox[{"z", "/", " ", "P"}]}], "-", " ", RowBox[{"4.0", "*", RowBox[{"z", "/", " ", "q"}]}], " ", "-", " ", RowBox[{"12.0", "/", " ", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{"zd", "*", "qrp"}], " ", "+", " ", "r2p2", " ", "+", " ", "q2"}], "]"}]}], "-", " ", RowBox[{"12.0", "*", RowBox[{"z", "/", " ", RowBox[{"Sqrt", "[", "r2", "]"}]}]}], " ", "-", " ", RowBox[{"12.0", "/", " ", RowBox[{"(", RowBox[{"r", "*", "p"}], ")"}]}], " ", "-", " ", RowBox[{"4.0", "*", RowBox[{"z", "/", " ", RowBox[{"(", RowBox[{"2", "*", "q"}], ")"}]}]}], "-", " ", RowBox[{"12.0", "*", RowBox[{"z", "/", " ", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{"vd", "*", "qr"}], " ", "+", " ", RowBox[{"4.0", "*", "q2"}], " ", "+", " ", "r2"}], "]"}]}]}], " ", "-", " ", RowBox[{"12.0", "/", " ", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{"vd", "*", "qrp"}], " ", "+", " ", "r2p2", " ", "+", " ", RowBox[{"4.0", "*", "q2"}]}], "]"}]}], "-", " ", RowBox[{"24.0", "/", " ", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{"zd", "*", "r2p"}], " ", "+", " ", "r2p2", " ", "+", " ", "r2"}], "]"}]}], "-", RowBox[{"12.0", "/", " ", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{"vd", "*", "qr"}], " ", "+", " ", RowBox[{"vd", "*", "qrp"}], " ", "+", " ", RowBox[{"2.0", "*", "r2p"}], " ", "+", " ", "r2p2", " ", "+", " ", RowBox[{"4.0", "*", "q2"}], " ", "+", " ", "r2"}], "]"}]}], "-", " ", RowBox[{"24.0", "/", " ", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{"vd", "*", "qr"}], " ", "+", " ", RowBox[{"vd", "*", "qrp"}], " ", "-", " ", RowBox[{"zd", "*", "r2p"}], " ", "+", " ", "r2p2", " ", "+", " ", RowBox[{"4.0", "*", "q2"}], " ", "+", " ", "r2"}], "]"}]}]}]}], " ", ";"}], "\[IndentingNewLine]"}]], "Input", CellChangeTimes->{ 3.5788956762396965`*^9, {3.5788957137421627`*^9, 3.578895724880582*^9}, { 3.578895806827526*^9, 3.578895904780098*^9}, {3.5788960148694916`*^9, 3.5788960197835*^9}, {3.5788961272676888`*^9, 3.578896146970524*^9}, { 3.5788962074362297`*^9, 3.578896248542302*^9}, {3.5788963748741245`*^9, 3.5788964212998056`*^9}, {3.5788976724220033`*^9, 3.5788976844028244`*^9}, {3.5788979771529384`*^9, 3.578897979212142*^9}, { 3.5789019299206524`*^9, 3.578901999605975*^9}, {3.5789024865140305`*^9, 3.578902487543632*^9}}], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"Vnn", " ", "=", RowBox[{ RowBox[{"6.0", "*", RowBox[{"z", "/", " ", RowBox[{"(", RowBox[{"r", "*", "p"}], ")"}]}]}], " ", "+", " ", RowBox[{"z", "*", RowBox[{"z", "/", " ", RowBox[{"(", RowBox[{"2", "*", "q"}], ")"}]}]}], "+", " ", RowBox[{"6.0", "*", RowBox[{"z", "/", " ", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{"vd", "*", "qrp"}], " ", "+", " ", "r2p2", " ", "+", " ", RowBox[{"4.0", "*", "q2"}]}], "]"}]}]}], " ", "+", " ", RowBox[{"6.0", "/", " ", RowBox[{"Sqrt", "[", RowBox[{"ad", "*", "r2p2"}], "]"}]}], "+", " ", RowBox[{"3.0", "/", " ", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{"ad", "*", "qrp"}], " ", "+", " ", RowBox[{"4.0", "*", "r2p2"}], " ", "+", " ", RowBox[{"4.0", "*", "q2"}]}], "]"}]}], "+", " ", RowBox[{"6.0", "/", " ", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{"ad", "*", "qrp"}], " ", "+", " ", RowBox[{"vd", "*", "r2p2"}], " ", "+", " ", RowBox[{"4.0", "*", "q2"}]}], "]"}]}]}]}], " ", ";"}], "\[IndentingNewLine]"}]], "Input", CellChangeTimes->{ 3.5788956762396965`*^9, {3.5788957137421627`*^9, 3.578895724880582*^9}, { 3.578895806827526*^9, 3.578895904780098*^9}, {3.5788960148694916`*^9, 3.5788960197835*^9}, {3.5788961272676888`*^9, 3.578896146970524*^9}, { 3.5788962074362297`*^9, 3.578896248542302*^9}, {3.5788963748741245`*^9, 3.5788964212998056`*^9}, {3.57889769338844*^9, 3.5788976974288473`*^9}, { 3.5788979843601513`*^9, 3.578897986466155*^9}, {3.578902015393203*^9, 3.578902039620045*^9}, {3.5789024960924473`*^9, 3.578902497106449*^9}}], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"func", " ", "=", " ", RowBox[{ "T", " ", "+", " ", "Vne", " ", "+", " ", "Vee", " ", "+", " ", "Vnn"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"t", " ", "=", " ", RowBox[{"FindMinimum", "[", RowBox[{"func", ",", RowBox[{"{", RowBox[{"P", ",", "0.25", ",", "0.27"}], "}"}], ",", RowBox[{"{", RowBox[{"Q", ",", "1.1", ",", "1.2"}], "}"}], ",", RowBox[{"{", RowBox[{"R", ",", "1.2", ",", "1.3"}], "}"}], ",", RowBox[{"{", RowBox[{"X", ",", "1.3", ",", "1.4"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Method", "\[Rule]", "Automatic"}], ",", RowBox[{"MaxIterations", "\[Rule]", "500"}]}], "}"}]}], "]"}]}], "\[IndentingNewLine]"}], "Input", CellChangeTimes->{ 3.5788956762396965`*^9, {3.5788957137421627`*^9, 3.578895724880582*^9}, { 3.578895806827526*^9, 3.578895904780098*^9}, {3.5788960148694916`*^9, 3.5788960197835*^9}, {3.5788961272676888`*^9, 3.578896146970524*^9}, { 3.5788962074362297`*^9, 3.578896248542302*^9}, {3.5788963748741245`*^9, 3.57889639532576*^9}, {3.5788978015434303`*^9, 3.5788978862671785`*^9}, { 3.5788979169212327`*^9, 3.5788979187152357`*^9}, {3.5789023036037087`*^9, 3.578902308065317*^9}, {3.578902360762209*^9, 3.57890236119901*^9}, { 3.578902644932308*^9, 3.578902662295139*^9}, {3.578903281959427*^9, 3.578903341645132*^9}, {3.5789033823612037`*^9, 3.5789034163380632`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"-", "79.6502549833623`"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", "\[Rule]", "0.26195329698878417`"}], ",", RowBox[{"Q", "\[Rule]", "1.1813491179927689`"}], ",", RowBox[{"R", "\[Rule]", "1.261694073789718`"}], ",", RowBox[{"X", "\[Rule]", "1.3567918035094888`"}]}], "}"}]}], "}"}]], "Output", CellChangeTimes->{ 3.578903430222088*^9, {3.5789304269399433`*^9, 3.578930454224391*^9}, { 3.578930537996538*^9, 3.5789305564669704`*^9}, 3.5789305967930417`*^9, { 3.578930629568699*^9, 3.578930760406129*^9}, 3.578930791247383*^9, 3.578930827907448*^9, 3.5789308964695683`*^9, 3.5789309357504373`*^9, { 3.5789309749845057`*^9, 3.578931060457056*^9}, {3.5789310943091154`*^9, 3.578931280261442*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{"u", " ", "=", " ", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"T", " ", "/.", "u"}], "\[IndentingNewLine]", RowBox[{"Vne", " ", "/.", "u"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"Vee", " ", "/.", "u"}], " ", "/.", "c"}], "\[IndentingNewLine]", RowBox[{"Vnn", " ", "/.", "u"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"ch", " ", "=", " ", RowBox[{ RowBox[{"X", "*", RowBox[{"(", RowBox[{"P", " ", "+", " ", "R"}], ")"}]}], " ", "/.", "u"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"dch", " ", "=", " ", RowBox[{"ch", "*", "0.529177"}]}], " ", StyleBox[ RowBox[{"(*", " ", RowBox[{"C", "-", RowBox[{"H", " ", "distance"}]}], " ", "*)"}], "Text"]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"cc", " ", "=", " ", RowBox[{ RowBox[{"2.", "*", RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q"}], ")"}]}], " ", "/.", "u"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"dcc", "=", RowBox[{"cc", "*", "0.529177"}]}], " ", StyleBox[ RowBox[{"(*", " ", RowBox[{"C", "-", RowBox[{"C", " ", "distance"}]}], " ", "*)"}], "Text"]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"Vne", " ", "+", " ", "Vee", " ", "+", " ", "Vnn"}], ")"}]}], "/", "T"}], " ", "/.", "u"}], " ", "/.", "c"}], " ", StyleBox[ RowBox[{"(*", " ", RowBox[{"Virial", " ", "Theorem"}], " ", "*)"}], "Text"]}], "\[IndentingNewLine]", RowBox[{ StyleBox[ RowBox[{"(*", RowBox[{ "atomization", " ", "energy", " ", "at", " ", "0", " ", "K", " ", "in", " ", RowBox[{"kcal", "/", "mol"}]}], " ", "*)"}], "Text"], RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"-", RowBox[{"t", "[", RowBox[{"[", "1", "]"}], "]"}]}], " ", "-", " ", RowBox[{"n", "*", "37.78899"}], " ", "-", RowBox[{ RowBox[{"(", " ", RowBox[{ RowBox[{"2", "*", "n"}], " ", "+", " ", "2"}], ")"}], "*", "0.50177"}]}], ")"}], "*", "627.50956", " "}]}]}], "Input", CellChangeTimes->{ 3.5788956762396965`*^9, {3.5788957137421627`*^9, 3.578895724880582*^9}, { 3.578895806827526*^9, 3.578895904780098*^9}, {3.5788960148694916`*^9, 3.5788960197835*^9}, {3.5788961272676888`*^9, 3.578896146970524*^9}, { 3.5788962074362297`*^9, 3.578896248542302*^9}, {3.5788963748741245`*^9, 3.57889639532576*^9}, 3.578896439614238*^9, {3.578902116762181*^9, 3.5789021273389993`*^9}, 3.5789025446709323`*^9, 3.5789025890062103`*^9, { 3.5789028686835012`*^9, 3.578902880882723*^9}, {3.5789029284316063`*^9, 3.578902967977676*^9}, {3.5789030373665977`*^9, 3.5789031226987476`*^9}, { 3.578903474307765*^9, 3.5789034887065907`*^9}, {3.578930301359723*^9, 3.5789303077401342`*^9}, {3.578930371263446*^9, 3.5789303949442873`*^9}, 3.578930448764382*^9}], Cell[BoxData["79.6501698556308`"], "Output", CellChangeTimes->{ 3.578903430222088*^9, {3.5789304269555435`*^9, 3.578930454239991*^9}, { 3.5789305380121384`*^9, 3.5789305564825706`*^9}, 3.578930596808642*^9, { 3.578930629584299*^9, 3.578930760421729*^9}, 3.5789307912629833`*^9, 3.5789308279230475`*^9, 3.578930896485168*^9, 3.578930935766037*^9, { 3.578930975000106*^9, 3.5789310604726562`*^9}, {3.5789310943247156`*^9, 3.5789312802770424`*^9}}], Cell[BoxData[ RowBox[{"-", "271.252926745166`"}]], "Output", CellChangeTimes->{ 3.578903430222088*^9, {3.5789304269555435`*^9, 3.578930454239991*^9}, { 3.5789305380121384`*^9, 3.5789305564825706`*^9}, 3.578930596808642*^9, { 3.578930629584299*^9, 3.578930760421729*^9}, 3.5789307912629833`*^9, 3.5789308279230475`*^9, 3.578930896485168*^9, 3.578930935766037*^9, { 3.578930975000106*^9, 3.5789310604726562`*^9}, {3.5789310943247156`*^9, 3.5789312802770424`*^9}}], Cell[BoxData["69.65069838184246`"], "Output", CellChangeTimes->{ 3.578903430222088*^9, {3.5789304269555435`*^9, 3.578930454239991*^9}, { 3.5789305380121384`*^9, 3.5789305564825706`*^9}, 3.578930596808642*^9, { 3.578930629584299*^9, 3.578930760421729*^9}, 3.5789307912629833`*^9, 3.5789308279230475`*^9, 3.578930896485168*^9, 3.578930935766037*^9, { 3.578930975000106*^9, 3.5789310604726562`*^9}, {3.5789310943247156`*^9, 3.5789312802770424`*^9}}], Cell[BoxData["42.30180352433044`"], "Output", CellChangeTimes->{ 3.578903430222088*^9, {3.5789304269555435`*^9, 3.578930454239991*^9}, { 3.5789305380121384`*^9, 3.5789305564825706`*^9}, 3.578930596808642*^9, { 3.578930629584299*^9, 3.578930760421729*^9}, 3.5789307912629833`*^9, 3.5789308279230475`*^9, 3.578930896485168*^9, 3.578930935766037*^9, { 3.578930975000106*^9, 3.5789310604726562`*^9}, {3.5789310943247156`*^9, 3.5789312802770424`*^9}}], Cell[BoxData["1.0939529349054957`"], "Output", CellChangeTimes->{ 3.578903430222088*^9, {3.5789304269555435`*^9, 3.578930454239991*^9}, { 3.5789305380121384`*^9, 3.5789305564825706`*^9}, 3.578930596808642*^9, { 3.578930629584299*^9, 3.578930760421729*^9}, 3.5789307912629833`*^9, 3.5789308279230475`*^9, 3.578930896485168*^9, 3.578930935766037*^9, { 3.578930975000106*^9, 3.5789310604726562`*^9}, {3.5789310943247156`*^9, 3.5789312802926426`*^9}}], Cell[BoxData["1.5275248841053866`"], "Output", CellChangeTimes->{ 3.578903430222088*^9, {3.5789304269555435`*^9, 3.578930454239991*^9}, { 3.5789305380121384`*^9, 3.5789305564825706`*^9}, 3.578930596808642*^9, { 3.578930629584299*^9, 3.578930760421729*^9}, 3.5789307912629833`*^9, 3.5789308279230475`*^9, 3.578930896485168*^9, 3.578930935766037*^9, { 3.578930975000106*^9, 3.5789310604726562`*^9}, {3.5789310943247156`*^9, 3.5789312802926426`*^9}}], Cell[BoxData["2.000001068770244`"], "Output", CellChangeTimes->{ 3.578903430222088*^9, {3.5789304269555435`*^9, 3.578930454239991*^9}, { 3.5789305380121384`*^9, 3.5789305564825706`*^9}, 3.578930596808642*^9, { 3.578930629584299*^9, 3.578930760421729*^9}, 3.5789307912629833`*^9, 3.5789308279230475`*^9, 3.578930896485168*^9, 3.578930935766037*^9, { 3.578930975000106*^9, 3.5789310604726562`*^9}, {3.5789310943247156`*^9, 3.5789312802926426`*^9}}], Cell[BoxData["666.1986514814874`"], "Output", CellChangeTimes->{ 3.578903430222088*^9, {3.5789304269555435`*^9, 3.578930454239991*^9}, { 3.5789305380121384`*^9, 3.5789305564825706`*^9}, 3.578930596808642*^9, { 3.578930629584299*^9, 3.578930760421729*^9}, 3.5789307912629833`*^9, 3.5789308279230475`*^9, 3.578930896485168*^9, 3.578930935766037*^9, { 3.578930975000106*^9, 3.5789310604726562`*^9}, {3.5789310943247156`*^9, 3.5789312802926426`*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"plot1", " ", "=", " ", RowBox[{ RowBox[{"Graphics", "[", RowBox[{"{", RowBox[{ RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q"}], ")"}]}], ",", "0"}], "}"}], ",", "P"}], "]"}], ",", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q"}], ")"}]}], ",", "0"}], "}"}], ",", "0.1"}], "]"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"P", " ", "+", " ", "Q"}], ",", "0"}], "}"}], ",", "P"}], "]"}], ",", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"P", " ", "+", " ", "Q"}], ",", "0"}], "}"}], ",", "0.1"}], "]"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{"0", ",", "0"}], "}"}], ",", "Q"}], "]"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q"}], ")"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{"P", " ", "+", " ", "R"}], ")"}], "/", "3."}]}], ",", RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"P", " ", "+", " ", "R"}], ")"}]}], "*", RowBox[{"Sqrt", "[", RowBox[{"8.", "/", " ", "9."}], "]"}]}]}], "}"}], ",", "R"}], "]"}], ",", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q"}], ")"}]}], "-", " ", RowBox[{"p", "*", RowBox[{ RowBox[{"(", RowBox[{"P", " ", "+", " ", "R"}], ")"}], "/", "3."}]}]}], ",", RowBox[{ RowBox[{"-", "p"}], "*", RowBox[{"(", RowBox[{"P", " ", "+", " ", "R"}], ")"}], "*", RowBox[{"Sqrt", "[", RowBox[{"8.", "/", " ", "9."}], "]"}]}]}], "}"}], ",", "0.1"}], "]"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"P", " ", "+", " ", "Q", " ", "+", RowBox[{ RowBox[{"(", RowBox[{"P", " ", "+", " ", "R"}], ")"}], "/", "3."}]}], ",", RowBox[{ RowBox[{"(", RowBox[{"P", " ", "+", " ", "R"}], ")"}], "*", RowBox[{"Sqrt", "[", RowBox[{"8.", "/", " ", "9."}], "]"}]}]}], "}"}], ",", "R"}], "]"}], ",", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q"}], ")"}], "+", " ", RowBox[{"p", "*", RowBox[{ RowBox[{"(", RowBox[{"P", " ", "+", " ", "R"}], ")"}], "/", "3."}]}]}], ",", RowBox[{"p", "*", RowBox[{"(", RowBox[{"P", " ", "+", " ", "R"}], ")"}], "*", RowBox[{"Sqrt", "[", RowBox[{"8.", "/", " ", "9."}], "]"}]}]}], "}"}], ",", "0.1"}], "]"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Dashing", "[", RowBox[{"{", RowBox[{"0.03", ",", "0.03"}], "}"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q"}], ")"}], "+", RowBox[{ RowBox[{"(", RowBox[{"P", " ", "+", " ", "R"}], ")"}], "/", "3."}]}], ",", RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"P", " ", "+", " ", "R"}], ")"}]}], "*", RowBox[{"Sqrt", "[", RowBox[{"2.", "/", " ", "9."}], "]"}]}]}], "}"}], ",", "R"}], "]"}], ",", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"P", " ", "+", " ", "Q", " ", "+", " ", RowBox[{"p", "*", RowBox[{ RowBox[{"(", RowBox[{"P", " ", "+", " ", "R"}], ")"}], "/", "3."}]}]}], ",", RowBox[{ RowBox[{"-", "p"}], "*", RowBox[{"(", RowBox[{"P", " ", "+", " ", "R"}], ")"}], "*", RowBox[{ RowBox[{"Sqrt", "[", "2.", "]"}], "/", "3."}]}]}], "}"}], ",", "0.1"}], "]"}], ",", RowBox[{"Circle", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q"}], ")"}]}], "-", RowBox[{ RowBox[{"(", RowBox[{"P", " ", "+", " ", "R"}], ")"}], "/", "3."}]}], ",", RowBox[{ RowBox[{"(", RowBox[{"P", " ", "+", " ", "R"}], ")"}], "*", RowBox[{"Sqrt", "[", RowBox[{"2.", "/", " ", "9."}], "]"}]}]}], "}"}], ",", "R"}], "]"}], ",", RowBox[{"Disk", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q"}], ")"}]}], "-", RowBox[{"p", "*", RowBox[{ RowBox[{"(", RowBox[{"P", " ", "+", " ", "R"}], ")"}], "/", "3."}]}]}], ",", RowBox[{"p", "*", RowBox[{"(", RowBox[{"P", " ", "+", " ", "R"}], ")"}], "*", RowBox[{ RowBox[{"Sqrt", "[", "2.", "]"}], "/", "3."}]}]}], "}"}], ",", "0.1"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Thickness", "[", "0.01", "]"}], ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q"}], ")"}]}], ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", " ", "+", " ", "Q"}], ",", "0"}], "}"}]}], "}"}], "]"}], ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"P", " ", "+", " ", "Q"}], ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", " ", "+", " ", "Q", " ", "+", " ", RowBox[{"p", "*", RowBox[{ RowBox[{"(", RowBox[{"P", " ", "+", " ", "R"}], ")"}], "/", "3."}]}]}], ",", RowBox[{"p", "*", RowBox[{"(", RowBox[{"P", " ", "+", " ", "R"}], ")"}], "*", RowBox[{"Sqrt", "[", RowBox[{"8.", "/", " ", "9."}], "]"}]}]}], "}"}]}], "}"}], "]"}], ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q"}], ")"}]}], ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q"}], ")"}]}], "-", " ", RowBox[{"p", "*", RowBox[{ RowBox[{"(", RowBox[{"P", " ", "+", " ", "R"}], ")"}], "/", "3."}]}]}], ",", RowBox[{ RowBox[{"-", "p"}], "*", RowBox[{"(", RowBox[{"P", " ", "+", " ", "R"}], ")"}], "*", RowBox[{"Sqrt", "[", RowBox[{"8.", "/", " ", "9."}], "]"}]}]}], "}"}]}], "}"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"{", RowBox[{ RowBox[{"Dashing", "[", RowBox[{"{", RowBox[{"0.03", ",", "0.03"}], "}"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"P", " ", "+", " ", "Q"}], ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"P", " ", "+", " ", "Q", " ", "+", " ", RowBox[{"p", "*", RowBox[{ RowBox[{"(", RowBox[{"P", " ", "+", " ", "R"}], ")"}], "/", "3."}]}]}], ",", RowBox[{ RowBox[{"-", "p"}], "*", RowBox[{"(", RowBox[{"P", " ", "+", " ", "R"}], ")"}], "*", RowBox[{"Sqrt", "[", RowBox[{"2.", "/", " ", "9."}], "]"}]}]}], "}"}]}], "}"}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q"}], ")"}]}], ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q"}], ")"}]}], " ", "-", " ", RowBox[{"p", "*", RowBox[{ RowBox[{"(", RowBox[{"P", " ", "+", " ", "R"}], ")"}], "/", "3."}]}]}], ",", RowBox[{"p", "*", RowBox[{"(", RowBox[{"P", " ", "+", " ", "R"}], ")"}], "*", RowBox[{"Sqrt", "[", RowBox[{"2.", "/", " ", "9."}], "]"}]}]}], "}"}]}], "}"}], "]"}]}], "}"}]}], "}"}]}], "}"}], "]"}], " ", "/.", "u"}]}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{"Show", "[", RowBox[{"plot1", ",", RowBox[{"{", RowBox[{ RowBox[{"AspectRatio", " ", "->", " ", "Automatic"}], ",", RowBox[{"PlotRange", " ", "->", " ", "All"}], ",", RowBox[{"Axes", " ", "->", " ", "True"}], ",", RowBox[{"GridLines", " ", "->", " ", "Automatic"}], ",", RowBox[{"Frame", " ", "->", " ", "True"}]}], "}"}]}], "]"}]}], "Input", CellChangeTimes->{ 3.5788956762396965`*^9, {3.5788957137421627`*^9, 3.578895724880582*^9}, { 3.578895806827526*^9, 3.578895904780098*^9}, {3.5788960148694916`*^9, 3.5788960197835*^9}, {3.5788961272676888`*^9, 3.578896146970524*^9}, { 3.5788962074362297`*^9, 3.578896248542302*^9}, {3.5788963748741245`*^9, 3.57889639532576*^9}, 3.578896439614238*^9, {3.578902116762181*^9, 3.5789021273389993`*^9}, 3.5789025446709323`*^9, 3.5789025890062103`*^9, 3.5789032247541265`*^9}], Cell[BoxData[ GraphicsBox[{CircleBox[{-1.443302414981553, 0}, 0.26195329698878417`], DiskBox[{-1.443302414981553, 0}, 0.1], CircleBox[{1.443302414981553, 0}, 0.26195329698878417`], DiskBox[{1.443302414981553, 0}, 0.1], CircleBox[{0, 0}, 1.1813491179927689`], CircleBox[{-1.9511848719077203`, -1.4365085173527103`}, 1.261694073789718], DiskBox[{-2.132393169685238, -1.9490429820157256`}, 0.1], CircleBox[{1.9511848719077203`, 1.4365085173527103`}, 1.261694073789718], DiskBox[{2.132393169685238, 1.9490429820157256`}, 0.1], {Dashing[{0.03, 0.03}], CircleBox[{1.9511848719077203`, -0.7182542586763552}, 1.261694073789718], DiskBox[{2.132393169685238, -0.9745214910078628}, 0.1], CircleBox[{-1.9511848719077203`, 0.7182542586763552}, 1.261694073789718], DiskBox[{-2.132393169685238, 0.9745214910078628}, 0.1]}, {Thickness[0.01], LineBox[{{-1.443302414981553, 0}, {1.443302414981553, 0}}], LineBox[{{1.443302414981553, 0}, {2.132393169685238, 1.9490429820157256`}}], LineBox[{{-1.443302414981553, 0}, {-2.132393169685238, -1.9490429820157256`}}], {Dashing[{0.03, 0.03}], LineBox[{{1.443302414981553, 0}, { 2.132393169685238, -0.9745214910078628}}], LineBox[{{-1.443302414981553, 0}, {-2.132393169685238, 0.9745214910078628}}]}}}, AspectRatio->Automatic, Axes->True, Frame->True, GridLines->Automatic, PlotRange->All]], "Output", CellChangeTimes->{ 3.578903430253288*^9, {3.578930426986744*^9, 3.5789304542711916`*^9}, { 3.578930538043339*^9, 3.578930556513771*^9}, 3.578930596839842*^9, { 3.5789306296154995`*^9, 3.5789307604529295`*^9}, 3.5789307912941837`*^9, 3.5789308279542475`*^9, 3.578930896516368*^9, 3.5789309357972374`*^9, { 3.5789309750313063`*^9, 3.578931060503856*^9}, {3.5789310943559155`*^9, 3.5789312803082423`*^9}}] }, Open ]] }, WindowSize->{872, 740}, WindowMargins->{{-1, Automatic}, {Automatic, -1}}, PrintingCopies->1, PrintingPageRange->{Automatic, Automatic}, PrivateNotebookOptions->{"VersionedStylesheet"->{"Default.nb"[8.] -> False}}, FrontEndVersion->"9.0 for Microsoft Windows (64-bit) (January 25, 2013)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[557, 20, 861, 17, 31, "Input"], Cell[1421, 39, 2429, 55, 172, "Input"], Cell[3853, 96, 2071, 53, 92, "Input"], Cell[5927, 151, 2276, 54, 92, "Input"], Cell[8206, 207, 3230, 80, 132, "Input"], Cell[11439, 289, 1768, 42, 92, "Input"], Cell[CellGroupData[{ Cell[13232, 335, 1464, 31, 92, "Input"], Cell[14699, 368, 799, 17, 31, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[15535, 390, 3106, 84, 252, "Input"], Cell[18644, 476, 465, 7, 31, "Output"], Cell[19112, 485, 482, 8, 31, "Output"], Cell[19597, 495, 466, 7, 31, "Output"], Cell[20066, 504, 466, 7, 31, "Output"], Cell[20535, 513, 467, 7, 31, "Output"], Cell[21005, 522, 467, 7, 31, "Output"], Cell[21475, 531, 466, 7, 31, "Output"], Cell[21944, 540, 466, 7, 31, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[22447, 552, 11962, 319, 412, "Input"], Cell[34412, 873, 1880, 37, 322, "Output"] }, Open ]] } ] *) (* End of internal cache information *)