(* 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[     44731,       1047]
NotebookOptionsPosition[     43501,       1007]
NotebookOutlinePosition[     43985,       1025]
CellTagsIndexPosition[     43942,       1022]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 StyleBox[
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{"C2H6", " ", "Molecule", " ", "after", " ", "Kimball"}], ",", " ", 
    RowBox[{"non", " ", "overlapping"}], ",", 
    RowBox[{
     RowBox[{"parameterized", " ", "with", " ", "G2"}], "//", 
     RowBox[{"6", "-", 
      RowBox[{
      "311", "g", " ", "\[IndentingNewLine]", "kin", " ", "param", " ", "1.0",
        " ", "sigma", " ", "varied"}]}]}], ",", " ", 
    RowBox[{
     RowBox[{"ethaneN1_new", ".", "nb"}], " ", "better", " ", "overall", " ", 
     "with", " ", "kin", " ", "params", " ", "and", " ", "sigma", " ", 
     "varied"}]}], " ", "*)"}], "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.578932625327005*^9, 
   3.5789326822047043`*^9}, {3.578932712624758*^9, 3.5789327663512526`*^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[{"c", " ", "=", " ", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"k1", " ", "->", " ", "1.0"}], ",", 
     RowBox[{"k2", " ", "->", " ", "1.0"}], ",", 
     RowBox[{"k3", " ", "->", " ", "1.0"}], ",", 
     RowBox[{"s1", " ", "->", " ", "0.273355"}], ",", 
     RowBox[{"s2", " ", "->", " ", "0.795"}], ",", 
     RowBox[{"s3", " ", "->", " ", "0.59"}]}], "}"}]}], 
  ";"}], "\[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.578904560584473*^9, 3.5789045711612916`*^9}, {
   3.5789046034221487`*^9, 3.5789046541846375`*^9}, 3.5789046872254953`*^9, {
   3.578904770233241*^9, 3.5789047765044527`*^9}, 3.5789048720234203`*^9, {
   3.578904916296298*^9, 3.5789049394623384`*^9}, {3.5789049867460213`*^9, 
   3.578904994842436*^9}, {3.578905066602562*^9, 3.578905103496627*^9}, {
   3.578905768057794*^9, 3.578905768556995*^9}, {3.578905805185859*^9, 
   3.5789058861656017`*^9}, {3.5789059186604586`*^9, 3.5789060027758064`*^9}, 
   3.5789060361442647`*^9, 3.5789060740211315`*^9, {3.578926941454946*^9, 
   3.578927025164693*^9}, {3.57892706903197*^9, 3.578927093414813*^9}, {
   3.57892712590967*^9, 3.578927145253704*^9}, {3.578927219353834*^9, 
   3.5789272206486363`*^9}, {3.5789272583539023`*^9, 3.57892729671437*^9}, {
   3.5789276082157173`*^9, 3.578927658791006*^9}, {3.578931850161643*^9, 
   3.5789318919697165`*^9}, {3.5789319243709736`*^9, 3.578932355992532*^9}, {
   3.578932397769405*^9, 3.5789325450648637`*^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.254", ",", "0.255"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"Q", ",", "1.176", ",", "1.177"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"R", ",", "1.257", ",", "1.258"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"X", ",", "1.358", ",", "1.36"}], "}"}], ",", 
    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}, {
   3.578927348709261*^9, 3.5789274057429614`*^9}, {3.5789274548830476`*^9, 
   3.578927458471054*^9}, {3.5789275258007727`*^9, 3.5789275845192757`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", "79.65021660708372`"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"P", "\[Rule]", "0.2544847392835131`"}], ",", 
     RowBox[{"Q", "\[Rule]", "1.1764241350529565`"}], ",", 
     RowBox[{"R", "\[Rule]", "1.2573110439748147`"}], ",", 
     RowBox[{"X", "\[Rule]", "1.3593698868015276`"}]}], "}"}]}], 
  "}"}]], "Output",
 CellChangeTimes->{
  3.578903430222088*^9, 3.578904578228104*^9, 3.5789046267441893`*^9, 
   3.5789046603622484`*^9, 3.5789046927947054`*^9, 3.5789047396727877`*^9, 
   3.5789047810908604`*^9, 3.578904877187029*^9, {3.578904920383505*^9, 
   3.5789049445323477`*^9}, 3.578904999382044*^9, {3.578905078536583*^9, 
   3.578905107693034*^9}, 3.578905772550602*^9, 3.5789058096162667`*^9, {
   3.578905841019122*^9, 3.57890585702475*^9}, 3.5789058907208095`*^9, {
   3.578905924682069*^9, 3.5789059489089117`*^9}, {3.57890598230857*^9, 
   3.5789060082670155`*^9}, 3.5789060404342723`*^9, 3.578906078389139*^9, 
   3.578926946774555*^9, {3.578926976867008*^9, 3.578927029579501*^9}, {
   3.5789270807475905`*^9, 3.57892709762682*^9}, {3.578927132336881*^9, 
   3.578927150635713*^9}, {3.5789272032858057`*^9, 3.5789272251882443`*^9}, {
   3.57892726245671*^9, 3.5789273016907787`*^9}, 3.578927412981374*^9, 
   3.57892747332228*^9, {3.5789275896984844`*^9, 3.5789276649530168`*^9}, {
   3.578931858008457*^9, 3.578931897242526*^9}, {3.578931930454984*^9, 
   3.5789325495264716`*^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.788991"}], " ", "-", 
     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.5789047231523585`*^9, 
   3.5789047343843784`*^9}, {3.578927185657775*^9, 3.5789271958913927`*^9}}],

Cell[BoxData["79.65023955900203`"], "Output",
 CellChangeTimes->{
  3.578903430222088*^9, 3.578904578228104*^9, 3.5789046267441893`*^9, 
   3.5789046603778486`*^9, 3.5789046928103056`*^9, 3.578904739688388*^9, 
   3.5789047811064606`*^9, 3.578904877202629*^9, {3.578904920383505*^9, 
   3.578904944547948*^9}, 3.578904999397644*^9, {3.578905078552183*^9, 
   3.578905107708634*^9}, 3.5789057725662017`*^9, 3.5789058096318665`*^9, {
   3.5789058410347223`*^9, 3.57890585704035*^9}, 3.578905890736409*^9, {
   3.578905924697669*^9, 3.578905948924512*^9}, {3.57890598232417*^9, 
   3.5789060082670155`*^9}, 3.5789060404342723`*^9, 3.578906078404739*^9, 
   3.5789269467901554`*^9, {3.5789269768826084`*^9, 3.578927029595101*^9}, {
   3.5789270807631903`*^9, 3.5789270976424203`*^9}, {3.578927132352481*^9, 
   3.5789271506513133`*^9}, {3.5789272032858057`*^9, 3.578927225203844*^9}, {
   3.57892726247231*^9, 3.578927301706379*^9}, 3.578927412981374*^9, 
   3.57892747333788*^9, {3.5789275896984844`*^9, 3.5789276649686165`*^9}, {
   3.578931858024057*^9, 3.5789318972581263`*^9}, {3.578931930470584*^9, 
   3.5789325495420713`*^9}}],

Cell[BoxData[
 RowBox[{"-", "276.27820074959703`"}]], "Output",
 CellChangeTimes->{
  3.578903430222088*^9, 3.578904578228104*^9, 3.5789046267441893`*^9, 
   3.5789046603778486`*^9, 3.5789046928103056`*^9, 3.578904739688388*^9, 
   3.5789047811064606`*^9, 3.578904877202629*^9, {3.578904920383505*^9, 
   3.578904944547948*^9}, 3.578904999397644*^9, {3.578905078552183*^9, 
   3.578905107708634*^9}, 3.5789057725662017`*^9, 3.5789058096318665`*^9, {
   3.5789058410347223`*^9, 3.57890585704035*^9}, 3.578905890736409*^9, {
   3.578905924697669*^9, 3.578905948924512*^9}, {3.57890598232417*^9, 
   3.5789060082670155`*^9}, 3.5789060404342723`*^9, 3.578906078404739*^9, 
   3.5789269467901554`*^9, {3.5789269768826084`*^9, 3.578927029595101*^9}, {
   3.5789270807631903`*^9, 3.5789270976424203`*^9}, {3.578927132352481*^9, 
   3.5789271506513133`*^9}, {3.5789272032858057`*^9, 3.578927225203844*^9}, {
   3.57892726247231*^9, 3.578927301706379*^9}, 3.578927412981374*^9, 
   3.57892747333788*^9, {3.5789275896984844`*^9, 3.5789276649686165`*^9}, {
   3.578931858024057*^9, 3.5789318972581263`*^9}, {3.578931930470584*^9, 
   3.5789325495420713`*^9}}],

Cell[BoxData["74.37359080001491`"], "Output",
 CellChangeTimes->{
  3.578903430222088*^9, 3.578904578228104*^9, 3.5789046267441893`*^9, 
   3.5789046603778486`*^9, 3.5789046928103056`*^9, 3.578904739688388*^9, 
   3.5789047811064606`*^9, 3.578904877202629*^9, {3.578904920383505*^9, 
   3.578904944547948*^9}, 3.578904999397644*^9, {3.578905078552183*^9, 
   3.578905107708634*^9}, 3.5789057725662017`*^9, 3.5789058096318665`*^9, {
   3.5789058410347223`*^9, 3.57890585704035*^9}, 3.578905890736409*^9, {
   3.578905924697669*^9, 3.578905948924512*^9}, {3.57890598232417*^9, 
   3.5789060082670155`*^9}, 3.5789060404342723`*^9, 3.578906078404739*^9, 
   3.5789269467901554`*^9, {3.5789269768826084`*^9, 3.578927029595101*^9}, {
   3.5789270807631903`*^9, 3.5789270976424203`*^9}, {3.578927132352481*^9, 
   3.5789271506513133`*^9}, {3.5789272032858057`*^9, 3.578927225203844*^9}, {
   3.57892726247231*^9, 3.578927301706379*^9}, 3.578927412981374*^9, 
   3.57892747333788*^9, {3.5789275896984844`*^9, 3.5789276649686165`*^9}, {
   3.578931858024057*^9, 3.5789318972581263`*^9}, {3.578931930470584*^9, 
   3.5789325495420713`*^9}}],

Cell[BoxData["42.60415378349641`"], "Output",
 CellChangeTimes->{
  3.578903430222088*^9, 3.578904578228104*^9, 3.5789046267441893`*^9, 
   3.5789046603778486`*^9, 3.5789046928103056`*^9, 3.578904739688388*^9, 
   3.5789047811064606`*^9, 3.578904877202629*^9, {3.578904920383505*^9, 
   3.578904944547948*^9}, 3.578904999397644*^9, {3.578905078552183*^9, 
   3.578905107708634*^9}, 3.5789057725662017`*^9, 3.5789058096318665`*^9, {
   3.5789058410347223`*^9, 3.57890585704035*^9}, 3.578905890736409*^9, {
   3.578905924697669*^9, 3.578905948924512*^9}, {3.57890598232417*^9, 
   3.5789060082670155`*^9}, 3.5789060404342723`*^9, 3.578906078404739*^9, 
   3.5789269467901554`*^9, {3.5789269768826084`*^9, 3.578927029595101*^9}, {
   3.5789270807631903`*^9, 3.5789270976424203`*^9}, {3.578927132352481*^9, 
   3.5789271506513133`*^9}, {3.5789272032858057`*^9, 3.578927225203844*^9}, {
   3.57892726247231*^9, 3.578927301706379*^9}, 3.578927412981374*^9, 
   3.57892747333788*^9, {3.5789275896984844`*^9, 3.5789276649686165`*^9}, {
   3.578931858024057*^9, 3.5789318972581263`*^9}, {3.578931930470584*^9, 
   3.5789325495576715`*^9}}],

Cell[BoxData["1.0875061824676495`"], "Output",
 CellChangeTimes->{
  3.578903430222088*^9, 3.578904578228104*^9, 3.5789046267441893`*^9, 
   3.5789046603778486`*^9, 3.5789046928103056`*^9, 3.578904739688388*^9, 
   3.5789047811064606`*^9, 3.578904877202629*^9, {3.578904920383505*^9, 
   3.578904944547948*^9}, 3.578904999397644*^9, {3.578905078552183*^9, 
   3.578905107708634*^9}, 3.5789057725662017`*^9, 3.5789058096318665`*^9, {
   3.5789058410347223`*^9, 3.57890585704035*^9}, 3.578905890736409*^9, {
   3.578905924697669*^9, 3.578905948924512*^9}, {3.57890598232417*^9, 
   3.5789060082670155`*^9}, 3.5789060404342723`*^9, 3.578906078404739*^9, 
   3.5789269467901554`*^9, {3.5789269768826084`*^9, 3.578927029595101*^9}, {
   3.5789270807631903`*^9, 3.5789270976424203`*^9}, {3.578927132352481*^9, 
   3.5789271506513133`*^9}, {3.5789272032858057`*^9, 3.578927225203844*^9}, {
   3.57892726247231*^9, 3.578927301706379*^9}, 3.578927412981374*^9, 
   3.57892747333788*^9, {3.5789275896984844`*^9, 3.5789276649686165`*^9}, {
   3.578931858024057*^9, 3.5789318972581263`*^9}, {3.578931930470584*^9, 
   3.5789325495576715`*^9}}],

Cell[BoxData["1.5144081307894999`"], "Output",
 CellChangeTimes->{
  3.578903430222088*^9, 3.578904578228104*^9, 3.5789046267441893`*^9, 
   3.5789046603778486`*^9, 3.5789046928103056`*^9, 3.578904739688388*^9, 
   3.5789047811064606`*^9, 3.578904877202629*^9, {3.578904920383505*^9, 
   3.578904944547948*^9}, 3.578904999397644*^9, {3.578905078552183*^9, 
   3.578905107708634*^9}, 3.5789057725662017`*^9, 3.5789058096318665`*^9, {
   3.5789058410347223`*^9, 3.57890585704035*^9}, 3.578905890736409*^9, {
   3.578905924697669*^9, 3.578905948924512*^9}, {3.57890598232417*^9, 
   3.5789060082670155`*^9}, 3.5789060404342723`*^9, 3.578906078404739*^9, 
   3.5789269467901554`*^9, {3.5789269768826084`*^9, 3.578927029595101*^9}, {
   3.5789270807631903`*^9, 3.5789270976424203`*^9}, {3.578927132352481*^9, 
   3.5789271506513133`*^9}, {3.5789272032858057`*^9, 3.578927225203844*^9}, {
   3.57892726247231*^9, 3.578927301706379*^9}, 3.578927412981374*^9, 
   3.57892747333788*^9, {3.5789275896984844`*^9, 3.5789276649686165`*^9}, {
   3.578931858024057*^9, 3.5789318972581263`*^9}, {3.578931930470584*^9, 
   3.5789325495576715`*^9}}],

Cell[BoxData["1.9999997118411892`"], "Output",
 CellChangeTimes->{
  3.578903430222088*^9, 3.578904578228104*^9, 3.5789046267441893`*^9, 
   3.5789046603778486`*^9, 3.5789046928103056`*^9, 3.578904739688388*^9, 
   3.5789047811064606`*^9, 3.578904877202629*^9, {3.578904920383505*^9, 
   3.578904944547948*^9}, 3.578904999397644*^9, {3.578905078552183*^9, 
   3.578905107708634*^9}, 3.5789057725662017`*^9, 3.5789058096318665`*^9, {
   3.5789058410347223`*^9, 3.57890585704035*^9}, 3.578905890736409*^9, {
   3.578905924697669*^9, 3.578905948924512*^9}, {3.57890598232417*^9, 
   3.5789060082670155`*^9}, 3.5789060404342723`*^9, 3.578906078404739*^9, 
   3.5789269467901554`*^9, {3.5789269768826084`*^9, 3.578927029595101*^9}, {
   3.5789270807631903`*^9, 3.5789270976424203`*^9}, {3.578927132352481*^9, 
   3.5789271506513133`*^9}, {3.5789272032858057`*^9, 3.578927225203844*^9}, {
   3.57892726247231*^9, 3.578927301706379*^9}, 3.578927412981374*^9, 
   3.57892747333788*^9, {3.5789275896984844`*^9, 3.5789276649686165`*^9}, {
   3.578931858024057*^9, 3.5789318972581263`*^9}, {3.578931930470584*^9, 
   3.5789325495576715`*^9}}],

Cell[BoxData["666.1733149806745`"], "Output",
 CellChangeTimes->{
  3.578903430222088*^9, 3.578904578228104*^9, 3.5789046267441893`*^9, 
   3.5789046603778486`*^9, 3.5789046928103056`*^9, 3.578904739688388*^9, 
   3.5789047811064606`*^9, 3.578904877202629*^9, {3.578904920383505*^9, 
   3.578904944547948*^9}, 3.578904999397644*^9, {3.578905078552183*^9, 
   3.578905107708634*^9}, 3.5789057725662017`*^9, 3.5789058096318665`*^9, {
   3.5789058410347223`*^9, 3.57890585704035*^9}, 3.578905890736409*^9, {
   3.578905924697669*^9, 3.578905948924512*^9}, {3.57890598232417*^9, 
   3.5789060082670155`*^9}, 3.5789060404342723`*^9, 3.578906078404739*^9, 
   3.5789269467901554`*^9, {3.5789269768826084`*^9, 3.578927029595101*^9}, {
   3.5789270807631903`*^9, 3.5789270976424203`*^9}, {3.578927132352481*^9, 
   3.5789271506513133`*^9}, {3.5789272032858057`*^9, 3.578927225203844*^9}, {
   3.57892726247231*^9, 3.578927301706379*^9}, 3.578927412981374*^9, 
   3.57892747333788*^9, {3.5789275896984844`*^9, 3.5789276649686165`*^9}, {
   3.578931858024057*^9, 3.5789318972581263`*^9}, {3.578931930470584*^9, 
   3.5789325495576715`*^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.4309088743364695`, 0}, 0.2544847392835131], 
   DiskBox[{-1.4309088743364695`, 0}, 0.1], 
   CircleBox[{1.4309088743364695`, 0}, 0.2544847392835131], 
   DiskBox[{1.4309088743364695`, 0}, 0.1], 
   CircleBox[{0, 0}, 1.1764241350529565`], 
   CircleBox[{-1.9348408020892454`, -1.4253347334815887`}, 
    1.2573110439748147`], 
   DiskBox[{-2.115938761921436, -1.9375571153071527`}, 0.1], 
   CircleBox[{1.9348408020892454`, 1.4253347334815887`}, 1.2573110439748147`],
    DiskBox[{2.115938761921436, 1.9375571153071527`}, 0.1], 
   {Dashing[{0.03, 0.03}], 
    CircleBox[{1.9348408020892454`, -0.7126673667407943}, 
     1.2573110439748147`], 
    DiskBox[{2.115938761921436, -0.9687785576535763}, 0.1], 
    CircleBox[{-1.9348408020892454`, 0.7126673667407943}, 
     1.2573110439748147`], 
    DiskBox[{-2.115938761921436, 0.9687785576535763}, 0.1]}, 
   {Thickness[0.01], 
    LineBox[{{-1.4309088743364695`, 0}, {1.4309088743364695`, 0}}], 
    LineBox[{{1.4309088743364695`, 0}, {2.115938761921436, 
      1.9375571153071527`}}], 
    LineBox[{{-1.4309088743364695`, 
       0}, {-2.115938761921436, -1.9375571153071527`}}], 
    {Dashing[{0.03, 0.03}], 
     LineBox[{{1.4309088743364695`, 0}, {
       2.115938761921436, -0.9687785576535763}}], 
     LineBox[{{-1.4309088743364695`, 0}, {-2.115938761921436, 
       0.9687785576535763}}]}}},
  AspectRatio->Automatic,
  Axes->True,
  Frame->True,
  GridLines->Automatic,
  PlotRange->All]], "Output",
 CellChangeTimes->{
  3.578903430253288*^9, 3.578904578259304*^9, 3.5789046267753897`*^9, 
   3.578904660409048*^9, 3.5789046928259053`*^9, 3.5789047397039876`*^9, 
   3.5789047811376605`*^9, 3.5789048772338295`*^9, {3.5789049204147053`*^9, 
   3.5789049445791473`*^9}, 3.578904999428844*^9, {3.5789050785833826`*^9, 
   3.578905107739834*^9}, 3.578905772597402*^9, 3.5789058096474667`*^9, {
   3.578905841050322*^9, 3.5789058570715504`*^9}, 3.5789058907676096`*^9, {
   3.578905924713269*^9, 3.5789059489557114`*^9}, {3.5789059823553705`*^9, 
   3.578906008298216*^9}, 3.578906040465472*^9, 3.5789060784359393`*^9, 
   3.5789269468213553`*^9, {3.578926976913808*^9, 3.578927029626301*^9}, {
   3.5789270807943907`*^9, 3.57892709767362*^9}, {3.5789271323836813`*^9, 
   3.578927150682513*^9}, {3.5789272033326063`*^9, 3.5789272252350445`*^9}, {
   3.5789272625035095`*^9, 3.578927301737579*^9}, 3.578927413012574*^9, 
   3.5789274733690805`*^9, {3.578927589729685*^9, 3.578927664999817*^9}, {
   3.578931858055257*^9, 3.5789318972893257`*^9}, {3.5789319305017843`*^9, 
   3.5789325495732718`*^9}}]
}, Open  ]]
},
WindowSize->{869, 740},
WindowMargins->{{Automatic, 306}, {Automatic, 6}},
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, 1119, 22, 70, "Input"],
Cell[1679, 44, 2524, 50, 112, "Input"],
Cell[4206, 96, 2071, 53, 92, "Input"],
Cell[6280, 151, 2276, 54, 92, "Input"],
Cell[8559, 207, 3230, 80, 132, "Input"],
Cell[11792, 289, 1768, 42, 92, "Input"],
Cell[CellGroupData[{
Cell[13585, 335, 1631, 33, 92, "Input"],
Cell[15219, 370, 1463, 26, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[16719, 401, 3083, 83, 252, "Input"],
Cell[19805, 486, 1130, 16, 31, "Output"],
Cell[20938, 504, 1148, 17, 31, "Output"],
Cell[22089, 523, 1130, 16, 31, "Output"],
Cell[23222, 541, 1130, 16, 31, "Output"],
Cell[24355, 559, 1131, 16, 31, "Output"],
Cell[25489, 577, 1131, 16, 31, "Output"],
Cell[26623, 595, 1131, 16, 31, "Output"],
Cell[27757, 613, 1130, 16, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[28924, 634, 11962, 319, 412, "Input"],
Cell[40889, 955, 2596, 49, 322, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)