(* 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[     42598,       1032]
NotebookOptionsPosition[     40608,        962]
NotebookOutlinePosition[     41093,        980]
CellTagsIndexPosition[     41050,        977]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["\<\
C2H6 Molecule after Kimball, non overlapping, parameterized with MP2/6-31G(d)\
\>", "Subsubtitle",
 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}, {3.5792337731768165`*^9, 3.579233777451224*^9}, 
   3.5792341967955604`*^9, 3.579234656987569*^9, {3.579237662522448*^9, 
   3.579237690805298*^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"}]}], ";", " ", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{"cos", " ", "and", " ", "sin", " ", "of", " ", "180"}], "-", 
     "109.5"}], " ", "*)"}], "\[IndentingNewLine]", " ", 
   RowBox[{"(*", 
    RowBox[{"Parameterlist", " ", "c"}], " ", "*)"}], "\[IndentingNewLine]", 
   RowBox[{"c", "=", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"k1", "->", "1.02669238"}], ",", 
      RowBox[{"k2", "->", "1.36467050"}], ",", 
      RowBox[{"k3", "->", "1.20088548"}], ",", 
      RowBox[{"s1", "->", "0.30102681"}], ",", 
      RowBox[{"s2", "->", "0.34863843"}], ",", 
      RowBox[{"s3", "->", "0.34858036"}]}], "}"}]}], ";"}]}]}], "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}, {3.5790943224634066`*^9, 3.579094402850005*^9}, {
   3.5792337813044305`*^9, 3.5792337913664484`*^9}, {3.5792339896271963`*^9, 
   3.5792339924820013`*^9}, 3.5792355216502876`*^9}],

Cell[BoxData[
 RowBox[{"(*", 
  RowBox[{"Terms", " ", "from", " ", 
   RowBox[{"ChemEdu", "/", 
    RowBox[{"Alkan", ".", "sys"}]}]}], "*)"}]], "Text", "Input",
 CellChangeTimes->{{3.5792338519413548`*^9, 3.5792338963078327`*^9}, {
   3.579233998924813*^9, 3.579234002029218*^9}, 3.579235530823104*^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.49412996948107`"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"P", "\[Rule]", "0.26224005195523936`"}], ",", 
     RowBox[{"Q", "\[Rule]", "1.2000710223961832`"}], ",", 
     RowBox[{"R", "\[Rule]", "1.261969408118396`"}], ",", 
     RowBox[{"X", "\[Rule]", "1.356989234003253`"}]}], "}"}]}], 
  "}"}]], "Output",
 CellChangeTimes->{3.5792377218337517`*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"u", " ", "=", " ", 
   RowBox[{"t", "[", 
    RowBox[{"[", "2", "]"}], "]"}]}], ";"}]], "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, 3.579101184396015*^9, {3.5792316082810135`*^9, 
   3.57923163459826*^9}, {3.5792316901987576`*^9, 3.579231853687045*^9}, 
   3.579231923980768*^9, {3.579232395803597*^9, 3.5792324440544815`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"Total", " ", "Energy"}], " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{"Etot", "=", 
   RowBox[{"t", "[", 
    RowBox[{"[", "1", "]"}], "]"}]}]}]], "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, 3.579101184396015*^9, {3.5792316082810135`*^9, 
   3.57923163459826*^9}, {3.5792316901987576`*^9, 3.579231853687045*^9}, 
   3.579231923980768*^9, {3.579232395803597*^9, 3.5792324440544815`*^9}, {
   3.579233677423848*^9, 3.5792336815734553`*^9}, 3.579235541368722*^9}],

Cell[BoxData[
 RowBox[{"-", "79.49412996948107`"}]], "Output",
 CellChangeTimes->{3.579237721849352*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"Kinetic", " ", "energy"}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{"T", " ", "/.", "u"}]}]], "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, 3.579101184396015*^9, {3.5792316082810135`*^9, 
   3.57923163459826*^9}, {3.5792316901987576`*^9, 3.579231853687045*^9}, 
   3.579231923980768*^9, 3.5792319734952555`*^9, {3.579233685333062*^9, 
   3.579233688671468*^9}, 3.5792355483575344`*^9}],

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

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{" ", 
  RowBox[{"(*", 
   RowBox[{"NUCLEUS", "-", 
    RowBox[{"ELECTRON", " ", "POTENTIAL", " ", "ENERGY"}]}], 
   StyleBox[
    RowBox[{" ", " "}]], "*)"}], "\[IndentingNewLine]", 
  RowBox[{"Vne", " ", "/.", "u"}]}]], "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, 3.579101184396015*^9, {3.5792316082810135`*^9, 
   3.57923163459826*^9}, {3.5792316901987576`*^9, 3.579231853687045*^9}, {
   3.579231923980768*^9, 3.579231930080379*^9}, 3.5792319760224595`*^9, {
   3.5792336918850737`*^9, 3.5792336945838785`*^9}, 3.5792355580451517`*^9}],

Cell[BoxData[
 RowBox[{"-", "270.4620043570094`"}]], "Output",
 CellChangeTimes->{3.579237721864952*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"ELECTRON", "-", 
    RowBox[{"ELECTRON", " ", "POTENTIAL", " ", "ENERGY"}]}], " ", "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"Vee", " ", "/.", "u"}], " ", "/.", "c"}]}]], "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, 3.579101184396015*^9, {3.5792316082810135`*^9, 
   3.57923163459826*^9}, {3.5792316901987576`*^9, 3.579231853687045*^9}, {
   3.579231923980768*^9, 3.5792319350255876`*^9}, 3.5792319781128635`*^9, {
   3.579233697766284*^9, 3.5792337003714886`*^9}, 3.579235564737563*^9}],

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

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"NUCLEUS", "-", 
    RowBox[{"NUCLEUS", " ", "POTENTIAL", " ", "ENERGY"}]}], " ", "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{"Vnn", " ", "/.", "u"}]}]], "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, 3.579101184396015*^9, {3.5792316082810135`*^9, 
   3.57923163459826*^9}, {3.5792316901987576`*^9, 3.579231853687045*^9}, {
   3.579231923980768*^9, 3.579231939689996*^9}, 3.5792319799224663`*^9, {
   3.579233703163893*^9, 3.5792337068142996`*^9}, 3.5792355762971835`*^9}],

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

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{" ", 
  RowBox[{"(*", 
   RowBox[{"C", "-", 
    RowBox[{"H", " ", "distance", "  ", "\[CapitalARing]"}]}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"ch", " ", "=", " ", 
     RowBox[{
      RowBox[{"X", "*", 
       RowBox[{"(", 
        RowBox[{"P", " ", "+", " ", "R"}], ")"}]}], " ", "/.", "u"}]}], ";"}],
    "\[IndentingNewLine]", 
   RowBox[{"dch", " ", "=", " ", 
    RowBox[{
    "ch", "*", "0.529177", 
     "                                             "}]}]}]}]], "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, 3.579101184396015*^9, {3.5792316082810135`*^9, 
   3.57923163459826*^9}, {3.5792316901987576`*^9, 3.579231853687045*^9}, {
   3.579231923980768*^9, 3.5792319444324045`*^9}, 3.57923198196607*^9, {
   3.579232330751483*^9, 3.579232335509491*^9}, {3.5792337100279055`*^9, 
   3.579233713132311*^9}, 3.579235583473196*^9}],

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

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{" ", 
  RowBox[{"(*", 
   RowBox[{"C", "-", 
    RowBox[{"C", " ", "distance", "  ", "\[CapitalARing]"}]}], " ", "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"cc", " ", "=", " ", 
     RowBox[{
      RowBox[{"2.", "*", 
       RowBox[{"(", 
        RowBox[{"P", " ", "+", " ", "Q"}], ")"}]}], " ", "/.", "u"}]}], ";"}],
    "\[IndentingNewLine]", 
   RowBox[{"dcc", "=", 
    RowBox[{
    "cc", "*", "0.529177", 
     "                                              "}]}]}]}]], "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, 3.579101184396015*^9, {3.5792316082810135`*^9, 
   3.57923163459826*^9}, {3.5792316901987576`*^9, 3.579231853687045*^9}, {
   3.579231923980768*^9, 3.5792319537300205`*^9}, 3.579231983994074*^9, 
   3.5792322853398027`*^9, 3.579232323949871*^9, {3.579233721228725*^9, 
   3.57923372403673*^9}, 3.579235589853607*^9}],

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

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"Virial", " ", "Theorem"}], " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{
     RowBox[{"-", 
      RowBox[{"(", 
       RowBox[{"Vne", " ", "+", " ", "Vee", " ", "+", " ", "Vnn"}], ")"}]}], 
     "/", "T"}], " ", "/.", "u"}], " ", "/.", "c", 
   "             "}]}]], "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, 3.579101184396015*^9, {3.5792316082810135`*^9, 
   3.57923163459826*^9}, {3.5792316901987576`*^9, 3.579231853687045*^9}, {
   3.579231923980768*^9, 3.579231986490078*^9}, {3.579233726953935*^9, 
   3.5792337305575414`*^9}, 3.579235598917223*^9}],

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

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
   "Atomization", " ", "Energy", " ", "at", " ", "0", " ", "K", " ", "in", 
    " ", 
    RowBox[{"kcal", "/", "mol"}]}], " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{
     RowBox[{"-", 
      RowBox[{"t", "[", 
       RowBox[{"[", "1", "]"}], "]"}]}], " ", "-", " ", 
     RowBox[{"n", "*", "37.78899"}], " ", "-", 
     RowBox[{
      RowBox[{"(", " ", 
       RowBox[{
        RowBox[{"2", "*", "n"}], " ", "+", " ", "2"}], ")"}], "*", 
      "0.50177"}]}], ")"}], "*", "627.5096", " "}]}]], "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, 3.579101184396015*^9, {3.5792316082810135`*^9, 
   3.57923163459826*^9}, {3.5792316901987576`*^9, 3.579231853687045*^9}, {
   3.579231923980768*^9, 3.5792319607344327`*^9}, {3.5792337327883453`*^9, 
   3.579233738778756*^9}, {3.579234893452784*^9, 3.579234901361998*^9}, 
   3.5792356093536415`*^9}],

