(* 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[     41325,        949]
NotebookOptionsPosition[     40057,        908]
NotebookOutlinePosition[     40542,        926]
CellTagsIndexPosition[     40499,        923]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 StyleBox[
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{"C2H6", " ", "Molecule", " ", "after", " ", "Kimball"}], ",", " ", 
    RowBox[{"non", " ", "overlapping"}], ",", 
    RowBox[{
     RowBox[{"parameterized", " ", "with", " ", "G3MP2"}], "//", 
     RowBox[{"6", "-", 
      RowBox[{"31", "g", 
       RowBox[{"(", "d", ")"}]}]}]}]}], " ", "*)"}], "Text"]], "Input",
 CellChangeTimes->{
  3.5788956762396965`*^9, {3.5788957137421627`*^9, 3.578895724880582*^9}, {
   3.578895806827526*^9, 3.578895904780098*^9}, {3.5788960148694916`*^9, 
   3.5788960197835*^9}, {3.5788961272676888`*^9, 3.578896146970524*^9}, {
   3.5788962074362297`*^9, 3.578896248542302*^9}, {3.5788963748741245`*^9, 
   3.578896407337781*^9}, 3.5789031386419754`*^9, {3.5789304190931296`*^9, 
   3.5789304212303333`*^9}, {3.5789317063449907`*^9, 3.5789317142698045`*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"Clear", "[", 
   RowBox[{"k1", ",", "k2", ",", "k3", ",", "s1", ",", "s2", ",", "s3"}], 
   "]"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{" ", 
  RowBox[{
   RowBox[{"n", "=", "2."}], ";", " ", 
   RowBox[{"cw", "=", 
    RowBox[{"1", "/", "3"}]}], ";", " ", 
   RowBox[{"sw", "=", 
    RowBox[{
     RowBox[{"Sqrt", "[", "2", "]"}], "/", "3"}]}], ";"}], " ", 
  StyleBox[
   RowBox[{"(*", 
    RowBox[{
     RowBox[{"cos", " ", "and", " ", "sin", " ", "of", " ", "180"}], "-", 
     "109.5"}], " ", "*)"}], "Text"], "\[IndentingNewLine]", 
  StyleBox[
   RowBox[{"(*", " ", 
    RowBox[{"Parameterlist", " ", "c"}], " ", "*)"}], 
   "Text"]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"c", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"k1", "->", "1.02027123"}], ",", 
     RowBox[{"k2", "->", "1.38902724"}], ",", 
     RowBox[{"k3", "->", "1.21474827"}], ",", 
     RowBox[{"s1", "->", "0.30256112"}], ",", 
     RowBox[{"s2", "->", "0.35724003"}], ",", 
     RowBox[{"s3", "->", "0.35319238"}]}], "}"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  StyleBox[
   RowBox[{"(*", " ", 
    RowBox[{"Terms", " ", "from", " ", 
     RowBox[{
      RowBox[{"ChemEdu", "/", "Kimball"}], "/", 
      RowBox[{"Alkan", ".", "sys"}]}]}], " ", "*)"}], "Text"]}]}], "Input",
 CellChangeTimes->{
  3.5788956762396965`*^9, {3.5788957137421627`*^9, 3.578895724880582*^9}, {
   3.578895806827526*^9, 3.578895904780098*^9}, {3.5788960148694916`*^9, 
   3.5788960197835*^9}, {3.5788961272676888`*^9, 3.578896146970524*^9}, {
   3.5788962074362297`*^9, 3.578896248542302*^9}, {3.5788963748741245`*^9, 
   3.578896407337781*^9}, {3.578897895970396*^9, 3.5788979299784555`*^9}, {
   3.578901683206219*^9, 3.5789017130022717`*^9}, 3.5789017480087333`*^9, 
   3.5789031559112062`*^9, {3.578930508668487*^9, 3.5789305512877617`*^9}, {
   3.578930588213027*^9, 3.5789305892270284`*^9}, {3.578930623203888*^9, 
   3.5789307850073724`*^9}, {3.5789308197642336`*^9, 
   3.5789308212306356`*^9}, {3.5789308916959596`*^9, 
   3.5789309168276043`*^9}, {3.578930968635295*^9, 3.578931053967445*^9}, {
   3.5789310902531085`*^9, 3.5789312758466344`*^9}, {3.57898059849006*^9, 
   3.57898062724091*^9}, {3.5789833698569264`*^9, 3.578983385831354*^9}, {
   3.578983821063319*^9, 3.5789838692362037`*^9}, {3.578984093892198*^9, 
   3.5789842193788185`*^9}, {3.578984259408489*^9, 3.5789842993913593`*^9}, 
   3.57898438517591*^9, {3.578984422584776*^9, 3.5789845278225603`*^9}, 
   3.578994102766264*^9, {3.578994136384323*^9, 3.5789944171536164`*^9}, {
   3.5789944872445393`*^9, 3.5789944880089407`*^9}, {3.578994526821809*^9, 
   3.5789946056019473`*^9}, {3.578994644555216*^9, 3.5789949234213057`*^9}, {
   3.5790013546913004`*^9, 3.5790013995725794`*^9}, {3.579001432987838*^9, 
   3.5790014462946615`*^9}, {3.5790031307327833`*^9, 3.579003231402541*^9}, {
   3.5790863790628405`*^9, 3.5790864565325766`*^9}, {3.579087975881403*^9, 
   3.5790880551589375`*^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[{"pr", "=", 
   RowBox[{"p", "*", "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}, {
   3.5789810992129383`*^9, 3.5789811049381485`*^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", "/", " ", "pr"}], " ", "-", " ", 
     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}, 3.578981172876268*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"Vnn", " ", "=", 
    RowBox[{
     RowBox[{"6.0", "*", 
      RowBox[{"z", "/", " ", 
       RowBox[{"(", 
        RowBox[{"r", "*", "p"}], ")"}]}]}], " ", "+", " ", 
     RowBox[{"z", "*", 
      RowBox[{"z", "/", " ", 
       RowBox[{"(", 
        RowBox[{"2", "*", "q"}], ")"}]}]}], "+", " ", 
     RowBox[{"6.0", "*", 
      RowBox[{"z", "/", " ", 
       RowBox[{"Sqrt", "[", 
        RowBox[{
         RowBox[{"vd", "*", "qrp"}], " ", "+", " ", "r2p2", " ", "+", " ", 
         RowBox[{"4.0", "*", "q2"}]}], "]"}]}]}], " ", "+", " ", 
     RowBox[{"6.0", "/", " ", 
      RowBox[{"Sqrt", "[", 
       RowBox[{"ad", "*", "r2p2"}], "]"}]}], "+", " ", 
     RowBox[{"3.0", "/", " ", 
      RowBox[{"Sqrt", "[", 
       RowBox[{
        RowBox[{"ad", "*", "qrp"}], " ", "+", " ", 
        RowBox[{"4.0", "*", "r2p2"}], " ", "+", " ", 
        RowBox[{"4.0", "*", "q2"}]}], "]"}]}], "+", " ", 
     RowBox[{"6.0", "/", " ", 
      RowBox[{"Sqrt", "[", 
       RowBox[{
        RowBox[{"ad", "*", "qrp"}], " ", "+", " ", 
        RowBox[{"vd", "*", "r2p2"}], " ", "+", " ", 
        RowBox[{"4.0", "*", "q2"}]}], "]"}]}]}]}], " ", ";"}], 
  "\[IndentingNewLine]"}]], "Input",
 CellChangeTimes->{
  3.5788956762396965`*^9, {3.5788957137421627`*^9, 3.578895724880582*^9}, {
   3.578895806827526*^9, 3.578895904780098*^9}, {3.5788960148694916`*^9, 
   3.5788960197835*^9}, {3.5788961272676888`*^9, 3.578896146970524*^9}, {
   3.5788962074362297`*^9, 3.578896248542302*^9}, {3.5788963748741245`*^9, 
   3.5788964212998056`*^9}, {3.57889769338844*^9, 3.5788976974288473`*^9}, {
   3.5788979843601513`*^9, 3.578897986466155*^9}, {3.578902015393203*^9, 
   3.578902039620045*^9}, {3.5789024960924473`*^9, 3.578902497106449*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"func", " ", "=", " ", 
   RowBox[{
   "T", " ", "+", " ", "Vne", " ", "+", " ", "Vee", " ", "+", " ", "Vnn"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"t", " ", "=", " ", 
  RowBox[{"FindMinimum", "[", 
   RowBox[{"func", ",", 
    RowBox[{"{", 
     RowBox[{"P", ",", "0.25", ",", "0.27"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"Q", ",", "1.1", ",", "1.2"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"R", ",", "1.2", ",", "1.3"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"X", ",", "1.3", ",", "1.4"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"Method", "\[Rule]", "Automatic"}], ",", 
      RowBox[{"MaxIterations", "\[Rule]", "500"}]}], "}"}]}], 
   "]"}]}], "\[IndentingNewLine]"}], "Input",
 CellChangeTimes->{
  3.5788956762396965`*^9, {3.5788957137421627`*^9, 3.578895724880582*^9}, {
   3.578895806827526*^9, 3.578895904780098*^9}, {3.5788960148694916`*^9, 
   3.5788960197835*^9}, {3.5788961272676888`*^9, 3.578896146970524*^9}, {
   3.5788962074362297`*^9, 3.578896248542302*^9}, {3.5788963748741245`*^9, 
   3.57889639532576*^9}, {3.5788978015434303`*^9, 3.5788978862671785`*^9}, {
   3.5788979169212327`*^9, 3.5788979187152357`*^9}, {3.5789023036037087`*^9, 
   3.578902308065317*^9}, {3.578902360762209*^9, 3.57890236119901*^9}, {
   3.578902644932308*^9, 3.578902662295139*^9}, {3.578903281959427*^9, 
   3.578903341645132*^9}, {3.5789033823612037`*^9, 3.5789034163380632`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", "79.64080110008578`"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"P", "\[Rule]", "0.2608822815444364`"}], ",", 
     RowBox[{"Q", "\[Rule]", "1.2187880874397026`"}], ",", 
     RowBox[{"R", "\[Rule]", "1.275633518434545`"}], ",", 
     RowBox[{"X", "\[Rule]", "1.3587096241370047`"}]}], "}"}]}], 
  "}"}]], "Output",
 CellChangeTimes->{
  3.578903430222088*^9, {3.5789304269399433`*^9, 3.578930454224391*^9}, {
   3.578930537996538*^9, 3.5789305564669704`*^9}, 3.5789305967930417`*^9, {
   3.578930629568699*^9, 3.578930760406129*^9}, 3.578930791247383*^9, 
   3.578930827907448*^9, 3.5789308964695683`*^9, 3.5789309357504373`*^9, {
   3.5789309749845057`*^9, 3.578931060457056*^9}, {3.5789310943091154`*^9, 
   3.578931280261442*^9}, {3.578980022393448*^9, 3.578980064653922*^9}, 
   3.5789804550322075`*^9, 3.5789807979920096`*^9, 3.578981298672488*^9, 
   3.5789836189493637`*^9, {3.578984132065465*^9, 3.5789842238092265`*^9}, {
   3.5789842631992955`*^9, 3.5789843051009693`*^9}, 3.5789843889979167`*^9, {
   3.5789844276859846`*^9, 3.578984537681778*^9}, 3.578994107758273*^9, {
   3.5789941432951355`*^9, 3.5789943264862576`*^9}, {3.578994377092746*^9, 
   3.578994420772823*^9}, 3.5789944932817497`*^9, {3.5789945313458166`*^9, 
   3.5789946092055535`*^9}, {3.5789946494692245`*^9, 
   3.5789946819484816`*^9}, {3.5789947144433384`*^9, 3.5789949267909117`*^9}, 
   3.579001460459486*^9, 3.5790032421371555`*^9, 3.579086474441408*^9, 
   3.5790880669606123`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"u", " ", "=", " ", 
   RowBox[{"t", "[", 
    RowBox[{"[", "2", "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"T", " ", "/.", "u"}], "\[IndentingNewLine]", 
 RowBox[{"Vne", " ", "/.", "u"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Vee", " ", "/.", "u"}], " ", "/.", "c"}], "\[IndentingNewLine]", 
 RowBox[{"Vnn", " ", "/.", "u"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"ch", " ", "=", " ", 
   RowBox[{
    RowBox[{"X", "*", 
     RowBox[{"(", 
      RowBox[{"P", " ", "+", " ", "R"}], ")"}]}], " ", "/.", "u"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"dch", " ", "=", " ", 
   RowBox[{"ch", "*", "0.529177"}]}], 
  "                                             ", 
  StyleBox[
   RowBox[{"(*", " ", 
    RowBox[{"C", "-", 
     RowBox[{"H", " ", "distance"}]}], " ", "*)"}], 
   "Text"]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"cc", " ", "=", " ", 
   RowBox[{
    RowBox[{"2.", "*", 
     RowBox[{"(", 
      RowBox[{"P", " ", "+", " ", "Q"}], ")"}]}], " ", "/.", "u"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"dcc", "=", 
   RowBox[{"cc", "*", "0.529177"}]}], 
  "                                               ", 
  StyleBox[
   RowBox[{"(*", " ", 
    RowBox[{"C", "-", 
     RowBox[{"C", " ", "distance"}]}], " ", "*)"}], 
   "Text"]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{
     RowBox[{"-", 
      RowBox[{"(", 
       RowBox[{"Vne", " ", "+", " ", "Vee", " ", "+", " ", "Vnn"}], ")"}]}], 
     "/", "T"}], " ", "/.", "u"}], " ", "/.", "c"}], "              ", 
  StyleBox[
   RowBox[{"(*", " ", 
    RowBox[{"Virial", " ", "Theorem"}], " ", "*)"}], 
   "Text"]}], "\[IndentingNewLine]", 
 RowBox[{
  StyleBox[
   RowBox[{"(*", 
    RowBox[{
    "atomization", " ", "energy", " ", "at", " ", "0", " ", "K", " ", "in", 
     " ", 
     RowBox[{"kcal", "/", "mol"}]}], " ", "*)"}], "Text"], 
  RowBox[{
   RowBox[{"(", 
    RowBox[{
     RowBox[{"-", 
      RowBox[{"t", "[", 
       RowBox[{"[", "1", "]"}], "]"}]}], " ", "-", " ", 
     RowBox[{"n", "*", "37.78899"}], " ", "-", 
     RowBox[{
      RowBox[{"(", " ", 
       RowBox[{
        RowBox[{"2", "*", "n"}], " ", "+", " ", "2"}], ")"}], "*", 
      "0.50177"}]}], ")"}], "*", "627.50956", " "}]}]}], "Input",
 CellChangeTimes->{
  3.5788956762396965`*^9, {3.5788957137421627`*^9, 3.578895724880582*^9}, {
   3.578895806827526*^9, 3.578895904780098*^9}, {3.5788960148694916`*^9, 
   3.5788960197835*^9}, {3.5788961272676888`*^9, 3.578896146970524*^9}, {
   3.5788962074362297`*^9, 3.578896248542302*^9}, {3.5788963748741245`*^9, 
   3.57889639532576*^9}, 3.578896439614238*^9, {3.578902116762181*^9, 
   3.5789021273389993`*^9}, 3.5789025446709323`*^9, 3.5789025890062103`*^9, {
   3.5789028686835012`*^9, 3.578902880882723*^9}, {3.5789029284316063`*^9, 
   3.578902967977676*^9}, {3.5789030373665977`*^9, 3.5789031226987476`*^9}, {
   3.578903474307765*^9, 3.5789034887065907`*^9}, {3.578930301359723*^9, 
   3.5789303077401342`*^9}, {3.578930371263446*^9, 3.5789303949442873`*^9}, 
   3.578930448764382*^9}],

