(* Content-type: application/vnd.wolfram.mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 8.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 157, 7] NotebookDataLength[ 43102, 1513] NotebookOptionsPosition[ 42276, 1484] NotebookOutlinePosition[ 42679, 1501] CellTagsIndexPosition[ 42636, 1498] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", " ", RowBox[{ RowBox[{"C2H6", " ", "Molek\[UDoubleDot]l", " ", "nach", " ", "Kimball"}], ",", "nicht\[UDoubleDot]berlappend", ",", RowBox[{ RowBox[{"parametrisiert", " ", "mit", " ", "G2"}], "//", RowBox[{"6", "-", RowBox[{"311", "g"}]}]}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"Clear", "[", RowBox[{"k1", ",", "k2", ",", "k3", ",", "s1", ",", "s2", ",", "s3"}], "]"}], ";", " ", RowBox[{"n", "=", "2."}], ";"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{"Terme", " ", "aus", " ", RowBox[{ RowBox[{"ChemEdu", "/", "Kimball"}], "/", RowBox[{"Alkan", ".", "sys"}]}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"T", "=", 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"}]}]}]}]}], ";"}], "\[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"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"Vee", "=", RowBox[{ RowBox[{"3.0", "*", "n", "*", RowBox[{"s1", "/", "P"}]}], "+", RowBox[{ RowBox[{"(", RowBox[{"n", " ", "-", " ", "1.0"}], ")"}], "*", "3.0", "*", RowBox[{"s2", "/", " ", "Q"}]}]}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"+", RowBox[{"(", RowBox[{"n", " ", "+", " ", "1.0"}], ")"}]}], "*", "6.0", "*", RowBox[{"s3", "/", " ", "R"}]}], " ", "+", " ", RowBox[{"8.0", "/", " ", RowBox[{"Sqrt", "[", "q2", "]"}]}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"+", " ", "24.0"}], "/", " ", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{"zd", "*", "qr"}], "+", "q2", " ", "+", " ", "r2"}], "]"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"+", " ", "24.0"}], "/", " ", RowBox[{"Sqrt", "[", "r2", "]"}]}], " ", "+", " ", RowBox[{"4.0", "/", RowBox[{"Sqrt", "[", RowBox[{"4.0", "*", "q2"}], "]"}]}]}], " ", "\[IndentingNewLine]", RowBox[{ RowBox[{"+", " ", "24.0"}], "/", " ", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{"vd", "*", "qr"}], " ", "+", " ", RowBox[{"4.0", "*", "q2"}], " ", "+", " ", "r2"}], "]"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ 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"}]}], "]"}]}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"+", " ", "24.0"}], "/", " ", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{"ad", "*", "qr"}], " ", "+", " ", RowBox[{"4.0", "*", "q2"}], " ", "+", " ", RowBox[{"vd", "*", "r2"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", 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"}]}]}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"-", " ", "4.0"}], "*", RowBox[{"z", "/", " ", RowBox[{"Sqrt", "[", "q2", "]"}]}]}], " ", "-", " ", RowBox[{"12.0", "/", " ", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{"zd", "*", "qrp"}], " ", "+", " ", "r2p2", " ", "+", " ", "q2"}], "]"}]}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"-", " ", "12.0"}], "*", RowBox[{"z", "/", " ", RowBox[{"Sqrt", "[", "r2", "]"}]}]}], " ", "-", " ", RowBox[{"12.0", "/", " ", RowBox[{"Sqrt", "[", "r2p2", "]"}]}], " ", "-", " ", RowBox[{"4.0", "*", RowBox[{"z", "/", " ", RowBox[{"Sqrt", "[", RowBox[{"4.0", "*", "q2"}], "]"}]}]}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"-", " ", "12.0"}], "*", RowBox[{"z", "/", " ", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{"vd", "*", "qr"}], " ", "+", " ", RowBox[{"4.0", "*", "q2"}], " ", "+", " ", "r2"}], "]"}]}]}], " ", "\[IndentingNewLine]", RowBox[{ RowBox[{"-", " ", "12.0"}], "/", " ", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{"vd", "*", "qrp"}], " ", "+", " ", "r2p2", " ", "+", " ", RowBox[{"4.0", "*", "q2"}]}], "]"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"-", " ", "24.0"}], "/", " ", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{"zd", "*", "r2p"}], " ", "+", " ", "r2p2", " ", "+", " ", "r2"}], "]"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"-", "12.0"}], "/", " ", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{"vd", "*", "qr"}], " ", "+", " ", RowBox[{"vd", "*", "qrp"}], " ", "+", " ", RowBox[{"2.0", "*", "r2p"}], " ", "+", " ", "r2p2", " ", "+", " ", RowBox[{"4.0", "*", "q2"}], " ", "+", " ", "r2"}], "]"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"-", " ", "24.0"}], "/", " ", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{"vd", "*", "qr"}], " ", "+", " ", RowBox[{"vd", "*", "qrp"}], " ", "-", " ", RowBox[{"zd", "*", "r2p"}], " ", "+", " ", "r2p2", " ", "+", " ", RowBox[{"4.0", "*", "q2"}], " ", "+", " ", "r2"}], "]"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"Vnn", " ", "=", RowBox[{ RowBox[{"6.0", "*", RowBox[{"z", "/", " ", RowBox[{"Sqrt", "[", "r2p2", "]"}]}]}], " ", "+", " ", RowBox[{"z", "*", RowBox[{"z", "/", " ", RowBox[{"Sqrt", "[", RowBox[{"4.0", "*", "q2"}], "]"}]}]}]}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"+", " ", "6.0"}], "*", RowBox[{"z", "/", " ", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{"vd", "*", "qrp"}], " ", "+", " ", "r2p2", " ", "+", " ", RowBox[{"4.0", "*", "q2"}]}], "]"}]}]}], " ", "+", " ", RowBox[{"6.0", "/", " ", RowBox[{"Sqrt", "[", RowBox[{"ad", "*", "r2p2"}], "]"}]}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{"+", " ", "3.0"}], "/", " ", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{"ad", "*", "qrp"}], " ", "+", " ", RowBox[{"4.0", "*", "r2p2"}], " ", "+", " ", RowBox[{"4.0", "*", "q2"}]}], "]"}]}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"+", " ", "6.0"}], "/", " ", RowBox[{"Sqrt", "[", RowBox[{ RowBox[{"ad", "*", "qrp"}], " ", "+", " ", RowBox[{"vd", "*", "r2p2"}], " ", "+", " ", RowBox[{"4.0", "*", "q2"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"func", " ", "=", " ", RowBox[{ "T", " ", "+", " ", "Vne", " ", "+", " ", "Vee", " ", "+", " ", "Vnn"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"(*", " ", RowBox[{"Parameterliste", " ", "c"}], " ", "*)"}], "\[IndentingNewLine]", 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[{ RowBox[{"func", "=", RowBox[{"func", " ", "/.", " ", "c"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"N", "[", RowBox[{ RowBox[{"t", " ", "=", " ", RowBox[{"FindMinimum", "[", RowBox[{"func", ",", RowBox[{"{", RowBox[{"P", ",", "0.26141206"}], "}"}], ",", RowBox[{"{", RowBox[{"Q", ",", "1.178240"}], "}"}], ",", RowBox[{"{", RowBox[{"R", ",", "1.2709725"}], "}"}], ",", RowBox[{"{", RowBox[{"X", ",", "1.3575745"}], "}"}]}], "]"}]}], ",", "10"}], "]"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"u", " ", "=", " ", RowBox[{"t", "[", RowBox[{"[", "2", "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"N", "[", RowBox[{ RowBox[{"Vne", " ", "/.", "u"}], ",", "10"}], "]"}], "\[IndentingNewLine]", RowBox[{"N", "[", RowBox[{ RowBox[{ RowBox[{"Vee", " ", "/.", "u"}], " ", "/.", "c"}], ",", "9"}], "]"}], "\[IndentingNewLine]", RowBox[{"N", "[", RowBox[{ RowBox[{"Vnn", " ", "/.", "u"}], ",", "8"}], "]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"ch", " ", "=", " ", RowBox[{ RowBox[{"X", "*", RowBox[{"(", RowBox[{"P", " ", "+", " ", "R"}], ")"}]}], " ", "/.", "u"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"dch", " ", "=", " ", RowBox[{"ch", "*", "0.529177"}]}], " ", RowBox[{"(*", " ", RowBox[{"C", "-", RowBox[{"H", " ", "Abstand"}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"cc", " ", "=", " ", RowBox[{ RowBox[{"2.", "*", RowBox[{"(", RowBox[{"P", " ", "+", " ", "Q"}], ")"}]}], " ", "/.", "u"}]}], ";"}], "\[IndentingNewLine]", RowBox[{"dcc", "=", RowBox[{"cc", "*", "0.529177"}]}], " ", RowBox[{"(*", " ", RowBox[{"C", "-", RowBox[{"C", " ", "Abstand"}]}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{"N", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{"-", RowBox[{"(", RowBox[{"Vne", " ", "+", " ", "Vee", " ", "+", " ", "Vnn"}], ")"}]}], "/", "T"}], " ", "/.", "u"}], " ", "/.", "c"}], ",", "9"}], "]"}], " ", RowBox[{"(*", " ", RowBox[{"Virial", " ", "Theorem"}], " ", "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"-", RowBox[{"t", "[", RowBox[{"[", "1", "]"}], "]"}]}], " ", "-", " ", RowBox[{"n", "*", "37.784301"}], " ", "-", " ", "n", " ", "-", " ", "1.0"}], ")"}], "*", "627.50956"}], " ", RowBox[{"(*", RowBox[{"Atomisierungsenergie", " ", "0", " ", "K", " ", "in", " ", RowBox[{"kcal", "/", "mol"}]}], " ", "*)"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"plot1", " ", "=", " ", 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."}], "]"}]}]}], "}"}]}], "}"}], "]"}]}], "}"}]}], "}"}]}], "}"}], "]"}]}], "\[IndentingNewLine]", RowBox[{"/.", "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"], Cell[OutputFormData["\<\ {-79.63122262216285, {P -> 0.2614120636865185, Q -> 1.178240088471064, R -> \ 1.2709725126384837, X -> 1.3575746431610034}}\ \>", "\<\ {-79.6312, {P -> 0.261412, Q -> 1.17824, R -> 1.27097, X -> 1.35757}}\ \>"], "Output"], Cell[OutputFormData["\<\ -271.1254026035281\ \>", "\<\ -271.125\ \>"], "Output"], Cell[OutputFormData["\<\ 69.66048820683633\ \>", "\<\ 69.6605\ \>"], "Output"], Cell[OutputFormData["\<\ 42.20246938858474\ \>", "\<\ 42.2025\ \>"], "Output"], Cell[OutputFormData["\<\ 1.1008609068628823\ \>", "\<\ 1.10086\ \>"], "Output"], Cell[OutputFormData["\<\ 1.523661613844586\ \>", "\<\ 1.52366\ \>"], "Output"], Cell[OutputFormData["\<\ 2.0000000029664085\ \>", "\<\ 2.\ \>"], "Output"], Cell[OutputFormData["\<\ 666.8045990603384\ \>", "\<\ 666.805\ \>"], "Output"], Cell[OutputFormData["\<\ Graphics[{Circle[{-1.4396521521575825, 0}, 0.2614120636865185], \ Disk[{-1.4396521521575825, 0}, 0.1], Circle[{1.4396521521575825, 0}, 0.2614120636865185], \ Disk[{1.4396521521575825, 0}, 0.1], Circle[{0, 0}, 1.178240088471064], Circle[{-1.950447010932583, \ -1.4447460337401115}, 1.2709725126384837], Disk[{-2.133094300287529, -1.9613505812130065}, 0.1], Circle[{1.950447010932583, 1.4447460337401115}, 1.2709725126384837], Disk[{2.133094300287529, 1.9613505812130065}, 0.1], {Dashing[{0.03, 0.03}], Circle[{1.950447010932583, -0.7223730168700557}, 1.2709725126384837], Disk[{2.133094300287529, -0.9806752906065033}, 0.1], Circle[{-1.950447010932583, 0.7223730168700557}, 1.2709725126384837], Disk[{-2.133094300287529, 0.9806752906065033}, 0.1]}, {Thickness[0.01], Line[{{-1.4396521521575825, 0}, {1.4396521521575825, 0}}], Line[{{1.4396521521575825, 0}, {2.133094300287529, 1.9613505812130065}}], Line[{{-1.4396521521575825, 0}, {-2.133094300287529, \ -1.9613505812130065}}], {Dashing[{0.03, 0.03}], Line[{{1.4396521521575825, 0}, {2.133094300287529, \ -0.9806752906065033}}], Line[{{-1.4396521521575825, 0}, {-2.133094300287529, \ 0.9806752906065033}}]}}}]\ \>", "\<\ -Graphics-\ \>"], "Output"], Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: .84302 MathPictureStart /Mabs { Mgmatrix idtransform Mtmatrix dtransform } bind def /Mabsadd { Mabs 3 -1 roll add 3 1 roll add exch } bind def %% Graphics %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10 scalefont setfont % Scaling calculations 0.5 0.14782 0.42151 0.14782 [ [.05654 -0.0125 -6 -9 ] [.05654 -0.0125 6 0 ] [.20436 -0.0125 -6 -9 ] [.20436 -0.0125 6 0 ] [.35218 -0.0125 -6 -9 ] [.35218 -0.0125 6 0 ] [.5 -0.0125 -3 -9 ] [.5 -0.0125 3 0 ] [.64782 -0.0125 -3 -9 ] [.64782 -0.0125 3 0 ] [.79564 -0.0125 -3 -9 ] [.79564 -0.0125 3 0 ] [.94346 -0.0125 -3 -9 ] [.94346 -0.0125 3 0 ] [ 0 0 -0.125 0 ] [-0.0125 .12587 -12 -4.5 ] [-0.0125 .12587 0 4.5 ] [-0.0125 .27369 -12 -4.5 ] [-0.0125 .27369 0 4.5 ] [-0.0125 .42151 -6 -4.5 ] [-0.0125 .42151 0 4.5 ] [-0.0125 .56933 -6 -4.5 ] [-0.0125 .56933 0 4.5 ] [-0.0125 .71715 -6 -4.5 ] [-0.0125 .71715 0 4.5 ] [ 0 0 -0.125 0 ] [ 0 .84302 .125 0 ] [ 1 0 .125 0 ] [ 0 0 0 0 ] [ 1 .84302 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 0 .5 r .25 Mabswid [ ] 0 setdash .05654 0 m .05654 .84302 L s .20436 0 m .20436 .84302 L s .35218 0 m .35218 .84302 L s .5 0 m .5 .84302 L s .64782 0 m .64782 .84302 L s .79564 0 m .79564 .84302 L s .94346 0 m .94346 .84302 L s 0 .12587 m 1 .12587 L s 0 .27369 m 1 .27369 L s 0 .42151 m 1 .42151 L s 0 .56933 m 1 .56933 L s 0 .71715 m 1 .71715 L s 0 g .05654 0 m .05654 .00625 L s [(-3)] .05654 -0.0125 0 1 Mshowa .20436 0 m .20436 .00625 L s [(-2)] .20436 -0.0125 0 1 Mshowa .35218 0 m .35218 .00625 L s [(-1)] .35218 -0.0125 0 1 Mshowa .5 0 m .5 .00625 L s [(0)] .5 -0.0125 0 1 Mshowa .64782 0 m .64782 .00625 L s [(1)] .64782 -0.0125 0 1 Mshowa .79564 0 m .79564 .00625 L s [(2)] .79564 -0.0125 0 1 Mshowa .94346 0 m .94346 .00625 L s [(3)] .94346 -0.0125 0 1 Mshowa .125 Mabswid .0861 0 m .0861 .00375 L s .11567 0 m .11567 .00375 L s .14523 0 m .14523 .00375 L s .1748 0 m .1748 .00375 L s .23392 0 m .23392 .00375 L s .26349 0 m .26349 .00375 L s .29305 0 m .29305 .00375 L s .32262 0 m .32262 .00375 L s .38174 0 m .38174 .00375 L s .41131 0 m .41131 .00375 L s .44087 0 m .44087 .00375 L s .47044 0 m .47044 .00375 L s .52956 0 m .52956 .00375 L s .55913 0 m .55913 .00375 L s .58869 0 m .58869 .00375 L s .61826 0 m .61826 .00375 L s .67738 0 m .67738 .00375 L s .70695 0 m .70695 .00375 L s .73651 0 m .73651 .00375 L s .76608 0 m .76608 .00375 L s .8252 0 m .8252 .00375 L s .85477 0 m .85477 .00375 L s .88433 0 m .88433 .00375 L s .9139 0 m .9139 .00375 L s .02698 0 m .02698 .00375 L s .97302 0 m .97302 .00375 L s .25 Mabswid 0 0 m 1 0 L s 0 .12587 m .00625 .12587 L s [(-2)] -0.0125 .12587 1 0 Mshowa 0 .27369 m .00625 .27369 L s [(-1)] -0.0125 .27369 1 0 Mshowa 0 .42151 m .00625 .42151 L s [(0)] -0.0125 .42151 1 0 Mshowa 0 .56933 m .00625 .56933 L s [(1)] -0.0125 .56933 1 0 Mshowa 0 .71715 m .00625 .71715 L s [(2)] -0.0125 .71715 1 0 Mshowa .125 Mabswid 0 .15543 m .00375 .15543 L s 0 .185 m .00375 .185 L s 0 .21456 m .00375 .21456 L s 0 .24413 m .00375 .24413 L s 0 .30325 m .00375 .30325 L s 0 .33282 m .00375 .33282 L s 0 .36238 m .00375 .36238 L s 0 .39195 m .00375 .39195 L s 0 .45107 m .00375 .45107 L s 0 .48064 m .00375 .48064 L s 0 .5102 m .00375 .5102 L s 0 .53977 m .00375 .53977 L s 0 .59889 m .00375 .59889 L s 0 .62846 m .00375 .62846 L s 0 .65802 m .00375 .65802 L s 0 .68759 m .00375 .68759 L s 0 .09631 m .00375 .09631 L s 0 .06674 m .00375 .06674 L s 0 .03718 m .00375 .03718 L s 0 .00761 m .00375 .00761 L s 0 .74671 m .00375 .74671 L s 0 .77628 m .00375 .77628 L s 0 .80584 m .00375 .80584 L s 0 .83541 m .00375 .83541 L s .25 Mabswid 0 0 m 0 .84302 L s .05654 .83677 m .05654 .84302 L s .20436 .83677 m .20436 .84302 L s .35218 .83677 m .35218 .84302 L s .5 .83677 m .5 .84302 L s .64782 .83677 m .64782 .84302 L s .79564 .83677 m .79564 .84302 L s .94346 .83677 m .94346 .84302 L s .125 Mabswid .0861 .83927 m .0861 .84302 L s .11567 .83927 m .11567 .84302 L s .14523 .83927 m .14523 .84302 L s .1748 .83927 m .1748 .84302 L s .23392 .83927 m .23392 .84302 L s .26349 .83927 m .26349 .84302 L s .29305 .83927 m .29305 .84302 L s .32262 .83927 m .32262 .84302 L s .38174 .83927 m .38174 .84302 L s .41131 .83927 m .41131 .84302 L s .44087 .83927 m .44087 .84302 L s .47044 .83927 m .47044 .84302 L s .52956 .83927 m .52956 .84302 L s .55913 .83927 m .55913 .84302 L s .58869 .83927 m .58869 .84302 L s .61826 .83927 m .61826 .84302 L s .67738 .83927 m .67738 .84302 L s .70695 .83927 m .70695 .84302 L s .73651 .83927 m .73651 .84302 L s .76608 .83927 m .76608 .84302 L s .8252 .83927 m .8252 .84302 L s .85477 .83927 m .85477 .84302 L s .88433 .83927 m .88433 .84302 L s .9139 .83927 m .9139 .84302 L s .02698 .83927 m .02698 .84302 L s .97302 .83927 m .97302 .84302 L s .25 Mabswid 0 .84302 m 1 .84302 L s .99375 .12587 m 1 .12587 L s .99375 .27369 m 1 .27369 L s .99375 .42151 m 1 .42151 L s .99375 .56933 m 1 .56933 L s .99375 .71715 m 1 .71715 L s .125 Mabswid .99625 .15543 m 1 .15543 L s .99625 .185 m 1 .185 L s .99625 .21456 m 1 .21456 L s .99625 .24413 m 1 .24413 L s .99625 .30325 m 1 .30325 L s .99625 .33282 m 1 .33282 L s .99625 .36238 m 1 .36238 L s .99625 .39195 m 1 .39195 L s .99625 .45107 m 1 .45107 L s .99625 .48064 m 1 .48064 L s .99625 .5102 m 1 .5102 L s .99625 .53977 m 1 .53977 L s .99625 .59889 m 1 .59889 L s .99625 .62846 m 1 .62846 L s .99625 .65802 m 1 .65802 L s .99625 .68759 m 1 .68759 L s .99625 .09631 m 1 .09631 L s .99625 .06674 m 1 .06674 L s .99625 .03718 m 1 .03718 L s .99625 .00761 m 1 .00761 L s .99625 .74671 m 1 .74671 L s .99625 .77628 m 1 .77628 L s .99625 .80584 m 1 .80584 L s .99625 .83541 m 1 .83541 L s .25 Mabswid 1 0 m 1 .84302 L s 0 .42151 m 1 .42151 L s .5 0 m .5 .84302 L s .5 Mabswid newpath .28719 .42151 .03864 0 365.73 arc s .28719 .42151 m .28719 .42151 .01478 0 365.73 arc F newpath .71281 .42151 .03864 0 365.73 arc s .71281 .42151 m .71281 .42151 .01478 0 365.73 arc F newpath .5 .42151 .17417 0 365.73 arc s newpath .21168 .20795 .18788 0 365.73 arc s .18469 .13158 m .18469 .13158 .01478 0 365.73 arc F newpath .78832 .63507 .18788 0 365.73 arc s .81531 .71144 m .81531 .71144 .01478 0 365.73 arc F [ .03 .03 ] 0 setdash newpath .78832 .31473 .18788 0 365.73 arc s .81531 .27655 m .81531 .27655 .01478 0 365.73 arc F newpath .21168 .52829 .18788 0 365.73 arc s .18469 .56647 m .18469 .56647 .01478 0 365.73 arc F .01 w [ ] 0 setdash .28719 .42151 m .71281 .42151 L s .71281 .42151 m .81531 .71144 L s .28719 .42151 m .18469 .13158 L s [ .03 .03 ] 0 setdash .71281 .42151 m .81531 .27655 L s .28719 .42151 m .18469 .56647 L s 0 0 m 1 0 L 1 .84302 L 0 .84302 L closepath clip newpath % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{406.813, 342.938}, ImageMargins->{{55, 0}, {0, 7}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["CompressedBitmap", "\<\ eJztncly3MgRhltqLuKmjaJ2idRKSaN9nbk7xuFwOMK++K4YO2Lm4LA9njuP 8wb2a/gNaSKbndlAf1VIoApokGaH2OpAFQr5/1WoJSsr8w9ff/nxr3/7+stP P3zd++3PX//x408//Gvv+7//fHRpfG40Ovfvo78/742K34ej0fTr8Ojf4X+K r9HoWfHf+eLyl+LXUvHrU/Frufj1QVM/67VPmu+zpl6fFLYv/0ma/JKLx/lG 49Llb+Hae779C2T9BNc+8+07kx/Pg3eUr9HTGgjmLfLapJAHSuYyZHqtqUuN UxvI7L227ahl1wMbV/Rc1sDlnWqj/rbcvM/rvZVrlO+4Uc++Lbl+Xy7+E4l/ FQiTCxXoF5zX1s7yneXrKd9q8X0u3Ga34VpKN5mSb7sqqySbhAdQ4gGUSPno yc3Lu1p8n/ewuXhZ59g8k3Vhsko5F4pc7wcq4avie02yjeXCNHFZE54PVHSR S2YW65C6palPBir/0+Jbsl2E1Cua+nCg8j9Uhq9A6jVNvT9Q+XeV4R1Ivamp dwYq/73iWwalm5B6R1NvLUL+4BB+F55M8hu6Gz3I32aAv6ktnFqIta5rA20/ NyDfPcj3ULm5OlAkOyrhHqQ+hlIuDRTJNchHmJ4o4qEisTfqEaTuQymbA0Vy WZHQTMJUGEOV/xLkIyQvIN/aQDFtKevPIPUbTR2q/JuQj5C8gnwrA8W0Aflo 7fAa8i0PFNM65HsB+d5AvqVFYArOu0hJ9A3kewv5xj0gaTMDI2XYS8j3Tsvu A0mbdrYK+V5Bvg/HLWs231H6dwOAsAL5XkO+j5V8RdoQ5Kd9nTcnSH7aW3oL +WZ35KbXaZOvb/nHkO9d6Hm2F7Rgqc9DPtL1VZ9XfBbT8SyX7pnK8in4jKXS tRHvEncvdkSMgUocquxq2XUV0tdMYfokb4nxpt//rGCI8vtbCm2bxcuOd5j9 yH915tr8YHK65Y8PuP3If6VyrTwZiZcdn/D0I//lYyZnr5Xms6ceQnzafAbB CeHScXuYvTZ9IWQhtVqkkJqBnhZfjuVHJM8LWpfVYpNtWlra0nPji+aOsFXK phZzEa6dh7JF/bJR5KLUISCm5y1rDW/JfTbkFH9vim/B9D2kvoRrfym+jze/ j34t673TbNPd/MlHplA0YNG1LbhGrUuUe5tNaiKuhIrXhKisV4pcpGgXHMIC qU3m2+OcRSK9eKSlJa1e8yaZQoTsx0k22o9b01TSU841U2kY9Dp2hz2uU+1n uGiEbt3ep9OFTlpco8EjruPvBx3JQAO7oBPjK8KeG911fZt2Q6lSCO2w1yGm NkW7QNL1rfSF+JoifgCptzSVrCNyIX6rvHoRx3fe4ohptU7Y70A+snHJysJS k/c4hYUrWrOPIfWupnaJ+J3Wex+ILymmp5B6H0ohizLWwAWNtxeKOL6/vqup hLNOZ9cIsbRq78yCJq5kadLl6ORF917r8yShix91+b9CN1lBlI2mSRBa1z4d KLAPWm2nEdi0xuoeepKAfZytsS6BebeXvVZkCyEgbgoXJ8C7P03LeBK8jgDS mxGRH7Vp08ZFCmLv3jbp0mjIy4X4k1Z5OmL/brjXFI+neYvAZvvgxx+qZK81 Xt2srhFAaa7UVbZprtPrVJ7XLi9X0/ys1ZcHnV2j8ry2eicTHe0F02wyKzrz l1BJpZLpZEHcvMrKe3/8ps3ee3g4mXHtZEL0ResrBVHc4C0+4OaqmzxIupyw NUJi3jtOOpLjQWbmczpgjUuwZNz7eIqAHR4engHrD1ijKSDBkbs+LAhOrqkf PW+hwLpqgFmANT8mUwU7Zwd7tUO8MqTR6tCLt/lRp57xlos/PH6rCbJ3JpwO mZevVwaO26swGBhuaeKk8ukPd/kaKWguZuUAnjdOZcGrIEuxg94aPAvNT9sO kAV5I0gLfJJY8GqaeZSftoUTx4IkJ2OXFkDH5Lx7Dbmx0z6FWLlEJrApLCwp lykskM0NeTXyskDWjBsdvgfNWRA7kZ3iygNITcGea2cpiFPaPB1tJZziZUXU fGTbckFT41620hB7bRLoectas4SYdu3N8xP1f1taHtnSnEzE5quLLMSstyVP KINBLDVGxtOE+KJiug6p21re7kARr6j8hJgsp8wi6hakki+lu05Z2mD32snl wb6udwQtGyupcQniM4fOscuTyU8EYadWRj7MPP7W5laUF0K1MHdgQcYImtuk sLGqNedlY1VR7kKq12fbwYz0gYfPEkBVEJzmVa55NWZxKqhkWjwRKXQt7hJu cq34pNFjx4OOU4tfxd/vim9J/Q5SRbcitfgQUvc09Y/FLzn58QXyFWboOQ4j ec8yUKVe0EolTzdUqWNtwaSifKhvsV0bzXZSLWszOCNt+AIkcyXQ6NwQ9TLn lSta0TxS7junJ89BJiIlfphqQgpt4Y/Q0JqW/uT6tRE5bZx/ew/4NaeEpCEt JbUYr0/D2aeNePhcCCny8tAyjtovkULmztU2U3w8bWby3BGvYdrQ4z2ey6P8 tM346Sk1MEuYdMvlLvjQ5YDxoCx2UHU5bHKKD20j0HBFasAU1VobZmJn18mI aF3rhg6eWE9PW2fEAQ2Tdcqead5cHDRwr7ChjYHge90SEBGkEaxbLU8/3nOu K30R4fU4QUTQvcEF0zgJ/sz0xw9acu5D6rI2frL7JKjUXoIuMIKWjwTVZEmu 6U1likAv6YPIHJRA03BC8724IXKb5u11xpOC2GvyGz8VkMst0qY+hZS5tpgi t1pe895GBzxyAYtrqW3l8wGeEbPs/TQLZnqx+J/aYC7nYfFQGO3BWBtzgGnz Qnm9u5k2+snxK1TRs1WsAe0JAX+qdlwsxalFV/BWSqnTOwMIyXrbzi+meLVo Am8qaBCeSPNYG9tBFRr6QyRbezsn7fUIkAPb5EUasfXzJRXJVi9W46M6/5SV zl8e43VklQ/bfL6gBdgM1PI8Mu7wND409NlSy88juSqIpfU+0tqpc0w7ve5d ePknmylIrBSL+FHn2DgHkrq5ZIqz4wqm4Mgfn1JW9VHlthxHnBacz8rxuqe+ rO1xr3JHeTo6fMS0UHum6GhgFo5krmLbOeksNPA83A0RtIwToxlx9UNr13Ul wowVdxVMYEkWVjoWfxeZC6/5Vpf80IrvnuImreAlZe+u5ovr6zzksFGyaPWk eyOjpi558YZ8uKocVLf3i483ZEJTXp7pU728dBnAQziQl+MWvicZaKifUT1W IbqkxBunZVsr6CZT4g1wkkLJI6VkA+7IRYnX4IfIMXMaa+3VtbKfK9oY24J7 iSsLyEh9Xi6ujBczqw26eBNpgjGLg4YA3r7YyxZNgR9ovdEd/mN93pBbK4qY XABe0/KMGeKNkHi7phS29iqtvJ6tNgF9Lus9NN2N87YRvZdaWR+87epb4OWt zTt5FfIF4yGafS51e3IbzThTQicR+/dUDNruz8XMtrYKCqO7o63HejXv/kEK 9ruKnYxjcmHfgXy7URY2lYW4yVFGFuQSBRDOxYK5XqVAxDtaE2YsRhOkFOzU H5m5KDlxplJoObzegg9zzEoGysaWjVdkGZMSsI/4IJNa77uRixmrk9tRZkj+ PjgidU+XfNgamsx5hQ8xSyPrlxTswegAcolmAbmOQBELe9oqyNmpMCOaBDsK 2NwkuzkfF8rPDc7Cu2TGZqo3FDv1Zjsqq9nvTEZZW1Ief3I3G+vEbHpIyznv kLXZgiaaUEgHE5wWmjLC1HdBzXnigtPLZLDTC9r7p51M9TIpg7pQQGrHTW16 I7NBV0Eco311TVpnyUZaprhmlMim51Nzo3MEXtptrpemiRal1c3iCiGxoz72 AobXYi7ttJdYEpbUJmbrQalEQ66T58F9noqVl6VWTYz0Drm0FWyi85UTarZe 7wyrkWfZW0yseCe8uRhttFdI40ctt7MtbVoi7S16/X8EQtTY1h5B8q6h2jgN Sd5WJlqXjNbZZslrTxq6Eum0nXvquLyLj77pJNWiLV3jdNJ+vZXi9UZU7VcV 5bj6LjTlss3yvpFlSpA3qbGK4uNA0RFbXvdcxJa1PJrR5FWH+R0c0zJO5JNz JOtBhkLtqRFDtkG6Vn5GV7yE5O+OD68v4tVyyWv6tAx8xI/KpTAjUq+HZI2f NEjhSBBt5GwzDQLEJRP2374El1+/LhLg+fL7OMz6CDYkbwMmMF021u5eaCk5 2Pnk4UNYWG/CB5kz17nLPn0sBMzW64iw7dVR2WDfG2tgpcxGcPiLH1WKsyEl XwixQR1IChu+M0rExbLKSbtmdm9KdGRjRbQnq01YSXEwn4OaMRUeIikl4DLN sO1MAc2wia5c0QaaHnilMkxp1h9/c+Or7bhQ4aRt8J5Ny03hFQBQJpDVNqTW ahMreno9phIzJaOXzVyhS1LOpgeYxVe7XG5chdCM5bk4Nit1RAdDqUvq5May RN7Fah3bKV4kyIqnZ6pZW0BM0vPt0DYpKKWUNS1PfpHLC39QpRSyg576giZk JFdKHO/whhJP5s1RBG2WTNs17SeuaGrQ/DuyEZTiQMdiBAfnCR2zzBtLYdfH Mg8Pmk0HD+fX9QspJJrliTnfodc6eBKkks8bjrz6EutzTRZpVLT5ywcB+2BL tvelc6GdXeIoJQi2rXSklM3iFxkY2EtGdohd8tE8jkQcsZlUkKTyglDkZzPU 3W1BgHcJnxss9VC2V2BejwnxDc3XJua7FzHdS146c7NgY73tkRILtzSVjOqy slDh2vuqNMdu+kxz/EWdoMVCt26Cpvq5WLBONzcLQTdK5nhrHGLh/jGPs/cW f8UnaLKdS5231BcbNKGI87KrqWTGk4sFWw/tQim5jfqtndse53W4977mG5V9 g7oaxDiUkwgwHzV7mQmgLXvz3kMdh03hiRRqpmbdSR1HmzM2QY7kVnLGn5sj umazijhHZqNNM2Ba/OR6kWzAewSleE+seZcoh7OzSupqJDtZJN+GO2gBRSeY cnFlXvy8XKWEJo8tpCtazHGINVJE0CnKLoNL2GSCnBF0w9rcwti6FNPiE2Hk 4JaO6HZJmPksfhpsAgf41ESSTK1hqnoaxL1hXTonSW7dj5BEhqVNo6aTStCm /8SPLaXivZafqZSQVjahfY7VftCMJI6bZ0f/gu4KbM4a3FTskgYaLPZzEFJ9 tUajmTmwNB6x1KdeZElTaUBxaYwzMkRGCCltYzSaOV8trVB2emmJkmufYrHw rbzi80AbPZmGDwExOZqjZ9Qhnn4eam0T4tY7fzbPWRROeakfRdF5d4m7jDpM Tue8wfoqOKXd0tosl0FBcpDGyr2zRuXT6/RcmoN+UuxPtI4Je4rtSTLeisfZ 5ijn+7c6vC0tkMLLKV/M3eJDi8yqOJPPwRNtr2QA06Vl3Uxw6Ig3NSu32n5m ZXqqdBAK6sGuZUBRfOSuDypaqX4kA/kofa9Nkl7+x72DmX4qLnHjfQZpIejN eadgqUc391ZxS9cKmIQXaOZDJx5SoFaHsRBo6zg6By3ZvEo42hgO4qwc22yO zjrhNuji/ulJD01rGEL3RtHRlNI6TI/znuLjV1XH3c+nYHqtmGhdZN2nF1Mb JXPc/TypxknFEUTnWA54cM6NRVdbNM5G/unjwSfiE0/v+bynyn3YsVk+5GRO HgywF6y3OPKX+iBSeu1rBYTx8nh9pQa0dypaHfeLNGr3pmDyIv9G76DnepDX j3h1p1XzLLTsWjCWnIOPF8oHqVH3teV7HSB2yYe02pVybcelMRcJXj6ezb+C o6BA5nHVS0+bI/heeqRpr5nk9dLQ4Wbat2FlJeoVh0WJtO6Ncv3HpTEFejB4 XGU6HpfA/Gt73a224cO7ndGcD68SJyVc2nNtM16O6jwJxSP8dseWWU0RR6Rp j68kuuTIGwPVy5ZIuNmELTLPf+i8Fo/cXJGqkQfoOmde8ejAq1LWJHWySiti FI2O4/r+qfjliAS8AqnfFt9y6feQuqap+WtW3gMx4/YOLHFHmfFrNHunNv9C Ofe6x69z4+flzXucufmALAY6jmjgJAFxJBKMc3IUPx7diBk56umNKXRTn0tG XbcVJ82FaSaSnxnvIu+Vtltaa9KTbyt2skO+BddI5xZkoXKqMy5LcOaVQalN zzOXtjehZLMIJ8MdQvxSyxsq4l3FdB1KNqv/+LGGk4SYDmNYCA06V3ZFU2kO QWv3l9rSqT8iqeqCpOdmgd5ZaQGVU2iV2pZUGm+JhVfK20li4Y7iDI4vkkom T6eHhb3ie7u4ch9SSSfbfKRx4ay8lXnQkWGfN8zJa5UqBae5hshfn95Dg/0i 7rIFe0+YNkJcOYXeXIK6oE8p55/7YOGN1vtJYsF7Yn5YLKS4eUhhwRM9UFmo OC/Nz0JT1yF9IH+r9T905F4vMt0hbxO9tYGzInokGVA1Ath5IOUFArRoymcA EwBWraObvJTJ6PLHMZfkRWF6rxKkYKoLOmxGLycHU5sw197pQDK6ydtdbydy koB90Go7jcCmNbb4pui3tvUCbGS01T1Al1H0GbY6u8I+X7Hik/KKfdJaWiSS 6f0nHckI8rVEErTJ9U+NUqKQkM0/xRLxYvqsUnnPeporhlxRa+hoCcWe6Q5T m/Dp8dVwZbD0hS7yqoziBu+5AMYVHfFxIyXOd/ywwhk6D7q4qrZjdNIy6UBN LrW9eRfvDZPvxFAugLSH/Ox0QiX1PHnQ9EL1RuNudOit7/2mXMCKXwsD9lzf UjLmOhUQ9xViXExq2naN7FxqDpjWT0XmTFyCU08C1rwm0iHqXV6IZr2XVouP tRw6ZbtQiG0Mlcjs0lLJbLsP2DNnvMsJuSwRyaKOQpg/UHri7YNGXrsW7tJ8 Z/NzIo+HbN9r/GZkQy6nmygsedaQ7JLNG7W5A5Ri+Usrv1yBtM0bJvkJ9aL0 Ohtj3cbRh1ZH/URR7xpi8QlCzBXnm3jd1va754TojW6cFaJ5lEyJqmazkEeL BFt3Io2GYW/wtPiZThFdXHc9jYJYMAE0YV7WFvDboiwB8U8pdZJvsoo6OFqP z1/7TfEtB4beQ2qR0PSoT+6llTeoNasTR9z7150d7iBuoRevNzIye0E8+ryG hDq8eUIReiF6vSEE3USJqdSb48ZfGfGl1DrvAENCZG6n32qrO5hF08a/QwqY 5oEYqs84BRDK/XUu58pe+b1uBIcqv9dRoHeP82TKX5GvR/m9DgyHKn/MCSEZ dbm8rY31TRiYpCJTR+dBSObaXU+7vUvPgl5x41vQB3Ww+hY3bgVg96Y42swt c9waw+7N5Rh1UfJ7+/XFyU92PXZvrsCiXco/acNs35HLIXQ38kfM3U4KhDpz S0vzLlkW3ZSKP/MNTHszcUdraRGIUqT2mo15FUm5IgalYCINCG0i5IqeMSRM ucKm9IEpHrne8uWKndMHJjs2TUi8u+JDQOK1APOaBg0BU9zu65GmdhlEMUV+ r1kXuQgbKqYtZZ2Q2L5el0EMU+SPuwS0fLYJe32gSMhdn6WSo6mhIjF3ORQN 2dxj0ab/4uZiO8r/HqSayy8y0hjCvOuGskUBic2q6D6kDqHVkGOxioTjUO0M QX5yDlbppyoxE4cl/65KSH29xRAl694hyC+8iq6MDpWYY6YXA5Vfxl/Zjp5d Y08/smVPS44hiC46P5GQtJhDkNBb3pmsPcna27DuLc8/hJP+2KtnzpOvrubj vnRJhXSW7yxft/nmjsGJ3vnXyYXRuf8Bv+IIag==\ \>"], ImageRangeCache->{{{87.125, 411.563}, {341.125, 67.8125}} -> {-5.57211, \ -1.61887, 0.0174104, 0.0174104}}], Cell[OutputFormData["\<\ Graphics[{Circle[{-1.4396521521575825, 0}, 0.2614120636865185], \ Disk[{-1.4396521521575825, 0}, 0.1], Circle[{1.4396521521575825, 0}, 0.2614120636865185], \ Disk[{1.4396521521575825, 0}, 0.1], Circle[{0, 0}, 1.178240088471064], Circle[{-1.950447010932583, \ -1.4447460337401115}, 1.2709725126384837], Disk[{-2.133094300287529, -1.9613505812130065}, 0.1], Circle[{1.950447010932583, 1.4447460337401115}, 1.2709725126384837], Disk[{2.133094300287529, 1.9613505812130065}, 0.1], {Dashing[{0.03, 0.03}], Circle[{1.950447010932583, -0.7223730168700557}, 1.2709725126384837], Disk[{2.133094300287529, -0.9806752906065033}, 0.1], Circle[{-1.950447010932583, 0.7223730168700557}, 1.2709725126384837], Disk[{-2.133094300287529, 0.9806752906065033}, 0.1]}, {Thickness[0.01], Line[{{-1.4396521521575825, 0}, {1.4396521521575825, 0}}], Line[{{1.4396521521575825, 0}, {2.133094300287529, 1.9613505812130065}}], Line[{{-1.4396521521575825, 0}, {-2.133094300287529, \ -1.9613505812130065}}], {Dashing[{0.03, 0.03}], Line[{{1.4396521521575825, 0}, {2.133094300287529, \ -0.9806752906065033}}], Line[{{-1.4396521521575825, 0}, {-2.133094300287529, \ 0.9806752906065033}}]}}}, {AspectRatio -> Automatic, PlotRange -> All, Axes -> True, GridLines -> \ Automatic, Frame -> True}]\ \>", "\<\ -Graphics-\ \>"], "Output"] }, Open ]] }, WindowSize->{808, 740}, WindowMargins->{{0, Automatic}, {Automatic, 0}}, PrintingCopies->1, PrintingPageRange->{Automatic, Automatic}, FrontEndVersion->"8.0 for Microsoft Windows (64-bit) (October 6, 2011)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[579, 22, 23326, 644, 1372, "Input"], Cell[23908, 668, 249, 6, 35, "Output"], Cell[24160, 676, 80, 4, 35, "Output"], Cell[24243, 682, 78, 4, 35, "Output"], Cell[24324, 688, 78, 4, 35, "Output"], Cell[24405, 694, 79, 4, 35, "Output"], Cell[24487, 700, 78, 4, 35, "Output"], Cell[24568, 706, 74, 4, 35, "Output"], Cell[24645, 712, 78, 4, 35, "Output"], Cell[24726, 718, 1275, 24, 35, "Output"], Cell[26004, 744, 14875, 709, 358, 6761, 573, "GraphicsData", "PostScript", \ "Graphics"], Cell[40882, 1455, 1378, 26, 35, "Output"] }, Open ]] } ] *) (* End of internal cache information *)