(* 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[     28542,        806]
NotebookOptionsPosition[     27376,        767]
NotebookOutlinePosition[     27862,        785]
CellTagsIndexPosition[     27819,        782]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{"(*", " ", 
  RowBox[{
   RowBox[{"C2H6", " ", "Molecule", " ", "after", " ", "Kimball"}], ",", 
   RowBox[{"non", " ", "overlapping"}], ",", 
   RowBox[{
    RowBox[{"parameterized", " ", "with", " ", "G2"}], "//", 
    RowBox[{"6", "-", 
     RowBox[{"311", "g"}]}]}]}], " ", "*)"}]], "Input",
 CellChangeTimes->{
  3.568565138949461*^9, {3.568565311002163*^9, 3.5685653384114113`*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"Clear", "[", 
    RowBox[{
    "k1", ",", "k2", ",", "k3", ",", "s1", ",", "s2", ",", "s3", ",", "Vee", 
     ",", "Vne", ",", "Vnn", ",", "P", ",", "Q", ",", "R", ",", "X", ",", 
     "func"}], "]"}], ";", " ", 
   RowBox[{"n", "=", "2."}], ";"}], "\[IndentingNewLine]", 
  RowBox[{"(*", " ", 
   RowBox[{"Term", " ", "from", " ", 
    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]"}], "Input",
 CellChangeTimes->{
  3.568565138949461*^9, {3.568565311002163*^9, 3.5685653384114113`*^9}, {
   3.5685654981712923`*^9, 3.5685655172969255`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"Vee", "=", 
    RowBox[{
     RowBox[{"3.0", "*", "n", "*", 
      RowBox[{"s1", "/", "P"}]}], "+", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"n", " ", "-", " ", "1.0"}], ")"}], "*", "3.0", "*", 
      RowBox[{"s2", "/", " ", "Q"}]}], "+", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"n", " ", "+", " ", "1.0"}], ")"}], "*", "6.0", "*", 
      RowBox[{"s3", "/", " ", "R"}]}], " ", "+", " ", 
     RowBox[{"8.0", "/", " ", 
      RowBox[{"Sqrt", "[", "q2", "]"}]}], "+", " ", 
     RowBox[{"24.0", "/", " ", 
      RowBox[{"Sqrt", "[", 
       RowBox[{
        RowBox[{"zd", "*", "qr"}], "+", "q2", " ", "+", " ", "r2"}], "]"}]}], 
     "+", " ", 
     RowBox[{"24.0", "/", " ", 
      RowBox[{"Sqrt", "[", "r2", "]"}]}], " ", "+", " ", 
     RowBox[{"4.0", "/", 
      RowBox[{"Sqrt", "[", 
       RowBox[{"4.0", "*", "q2"}], "]"}]}], " ", "+", " ", 
     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"}]}], "]"}]}]}]}], ";"}], 
  "\[IndentingNewLine]"}]], "Input",
 CellChangeTimes->{
  3.568565138949461*^9, {3.568565311002163*^9, 3.568565344401822*^9}, {
   3.5685656015838737`*^9, 3.568565613767495*^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", "/", " ", 
       RowBox[{"Sqrt", "[", "q2", "]"}]}]}], " ", "-", " ", 
     RowBox[{"12.0", "/", " ", 
      RowBox[{"Sqrt", "[", 
       RowBox[{
        RowBox[{"zd", "*", "qrp"}], " ", "+", " ", "r2p2", " ", "+", " ", 
        "q2"}], "]"}]}], "-", " ", 
     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"}], "]"}]}]}], "-", " ", 
     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.568565138949461*^9, {3.568565311002163*^9, 3.568565349737031*^9}, {
   3.5685655577789965`*^9, 3.568565584923044*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"Vnn", " ", "=", 
    RowBox[{
     RowBox[{"6.0", "*", 
      RowBox[{"z", "/", " ", 
       RowBox[{"Sqrt", "[", "r2p2", "]"}]}]}], " ", "+", " ", 
     RowBox[{"z", "*", 
      RowBox[{"z", "/", " ", 
       RowBox[{"Sqrt", "[", 
        RowBox[{"4.0", "*", "q2"}], "]"}]}]}], "+", " ", 
     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.568565138949461*^9, {3.568565311002163*^9, 3.568565349737031*^9}, {
   3.568565627433119*^9, 3.5685656328619285`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  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[{
   RowBox[{"func", "=", 
    RowBox[{"func", " ", "/.", " ", "c"}]}], ";"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{"N", "[", 
  RowBox[{
   RowBox[{"t", " ", "=", " ", 
    RowBox[{"FindMinimum", "[", 
     RowBox[{"func", ",", 
      RowBox[{"{", 
       RowBox[{"P", ",", "0.25"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"Q", ",", "1."}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"R", ",", "1.2"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"X", ",", "1.35"}], "}"}]}], "]"}]}], ",", "10"}], 
  "]"}], "\[IndentingNewLine]"}], "Input",
 CellChangeTimes->{3.568565138949461*^9, 3.568565311002163*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", "79.63122262216285`"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"P", "\[Rule]", "0.261412063549814`"}], ",", 
     RowBox[{"Q", "\[Rule]", "1.178240111695778`"}], ",", 
     RowBox[{"R", "\[Rule]", "1.270972496347005`"}], ",", 
     RowBox[{"X", "\[Rule]", "1.3575746415402064`"}]}], "}"}]}], 
  "}"}]], "Output",
 CellChangeTimes->{3.5685656629231815`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 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[{
  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[{
  RowBox[{"dcc", "=", 
   RowBox[{"cc", "*", "0.529177"}]}], "  ", 
  RowBox[{"(*", " ", 
   RowBox[{"C", "-", 
    RowBox[{"C", " ", "Abstand"}]}], " ", "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  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[{"-", 
      RowBox[{"t", "[", 
       RowBox[{"[", "1", "]"}], "]"}]}], " ", "-", " ", 
     RowBox[{"n", "*", "37.784301"}], " ", "-", " ", "n", " ", "-", " ", 
     "1.0"}], ")"}], "*", "627.50956"}], " ", 
  RowBox[{"(*", 
   RowBox[{"Atomisierungsenergie", " ", "0", " ", "K", " ", "in", " ", 
    RowBox[{"kcal", "/", "mol"}]}], " ", 
   "*)"}]}], "\[IndentingNewLine]"}], "Input",
 CellChangeTimes->{3.568565138949461*^9, 3.568565311002163*^9}],

Cell[BoxData[
 RowBox[{"-", "271.1254029126776`"}]], "Output",
 CellChangeTimes->{3.5685656629387817`*^9}],

Cell[BoxData["69.66048829986788`"], "Output",
 CellChangeTimes->{3.5685656629387817`*^9}],

Cell[BoxData["42.202469363840784`"], "Output",
 CellChangeTimes->{3.5685656629387817`*^9}],

Cell[BoxData["1.100860893746612`"], "Output",
 CellChangeTimes->{3.5685656629387817`*^9}],

Cell[BoxData["1.5236616382798733`"], "Output",
 CellChangeTimes->{3.5685656629387817`*^9}],

Cell[BoxData["1.9999999999416902`"], "Output",
 CellChangeTimes->{3.5685656629387817`*^9}],

Cell[BoxData["666.8045990603384`"], "Output",
 CellChangeTimes->{3.5685656629387817`*^9}]
}, Open  ]],

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]"}]], "Input",
 CellChangeTimes->{{3.568565138949461*^9, 3.568565161803501*^9}, {
  3.568565383589091*^9, 3.5685653861630955`*^9}, {3.5685654285327697`*^9, 
  3.568565432042776*^9}, {3.5685654814324627`*^9, 3.5685654849424686`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Show", "[", 
  RowBox[{"plot1", ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"AspectRatio", " ", "->", " ", "Automatic"}], ",", 
     RowBox[{"PlotRange", " ", "->", " ", "All"}], ",", 
     RowBox[{"Axes", " ", "->", " ", "True"}], ",", 
     RowBox[{"GridLines", " ", "->", " ", "Automatic"}], ",", 
     RowBox[{"Frame", " ", "->", " ", "True"}]}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{
  3.568565138949461*^9, 3.5685651739871225`*^9, {3.568565395803912*^9, 
   3.5685653987523174`*^9}}],

Cell[BoxData[
 GraphicsBox[{CircleBox[{-1.439652175245592, 0}, 0.261412063549814], 
   DiskBox[{-1.439652175245592, 0}, 0.1], 
   CircleBox[{1.439652175245592, 0}, 0.261412063549814], 
   DiskBox[{1.439652175245592, 0}, 0.1], CircleBox[{0, 0}, 1.178240111695778],
    CircleBox[{-1.9504470285445317`, -1.444746018251472}, 1.270972496347005], 
   DiskBox[{-2.1330943151134822`, -1.9613505578443824`}, 0.1], 
   CircleBox[{1.9504470285445317`, 1.444746018251472}, 1.270972496347005], 
   DiskBox[{2.1330943151134822`, 1.9613505578443824`}, 0.1], 
   {Dashing[{0.03, 0.03}], 
    CircleBox[{1.9504470285445317`, -0.722373009125736}, 1.270972496347005], 
    DiskBox[{2.1330943151134822`, -0.9806752789221912}, 0.1], 
    CircleBox[{-1.9504470285445317`, 0.722373009125736}, 1.270972496347005], 
    DiskBox[{-2.1330943151134822`, 0.9806752789221912}, 0.1]}, 
   {Thickness[0.01], 
    LineBox[{{-1.439652175245592, 0}, {1.439652175245592, 0}}], 
    LineBox[{{1.439652175245592, 0}, {2.1330943151134822`, 
      1.9613505578443824`}}], 
    LineBox[{{-1.439652175245592, 
       0}, {-2.1330943151134822`, -1.9613505578443824`}}], 
    {Dashing[{0.03, 0.03}], 
     LineBox[{{1.439652175245592, 0}, {
       2.1330943151134822`, -0.9806752789221912}}], 
     LineBox[{{-1.439652175245592, 0}, {-2.1330943151134822`, 
       0.9806752789221912}}]}}},
  AspectRatio->Automatic,
  Axes->True,
  Frame->True,
  GridLines->Automatic,
  PlotRange->All]], "Output",
 CellChangeTimes->{3.5685656629543815`*^9}]
}, Open  ]]
},
WindowSize->{1163, 740},
WindowMargins->{{327, Automatic}, {Automatic, 0}},
PrintingCopies->1,
PrintingPageRange->{Automatic, Automatic},
PrivateNotebookOptions->{"VersionedStylesheet"->{"Default.nb"[8.] -> False}},
FrontEndVersion->"9.0 for Microsoft Windows (64-bit) (November 20, 2012)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[557, 20, 411, 10, 31, "Input"],
Cell[971, 32, 2164, 59, 132, "Input"],
Cell[3138, 93, 1886, 50, 92, "Input"],
Cell[5027, 145, 2843, 74, 112, "Input"],
Cell[7873, 221, 1353, 36, 72, "Input"],
Cell[CellGroupData[{
Cell[9251, 261, 1419, 39, 152, "Input"],
Cell[10673, 302, 422, 11, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[11132, 318, 2225, 68, 232, "Input"],
Cell[13360, 388, 106, 2, 31, "Output"],
Cell[13469, 392, 89, 1, 31, "Output"],
Cell[13561, 395, 90, 1, 31, "Output"],
Cell[13654, 398, 89, 1, 31, "Output"],
Cell[13746, 401, 90, 1, 31, "Output"],
Cell[13839, 404, 90, 1, 31, "Output"],
Cell[13932, 407, 89, 1, 31, "Output"]
}, Open  ]],
Cell[14036, 411, 11271, 305, 272, "Input"],
Cell[CellGroupData[{
Cell[25332, 720, 526, 12, 31, "Input"],
Cell[25861, 734, 1499, 30, 355, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)