Cell[BoxData["79.64067283404866`"], "Output",
 CellChangeTimes->{
  3.578903430222088*^9, {3.5789304269555435`*^9, 3.578930454239991*^9}, {
   3.5789305380121384`*^9, 3.5789305564825706`*^9}, 3.578930596808642*^9, {
   3.578930629584299*^9, 3.578930760421729*^9}, 3.5789307912629833`*^9, 
   3.5789308279230475`*^9, 3.578930896485168*^9, 3.578930935766037*^9, {
   3.578930975000106*^9, 3.5789310604726562`*^9}, {3.5789310943247156`*^9, 
   3.5789312802770424`*^9}, {3.578980022393448*^9, 3.5789800646695223`*^9}, 
   3.5789804550478077`*^9, 3.5789807980076094`*^9, 3.5789812986880884`*^9, 
   3.578983618964964*^9, {3.578984132081065*^9, 3.5789842238248262`*^9}, {
   3.5789842632148952`*^9, 3.5789843051165695`*^9}, 3.5789843890135164`*^9, {
   3.5789844277015843`*^9, 3.5789845376973777`*^9}, 3.578994107773873*^9, {
   3.5789941433107357`*^9, 3.5789943265018573`*^9}, {3.578994377108346*^9, 
   3.578994420772823*^9}, 3.57899449329735*^9, {3.578994531361417*^9, 
   3.5789946092211533`*^9}, {3.5789946494692245`*^9, 
   3.5789946819484816`*^9}, {3.5789947144433384`*^9, 3.578994926806512*^9}, 
   3.579001460475086*^9, 3.5790032421481557`*^9, 3.5790864744570084`*^9, 
   3.5790880669796133`*^9}],

Cell[BoxData[
 RowBox[{"-", "269.8770560419033`"}]], "Output",
 CellChangeTimes->{
  3.578903430222088*^9, {3.5789304269555435`*^9, 3.578930454239991*^9}, {
   3.5789305380121384`*^9, 3.5789305564825706`*^9}, 3.578930596808642*^9, {
   3.578930629584299*^9, 3.578930760421729*^9}, 3.5789307912629833`*^9, 
   3.5789308279230475`*^9, 3.578930896485168*^9, 3.578930935766037*^9, {
   3.578930975000106*^9, 3.5789310604726562`*^9}, {3.5789310943247156`*^9, 
   3.5789312802770424`*^9}, {3.578980022393448*^9, 3.5789800646695223`*^9}, 
   3.5789804550478077`*^9, 3.5789807980076094`*^9, 3.5789812986880884`*^9, 
   3.578983618964964*^9, {3.578984132081065*^9, 3.5789842238248262`*^9}, {
   3.5789842632148952`*^9, 3.5789843051165695`*^9}, 3.5789843890135164`*^9, {
   3.5789844277015843`*^9, 3.5789845376973777`*^9}, 3.578994107773873*^9, {
   3.5789941433107357`*^9, 3.5789943265018573`*^9}, {3.578994377108346*^9, 
   3.578994420772823*^9}, 3.57899449329735*^9, {3.578994531361417*^9, 
   3.5789946092211533`*^9}, {3.5789946494692245`*^9, 
   3.5789946819484816`*^9}, {3.5789947144433384`*^9, 3.578994926806512*^9}, 
   3.579001460475086*^9, 3.5790032421481557`*^9, 3.5790864744570084`*^9, 
   3.5790880669806137`*^9}],

Cell[BoxData["68.98746575274905`"], "Output",
 CellChangeTimes->{
  3.578903430222088*^9, {3.5789304269555435`*^9, 3.578930454239991*^9}, {
   3.5789305380121384`*^9, 3.5789305564825706`*^9}, 3.578930596808642*^9, {
   3.578930629584299*^9, 3.578930760421729*^9}, 3.5789307912629833`*^9, 
   3.5789308279230475`*^9, 3.578930896485168*^9, 3.578930935766037*^9, {
   3.578930975000106*^9, 3.5789310604726562`*^9}, {3.5789310943247156`*^9, 
   3.5789312802770424`*^9}, {3.578980022393448*^9, 3.5789800646695223`*^9}, 
   3.5789804550478077`*^9, 3.5789807980076094`*^9, 3.5789812986880884`*^9, 
   3.578983618964964*^9, {3.578984132081065*^9, 3.5789842238248262`*^9}, {
   3.5789842632148952`*^9, 3.5789843051165695`*^9}, 3.5789843890135164`*^9, {
   3.5789844277015843`*^9, 3.5789845376973777`*^9}, 3.578994107773873*^9, {
   3.5789941433107357`*^9, 3.5789943265018573`*^9}, {3.578994377108346*^9, 
   3.578994420772823*^9}, 3.57899449329735*^9, {3.578994531361417*^9, 
   3.5789946092211533`*^9}, {3.5789946494692245`*^9, 
   3.5789946819484816`*^9}, {3.5789947144433384`*^9, 3.578994926806512*^9}, 
   3.579001460475086*^9, 3.5790032421481557`*^9, 3.5790864744570084`*^9, 
   3.579088066981614*^9}],

Cell[BoxData["41.60811635501984`"], "Output",
 CellChangeTimes->{
  3.578903430222088*^9, {3.5789304269555435`*^9, 3.578930454239991*^9}, {
   3.5789305380121384`*^9, 3.5789305564825706`*^9}, 3.578930596808642*^9, {
   3.578930629584299*^9, 3.578930760421729*^9}, 3.5789307912629833`*^9, 
   3.5789308279230475`*^9, 3.578930896485168*^9, 3.578930935766037*^9, {
   3.578930975000106*^9, 3.5789310604726562`*^9}, {3.5789310943247156`*^9, 
   3.5789312802770424`*^9}, {3.578980022393448*^9, 3.5789800646695223`*^9}, 
   3.5789804550478077`*^9, 3.5789807980076094`*^9, 3.5789812986880884`*^9, 
   3.578983618964964*^9, {3.578984132081065*^9, 3.5789842238248262`*^9}, {
   3.5789842632148952`*^9, 3.5789843051165695`*^9}, 3.5789843890135164`*^9, {
   3.5789844277015843`*^9, 3.5789845376973777`*^9}, 3.578994107773873*^9, {
   3.5789941433107357`*^9, 3.5789943265018573`*^9}, {3.578994377108346*^9, 
   3.578994420772823*^9}, 3.57899449329735*^9, {3.578994531361417*^9, 
   3.5789946092211533`*^9}, {3.5789946494692245`*^9, 
   3.5789946819484816`*^9}, {3.5789947144433384`*^9, 3.578994926806512*^9}, 
   3.579001460475086*^9, 3.5790032421481557`*^9, 3.5790864744570084`*^9, 
   3.579088066982614*^9}],

Cell[BoxData["1.1047516070305332`"], "Output",
 CellChangeTimes->{
  3.578903430222088*^9, {3.5789304269555435`*^9, 3.578930454239991*^9}, {
   3.5789305380121384`*^9, 3.5789305564825706`*^9}, 3.578930596808642*^9, {
   3.578930629584299*^9, 3.578930760421729*^9}, 3.5789307912629833`*^9, 
   3.5789308279230475`*^9, 3.578930896485168*^9, 3.578930935766037*^9, {
   3.578930975000106*^9, 3.5789310604726562`*^9}, {3.5789310943247156`*^9, 
   3.5789312802770424`*^9}, {3.578980022393448*^9, 3.5789800646695223`*^9}, 
   3.5789804550478077`*^9, 3.5789807980076094`*^9, 3.5789812986880884`*^9, 
   3.578983618964964*^9, {3.578984132081065*^9, 3.5789842238248262`*^9}, {
   3.5789842632148952`*^9, 3.5789843051165695`*^9}, 3.5789843890135164`*^9, {
   3.5789844277015843`*^9, 3.5789845376973777`*^9}, 3.578994107773873*^9, {
   3.5789941433107357`*^9, 3.5789943265018573`*^9}, {3.578994377108346*^9, 
   3.578994420772823*^9}, 3.57899449329735*^9, {3.578994531361417*^9, 
   3.5789946092211533`*^9}, {3.5789946494692245`*^9, 
   3.5789946819484816`*^9}, {3.5789947144433384`*^9, 3.578994926806512*^9}, 
   3.579001460475086*^9, 3.5790032421481557`*^9, 3.5790864744570084`*^9, 
   3.579088066984614*^9}],

Cell[BoxData["1.5660150536958393`"], "Output",
 CellChangeTimes->{
  3.578903430222088*^9, {3.5789304269555435`*^9, 3.578930454239991*^9}, {
   3.5789305380121384`*^9, 3.5789305564825706`*^9}, 3.578930596808642*^9, {
   3.578930629584299*^9, 3.578930760421729*^9}, 3.5789307912629833`*^9, 
   3.5789308279230475`*^9, 3.578930896485168*^9, 3.578930935766037*^9, {
   3.578930975000106*^9, 3.5789310604726562`*^9}, {3.5789310943247156`*^9, 
   3.5789312802770424`*^9}, {3.578980022393448*^9, 3.5789800646695223`*^9}, 
   3.5789804550478077`*^9, 3.5789807980076094`*^9, 3.5789812986880884`*^9, 
   3.578983618964964*^9, {3.578984132081065*^9, 3.5789842238248262`*^9}, {
   3.5789842632148952`*^9, 3.5789843051165695`*^9}, 3.5789843890135164`*^9, {
   3.5789844277015843`*^9, 3.5789845376973777`*^9}, 3.578994107773873*^9, {
   3.5789941433107357`*^9, 3.5789943265018573`*^9}, {3.578994377108346*^9, 
   3.578994420772823*^9}, 3.57899449329735*^9, {3.578994531361417*^9, 
   3.5789946092211533`*^9}, {3.5789946494692245`*^9, 
   3.5789946819484816`*^9}, {3.5789947144433384`*^9, 3.578994926806512*^9}, 
   3.579001460475086*^9, 3.5790032421481557`*^9, 3.5790864744570084`*^9, 
   3.579088066985614*^9}],

Cell[BoxData["2.0000016105594365`"], "Output",
 CellChangeTimes->{
  3.578903430222088*^9, {3.5789304269555435`*^9, 3.578930454239991*^9}, {
   3.5789305380121384`*^9, 3.5789305564825706`*^9}, 3.578930596808642*^9, {
   3.578930629584299*^9, 3.578930760421729*^9}, 3.5789307912629833`*^9, 
   3.5789308279230475`*^9, 3.578930896485168*^9, 3.578930935766037*^9, {
   3.578930975000106*^9, 3.5789310604726562`*^9}, {3.5789310943247156`*^9, 
   3.5789312802770424`*^9}, {3.578980022393448*^9, 3.5789800646695223`*^9}, 
   3.5789804550478077`*^9, 3.5789807980076094`*^9, 3.5789812986880884`*^9, 
   3.578983618964964*^9, {3.578984132081065*^9, 3.5789842238248262`*^9}, {
   3.5789842632148952`*^9, 3.5789843051165695`*^9}, 3.5789843890135164`*^9, {
   3.5789844277015843`*^9, 3.5789845376973777`*^9}, 3.578994107773873*^9, {
   3.5789941433107357`*^9, 3.5789943265018573`*^9}, {3.578994377108346*^9, 
   3.578994420772823*^9}, 3.57899449329735*^9, {3.578994531361417*^9, 
   3.5789946092211533`*^9}, {3.5789946494692245`*^9, 
   3.5789946819484816`*^9}, {3.5789947144433384`*^9, 3.578994926806512*^9}, 
   3.579001460475086*^9, 3.5790032421481557`*^9, 3.5790864744570084`*^9, 
   3.579088066986614*^9}],

Cell[BoxData["660.2662493463464`"], "Output",
 CellChangeTimes->{
  3.578903430222088*^9, {3.5789304269555435`*^9, 3.578930454239991*^9}, {
   3.5789305380121384`*^9, 3.5789305564825706`*^9}, 3.578930596808642*^9, {
   3.578930629584299*^9, 3.578930760421729*^9}, 3.5789307912629833`*^9, 
   3.5789308279230475`*^9, 3.578930896485168*^9, 3.578930935766037*^9, {
   3.578930975000106*^9, 3.5789310604726562`*^9}, {3.5789310943247156`*^9, 
   3.5789312802770424`*^9}, {3.578980022393448*^9, 3.5789800646695223`*^9}, 
   3.5789804550478077`*^9, 3.5789807980076094`*^9, 3.5789812986880884`*^9, 
   3.578983618964964*^9, {3.578984132081065*^9, 3.5789842238248262`*^9}, {
   3.5789842632148952`*^9, 3.5789843051165695`*^9}, 3.5789843890135164`*^9, {
   3.5789844277015843`*^9, 3.5789845376973777`*^9}, 3.578994107773873*^9, {
   3.5789941433107357`*^9, 3.5789943265018573`*^9}, {3.578994377108346*^9, 
   3.578994420772823*^9}, 3.57899449329735*^9, {3.578994531361417*^9, 
   3.5789946092211533`*^9}, {3.5789946494692245`*^9, 
   3.5789946819484816`*^9}, {3.5789947144433384`*^9, 3.578994926806512*^9}, 
   3.579001460475086*^9, 3.5790032421481557`*^9, 3.5790864744570084`*^9, 
   3.579088066987614*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"plot1", " ", "=", " ", 
    RowBox[{
     RowBox[{"Graphics", "[", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"-", "q"}], ",", "0"}], "}"}], ",", "P"}], "]"}], ",", 
        RowBox[{"Disk", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"-", "q"}], ",", "0"}], "}"}], ",", "0.1"}], "]"}], ",", 
        "\[IndentingNewLine]", 
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{"q", ",", "0"}], "}"}], ",", "P"}], "]"}], ",", 
        RowBox[{"Disk", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{"q", ",", "0"}], "}"}], ",", "0.1"}], "]"}], ",", 
        "\[IndentingNewLine]", 
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{"0", ",", "0"}], "}"}], ",", "Q"}], "]"}], ",", 
        "\[IndentingNewLine]", 
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{
             RowBox[{"-", "q"}], "-", 
             RowBox[{"r", "*", "cw"}]}], ",", 
            RowBox[{
             RowBox[{"-", "r"}], "*", "2", "*", "sw"}]}], "}"}], ",", "R"}], 
         "]"}], ",", 
        RowBox[{"Disk", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{
             RowBox[{"-", "q"}], "-", " ", 
             RowBox[{"pr", "*", "cw"}]}], ",", 
            RowBox[{
             RowBox[{"-", "pr"}], "*", "2", "*", "sw"}]}], "}"}], ",", 
          "0.1"}], "]"}], ",", 
        RowBox[{"Circle", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"q", " ", "+", 
             RowBox[{"r", "*", "cw"}]}], ",", 
            RowBox[{"r", "*", "2", "*", "sw"}]}], "}"}], ",", "R"}], "]"}], 
        ",", 
        RowBox[{"Disk", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"q", "+", " ", 
             RowBox[{"pr", "*", "cw"}]}], ",", 
            RowBox[{"pr", "*", "2", "*", "sw"}]}], "}"}], ",", "0.1"}], "]"}],
         ",", "\[IndentingNewLine]", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"Dashing", "[", 
           RowBox[{"{", 
            RowBox[{"0.02", ",", "0.02"}], "}"}], "]"}], ",", 
          "\[IndentingNewLine]", 
          RowBox[{"Circle", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{"q", "+", 
               RowBox[{"r", "*", "cw"}]}], ",", 
              RowBox[{
               RowBox[{"-", "r"}], "*", "sw"}]}], "}"}], ",", "R"}], "]"}], 
          ",", 
          RowBox[{"Disk", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{"q", " ", "+", " ", 
               RowBox[{"pr", "*", "cw"}]}], ",", 
              RowBox[{
               RowBox[{"-", "pr"}], "*", "sw"}]}], "}"}], ",", "0.1"}], "]"}],
           ",", 
          RowBox[{"Circle", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{
               RowBox[{"-", "q"}], "-", 
               RowBox[{"r", "*", "cw"}]}], ",", 
              RowBox[{"r", "*", "sw"}]}], "}"}], ",", "R"}], "]"}], ",", 
          RowBox[{"Disk", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{
               RowBox[{"-", "q"}], "-", 
               RowBox[{"pr", "*", "cw"}]}], ",", 
              RowBox[{"pr", "*", "sw"}]}], "}"}], ",", "0.1"}], "]"}]}], 
         "}"}], ",", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"Thickness", "[", "0.01", "]"}], ",", 
          RowBox[{"Line", "[", 
           RowBox[{"{", 
            RowBox[{
             RowBox[{"{", 
              RowBox[{
               RowBox[{"-", "q"}], ",", "0"}], "}"}], ",", 
             RowBox[{"{", 
              RowBox[{"q", ",", "0"}], "}"}]}], "}"}], "]"}], ",", 
          RowBox[{"Line", "[", 
           RowBox[{"{", 
            RowBox[{
             RowBox[{"{", 
              RowBox[{"q", ",", "0"}], "}"}], ",", 
             RowBox[{"{", 
              RowBox[{
               RowBox[{"q", " ", "+", " ", 
                RowBox[{"pr", "*", "cw"}]}], ",", 
               RowBox[{"pr", "*", "2", "*", "sw"}]}], "}"}]}], "}"}], "]"}], 
          ",", 
          RowBox[{"Line", "[", 
           RowBox[{"{", 
            RowBox[{
             RowBox[{"{", 
              RowBox[{
               RowBox[{"-", "q"}], ",", "0"}], "}"}], ",", 
             RowBox[{"{", 
              RowBox[{
               RowBox[{
                RowBox[{"-", "q"}], "-", " ", 
                RowBox[{"pr", "/", "3."}]}], ",", 
               RowBox[{
                RowBox[{"-", "pr"}], "*", "2", "*", "sw"}]}], "}"}]}], "}"}], 
           "]"}], ",", "\[IndentingNewLine]", 
          RowBox[{"{", 
           RowBox[{
            RowBox[{"Dashing", "[", 
             RowBox[{"{", 
              RowBox[{"0.01", ",", "0.01"}], "}"}], "]"}], ",", 
            "\[IndentingNewLine]", 
            RowBox[{"Line", "[", 
             RowBox[{"{", 
              RowBox[{
               RowBox[{"{", 
                RowBox[{"q", ",", "0"}], "}"}], ",", 
               RowBox[{"{", 
                RowBox[{
                 RowBox[{"q", " ", "+", " ", 
                  RowBox[{"pr", "*", "cw"}]}], ",", 
                 RowBox[{
                  RowBox[{"-", "pr"}], "*", "sw"}]}], "}"}]}], "}"}], "]"}], 
            ",", "\[IndentingNewLine]", 
            RowBox[{"Line", "[", 
             RowBox[{"{", 
              RowBox[{
               RowBox[{"{", 
                RowBox[{
                 RowBox[{"-", "q"}], ",", "0"}], "}"}], ",", 
               RowBox[{"{", 
                RowBox[{
                 RowBox[{
                  RowBox[{"-", "q"}], " ", "-", " ", 
                  RowBox[{"pr", "*", "cw"}]}], ",", 
                 RowBox[{"pr", "*", "sw"}]}], "}"}]}], "}"}], "]"}]}], 
           "}"}]}], "}"}]}], "}"}], "]"}], " ", "/.", "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, 3.578979713185305*^9, {3.5789799418037066`*^9, 
   3.578979956342932*^9}, {3.5789800116606293`*^9, 3.578980061502717*^9}, {
   3.578980173900914*^9, 3.5789802087045755`*^9}, {3.5789802638886724`*^9, 
   3.578980291485121*^9}, {3.5789803441440134`*^9, 3.5789803646424494`*^9}, {
   3.5789806842590103`*^9, 3.5789806925426245`*^9}, {3.5789807304262915`*^9, 
   3.5789807334994965`*^9}, {3.5789807825235825`*^9, 
   3.5789807838651853`*^9}, {3.578981252040807*^9, 3.5789812780148525`*^9}, {
   3.5789834407746506`*^9, 3.57898345180387*^9}, {3.57898349161514*^9, 
   3.578983600650532*^9}}],

Cell[BoxData[
 GraphicsBox[{CircleBox[{-1.4796703689841388`, 0}, 0.2608822815444364], 
   DiskBox[{-1.4796703689841388`, 0}, 0.1], 
   CircleBox[{1.4796703689841388`, 0}, 0.2608822815444364], 
   DiskBox[{1.4796703689841388`, 0}, 0.1], 
   CircleBox[{0, 0}, 1.2187880874397026`], 
   CircleBox[{-1.9918423023104659`, -1.4486409887538807`}, 1.275633518434545],
    DiskBox[{-2.1755633040074756`, -1.968282453339244}, 0.1], 
   CircleBox[{1.9918423023104659`, 1.4486409887538807`}, 1.275633518434545], 
   DiskBox[{2.1755633040074756`, 1.968282453339244}, 0.1], 
   {Dashing[{0.02, 0.02}], 
    CircleBox[{1.9918423023104659`, -0.7243204943769403}, 1.275633518434545], 
    DiskBox[{2.1755633040074756`, -0.984141226669622}, 0.1], 
    CircleBox[{-1.9918423023104659`, 0.7243204943769403}, 1.275633518434545], 
    DiskBox[{-2.1755633040074756`, 0.984141226669622}, 0.1]}, 
   {Thickness[0.01], 
    LineBox[{{-1.4796703689841388`, 0}, {1.4796703689841388`, 0}}], 
    LineBox[{{1.4796703689841388`, 0}, {2.1755633040074756`, 
      1.968282453339244}}], 
    LineBox[{{-1.4796703689841388`, 
       0}, {-2.1755633040074756`, -1.968282453339244}}], 
    {Dashing[{0.01, 0.01}], 
     LineBox[{{1.4796703689841388`, 0}, {
       2.1755633040074756`, -0.984141226669622}}], 
     LineBox[{{-1.4796703689841388`, 0}, {-2.1755633040074756`, 
       0.984141226669622}}]}}},
  AspectRatio->Automatic,
  Axes->True,
  Frame->True,
  GridLines->Automatic,
  PlotRange->All]], "Output",
 CellChangeTimes->{
  3.578903430253288*^9, {3.578930426986744*^9, 3.5789304542711916`*^9}, {
   3.578930538043339*^9, 3.578930556513771*^9}, 3.578930596839842*^9, {
   3.5789306296154995`*^9, 3.5789307604529295`*^9}, 3.5789307912941837`*^9, 
   3.5789308279542475`*^9, 3.578930896516368*^9, 3.5789309357972374`*^9, {
   3.5789309750313063`*^9, 3.578931060503856*^9}, {3.5789310943559155`*^9, 
   3.5789312803082423`*^9}, {3.5789800224402485`*^9, 3.5789800646851225`*^9}, 
   3.578980455063408*^9, 3.57898079803881*^9, 3.578981298719289*^9, 
   3.578983618980564*^9, {3.5789841321122656`*^9, 3.5789842238404264`*^9}, {
   3.5789842632304955`*^9, 3.5789843051321692`*^9}, 3.578984389044717*^9, {
   3.5789844277171845`*^9, 3.5789845377285776`*^9}, 3.578994107805073*^9, {
   3.578994143341935*^9, 3.578994326533057*^9}, {3.578994377139546*^9, 
   3.578994420804023*^9}, 3.5789944933285503`*^9, {3.578994531392617*^9, 
   3.5789946092523537`*^9}, {3.5789946495004244`*^9, 
   3.5789946819796815`*^9}, {3.578994714474539*^9, 3.5789949268377113`*^9}, 
   3.579001460521886*^9, 3.5790032421771574`*^9, 3.5790864744882083`*^9, 
   3.579088067007615*^9}]
}, Open  ]],

