(* 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[     43559,       1010]
NotebookOptionsPosition[     41793,        947]
NotebookOutlinePosition[     42275,        965]
CellTagsIndexPosition[     42232,        962]
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[{
   "z", ",", "\[Sigma]", ",", "k", ",", "nc", ",", "R1", ",", "R2", ",", "p", 
    ",", "Ekin", ",", "Epot", ",", "vne", ",", "vnn", ",", "vee", ",", "vir", 
    ",", "ow", ",", "cn", ",", "n", ",", "w", ",", "c", ",", "c2", ",", "ec", 
    ",", "\[IndentingNewLine]", "dwnr", ",", "dwr", ",", "dnr", ",", "d", ",",
     "d1", ",", "d2", ",", "pi", ",", "i", ",", "j", ",", "t", ",", 
    "\[Lambda]", ",", "cs", ",", "ss", ",", "wi"}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"bohr", "=", "0.529177"}], ";"}], " ", "\[IndentingNewLine]", " ", 
  RowBox[{"(*", 
   RowBox[{"Parameterlist", " ", "c"}], " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{"(*", 
   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"}]}], "}"}]}], ";"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"k", "=", 
    RowBox[{"{", 
     RowBox[{
     "1.02669238", ",", "1.02669238", ",", "1.36467050", ",", "1.20088548", 
      ",", "1.20088548", ",", "1.20088548", ",", "1.20088548", ",", 
      "1.20088548", ",", "1.20088548"}], "}"}]}], ";"}], " "}], "\n", 
 RowBox[{
  RowBox[{"\[Sigma]", "=", 
   RowBox[{"{", 
    RowBox[{
    "0.30102681", ",", "0.30102681", ",", "0.34863843", ",", "0.34858036", 
     ",", "0.34858036", ",", "0.34858036", ",", "0.34858036", ",", 
     "0.34858036", ",", "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, {
   3.579254526481946*^9, 3.579254528197949*^9}, {3.579256065565049*^9, 
   3.5792562478981695`*^9}, {3.579256291921447*^9, 3.579256298769859*^9}, {
   3.5792570670400085`*^9, 3.5792571152284927`*^9}, {3.579257266377158*^9, 
   3.5792572671259594`*^9}, {3.5792573858077683`*^9, 
   3.5792574224522324`*^9}, {3.579258705944487*^9, 3.579258744117754*^9}, {
   3.5792587970174465`*^9, 3.5792588039126587`*^9}, 3.5792608573214664`*^9, {
   3.579260900939143*^9, 3.579260918489174*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"Terms", " ", "analogous", " ", "to", " ", "CH4_mat"}], " ", 
   "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"z", "=", "6.0"}], ";", " ", 
    RowBox[{"nc", "=", "9"}], ";", " ", 
    RowBox[{"n", "=", "2"}], ";", 
    RowBox[{"cs", "=", 
     RowBox[{"Cos", "[", "wi", "]"}]}], ";", " ", 
    RowBox[{"ss", "=", 
     RowBox[{"Sin", "[", "wi", "]"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"dwr", "=", 
     RowBox[{"Table", "[", 
      RowBox[{"0", ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "nc"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"j", ",", "1", ",", "nc"}], "}"}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"dnr", "=", "dwr"}], ";", 
    RowBox[{"dwnr", "=", "dwr"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"c2", "=", 
     RowBox[{"Table", "[", 
      RowBox[{"0", ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "nc"}], "}"}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"ow", "=", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"-", "2"}], ",", 
       RowBox[{"-", "2"}], ",", 
       RowBox[{"-", "2"}], ",", 
       RowBox[{"-", "2"}], ",", 
       RowBox[{"-", "2"}], ",", 
       RowBox[{"-", "2"}], ",", 
       RowBox[{"-", "2"}], ",", 
       RowBox[{"-", "2"}], ",", 
       RowBox[{"-", "2"}]}], "}"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"cn", "=", 
     RowBox[{"{", 
      RowBox[{
      "z", ",", "z", ",", "0", ",", "1", ",", "1", ",", "1", ",", "1", ",", 
       "1", ",", "1"}], "}"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"c", "=", 
     RowBox[{"{", 
      RowBox[{
      "R1", ",", "R1", ",", "R2", ",", "R3", ",", "R3", ",", "R3", ",", "R3", 
       ",", "R3", ",", "R3"}], "}"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"ec", "=", 
     RowBox[{"{", 
      RowBox[{
      "0", ",", "0", ",", "0", ",", "p", ",", "p", ",", "p", ",", "p", ",", 
       "p", ",", "p"}], "}"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"d1", "=", 
     RowBox[{"R1", "+", "R2"}]}], ";", " ", 
    RowBox[{"d2", "=", 
     RowBox[{"d1", "+", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"R1", "+", "R3"}], ")"}], "*", "cs"}]}]}], ";", " ", 
    RowBox[{"d3", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"R1", "+", "R3"}], ")"}], "*", "ss"}]}], ";", " ", 
    RowBox[{"d4", "=", " ", 
     RowBox[{"d3", "/", "2"}]}], ";", " ", 
    RowBox[{"d5", "=", " ", 
     RowBox[{"d3", "*", 
      RowBox[{
       RowBox[{"Sqrt", "[", "3", "]"}], "/", "2"}]}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"d6", "=", 
     RowBox[{"d1", "+", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"R1", "+", "R3", "+", "p"}], ")"}], "*", "cs"}]}]}], ";", " ", 
    RowBox[{"d7", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"R1", "+", "R3", "+", "p"}], ")"}], "*", "ss"}]}], ";", " ", 
    RowBox[{"d8", "=", 
     RowBox[{"d7", "/", "2"}]}], ";", " ", 
    RowBox[{"d9", "=", 
     RowBox[{"d7", "*", 
      RowBox[{
       RowBox[{"Sqrt", "[", "3", "]"}], "/", "2"}]}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"w", " ", "=", " ", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"d1", ",", " ", "0", ",", " ", "0"}], "}"}], ",", " ", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "d1"}], ",", " ", "0", ",", "0"}], "}"}], ",", " ", 
       RowBox[{"{", 
        RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", "\[IndentingNewLine]", 
       RowBox[{"{", 
        RowBox[{"d2", ",", " ", "d3", ",", " ", "0"}], "}"}], ",", " ", 
       RowBox[{"{", 
        RowBox[{"d2", ",", " ", 
         RowBox[{"-", "d4"}], ",", " ", "d5"}], "}"}], ",", " ", 
       RowBox[{"{", 
        RowBox[{"d2", ",", " ", 
         RowBox[{"-", "d4"}], ",", " ", 
         RowBox[{"-", "d5"}]}], "}"}], ",", "\[IndentingNewLine]", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "d2"}], ",", 
         RowBox[{"-", "d3"}], ",", "0"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "d2"}], ",", "d4", ",", "d5"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "d2"}], ",", "d4", ",", 
         RowBox[{"-", "d5"}]}], "}"}]}], "}"}]}], ";"}], " ", "\n", 
   RowBox[{
    RowBox[{"n", " ", "=", " ", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"d1", ",", " ", "0", ",", " ", "0"}], "}"}], ",", " ", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "d1"}], ",", " ", "0", ",", "0"}], "}"}], ",", " ", 
       RowBox[{"{", 
        RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", "\[IndentingNewLine]", 
       RowBox[{"{", 
        RowBox[{"d6", ",", " ", "d7", ",", " ", "0"}], "}"}], ",", " ", 
       RowBox[{"{", 
        RowBox[{"d6", ",", " ", 
         RowBox[{"-", "d8"}], ",", " ", "d9"}], "}"}], ",", " ", 
       RowBox[{"{", 
        RowBox[{"d6", ",", " ", 
         RowBox[{"-", "d8"}], ",", " ", 
         RowBox[{"-", "d9"}]}], "}"}], ",", "\[IndentingNewLine]", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "d6"}], ",", 
         RowBox[{"-", "d7"}], ",", "0"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "d6"}], ",", "d8", ",", "d9"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "d6"}], ",", "d8", ",", 
         RowBox[{"-", "d9"}]}], "}"}]}], "}"}]}], ";"}], " ", 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"For", "[", 
     RowBox[{
      RowBox[{"i", " ", "=", " ", "1"}], ",", " ", 
      RowBox[{"i", " ", "\[LessEqual]", " ", "nc"}], ",", " ", 
      RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{
        RowBox[{"c2", "[", 
         RowBox[{"[", "i", "]"}], "]"}], " ", "=", " ", 
        RowBox[{
         RowBox[{"k", "[", 
          RowBox[{"[", "i", "]"}], "]"}], "/", 
         RowBox[{
          RowBox[{"c", "[", 
           RowBox[{"[", "i", "]"}], "]"}], "^", "2"}]}]}], ";", 
       "\[IndentingNewLine]", 
       RowBox[{"For", "[", 
        RowBox[{
         RowBox[{"j", " ", "=", " ", "1"}], ",", " ", 
         RowBox[{"j", " ", "\[LessEqual]", " ", "nc"}], ",", " ", 
         RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", 
         RowBox[{
          RowBox[{"d", " ", "=", " ", 
           RowBox[{"Norm", "[", 
            RowBox[{
             RowBox[{"w", "[", 
              RowBox[{"[", "i", "]"}], "]"}], " ", "-", " ", 
             RowBox[{"w", "[", 
              RowBox[{"[", "j", "]"}], "]"}]}], "]"}]}], ";", 
          "\[IndentingNewLine]", 
          RowBox[{
           RowBox[{"dwr", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "j"}], "]"}], "]"}], " ", "=", " ", 
           RowBox[{"If", "[", 
            RowBox[{
             RowBox[{"i", " ", "\[NotEqual]", " ", "j"}], ",", " ", 
             RowBox[{"1", "/", "d"}], ",", 
             RowBox[{"1.5", "*", 
              RowBox[{"(", 
               RowBox[{
                RowBox[{"\[Sigma]", "[", 
                 RowBox[{"[", "i", "]"}], "]"}], "/", 
                RowBox[{"c", "[", 
                 RowBox[{"[", "i", "]"}], "]"}]}], ")"}]}]}], "]"}]}], ";", 
          "\[IndentingNewLine]", 
          RowBox[{"d", "=", " ", 
           RowBox[{"Norm", "[", 
            RowBox[{
             RowBox[{"n", "[", 
              RowBox[{"[", "i", "]"}], "]"}], " ", "-", " ", 
             RowBox[{"n", "[", 
              RowBox[{"[", "j", "]"}], "]"}]}], "]"}]}], ";", 
          "\[IndentingNewLine]", 
          RowBox[{
           RowBox[{"dnr", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "j"}], "]"}], "]"}], " ", "=", " ", 
           RowBox[{"If", "[", 
            RowBox[{
             RowBox[{"i", " ", "\[NotEqual]", " ", "j"}], ",", " ", 
             RowBox[{"1", "/", "d"}], ",", " ", "0"}], "]"}]}], ";", 
          "\[IndentingNewLine]", 
          RowBox[{"d", "=", " ", 
           RowBox[{"Norm", "[", 
            RowBox[{
             RowBox[{"w", "[", 
              RowBox[{"[", "i", "]"}], "]"}], " ", "-", " ", 
             RowBox[{"n", "[", 
              RowBox[{"[", "j", "]"}], "]"}]}], "]"}]}], ";", 
          "\[IndentingNewLine]", 
          RowBox[{
           RowBox[{"dwnr", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "j"}], "]"}], "]"}], " ", "=", " ", 
           RowBox[{"If", "[", 
            RowBox[{
             RowBox[{"i", " ", "\[NotEqual]", " ", "j"}], ",", " ", 
             RowBox[{"1", "/", "d"}], ",", " ", 
             RowBox[{
              RowBox[{"(", 
               RowBox[{"1.5", " ", "-", " ", 
                RowBox[{"0.5", "*", 
                 RowBox[{
                  RowBox[{"(", 
                   RowBox[{
                    RowBox[{"ec", "[", 
                    RowBox[{"[", "i", "]"}], "]"}], "/", 
                    RowBox[{"c", "[", 
                    RowBox[{"[", "i", "]"}], "]"}]}], ")"}], "^", "2"}]}]}], 
               ")"}], "/", 
              RowBox[{"c", "[", 
               RowBox[{"[", "i", "]"}], "]"}]}]}], "]"}]}]}]}], "]"}]}]}], 
     "]"}], ";"}], " ", "\n", 
   RowBox[{
    RowBox[{"Ekin", " ", "=", " ", 
     RowBox[{
      RowBox[{"-", 
       RowBox[{"(", 
        RowBox[{"9", "/", "8"}], ")"}]}], "*", 
      RowBox[{"ow", " ", ".", " ", "c2"}]}]}], ";"}], " ", "\n", 
   RowBox[{
    RowBox[{"vee", " ", "=", " ", 
     RowBox[{"0.5", "*", 
      RowBox[{"ow", " ", ".", " ", 
       RowBox[{"(", 
        RowBox[{"ow", " ", ".", " ", "dwr"}], ")"}]}]}]}], ";"}], " ", "\n", 
   RowBox[{
    RowBox[{"vnn", " ", "=", " ", 
     RowBox[{"0.5", "*", 
      RowBox[{"cn", " ", ".", " ", 
       RowBox[{"(", 
        RowBox[{"cn", " ", ".", " ", "dnr"}], ")"}]}]}]}], ";"}], " ", "\n", 
   RowBox[{
    RowBox[{"vne", " ", "=", " ", 
     RowBox[{"cn", " ", ".", " ", 
      RowBox[{"(", 
       RowBox[{"ow", " ", ".", " ", "dwnr"}], ")"}]}]}], ";"}], " ", "\n", 
   RowBox[{
    RowBox[{"Epot", " ", "=", " ", 
     RowBox[{"vne", " ", "+", " ", "vee", " ", "+", " ", "vnn"}]}], ";"}], 
   " ", "\n", 
   RowBox[{
    RowBox[{"func", " ", "=", " ", 
     RowBox[{"Ekin", " ", "+", " ", "Epot"}]}], ";"}], " "}]}]], "Text", \
"Input",
 CellChangeTimes->{{3.5792338519413548`*^9, 3.5792338963078327`*^9}, {
   3.579233998924813*^9, 3.579234002029218*^9}, 3.579235530823104*^9, {
   3.579254579459639*^9, 3.579254604497683*^9}, {3.5792546701893983`*^9, 
   3.5792549883675575`*^9}, {3.579255040643249*^9, 3.5792550505492663`*^9}, {
   3.579255231415984*^9, 3.5792553234093456`*^9}, {3.579255363563816*^9, 
   3.5792558209098196`*^9}, {3.5792558541534777`*^9, 3.579255854403078*^9}, {
   3.5792559106099772`*^9, 3.57925596911008*^9}, 3.5792560056141443`*^9, {
   3.5792566762281218`*^9, 3.579256679098527*^9}, 3.5792567540254583`*^9, {
   3.5792579024650755`*^9, 3.5792579089702873`*^9}, {3.579258007921261*^9, 
   3.57925801311607*^9}, {3.5792585291025763`*^9, 3.579258546933408*^9}, {
   3.5792589863081794`*^9, 3.579258990551387*^9}, {3.57925925426985*^9, 
   3.579259268918276*^9}, {3.5792608668842835`*^9, 3.579260889363923*^9}, {
   3.579260931390397*^9, 3.579260934323202*^9}, {3.5792785711210957`*^9, 
   3.5792785751615033`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"vir", "=", 
   RowBox[{
    RowBox[{"-", "Epot"}], "/", "Ekin"}]}], ";"}]], "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}, {3.579256422977277*^9, 
   3.5792564375789022`*^9}, 3.579256757613465*^9, {3.5792569515842056`*^9, 
   3.5792569528634076`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"func", " ", "=", 
   RowBox[{"Ekin", "+", "Epot"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"t", " ", "=", " ", 
  RowBox[{"NMinimize", "[", 
   RowBox[{
    RowBox[{"{", "func", "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{"R1", ",", " ", "0.25", ",", " ", "0.27"}], "}"}], ",", " ", 
      RowBox[{"{", 
       RowBox[{"R2", ",", " ", "1.19", ",", " ", "1.22"}], "}"}], ",", " ", 
      "\n", "   ", 
      RowBox[{"{", 
       RowBox[{"R3", ",", " ", "1.24", ",", " ", "1.28"}], "}"}], ",", " ", 
      RowBox[{"{", 
       RowBox[{"p", ",", " ", "0.53", ",", "0.56"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"wi", ",", "1.2", ",", "1.25"}], "}"}]}], "}"}]}], 
   "]"}]}], "\[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}, {3.579256407065249*^9, 
   3.5792564083132515`*^9}, {3.5792568138671637`*^9, 3.579256890338498*^9}, {
   3.57925694283259*^9, 3.5792569442365923`*^9}, {3.5792572545055375`*^9, 
   3.579257255191939*^9}, 3.579257291852003*^9, {3.579260811769387*^9, 
   3.579260834217826*^9}, {3.5792609483476267`*^9, 3.5792609527312346`*^9}, {
   3.5792609956469097`*^9, 3.5792610021521215`*^9}, 3.5792610533046107`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", "79.4941838100811`"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"R1", "\[Rule]", "0.2622399840848178`"}], ",", 
     RowBox[{"R2", "\[Rule]", "1.2000052975561055`"}], ",", 
     RowBox[{"R3", "\[Rule]", "1.261980036069545`"}], ",", 
     RowBox[{"p", "\[Rule]", "0.5440916298444888`"}], ",", 
     RowBox[{"wi", "\[Rule]", "1.2270391929835192`"}]}], "}"}]}], 
  "}"}]], "Output",
 CellChangeTimes->{3.579261065956233*^9, 3.5792785816667147`*^9, 
  3.5792787906758814`*^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.4941838100811`"}]], "Output",
 CellChangeTimes->{3.579261065971833*^9, 3.5792785816823144`*^9, 
  3.5792787906914816`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"Kinetic", " ", "energy"}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{"Ekin", " ", "/.", "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, {3.579256527825061*^9, 
   3.5792565289950633`*^9}}],