Cell[BoxData["568.2287489370774`"], "Output",
 CellChangeTimes->{3.579237721942952*^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[{
  RowBox[{"Show", "[", 
   RowBox[{"plot1", ",", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"AspectRatio", " ", "->", " ", "Automatic"}], ",", 
      RowBox[{"PlotRange", " ", "->", " ", "All"}], ",", 
      RowBox[{"Axes", " ", "->", " ", "True"}], ",", 
      RowBox[{"GridLines", " ", "->", " ", "Automatic"}], ",", 
      RowBox[{"Frame", " ", "->", " ", "True"}]}], "}"}]}], "]"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"Projection", " ", "on", " ", "xy", " ", "plane"}], " ", 
   "*)"}]}]}], "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}, {3.5792313709889965`*^9, 3.57923138446742*^9}, {
   3.579233747530371*^9, 3.579233755969986*^9}, 3.5792356203360605`*^9}],

Cell[BoxData[
 GraphicsBox[{CircleBox[{-1.4623110743514225`, 0}, 0.26224005195523936`], 
   DiskBox[{-1.4623110743514225`, 0}, 0.1], 
   CircleBox[{1.4623110743514225`, 0}, 0.26224005195523936`], 
   DiskBox[{1.4623110743514225`, 0}, 0.1], 
   CircleBox[{0, 0}, 1.2000710223961832`], 
   CircleBox[{-1.9703808943759675`, -1.4370384602223387`}, 1.261969408118396],
    DiskBox[{-2.1517563502467008`, -1.9500457193703256`}, 0.1], 
   CircleBox[{1.9703808943759675`, 1.4370384602223387`}, 1.261969408118396], 
   DiskBox[{2.1517563502467008`, 1.9500457193703256`}, 0.1], 
   {Dashing[{0.02, 0.02}], 
    CircleBox[{1.9703808943759675`, -0.7185192301111694}, 1.261969408118396], 
    DiskBox[{2.1517563502467008`, -0.9750228596851628}, 0.1], 
    CircleBox[{-1.9703808943759675`, 0.7185192301111694}, 1.261969408118396], 
    DiskBox[{-2.1517563502467008`, 0.9750228596851628}, 0.1]}, 
   {Thickness[0.01], 
    LineBox[{{-1.4623110743514225`, 0}, {1.4623110743514225`, 0}}], 
    LineBox[{{1.4623110743514225`, 0}, {2.1517563502467008`, 
      1.9500457193703256`}}], 
    LineBox[{{-1.4623110743514225`, 
       0}, {-2.1517563502467008`, -1.9500457193703256`}}], 
    {Dashing[{0.01, 0.01}], 
     LineBox[{{1.4623110743514225`, 0}, {
       2.1517563502467008`, -0.9750228596851628}}], 
     LineBox[{{-1.4623110743514225`, 0}, {-2.1517563502467008`, 
       0.9750228596851628}}]}}},
  AspectRatio->Automatic,
  Axes->True,
  Frame->True,
  GridLines->Automatic,
  PlotRange->All]], "Output",
 CellChangeTimes->{3.5792377219585524`*^9}]
}, Open  ]],

Cell[BoxData["\[IndentingNewLine]"], "Input",
 CellChangeTimes->{3.5789832731211567`*^9}]
}, Open  ]]
},
WindowSize->{872, 739},
WindowMargins->{{471, Automatic}, {Automatic, -3}},
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[CellGroupData[{
Cell[579, 22, 735, 12, 33, "Subsubtitle"],
Cell[1317, 36, 2870, 53, 112, "Input"],
Cell[4190, 91, 302, 6, 31, "Text"],
Cell[4495, 99, 2188, 56, 92, "Input"],
Cell[6686, 157, 2276, 54, 92, "Input"],
Cell[8965, 213, 3198, 78, 132, "Input"],
Cell[12166, 293, 1768, 42, 92, "Input"],
Cell[CellGroupData[{
Cell[13959, 339, 1464, 31, 92, "Input"],
Cell[15426, 372, 424, 11, 31, "Output"]
}, Open  ]],
Cell[15865, 386, 1126, 18, 31, "Input"],
Cell[CellGroupData[{
Cell[17016, 408, 1279, 21, 52, "Input"],
Cell[18298, 431, 104, 2, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[18439, 438, 1201, 19, 52, "Input"],
Cell[19643, 459, 87, 1, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[19767, 465, 1325, 22, 52, "Input"],
Cell[21095, 489, 104, 2, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[21236, 496, 1318, 22, 52, "Input"],
Cell[22557, 520, 89, 1, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[22683, 526, 1286, 21, 52, "Input"],
Cell[23972, 549, 87, 1, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[24096, 555, 1654, 33, 72, "Input"],
Cell[25753, 590, 88, 1, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[25878, 596, 1649, 33, 72, "Input"],
Cell[27530, 631, 88, 1, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[27655, 637, 1400, 26, 52, "Input"],
Cell[29058, 665, 88, 1, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[29183, 671, 1687, 34, 52, "Input"],
Cell[30873, 707, 87, 1, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[30997, 713, 7951, 209, 372, "Input"],
Cell[38951, 924, 1537, 31, 320, "Output"]
}, Open  ]],
Cell[40503, 958, 89, 1, 52, "Input"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)