Cell[BoxData["\[IndentingNewLine]"], "Input",
 CellChangeTimes->{3.5789832731211567`*^9}]
},
WindowSize->{872, 739},
WindowMargins->{{Automatic, 249}, {Automatic, 31}},
PrintingCopies->1,
PrintingPageRange->{Automatic, Automatic},
PrivateNotebookOptions->{"VersionedStylesheet"->{"Default.nb"[8.] -> False}},
FrontEndVersion->"9.0 for Microsoft Windows (64-bit) (January 25, 2013)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[557, 20, 861, 17, 31, "Input"],
Cell[1421, 39, 2979, 63, 132, "Input"],
Cell[4403, 104, 2188, 56, 92, "Input"],
Cell[6594, 162, 2276, 54, 92, "Input"],
Cell[8873, 218, 3198, 78, 132, "Input"],
Cell[12074, 298, 1768, 42, 92, "Input"],
Cell[CellGroupData[{
Cell[13867, 344, 1464, 31, 92, "Input"],
Cell[15334, 377, 1533, 27, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[16904, 409, 3106, 84, 252, "Input"],
Cell[20013, 495, 1199, 17, 31, "Output"],
Cell[21215, 514, 1216, 18, 31, "Output"],
Cell[22434, 534, 1197, 17, 31, "Output"],
Cell[23634, 553, 1197, 17, 31, "Output"],
Cell[24834, 572, 1198, 17, 31, "Output"],
Cell[26035, 591, 1198, 17, 31, "Output"],
Cell[27236, 610, 1198, 17, 31, "Output"],
Cell[28437, 629, 1197, 17, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[29671, 651, 7650, 202, 332, "Input"],
Cell[37324, 855, 2625, 47, 320, "Output"]
}, Open  ]],
Cell[39964, 905, 89, 1, 52, "Input"]
}
]
*)

(* End of internal cache information *)