Cell[BoxData["79.49419683114962`"], "Output",
 CellChangeTimes->{3.5792610659874334`*^9, 3.5792785816979146`*^9, 
  3.5792787906914816`*^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, {
   3.579256537949479*^9, 3.5792565381210794`*^9}}],

Cell[BoxData[
 RowBox[{"-", "270.426103019044`"}]], "Output",
 CellChangeTimes->{3.579261066003033*^9, 3.5792785817135143`*^9, 
  3.579278790707082*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"ELECTRON", "-", 
    RowBox[{"ELECTRON", " ", "POTENTIAL", " ", "ENERGY"}]}], " ", "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{"vee", " ", "/.", "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.5792319350255876`*^9}, 3.5792319781128635`*^9, {
   3.579233697766284*^9, 3.5792337003714886`*^9}, 3.579235564737563*^9, {
   3.5792565420990863`*^9, 3.579256542270686*^9}, 3.5792566391000566`*^9}],

Cell[BoxData["69.40599530071324`"], "Output",
 CellChangeTimes->{3.579261066018633*^9, 3.5792785817291145`*^9, 
  3.579278790707082*^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, {
   3.579256545499892*^9, 3.5792565456870923`*^9}}],

Cell[BoxData["42.031727077100044`"], "Output",
 CellChangeTimes->{3.579261066034233*^9, 3.5792785817291145`*^9, 
  3.5792787907226815`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{" ", 
  RowBox[{"(*", 
   RowBox[{"C", "-", 
    RowBox[{"H", " ", "distance", "  ", "\[CapitalARing]"}]}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"ch", " ", "=", 
     RowBox[{
      RowBox[{"R1", "+", "R3", "+", "p"}], "/.", "u"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{"dch", " ", "=", " ", 
    RowBox[{
    "ch", "*", "bohr", 
     "                                             "}]}]}]}]], "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, {3.5792565803815536`*^9, 
   3.5792566233284287`*^9}, {3.579257470812318*^9, 3.5792574914355536`*^9}}],

