(* 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[     44855,       1066]
NotebookOptionsPosition[     42819,        995]
NotebookOutlinePosition[     43303,       1013]
CellTagsIndexPosition[     43260,       1010]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["\<\
C2H6 Molecule after Kimball, non overlapping, parameterized with MP2/6-31G(d)\
\[AliasDelimiter]\[AliasDelimiter]\
\>", "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}, {3.57924170619635*^9, 3.579241706555151*^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[{
  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, {3.579240085619104*^9, 
   3.579240093216317*^9}, {3.579240186317281*^9, 3.579240221245742*^9}, {
   3.5792416408634357`*^9, 3.5792416835139103`*^9}, 3.5792418843954635`*^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}, {
   3.579239787767781*^9, 3.579239823413843*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"Vir", "=", 
    RowBox[{
     RowBox[{"-", 
      RowBox[{"(", 
       RowBox[{"Vee", "+", "Vne", "+", "Vnn"}], ")"}]}], "/", "T"}]}], ";"}], 
  "\[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}, {
   3.579239787767781*^9, 3.579239823413843*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"func", " ", "=", " ", 
   RowBox[{
   "T", " ", "+", " ", "Vne", " ", "+", " ", "Vee", " ", "+", " ", "Vnn"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"t", " ", "=", " ", 
  RowBox[{"NMinimize", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"func", ",", 
      RowBox[{"Vir", "\[Equal]", "2.0"}]}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{"P", ",", " ", "0.25", ",", " ", "0.27"}], "}"}], ",", " ", 
      RowBox[{"{", 
       RowBox[{"Q", ",", " ", "1.19", ",", " ", "1.23"}], "}"}], ",", " ", 
      "\n", "   ", 
      RowBox[{"{", 
       RowBox[{"R", ",", " ", "1.24", ",", " ", "1.28"}], "}"}], ",", " ", 
      RowBox[{"{", 
       RowBox[{"X", ",", " ", "1.34", ",", " ", "1.37"}], "}"}]}], "}"}]}], 
   "]"}]}], "\[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.579238017719472*^9, 3.5792380215726786`*^9}, {3.579238063037551*^9, 
   3.5792381365136805`*^9}, {3.5792383627608776`*^9, 
   3.5792384035705495`*^9}, {3.579238444427021*^9, 3.579238476875078*^9}, {
   3.5792385863872705`*^9, 3.579238636338558*^9}, {3.579238680533436*^9, 
   3.5792387403283405`*^9}, {3.579238925812667*^9, 3.5792389587287245`*^9}, {
   3.5792390123772187`*^9, 3.5792390173224277`*^9}, {3.579239059130501*^9, 
   3.579239067835316*^9}, {3.579239175990306*^9, 3.579239176988708*^9}, {
   3.5792396087730665`*^9, 3.5792396113470707`*^9}, {3.5792396550583477`*^9, 
   3.579239678083988*^9}, {3.579239718129258*^9, 3.5792397203444624`*^9}, {
   3.57923983269586*^9, 3.579239865268717*^9}, {3.579240600888009*^9, 
   3.579240673256536*^9}, {3.579240710556202*^9, 3.579240711445403*^9}, 
   3.579241293248025*^9, {3.5792413542909327`*^9, 3.579241380077778*^9}, {
   3.579241480526354*^9, 3.5792415281532373`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", "79.4941299695062`"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"P", "\[Rule]", "0.2622400004544331`"}], ",", 
     RowBox[{"Q", "\[Rule]", "1.2000677311822714`"}], ",", 
     RowBox[{"R", "\[Rule]", "1.2619699087119165`"}], ",", 
     RowBox[{"X", "\[Rule]", "1.356988625309145`"}]}], "}"}]}], 
  "}"}]], "Output",
 CellChangeTimes->{3.579241969540413*^9, 3.5792656200135965`*^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.4941299695062`"}]], "Output",
 CellChangeTimes->{3.5792419695872135`*^9, 3.579265620060396*^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.4941299695062`"], "Output",
 CellChangeTimes->{3.5792419696028132`*^9, 3.5792656200759964`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{" ", 
  RowBox[{"(*", 
   RowBox[{"NUCLEUS", "-", 
    RowBox[{"ELECTRON", " ", "POTENTIAL", " ", "ENERGY"}]}], "  ", "*)"}], 
  "\[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.462120556669`"}]], "Output",
 CellChangeTimes->{3.579241969618413*^9, 3.5792656200915966`*^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.42997476709529`"], "Output",
 CellChangeTimes->{3.579241969634013*^9, 3.5792656201071963`*^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.0438858505613`"], "Output",
 CellChangeTimes->{3.579241969634013*^9, 3.5792656201071963`*^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.0945155798166053`"], "Output",
 CellChangeTimes->{3.5792419696496134`*^9, 3.579265620122796*^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.547639237008633`"], "Output",
 CellChangeTimes->{3.579241969665213*^9, 3.5792656201383963`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"Virial", " ", "Theorem"}], " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"Vir", "/.", "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, {3.579240004389761*^9, 
   3.579240005544163*^9}}],

Cell[BoxData["2.`"], "Output",
 CellChangeTimes->{3.5792419696808133`*^9, 3.5792656201539965`*^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.2287489528524`"], "Output",
 CellChangeTimes->{3.579241969696413*^9, 3.5792656201539965`*^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.4623077316367046`, 0}, 0.2622400004544331], 
   DiskBox[{-1.4623077316367046`, 0}, 0.1], 
   CircleBox[{1.4623077316367046`, 0}, 0.2622400004544331], 
   DiskBox[{1.4623077316367046`, 0}, 0.1], 
   CircleBox[{0, 0}, 1.2000677311822714`], 
   CircleBox[{-1.970377701358821, -1.43703888363101}, 1.2619699087119165`], 
   DiskBox[{-2.151752901410778, -1.9500454192142325`}, 0.1], 
   CircleBox[{1.970377701358821, 1.43703888363101}, 1.2619699087119165`], 
   DiskBox[{2.151752901410778, 1.9500454192142325`}, 0.1], 
   {Dashing[{0.02, 0.02}], 
    CircleBox[{1.970377701358821, -0.718519441815505}, 1.2619699087119165`], 
    DiskBox[{2.151752901410778, -0.9750227096071162}, 0.1], 
    CircleBox[{-1.970377701358821, 0.718519441815505}, 1.2619699087119165`], 
    DiskBox[{-2.151752901410778, 0.9750227096071162}, 0.1]}, 
   {Thickness[0.01], 
    LineBox[{{-1.4623077316367046`, 0}, {1.4623077316367046`, 0}}], 
    LineBox[{{1.4623077316367046`, 0}, {2.151752901410778, 
      1.9500454192142325`}}], 
    LineBox[{{-1.4623077316367046`, 
       0}, {-2.151752901410778, -1.9500454192142325`}}], 
    {Dashing[{0.01, 0.01}], 
     LineBox[{{1.4623077316367046`, 0}, {
       2.151752901410778, -0.9750227096071162}}], 
     LineBox[{{-1.4623077316367046`, 0}, {-2.151752901410778, 
       0.9750227096071162}}]}}},
  AspectRatio->Automatic,
  Axes->True,
  Frame->True,
  GridLines->Automatic,
  PlotRange->All]], "Output",
 CellChangeTimes->{3.5792419697120132`*^9, 3.579265620185197*^9}]
}, Open  ]],

Cell[BoxData["\[IndentingNewLine]"], "Input",
 CellChangeTimes->{3.5789832731211567`*^9}]
}, Open  ]]
},
WindowSize->{891, 810},
WindowMargins->{{11, Automatic}, {Automatic, -1}},
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, 816, 13, 33, "Subsubtitle"],
Cell[1398, 37, 3041, 56, 112, "Input"],
Cell[4442, 95, 302, 6, 31, "Text"],
Cell[4747, 103, 2188, 56, 92, "Input"],
Cell[6938, 161, 2276, 54, 92, "Input"],
Cell[9217, 217, 3198, 78, 132, "Input"],
Cell[12418, 297, 1818, 43, 92, "Input"],
Cell[14239, 342, 822, 17, 52, "Input"],
Cell[CellGroupData[{
Cell[15086, 363, 2466, 46, 92, "Input"],
Cell[17555, 411, 445, 11, 31, "Output"]
}, Open  ]],
Cell[18015, 425, 1126, 18, 31, "Input"],
Cell[CellGroupData[{
Cell[19166, 447, 1279, 21, 52, "Input"],
Cell[20448, 470, 127, 2, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[20612, 477, 1201, 19, 52, "Input"],
Cell[21816, 498, 112, 1, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[21965, 504, 1295, 21, 52, "Input"],
Cell[23263, 527, 127, 2, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[23427, 534, 1318, 22, 52, "Input"],
Cell[24748, 558, 111, 1, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[24896, 564, 1286, 21, 52, "Input"],
Cell[26185, 587, 110, 1, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[26332, 593, 1654, 33, 72, "Input"],
Cell[27989, 628, 112, 1, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[28138, 634, 1649, 33, 72, "Input"],
Cell[29790, 669, 111, 1, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[29938, 675, 1300, 21, 52, "Input"],
Cell[31241, 698, 98, 1, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[31376, 704, 1687, 34, 52, "Input"],
Cell[33066, 740, 111, 1, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[33214, 746, 7951, 209, 372, "Input"],
Cell[41168, 957, 1531, 31, 320, "Output"]
}, Open  ]],
Cell[42714, 991, 89, 1, 52, "Input"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)