Cell[BoxData["1.0945029540114424`"], "Output",
 CellChangeTimes->{3.5792610660498333`*^9, 3.5792785817447147`*^9, 
  3.5792787907226815`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{" ", 
  RowBox[{"(*", 
   RowBox[{"C", "-", 
    RowBox[{"C", " ", "distance", "  ", "\[CapitalARing]"}]}], " ", "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"cc", " ", "=", "  ", 
     RowBox[{
      RowBox[{"2", "*", "d1"}], "  ", "/.", "u"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{"dcc", "=", 
    RowBox[{
    "cc", "*", "bohr", 
     "                                              "}]}]}]}]], "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, {3.5792566088672037`*^9, 
   3.5792566162616167`*^9}}],

Cell[BoxData["1.5475731428057977`"], "Output",
 CellChangeTimes->{3.5792610660654335`*^9, 3.579278581760315*^9, 
  3.5792787907382812`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"Virial", " ", "Theorem"}], " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{"vir", "/.", "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.579231986490078*^9}, {3.579233726953935*^9, 
   3.5792337305575414`*^9}, 3.579235598917223*^9, {3.579240004389761*^9, 
   3.579240005544163*^9}, {3.5792566498484755`*^9, 3.579256654731284*^9}}],

Cell[BoxData["1.9999998362010167`"], "Output",
 CellChangeTimes->{3.5792610660654335`*^9, 3.579278581775915*^9, 
  3.5792787907538815`*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
   "Atomization", " ", "Energy", " ", "at", " ", "0", " ", "K", " ", "in", 
    " ", 
    RowBox[{"kcal", "/", "mol"}]}], " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "Etot"}], " ", "-", " ", 
      RowBox[{"n", "*", "37.78899"}], " ", "-", 
      RowBox[{
       RowBox[{"(", " ", 
        RowBox[{
         RowBox[{"2", "*", "n"}], " ", "+", " ", "2"}], ")"}], "*", 
       "0.50177"}]}], ")"}], "*", "627.5096"}], " ", "*)"}], 
  "\[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.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, {3.5792569813646574`*^9, 3.579256983189861*^9}, 
   3.5792575835225153`*^9, {3.579259029379855*^9, 3.5792590354014654`*^9}, 
   3.579259348790416*^9}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"plot1", "=", 
   RowBox[{
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"Graphics3D", "[", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"Opacity", "[", "0.6", "]"}], ",", 
         RowBox[{"Sphere", "[", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{"w", "[", 
             RowBox[{"[", "i", "]"}], "]"}], "}"}], ",", 
           RowBox[{"c", "[", 
            RowBox[{"[", "i", "]"}], "]"}]}], "]"}], ",", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"Darker", "[", 
            RowBox[{"Red", ",", "1.0"}], "]"}], ",", 
           RowBox[{"If", "[", 
            RowBox[{
             RowBox[{"i", "\[NotEqual]", "3"}], ",", 
             RowBox[{"Sphere", "[", 
              RowBox[{
               RowBox[{"{", 
                RowBox[{"n", "[", 
                 RowBox[{"[", "i", "]"}], "]"}], "}"}], ",", "0.08"}], 
              "]"}]}], "]"}]}], "}"}]}], "}"}], " ", "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "nc"}], "}"}]}], "]"}], "   ", "/.", " ", 
    "u"}]}], "  ", ";"}], "\n", 
 RowBox[{"Show", "[", 
  RowBox[{"plot1", ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"AspectRatio", " ", "\[Rule]", " ", "Automatic"}], ",", 
     RowBox[{"Axes", " ", "->", " ", "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.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, {3.5792569813646574`*^9, 3.579256983189861*^9}, 
   3.5792575835225153`*^9, 3.579257631055799*^9, {3.5792787364189863`*^9, 
   3.579278783484269*^9}}],

Cell[BoxData[
 Graphics3DBox[{
   {Opacity[0.6], 
    SphereBox[{{1.4622452816409233`, 0, 0}}, 0.2622399840848178], 
    {RGBColor[0., 0., 0.], SphereBox[{{1.4622452816409233`, 0, 0}}, 0.08]}}, 
   {Opacity[0.6], 
    SphereBox[{{-1.4622452816409233`, 0, 0}}, 0.2622399840848178], 
    {RGBColor[0., 0., 0.], SphereBox[{{-1.4622452816409233`, 0, 0}}, 0.08]}}, 
   {Opacity[0.6], SphereBox[{{0, 0, 0}}, 1.2000052975561055`], 
    {RGBColor[0., 0., 0.]}}, 
   {Opacity[0.6], 
    SphereBox[{{1.975948255396272, 1.4350456176004571`, 0}}, 
     1.261980036069545], 
    {RGBColor[0., 0., 0.], 
     SphereBox[{{2.1593217084087852`, 1.9473051986695669`, 0}}, 0.08]}}, 
   {Opacity[0.6], 
    SphereBox[{{1.975948255396272, -0.7175228088002286, 1.242785960431525}}, 
     1.261980036069545], 
    {RGBColor[0., 0., 0.], 
     SphereBox[{{2.1593217084087852`, -0.9736525993347834, 
      1.6864157709693481`}}, 0.08]}}, 
   {Opacity[0.6], 
    SphereBox[{{1.975948255396272, -0.7175228088002286, -1.242785960431525}}, 
     1.261980036069545], 
    {RGBColor[0., 0., 0.], 
     SphereBox[{{
      2.1593217084087852`, -0.9736525993347834, -1.6864157709693481`}}, 
      0.08]}}, 
   {Opacity[0.6], 
    SphereBox[{{-1.975948255396272, -1.4350456176004571`, 0}}, 
     1.261980036069545], 
    {RGBColor[0., 0., 0.], 
     SphereBox[{{-2.1593217084087852`, -1.9473051986695669`, 0}}, 0.08]}}, 
   {Opacity[0.6], 
    SphereBox[{{-1.975948255396272, 0.7175228088002286, 1.242785960431525}}, 
     1.261980036069545], 
    {RGBColor[0., 0., 0.], 
     SphereBox[{{-2.1593217084087852`, 0.9736525993347834, 
      1.6864157709693481`}}, 0.08]}}, 
   {Opacity[0.6], 
    SphereBox[{{-1.975948255396272, 0.7175228088002286, -1.242785960431525}}, 
     1.261980036069545], 
    {RGBColor[0., 0., 0.], 
     SphereBox[{{-2.1593217084087852`, 
      0.9736525993347834, -1.6864157709693481`}}, 0.08]}}},
  AspectRatio->Automatic,
  AutomaticImageSize->True,
  Axes->True,
  ImageSize->{320.7506898952323, 298.256308784256},
  ViewPoint->{-0.6538785064850376, -2.8134307880074383`, 1.7626826429759674`},
  ViewVertical->{-0.14343321182259303`, -0.07033009218808299, 
   1.2770970335871166`}]], "Output",
 CellChangeTimes->{3.5792610660810337`*^9, 3.579278581791515*^9, 
  3.5792787907694817`*^9}]
}, Open  ]]
}, Open  ]]
},
WindowSize->{808, 810},
WindowMargins->{{7, 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, 4238, 77, 212, "Input"],
Cell[5639, 116, 11335, 303, 672, "Text"],
Cell[16977, 421, 843, 15, 31, "Input"],
Cell[CellGroupData[{
Cell[17845, 440, 2822, 49, 92, "Input"],
Cell[20670, 491, 536, 13, 52, "Output"]
}, Open  ]],
Cell[21221, 507, 1126, 18, 31, "Input"],
Cell[CellGroupData[{
Cell[22372, 529, 1279, 21, 52, "Input"],
Cell[23654, 552, 154, 3, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[23845, 560, 1256, 20, 52, "Input"],
Cell[25104, 582, 140, 2, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[25281, 589, 1347, 22, 52, "Input"],
Cell[26631, 613, 152, 3, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[26820, 621, 1369, 22, 52, "Input"],
Cell[28192, 645, 136, 2, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[28365, 652, 1338, 22, 52, "Input"],
Cell[29706, 676, 139, 2, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[29882, 683, 1687, 32, 72, "Input"],
Cell[31572, 717, 141, 2, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[31750, 724, 1634, 32, 72, "Input"],
Cell[33387, 758, 139, 2, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[33563, 765, 1312, 20, 52, "Input"],
Cell[34878, 787, 139, 2, 31, "Output"]
}, Open  ]],
Cell[35032, 792, 1841, 36, 72, "Input"],
Cell[CellGroupData[{
Cell[36898, 832, 2585, 56, 92, "Input"],
Cell[39486, 890, 2279, 53, 313, "Output"]
}, Open  ]]
}, Open  ]]
}
]
*)

(* End of internal cache information *)