(* 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[    135842,       3366]
NotebookOptionsPosition[    130688,       3205]
NotebookOutlinePosition[    131223,       3226]
CellTagsIndexPosition[    131180,       3223]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{
     RowBox[{
     "Kimball", " ", "model", " ", "and", " ", "energy", " ", "derived", " ", 
      "from", " ", 
      RowBox[{"MP2", "/", "6"}]}], "-", 
     RowBox[{"31", "G", 
      RowBox[{"(", "d", ")"}], " ", "data", " ", "source"}]}], ",", " ", 
    RowBox[{"ES", " ", "04", " ", "July", " ", "2013"}]}], " ", "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"Clear", "[", 
    RowBox[{"a", ",", "u", ",", "u1", ",", "u2", ",", "t"}], "]"}], 
   ";"}]}]], "Input",
 CellChangeTimes->{{3.568278170270548*^9, 3.568278174310955*^9}, {
   3.5682782109554195`*^9, 3.5682782183342323`*^9}, 3.568284752749406*^9, 
   3.5682847957898817`*^9, {3.5682872948362713`*^9, 3.568287369232802*^9}, 
   3.5682960154833097`*^9, {3.568305453232698*^9, 3.5683054626863146`*^9}, {
   3.568307776310778*^9, 3.56830777732478*^9}, {3.5804443791490765`*^9, 
   3.58044438139548*^9}, {3.5804469366841097`*^9, 3.5804469795061855`*^9}, {
   3.580459588972727*^9, 3.580459618113578*^9}, {3.5805250613979683`*^9, 
   3.5805250682307806`*^9}, {3.580533517520193*^9, 3.5805335275198107`*^9}, {
   3.580623964920255*^9, 3.5806239657314568`*^9}, {3.5819501929229097`*^9, 
   3.581950196339316*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"t", "=", 
   RowBox[{"1.889726", "*", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{"0.0000000009", ",", "0.0000000019", ",", 
        RowBox[{"-", "0.0793743911"}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"-", "0.0000000014"}], ",", 
        RowBox[{"-", "0.9398663736"}], ",", "0.3093075561"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"-", "0.8139481308"}], ",", "0.4699331753", ",", 
        "0.3093075629"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"0.8139481313", ",", "0.4699331772", ",", "0.3093075625"}], 
       "}"}]}], "}"}]}]}], ";"}]], "Input",
 CellChangeTimes->{{3.567862912170579*^9, 3.567862949329844*^9}, {
   3.5678630095303497`*^9, 3.5678630107471523`*^9}, {3.5679573013548217`*^9, 
   3.5679573132576427`*^9}, 3.568114367533684*^9, {3.5682782383334675`*^9, 
   3.5682782435438766`*^9}, 3.5682847420789876`*^9, 3.568284795821082*^9, {
   3.568287382898426*^9, 3.5682873928356433`*^9}, 3.5682960154833097`*^9, 
   3.5683055034491863`*^9, {3.5804444043275204`*^9, 3.580444406558325*^9}, {
   3.5804469063576565`*^9, 3.58044691935248*^9}, 3.580447101826*^9, 
   3.58048985172836*^9, {3.5805250473735437`*^9, 3.58052505082115*^9}, {
   3.5805257978758163`*^9, 3.5805257998726196`*^9}, 3.580533417945218*^9, {
   3.5805334506428757`*^9, 3.5805334825293317`*^9}, {3.5819183557945037`*^9, 
   3.581918381690549*^9}, {3.58191857001408*^9, 3.581918605800543*^9}, 
   3.5819193462275963`*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"u", "=", 
  RowBox[{"ChemicalData", "\[NoBreak]", "[", 
   RowBox[{"\"\<Ammonia\>\"", ",", "\"\<VertexTypes\>\""}], "]"}]}]], "Input",
 CellChangeTimes->{{3.567869377256353*^9, 3.5678694101568108`*^9}, {
   3.5678695244894114`*^9, 3.56786953511303*^9}, 3.56786988243204*^9, 
   3.5678699252853155`*^9, {3.567870069179968*^9, 3.5678700766835814`*^9}, {
   3.5678703216664114`*^9, 3.567870325738019*^9}, 3.5679271251072817`*^9, {
   3.567942840466694*^9, 3.567942854381918*^9}, {3.567942932272855*^9, 
   3.567942950509287*^9}, {3.567943087867528*^9, 3.567943093733138*^9}, 
   3.567943185211699*^9, 3.5679432253505697`*^9, 3.5679433277647495`*^9, {
   3.5679434430645523`*^9, 3.5679434440941534`*^9}, {3.568278260766307*^9, 
   3.5682782624355097`*^9}, 3.568284752702606*^9, 3.568284795852282*^9, 
   3.56829601549891*^9, 3.5683055034647865`*^9, {3.5804444900184712`*^9, 
   3.5804444924208755`*^9}, {3.5804469928442087`*^9, 3.580447008101035*^9}, {
   3.5805253295468397`*^9, 3.580525331153642*^9}, {3.5819184024697857`*^9, 
   3.5819184121574025`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"\<\"N\"\>", ",", "\<\"H\"\>", ",", "\<\"H\"\>", ",", "\<\"H\"\>"}],
   "}"}]], "Output",
 CellChangeTimes->{
  3.581932641399249*^9, 3.5819327491174383`*^9, 3.5819327949035187`*^9, {
   3.5819329248049464`*^9, 3.581932952838196*^9}, 3.581950500742651*^9, 
   3.581997263261365*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"n", "=", 
  RowBox[{"Length", "[", "u", "]"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"num1", "=", 
   RowBox[{"Table", "[", 
    RowBox[{"i", ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"u1", "=", 
  RowBox[{"Table", "[", 
   RowBox[{
    RowBox[{"If", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"u", "[", 
        RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\<C\>\""}], ",", 
      "0.255", ",", 
      RowBox[{"If", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"u", "[", 
          RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\<H\>\""}], ",", 
        "0.12", ",", 
        RowBox[{"If", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"u", "[", 
            RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\<N\>\""}], ",", 
          "0.24", ",", 
          RowBox[{"If", "[", 
           RowBox[{
            RowBox[{
             RowBox[{"u", "[", 
              RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\<O\>\""}], ",", 
            "0.18", ",", 
            RowBox[{"If", "[", 
             RowBox[{
              RowBox[{
               RowBox[{"u", "[", 
                RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\<Cl\>\""}], 
              ",", "0.5"}], "]"}]}], "]"}]}], "]"}]}], "]"}]}], "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], 
   "]"}]}], "\[IndentingNewLine]", 
 RowBox[{"u2", "=", 
  RowBox[{"Table", "[", 
   RowBox[{
    RowBox[{"If", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"u", "[", 
        RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\<C\>\""}], ",", 
      RowBox[{"(*", "0.25618", "*)"}], "0.2562", ",", 
      RowBox[{"If", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"u", "[", 
          RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\<H\>\""}], ",", 
        "0.6", ",", 
        RowBox[{"If", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"u", "[", 
            RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\<N\>\""}], ",", 
          "0.21522", ",", 
          RowBox[{"If", "[", 
           RowBox[{
            RowBox[{
             RowBox[{"u", "[", 
              RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\<O\>\""}], ",", 
            StyleBox["0.1792",
             FontColor->RGBColor[1, 0, 0]], ",", 
            RowBox[{"If", "[", 
             RowBox[{
              RowBox[{
               RowBox[{"u", "[", 
                RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\<Cl\>\""}], 
              ",", "0.48698"}], "]"}]}], "]"}]}], "]"}]}], "]"}]}], "]"}], 
    ",", 
    RowBox[{"{", 
     RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], 
   "]"}]}], "\[IndentingNewLine]", 
 RowBox[{"cn", "=", 
  RowBox[{"Table", "[", 
   RowBox[{
    RowBox[{"If", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"u", "[", 
        RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\<C\>\""}], ",", "6.0",
       ",", 
      RowBox[{"If", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"u", "[", 
          RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\<H\>\""}], ",", 
        "1.0", ",", 
        RowBox[{"If", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"u", "[", 
            RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\<N\>\""}], ",", 
          "7.0", ",", 
          RowBox[{"If", "[", 
           RowBox[{
            RowBox[{
             RowBox[{"u", "[", 
              RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\<O\>\""}], ",", 
            "8.0", ",", 
            RowBox[{"If", "[", 
             RowBox[{
              RowBox[{
               RowBox[{"u", "[", 
                RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\<Cl\>\""}], 
              ",", "7.0"}], "]"}]}], "]"}]}], "]"}]}], "]"}]}], "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], 
   "]"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"bc", "=", "0"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"For", "[", 
   RowBox[{
    RowBox[{"i", "=", "1"}], ",", " ", 
    RowBox[{"i", "<", 
     RowBox[{"n", "+", "1"}]}], ",", 
    RowBox[{"i", "++"}], ",", 
    RowBox[{"If", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"u", "[", 
        RowBox[{"[", "i", "]"}], "]"}], "\[NotEqual]", "\"\<H\>\""}], ",", 
      " ", 
      RowBox[{"bc", "=", 
       RowBox[{"bc", "+", "1"}]}]}], "]"}]}], "]"}], 
  ";"}], "\[IndentingNewLine]", "bc", "\[IndentingNewLine]", 
 RowBox[{"bh", "=", 
  RowBox[{"n", "-", "bc"}]}]}], "Input",
 CellChangeTimes->{{3.568277663590458*^9, 3.568277712652544*^9}, 
   3.5682785838808746`*^9, 3.580447125101241*^9, {3.580447338150815*^9, 
   3.580447473683853*^9}, {3.580447529251151*^9, 3.5804475516215897`*^9}, {
   3.580448220971966*^9, 3.580448222141968*^9}, {3.5804483597810097`*^9, 
   3.5804483678618236`*^9}, {3.5804484008402815`*^9, 
   3.5804484262371264`*^9}, {3.5804484577023816`*^9, 
   3.5804484578895817`*^9}, {3.5804485097752733`*^9, 3.5804485333157144`*^9}, 
   3.58044856532697*^9, {3.580449794250329*^9, 3.5804497945467296`*^9}, {
   3.580449861034046*^9, 3.580449943339791*^9}, {3.580450110182084*^9, 
   3.5804501195421004`*^9}, {3.5804595061521816`*^9, 3.580459540721842*^9}, 
   3.5804625847395887`*^9, {3.5804626281388645`*^9, 3.580462631758071*^9}, {
   3.5804626688549366`*^9, 3.5804627184318237`*^9}, {3.5804763932858744`*^9, 
   3.5804763978722825`*^9}, 3.5804765357297244`*^9, 3.5804765911254215`*^9, {
   3.580476629189489*^9, 3.5804766297510896`*^9}, {3.580476668641958*^9, 
   3.5804766699055605`*^9}, {3.580477336588331*^9, 3.5804773389439354`*^9}, {
   3.58047738137601*^9, 3.580477410969262*^9}, 3.580623917574172*^9, {
   3.5819203240109224`*^9, 3.581920324900124*^9}}],

Cell[BoxData["4"], "Output",
 CellChangeTimes->{
  3.5819326414304485`*^9, 3.581932749133038*^9, 3.5819327949191184`*^9, {
   3.5819329248205466`*^9, 3.581932952853796*^9}, 3.5819505007582507`*^9, 
   3.581997263276965*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"0.24`", ",", "0.12`", ",", "0.12`", ",", "0.12`"}], 
  "}"}]], "Output",
 CellChangeTimes->{
  3.5819326414304485`*^9, 3.581932749133038*^9, 3.5819327949191184`*^9, {
   3.5819329248205466`*^9, 3.581932952853796*^9}, 3.5819505007582507`*^9, 
   3.581997263276965*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"0.21522`", ",", "0.6`", ",", "0.6`", ",", "0.6`"}], 
  "}"}]], "Output",
 CellChangeTimes->{
  3.5819326414304485`*^9, 3.581932749133038*^9, 3.5819327949191184`*^9, {
   3.5819329248205466`*^9, 3.581932952853796*^9}, 3.5819505007582507`*^9, 
   3.581997263276965*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"7.`", ",", "1.`", ",", "1.`", ",", "1.`"}], "}"}]], "Output",
 CellChangeTimes->{
  3.5819326414304485`*^9, 3.581932749133038*^9, 3.5819327949191184`*^9, {
   3.5819329248205466`*^9, 3.581932952853796*^9}, 3.5819505007582507`*^9, 
   3.581997263276965*^9}],

Cell[BoxData["1"], "Output",
 CellChangeTimes->{
  3.5819326414304485`*^9, 3.581932749133038*^9, 3.5819327949191184`*^9, {
   3.5819329248205466`*^9, 3.581932952853796*^9}, 3.5819505007582507`*^9, 
   3.5819972632925653`*^9}],

Cell[BoxData["3"], "Output",
 CellChangeTimes->{
  3.5819326414304485`*^9, 3.581932749133038*^9, 3.5819327949191184`*^9, {
   3.5819329248205466`*^9, 3.581932952853796*^9}, 3.5819505007582507`*^9, 
   3.5819972632925653`*^9}]
}, Open  ]],

Cell[BoxData[{
 RowBox[{
  RowBox[{"Plot1", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"Graphics3D", "[", 
      RowBox[{"{", 
       RowBox[{"Sphere", "[", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{"t", "[", 
           RowBox[{"[", "i", "]"}], "]"}], "}"}], ",", 
         RowBox[{"u1", "[", 
          RowBox[{"[", "i", "]"}], "]"}]}], "]"}], "}"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Plot1o", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"Graphics3D", "[", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"Opacity", "[", "0.8", "]"}], ",", 
        RowBox[{"If", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"u", "[", 
            RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\<N\>\""}], ",", 
          RowBox[{"Lighter", "[", "Blue", "]"}], ",", 
          RowBox[{"If", "[", 
           RowBox[{
            RowBox[{
             RowBox[{"u", "[", 
              RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\<H\>\""}], ",", 
            RowBox[{"Lighter", "[", "White", "]"}], ",", 
            RowBox[{"If", "[", 
             RowBox[{
              RowBox[{
               RowBox[{"u", "[", 
                RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\<O\>\""}], 
              ",", 
              RowBox[{"Lighter", "[", 
               RowBox[{"Red", ",", "1"}], "]"}], ",", 
              RowBox[{"If", "[", 
               RowBox[{
                RowBox[{
                 RowBox[{"u", "[", 
                  RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\<Cl\>\""}], 
                ",", 
                RowBox[{"Lighter", "[", "Green", "]"}], ",", 
                RowBox[{"If", "[", 
                 RowBox[{
                  RowBox[{
                   RowBox[{"u", "[", 
                    RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\<B\>\""}],
                   ",", 
                  RowBox[{"Lighter", "[", "Orange", "]"}], ",", 
                  RowBox[{"Darker", "[", "Brown", "]"}]}], "]"}]}], "]"}]}], 
             "]"}]}], "]"}]}], "]"}], ",", 
        RowBox[{"Sphere", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{"t", "[", 
            RowBox[{"[", "i", "]"}], "]"}], "}"}], ",", 
          RowBox[{"u1", "[", 
           RowBox[{"[", "i", "]"}], "]"}]}], "]"}]}], "}"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]"}], "Input",
 CellChangeTimes->{{3.5679440895764875`*^9, 3.5679441000441055`*^9}, {
   3.5679441319773617`*^9, 3.5679441523041973`*^9}, {3.567944240210352*^9, 
   3.5679442813476243`*^9}, {3.567944323966899*^9, 3.567944345931738*^9}, {
   3.567944397567828*^9, 3.567944456941533*^9}, {3.567944491557994*^9, 
   3.5679445468444905`*^9}, {3.5679446505690727`*^9, 3.567944683875131*^9}, {
   3.567944736244423*^9, 3.5679448177233667`*^9}, {3.567945549707852*^9, 
   3.5679455499730525`*^9}, {3.5679456039959474`*^9, 3.567945606195551*^9}, {
   3.567945704460124*^9, 3.5679457164877453`*^9}, {3.5679457857362666`*^9, 
   3.5679458050023003`*^9}, {3.567946353670864*^9, 3.5679464127169676`*^9}, {
   3.5679466013368993`*^9, 3.5679466311173515`*^9}, {3.5679467590999765`*^9, 
   3.567946915739851*^9}, {3.5679469466123056`*^9, 3.567946946815106*^9}, {
   3.5679471385082426`*^9, 3.5679471497870626`*^9}, {3.5679483907224417`*^9, 
   3.5679484036080647`*^9}, {3.567948441937332*^9, 3.5679484421869326`*^9}, {
   3.567948496740228*^9, 3.56794849772303*^9}, {3.567948847912445*^9, 
   3.5679488949933276`*^9}, {3.568011728877214*^9, 3.5680117308213253`*^9}, {
   3.568011766325356*^9, 3.568011773436763*^9}, {3.5681173314452896`*^9, 
   3.5681173370612993`*^9}, {3.568117919800323*^9, 3.5681179678952074`*^9}, {
   3.56811809459863*^9, 3.568118105409449*^9}, {3.5681181818963833`*^9, 
   3.5681182389144835`*^9}, {3.568118274888147*^9, 3.5681182883977704`*^9}, {
   3.568118518513775*^9, 3.56811852739019*^9}, {3.5681188250543127`*^9, 
   3.568118862634779*^9}, {3.568118972396572*^9, 3.5681189862493963`*^9}, {
   3.568119020694257*^9, 3.568119047136303*^9}, {3.5681190903015785`*^9, 
   3.568119136930061*^9}, {3.5681191671005135`*^9, 3.568119202808976*^9}, {
   3.5681194788046613`*^9, 3.568119487603077*^9}, {3.5681195353235607`*^9, 
   3.5681195410331707`*^9}, 3.56811959201406*^9, {3.5681196426985493`*^9, 
   3.568119659702579*^9}, {3.568119726813897*^9, 3.56811974027672*^9}, {
   3.5682659184078283`*^9, 3.5682659267694435`*^9}, {3.568265972477524*^9, 
   3.568266002725977*^9}, 3.568277802649102*^9, {3.5804449976745625`*^9, 
   3.580445016550596*^9}, {3.5804471386732645`*^9, 3.580447144570075*^9}, {
   3.5804597731466503`*^9, 3.580459773926652*^9}, {3.580459890942457*^9, 
   3.5804599070572853`*^9}, {3.5804612382856236`*^9, 
   3.5804612398300266`*^9}, {3.580461274914488*^9, 3.580461275585289*^9}, {
   3.5804613451302114`*^9, 3.580461345504612*^9}, {3.580465262094291*^9, 
   3.5804652627026925`*^9}, {3.580465328893608*^9, 3.5804653540408525`*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"Plot2", "=", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"Graphics3D", "[", 
       RowBox[{"Text", "[", 
        RowBox[{"i", ",", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{
            RowBox[{"t", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "1"}], "]"}], "]"}], "-", "0.5"}], ",", 
           RowBox[{"t", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "2"}], "]"}], "]"}], ",", 
           RowBox[{"t", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "3"}], "]"}], "]"}]}], "}"}], " ", ",", 
         RowBox[{"{", 
          RowBox[{"1", ",", "0"}], "}"}]}], "]"}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}], 
  "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{"Distance", " ", "Matrix", " ", "Bohr"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"dc", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"Norm", "[", 
      RowBox[{
       RowBox[{"t", "[", 
        RowBox[{"[", "i", "]"}], "]"}], "-", 
       RowBox[{"t", "[", 
        RowBox[{"[", "j", "]"}], "]"}]}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "n"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"j", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.5679440895764875`*^9, 3.5679441000441055`*^9}, {
   3.5679441319773617`*^9, 3.5679441523041973`*^9}, {3.567944240210352*^9, 
   3.5679442813476243`*^9}, {3.567944323966899*^9, 3.567944345931738*^9}, {
   3.567944397567828*^9, 3.567944456941533*^9}, {3.567944491557994*^9, 
   3.5679445468444905`*^9}, {3.5679446505690727`*^9, 3.567944683875131*^9}, {
   3.567944736244423*^9, 3.5679448177233667`*^9}, {3.567945549707852*^9, 
   3.5679455499730525`*^9}, {3.5679456039959474`*^9, 3.567945606195551*^9}, {
   3.567945704460124*^9, 3.5679457164877453`*^9}, {3.5679457857362666`*^9, 
   3.5679458050023003`*^9}, {3.567946353670864*^9, 3.5679464127169676`*^9}, {
   3.5679466013368993`*^9, 3.5679466311173515`*^9}, {3.5679467590999765`*^9, 
   3.567946915739851*^9}, {3.5679469466123056`*^9, 3.567946946815106*^9}, {
   3.5679471385082426`*^9, 3.5679471497870626`*^9}, {3.5679483907224417`*^9, 
   3.5679484036080647`*^9}, {3.567948441937332*^9, 3.5679484421869326`*^9}, {
   3.567948496740228*^9, 3.56794849772303*^9}, {3.567948847912445*^9, 
   3.5679488949933276`*^9}, {3.568011728877214*^9, 3.5680117308213253`*^9}, {
   3.568011766325356*^9, 3.568011773436763*^9}, {3.5681173314452896`*^9, 
   3.5681173370612993`*^9}, {3.568117919800323*^9, 3.5681179678952074`*^9}, {
   3.56811809459863*^9, 3.568118105409449*^9}, {3.5681181818963833`*^9, 
   3.5681182389144835`*^9}, {3.568118274888147*^9, 3.5681182883977704`*^9}, {
   3.568118518513775*^9, 3.56811852739019*^9}, {3.5681188250543127`*^9, 
   3.568118862634779*^9}, {3.568118972396572*^9, 3.5681189862493963`*^9}, {
   3.568119020694257*^9, 3.568119047136303*^9}, {3.5681190903015785`*^9, 
   3.568119136930061*^9}, {3.5681191671005135`*^9, 3.568119202808976*^9}, {
   3.5681194788046613`*^9, 3.568119487603077*^9}, {3.5681195353235607`*^9, 
   3.5681195410331707`*^9}, 3.56811959201406*^9, {3.5681196426985493`*^9, 
   3.568119659702579*^9}, {3.568119726813897*^9, 3.56811974027672*^9}, {
   3.5681199008322024`*^9, 3.568119920691037*^9}, {3.5681199536070952`*^9, 
   3.568119958521104*^9}, {3.568120000641178*^9, 3.568120067206495*^9}, {
   3.5681211528280015`*^9, 3.5681211643408213`*^9}, 3.5681213201694956`*^9, {
   3.568121503313817*^9, 3.568121508961027*^9}, {3.568266015439999*^9, 
   3.5682660461096525`*^9}, 3.568277806081108*^9, {3.5804471990141706`*^9, 
   3.5804472170010023`*^9}, {3.5804476007616763`*^9, 3.5804476229917154`*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"x", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"t", "[", 
      RowBox[{"[", 
       RowBox[{"i", ",", "1"}], "]"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"y", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"t", "[", 
      RowBox[{"[", 
       RowBox[{"i", ",", "2"}], "]"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"z", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"t", "[", 
      RowBox[{"[", 
       RowBox[{"i", ",", "3"}], "]"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.5678630319943895`*^9, 3.5678630956269007`*^9}, {
  3.5678632280087337`*^9, 3.567863228336334*^9}, {3.5678632621727934`*^9, 
  3.5678633509213495`*^9}, {3.5678644293824434`*^9, 3.567864451908883*^9}, {
  3.567864923404111*^9, 3.56786494544695*^9}, {3.5678650143210707`*^9, 
  3.5678650845367947`*^9}, {3.567869231084096*^9, 3.567869278118179*^9}, {
  3.5682778102931156`*^9, 3.5682778156907253`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"(*", " ", 
    RowBox[{"Distance", " ", "matrix", " ", "Bohr"}], " ", "*)"}], "Text"], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"Dij", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"If", "[", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{
            RowBox[{"dc", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "j"}], "]"}], "]"}], "<", "3.0"}], "||", 
           RowBox[{"(", 
            RowBox[{
             RowBox[{
              RowBox[{"u", "[", 
               RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\<Cl\>\""}], 
             " ", "&&", " ", 
             RowBox[{
              RowBox[{"dc", "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", "j"}], "]"}], "]"}], "<", "3.24"}]}], 
            ")"}]}], ")"}], ",", 
         RowBox[{"dc", "[", 
          RowBox[{"[", 
           RowBox[{"i", ",", "j"}], "]"}], "]"}], ",", "0"}], "]"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "n"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"j", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{"MatrixForm", "[", 
    RowBox[{"Prepend", "[", 
     RowBox[{"Dij", ",", "num1"}], "]"}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.567863535157673*^9, 3.5678635412104836`*^9}, {
   3.5678638182825704`*^9, 3.567863827283786*^9}, {3.5679503910983553`*^9, 
   3.5679504266352177`*^9}, {3.5679504653700857`*^9, 3.5679505048069553`*^9}, 
   3.567950534915008*^9, {3.567951107623214*^9, 3.567951136046464*^9}, 
   3.567951302545556*^9, {3.5679562986474605`*^9, 3.567956305355472*^9}, {
   3.5679563590975666`*^9, 3.567956370688387*^9}, {3.567956414930065*^9, 
   3.567956430467692*^9}, {3.567956484740187*^9, 3.567956500667815*^9}, {
   3.5679565456270943`*^9, 3.5679565868267665`*^9}, {3.5679569267045636`*^9, 
   3.5679569448941956`*^9}, {3.567956976562251*^9, 3.567957000711094*^9}, {
   3.5679570348595533`*^9, 3.5679570632672033`*^9}, {3.5679573263928657`*^9, 
   3.5679573286548696`*^9}, {3.5679573611029263`*^9, 
   3.5679573637393312`*^9}, {3.5679592820903006`*^9, 3.567959317127962*^9}, {
   3.567959397655304*^9, 3.5679594004945087`*^9}, {3.5679594704762316`*^9, 
   3.5679594724574347`*^9}, {3.567959655898157*^9, 3.567959673073787*^9}, 
   3.5681143121223865`*^9, 3.568121866388855*^9, 3.568220668221574*^9, {
   3.568220718968463*^9, 3.568220723742071*^9}, {3.568220862051914*^9, 
   3.5682208641111183`*^9}, {3.568220906761593*^9, 3.5682209083527956`*^9}, {
   3.5682209451844606`*^9, 3.5682209651680956`*^9}, {3.568221147782016*^9, 
   3.5682211695908546`*^9}, {3.568221237279373*^9, 3.568221253612602*^9}, {
   3.568221305420293*^9, 3.568221306247094*^9}, {3.5682215161142635`*^9, 
   3.568221546206716*^9}, 3.56822158247678*^9, {3.568277819185131*^9, 
   3.5682778284359474`*^9}, {3.568286391407484*^9, 3.568286398053096*^9}, {
   3.5683070136254387`*^9, 3.5683070186174474`*^9}, 3.5804476741754055`*^9}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"1", "2", "3", "4"},
     {"0.`", "1.9219753314099903`", "1.9219752850062577`", 
      "1.9219752844103315`"},
     {"1.9219753314099903`", "0.`", "0", "0"},
     {"1.9219752850062577`", "0", "0.`", "0"},
     {"1.9219752844103315`", "0", "0", "0.`"}
    },
    GridBoxAlignment->{
     "Columns" -> {{Center}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
      "RowsIndexed" -> {}},
    GridBoxSpacings->{"Columns" -> {
        Offset[0.27999999999999997`], {
         Offset[0.7]}, 
        Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
        Offset[0.2], {
         Offset[0.4]}, 
        Offset[0.2]}, "RowsIndexed" -> {}}], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$]]]], "Output",
 CellChangeTimes->{
  3.581932641477249*^9, 3.581932749164238*^9, 3.5819327949503183`*^9, {
   3.5819329248517466`*^9, 3.581932952884996*^9}, 3.5819505007894506`*^9, 
   3.5819972633081656`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"For", "[", 
  RowBox[{
   RowBox[{"i", "=", "1"}], " ", ",", 
   RowBox[{"i", "<", 
    RowBox[{"n", "+", "1"}]}], ",", 
   RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
   RowBox[{"For", "[", 
    RowBox[{
     RowBox[{"j", "=", "1"}], ",", " ", 
     RowBox[{"j", "<", 
      RowBox[{"bc", "+", "1"}]}], ",", 
     RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", 
     RowBox[{"If", "[", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{
         RowBox[{
          RowBox[{"Dij", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "j"}], "]"}], "]"}], "\[NotEqual]", "0"}], "  ",
          "&&", " ", 
         RowBox[{
          RowBox[{"u", "[", 
           RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", " ", "\"\<H\>\""}]}], 
        ")"}], ",", 
       RowBox[{"If", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"u", "[", 
           RowBox[{"[", "j", "]"}], "]"}], "\[Equal]", "\"\<C\>\""}], ",", 
         RowBox[{
          RowBox[{"u2", "[", 
           RowBox[{"[", "i", "]"}], "]"}], "=", "0.59"}], ",", 
         RowBox[{"If", "[", 
          RowBox[{
           RowBox[{
            RowBox[{"u", "[", 
             RowBox[{"[", "j", "]"}], "]"}], "\[Equal]", "\"\<F\>\""}], ",", 
           RowBox[{
            RowBox[{"u2", "[", 
             RowBox[{"[", "i", "]"}], "]"}], "=", "0.5"}], ",", 
           RowBox[{"If", "[", 
            RowBox[{
             RowBox[{
              RowBox[{"u", "[", 
               RowBox[{"[", "j", "]"}], "]"}], "\[Equal]", "\"\<N\>\""}], ",", 
             RowBox[{
              RowBox[{"u2", "[", 
               RowBox[{"[", "i", "]"}], "]"}], "=", "0.5541"}], ",", 
             RowBox[{"If", "[", 
              RowBox[{
               RowBox[{
                RowBox[{"u", "[", 
                 RowBox[{"[", "j", "]"}], "]"}], "\[Equal]", "\"\<O\>\""}], 
               ",", 
               RowBox[{
                RowBox[{"u2", "[", 
                 RowBox[{"[", "i", "]"}], "]"}], "=", 
                StyleBox["0.52",
                 FontColor->RGBColor[1, 0, 0]]}]}], "]"}]}], "]"}]}], "]"}]}],
         "]"}]}], "]"}]}], "]"}]}], 
  "]"}], "\[IndentingNewLine]", "u2"}], "Input",
 CellChangeTimes->{{3.5679585418846006`*^9, 3.5679585738334565`*^9}, {
   3.567959075187137*^9, 3.5679590758579383`*^9}, {3.5681145227851562`*^9, 
   3.5681146476945763`*^9}, {3.568280038115429*^9, 3.5682801254131823`*^9}, {
   3.5682851622189255`*^9, 3.5682851860713673`*^9}, 3.5682852194866257`*^9, {
   3.568296148255143*^9, 3.5682961775675945`*^9}, {3.568305693114319*^9, 
   3.5683057123647532`*^9}, {3.5683062281640587`*^9, 3.568306232532067*^9}, {
   3.568306348658671*^9, 3.5683063634318967`*^9}, {3.5683064594344654`*^9, 
   3.5683064674060793`*^9}, {3.5683072458942466`*^9, 3.568307254973462*^9}, {
   3.56830729910594*^9, 3.56830731050956*^9}, {3.568307613368492*^9, 
   3.568307617174899*^9}, 3.580447736825115*^9, 3.5804627951995583`*^9, {
   3.5804628372416325`*^9, 3.5804628374600325`*^9}, 3.580476792802576*^9, 
   3.580622794013399*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"0.21522`", ",", "0.5541`", ",", "0.5541`", ",", "0.5541`"}], 
  "}"}]], "Output",
 CellChangeTimes->{
  3.581932641477249*^9, 3.581932749164238*^9, 3.5819327949659185`*^9, {
   3.5819329248517466`*^9, 3.5819329529005957`*^9}, 3.581950500805051*^9, 
   3.5819972633081656`*^9}]
}, Open  ]],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"Fij", "=", 
    RowBox[{"Dij", "-", "u2"}]}], ";"}], "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"MatrixForm", "[", "Fij", "]"}], ";"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Gij", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"If", "[", 
      RowBox[{
       RowBox[{"i", "\[GreaterEqual]", "j"}], ",", 
       RowBox[{"If", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"Fij", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "j"}], "]"}], "]"}], ">", "0"}], ",", 
         RowBox[{"Fij", "[", 
          RowBox[{"[", 
           RowBox[{"i", ",", "j"}], "]"}], "]"}], ",", "0"}], "]"}], ",", 
       "\"\< \>\""}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "n"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"j", ",", "1", ",", "n"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"MatrixForm", "[", "Gij", "]"}], ";"}], "*)"}]}]}], "Input",
 CellChangeTimes->{{3.567958596890297*^9, 3.567958603270708*^9}, {
   3.5679587604565845`*^9, 3.5679587707682023`*^9}, {3.5679588272871017`*^9, 
   3.567958828254304*^9}, {3.567959580066424*^9, 3.5679595814704266`*^9}, {
   3.567959624994503*^9, 3.5679596300177116`*^9}, {3.5679597079086485`*^9, 
   3.567959902877791*^9}, {3.5679599961659546`*^9, 3.5679600446196404`*^9}, {
   3.5679609611992497`*^9, 3.567960974006872*^9}, 3.567961975154231*^9, {
   3.5679622637079377`*^9, 3.5679622945179915`*^9}, 3.5681142413606625`*^9, 
   3.568121878697276*^9, {3.5681219332193723`*^9, 3.56812197189184*^9}, {
   3.5682778370003624`*^9, 3.5682778397927675`*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"Hij", "=", 
   RowBox[{"Transpose", "[", "Gij", "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"MatrixForm", "[", "Hij", "]"}], ";"}], "*)"}]}]}], "Input",
 CellChangeTimes->{{3.5679601758938704`*^9, 3.5679601965327067`*^9}, {
   3.567960427147912*^9, 3.567960432217921*^9}, {3.5679605764089737`*^9, 
   3.5679605968294096`*^9}, {3.5679607389456596`*^9, 3.5679607643581038`*^9}, 
   3.5681142501278777`*^9, 3.5681218877140923`*^9}],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"Kij", "=", 
    RowBox[{"Hij", "-", "u2"}]}], ";"}], "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"MatrixForm", "[", "Kij", "]"}], ";"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Lij", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"If", "[", 
      RowBox[{
       RowBox[{"i", "\[LessEqual]", "j"}], ",", 
       RowBox[{"If", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"Kij", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "j"}], "]"}], "]"}], ">", "0"}], ",", 
         RowBox[{"Kij", "[", 
          RowBox[{"[", 
           RowBox[{"i", ",", "j"}], "]"}], "]"}], ",", "0"}], "]"}], ",", 
       "\"\< \>\""}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "n"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"j", ",", "1", ",", "n"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"MatrixForm", "[", "Lij", "]"}], ";"}], "*)"}]}]}], "Input",
 CellChangeTimes->{{3.5679602339571724`*^9, 3.567960257419614*^9}, {
   3.567960343406965*^9, 3.567960353500182*^9}, {3.5679604155258913`*^9, 
   3.567960421422702*^9}, {3.567960854323462*^9, 3.567960874197897*^9}, {
   3.5679611596939983`*^9, 3.567961170270817*^9}, {3.5679614000124207`*^9, 
   3.567961407016833*^9}, {3.5679615066074076`*^9, 3.567961509649413*^9}, 
   3.5679616124379935`*^9, 3.5679616766477065`*^9, {3.5679617424486217`*^9, 
   3.5679617456310277`*^9}, 3.567961801276325*^9, {3.5679618582788253`*^9, 
   3.5679618970760937`*^9}, {3.5679620279135237`*^9, 3.567962030113127*^9}, 
   3.567962107894864*^9, {3.567962149468937*^9, 3.5679621572533503`*^9}, 
   3.5679623389156694`*^9, {3.567962645815009*^9, 3.567962648888214*^9}, {
   3.567962691647889*^9, 3.567962693020692*^9}, {3.568048161486014*^9, 
   3.568048162827616*^9}, 3.5680482096900983`*^9, 3.5681142588638926`*^9, 
   3.568121899632513*^9, 3.5681219958534822`*^9, {3.568277844675576*^9, 
   3.5682778463603787`*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"Lij", "=", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"If", "[", 
       RowBox[{
        RowBox[{"i", "\[LessEqual]", "j"}], ",", 
        RowBox[{"If", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"Lij", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "j"}], "]"}], "]"}], ">", "0"}], ",", 
          RowBox[{"If", "[", 
           RowBox[{
            RowBox[{
             RowBox[{"u", "[", 
              RowBox[{"[", "j", "]"}], "]"}], "\[NotEqual]", "\"\<H\>\""}], 
            ",", 
            RowBox[{"0.5", "*", 
             RowBox[{"Lij", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "j"}], "]"}], "]"}]}], ",", 
            RowBox[{"Lij", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "j"}], "]"}], "]"}]}], "]"}], ",", "0"}], 
         "]"}], ",", "\"\< \>\""}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "n"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"j", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}], 
  "  "}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Qij", "=", "Lij"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"MatrixForm", "[", "Lij", "]"}], ";"}], "*)"}]}]}], "Input",
 CellChangeTimes->{{3.5679627617388124`*^9, 3.567962837461345*^9}, {
   3.567963118121438*^9, 3.567963118698639*^9}, {3.5679631520670977`*^9, 
   3.5679631541731014`*^9}, {3.5679631898191643`*^9, 
   3.5679632154656096`*^9}, {3.5679635823158536`*^9, 3.56796361468591*^9}, {
   3.567963691516045*^9, 3.567963693590849*^9}, {3.567963737957327*^9, 
   3.567963807299449*^9}, {3.568048114904332*^9, 3.5680481161835346`*^9}, 
   3.5680482221389203`*^9, {3.5680485003498087`*^9, 3.5680485123618298`*^9}, 
   3.5681142659931054`*^9, 3.5681219084621286`*^9, {3.568277848107582*^9, 
   3.5682778496519847`*^9}, {3.568285867917365*^9, 3.568285868634966*^9}, {
   3.568297947687104*^9, 3.568297949137906*^9}, 3.5683030673497076`*^9, {
   3.5683057620820403`*^9, 3.5683057631428423`*^9}, 3.5804497010089655`*^9}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"Nij", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"If", "[", 
      RowBox[{
       RowBox[{"i", "\[LessEqual]", "j"}], ",", 
       RowBox[{"If", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"Lij", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "j"}], "]"}], "]"}], ">", "0"}], ",", 
         RowBox[{"2.25", "/", 
          RowBox[{
           RowBox[{"Lij", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "j"}], "]"}], "]"}], "^", "2"}]}], ",", "0"}], 
        "]"}], ",", "\"\< \>\""}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "n"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"j", ",", "1", ",", "n"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"MatrixForm", "[", "Nij", "]"}], ";"}], "*)"}]}]}], "Input",
 CellChangeTimes->{{3.5679628944794455`*^9, 3.5679629710755796`*^9}, {
   3.567963889963994*^9, 3.567963933800071*^9}, {3.5680478391738477`*^9, 
   3.5680478414358516`*^9}, 3.568048231171336*^9, {3.5680482823550262`*^9, 
   3.56804828449223*^9}, 3.5680483349271183`*^9, 3.568114273059918*^9, 
   3.5681219175101447`*^9, {3.568277851305588*^9, 3.568277857295998*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{"Total", " ", "kinetic", " ", "energy", " ", "in", " ", "Eh"}], 
    ",", " ", 
    RowBox[{"=", " ", 
     RowBox[{"|", "Etot", "|"}]}]}], " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"Tval", "=", 
    RowBox[{"Sum", "[", 
     RowBox[{
      RowBox[{"If", "[", 
       RowBox[{
        RowBox[{"i", "\[LessEqual]", "j"}], ",", 
        RowBox[{"Nij", "[", 
         RowBox[{"[", 
          RowBox[{"i", ",", "j"}], "]"}], "]"}], ",", "0"}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "n"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"j", ",", "1", ",", "n"}], "}"}]}], "]"}]}], "\n", 
   RowBox[{"Tcor", "=", 
    RowBox[{"Sum", "[", 
     RowBox[{
      RowBox[{"If", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"u", "[", 
          RowBox[{"[", "i", "]"}], "]"}], "\[NotEqual]", "\"\<H\>\""}], ",", 
        RowBox[{"2.25", "/", 
         RowBox[{
          RowBox[{"u2", "[", 
           RowBox[{"[", "i", "]"}], "]"}], "^", "2"}]}], ",", "0"}], "]"}], 
      ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], "\n", 
   RowBox[{"Ttot", "=", 
    RowBox[{"Tval", "+", "Tcor"}]}], "  ", 
   RowBox[{"(*", " ", 
    RowBox[{"plus", " ", "\[Pi]", " ", "and", " ", "lone", " ", 
     RowBox[{"pairs", ":", " ", 
      RowBox[{"+", " ", "8.7616"}]}]}], " ", "*)"}], "\n", 
   RowBox[{"TvalH", "=", 
    RowBox[{"Sum", "[", 
     RowBox[{
      RowBox[{"If", "[", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"i", "\[LessEqual]", "j"}], " ", "&&", " ", 
          RowBox[{
           RowBox[{"u", "[", 
            RowBox[{"[", "j", "]"}], "]"}], "\[Equal]", " ", "\"\<H\>\""}]}], 
         ")"}], ",", 
        RowBox[{"Nij", "[", 
         RowBox[{"[", 
          RowBox[{"i", ",", "j"}], "]"}], "]"}], ",", "0"}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "n"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"j", ",", "1", ",", "n"}], "}"}]}], "]"}]}], "\n", 
   RowBox[{"TvalC", "=", 
    RowBox[{"Tval", "-", "TvalH"}]}]}]}]], "Input",
 CellChangeTimes->{{3.567963992456174*^9, 3.567964003890994*^9}, {
   3.5679640974287586`*^9, 3.5679642038209453`*^9}, {3.56796433180357*^9, 
   3.567964353487608*^9}, {3.567964403641696*^9, 3.567964453046983*^9}, {
   3.5679646160984693`*^9, 3.5679646168316708`*^9}, {3.567964688747797*^9, 
   3.567964762255126*^9}, {3.568277859480002*^9, 3.568277863598409*^9}, {
   3.568286113290196*^9, 3.5682861137113967`*^9}, {3.5682870401502237`*^9, 
   3.5682870902263117`*^9}, {3.568297787256422*^9, 3.5682977875216227`*^9}, 
   3.568305788461687*^9, 3.580447849972114*^9}],

Cell[BoxData["5.080481444811868`"], "Output",
 CellChangeTimes->{
  3.581932641508449*^9, 3.581932749195438*^9, 3.5819327949971185`*^9, {
   3.581932924882947*^9, 3.581932952916196*^9}, 3.581950500820651*^9, 
   3.581997263339365*^9}],

Cell[BoxData["48.57549825442976`"], "Output",
 CellChangeTimes->{
  3.581932641508449*^9, 3.581932749195438*^9, 3.5819327949971185`*^9, {
   3.581932924882947*^9, 3.581932952916196*^9}, 3.581950500820651*^9, 
   3.581997263339365*^9}],

Cell[BoxData["53.65597969924163`"], "Output",
 CellChangeTimes->{
  3.581932641508449*^9, 3.581932749195438*^9, 3.5819327949971185`*^9, {
   3.581932924882947*^9, 3.581932952916196*^9}, 3.581950500820651*^9, 
   3.581997263339365*^9}],

Cell[BoxData["5.080481444811868`"], "Output",
 CellChangeTimes->{
  3.581932641508449*^9, 3.581932749195438*^9, 3.5819327949971185`*^9, {
   3.581932924882947*^9, 3.581932952916196*^9}, 3.581950500820651*^9, 
   3.581997263339365*^9}],

Cell[BoxData["0.`"], "Output",
 CellChangeTimes->{
  3.581932641508449*^9, 3.581932749195438*^9, 3.5819327949971185`*^9, {
   3.581932924882947*^9, 3.581932952916196*^9}, 3.581950500820651*^9, 
   3.581997263339365*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"Mij", "=", 
     RowBox[{"Prepend", "[", 
      RowBox[{"Qij", ",", "u"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Oij", "=", 
     RowBox[{"Transpose", "[", "Mij", "]"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Pij", "=", 
     RowBox[{"Prepend", "[", 
      RowBox[{"Oij", ",", 
       RowBox[{"Prepend", "[", 
        RowBox[{"u", ",", "\"\< \>\""}], "]"}]}], "]"}]}], ";"}]}]}]], "Input",
 CellChangeTimes->{{3.568045828002715*^9, 3.56804588761042*^9}, 
   3.5680459828797874`*^9, {3.568046033938677*^9, 3.5680460518787084`*^9}, {
   3.568046989190755*^9, 3.5680470284404235`*^9}, {3.568047061028881*^9, 
   3.568047085567724*^9}, {3.5680471562982483`*^9, 3.568047169199471*^9}, {
   3.568047223393966*^9, 3.568047364059413*^9}, {3.568047417286707*^9, 
   3.5680474225127163`*^9}, {3.5680474583459787`*^9, 3.568047487861231*^9}, {
   3.5680475345521126`*^9, 3.5680475468605347`*^9}, {3.5680479168775845`*^9, 
   3.5680479311516094`*^9}, {3.568048066918648*^9, 3.568048067667449*^9}, {
   3.5680483816336*^9, 3.5680483930216208`*^9}, {3.568048424159275*^9, 
   3.568048450211321*^9}, {3.5680486326224413`*^9, 3.5680486370996494`*^9}, {
   3.5680487358166227`*^9, 3.5680487967191296`*^9}, 3.568048834377596*^9, 
   3.568048984418659*^9, {3.5680490255559316`*^9, 3.5680490790796256`*^9}, {
   3.5680491274553103`*^9, 3.5680491280949116`*^9}, {3.5680491782333994`*^9, 
   3.568049231944294*^9}, {3.568049361533722*^9, 3.5680493691933355`*^9}, {
   3.568049422171028*^9, 3.568049427490638*^9}, {3.568049476427923*^9, 
   3.5680495489368505`*^9}, {3.568049620634577*^9, 3.5680496274361887`*^9}, {
   3.5680496655314555`*^9, 3.568049693876705*^9}, 3.568049746573598*^9, {
   3.568049846460573*^9, 3.568049932759925*^9}, {3.568050009059659*^9, 
   3.5680500167348723`*^9}, {3.5680501572755194`*^9, 3.568050179911159*^9}, 
   3.5681087021497335`*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"Eij", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"If", "[", 
      RowBox[{
       RowBox[{"i", "\[GreaterEqual]", "j"}], ",", 
       RowBox[{"Dij", "[", 
        RowBox[{"[", 
         RowBox[{"i", ",", "j"}], "]"}], "]"}], ",", "0"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", 
       RowBox[{"bc", "-", "1"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"j", ",", "1", ",", 
       RowBox[{"bc", "-", "1"}]}], "}"}]}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.568045828002715*^9, 3.56804588761042*^9}, 
   3.5680459828797874`*^9, {3.568046033938677*^9, 3.5680460518787084`*^9}, {
   3.568046989190755*^9, 3.5680470284404235`*^9}, {3.568047061028881*^9, 
   3.568047085567724*^9}, {3.5680471562982483`*^9, 3.568047169199471*^9}, {
   3.568047223393966*^9, 3.568047364059413*^9}, {3.568047417286707*^9, 
   3.5680474225127163`*^9}, {3.5680474583459787`*^9, 3.568047487861231*^9}, {
   3.5680475345521126`*^9, 3.5680475468605347`*^9}, {3.5680479168775845`*^9, 
   3.5680479311516094`*^9}, {3.568048066918648*^9, 3.568048067667449*^9}, {
   3.5680483816336*^9, 3.5680483930216208`*^9}, {3.568048424159275*^9, 
   3.568048450211321*^9}, {3.5680486326224413`*^9, 3.5680486370996494`*^9}, {
   3.5680487358166227`*^9, 3.5680487967191296`*^9}, 3.568048834377596*^9, 
   3.568048984418659*^9, {3.5680490255559316`*^9, 3.5680490790796256`*^9}, {
   3.5680491274553103`*^9, 3.5680491280949116`*^9}, {3.5680491782333994`*^9, 
   3.568049231944294*^9}, {3.568049361533722*^9, 3.5680493691933355`*^9}, {
   3.568049422171028*^9, 3.568049427490638*^9}, {3.568049476427923*^9, 
   3.5680495489368505`*^9}, {3.568049620634577*^9, 3.5680496274361887`*^9}, {
   3.5680496655314555`*^9, 3.568049693876705*^9}, 3.568049746573598*^9, {
   3.568049846460573*^9, 3.568049932759925*^9}, {3.568050009059659*^9, 
   3.5680500167348723`*^9}, {3.5680501572755194`*^9, 3.568050179911159*^9}, {
   3.5681087021497335`*^9, 3.5681087223049684`*^9}, 3.568114284401138*^9, 
   3.5681220092695055`*^9, {3.5804479260534477`*^9, 3.580447971730328*^9}, {
   3.580449452781329*^9, 3.5804494665249534`*^9}, 3.5804629573930435`*^9, 
   3.5804651851393557`*^9}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"nc", "=", "0"}], ";"}], "  "}], "\n", 
 RowBox[{"For", "[", 
  RowBox[{
   RowBox[{"i", "=", "1"}], ",", 
   RowBox[{"i", "<", 
    RowBox[{"bc", "+", "1"}]}], ",", 
   RowBox[{"i", "++"}], ",", 
   RowBox[{"If", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"u", "[", 
       RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\<C\>\""}], ",", 
     RowBox[{"nc", "=", 
      RowBox[{"nc", "+", "1"}]}], ",", "nc"}], "]"}]}], 
  "]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"cnt", "=", 
   RowBox[{"Table", "[", 
    RowBox[{"0", ",", 
     RowBox[{"{", "nc", "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", "nc", "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"no", "=", "0"}], ";"}], "  "}], "\n", 
 RowBox[{"For", "[", 
  RowBox[{
   RowBox[{"i", "=", "1"}], ",", 
   RowBox[{"i", "<", 
    RowBox[{"bc", "+", "1"}]}], ",", 
   RowBox[{"i", "++"}], ",", 
   RowBox[{"If", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"u", "[", 
       RowBox[{"[", "i", "]"}], "]"}], "\[Equal]", "\"\<O\>\""}], ",", 
     RowBox[{"no", "=", 
      RowBox[{"no", "+", "1"}]}], ",", "no"}], "]"}]}], 
  "]"}], "\[IndentingNewLine]", "no"}], "Input",
 CellChangeTimes->{{3.5804503559293156`*^9, 3.5804503572241178`*^9}, {
   3.580450419593027*^9, 3.5804504349902544`*^9}, 3.5804645273330007`*^9, 
   3.580465176465741*^9}],

Cell[BoxData["0"], "Output",
 CellChangeTimes->{
  3.581932641524049*^9, 3.581932749226638*^9, 3.5819327950127187`*^9, {
   3.5819329248985467`*^9, 3.5819329529473963`*^9}, 3.5819505008518505`*^9, 
   3.581997263354965*^9}],

Cell[BoxData["0"], "Output",
 CellChangeTimes->{
  3.581932641524049*^9, 3.581932749226638*^9, 3.5819327950127187`*^9, {
   3.5819329248985467`*^9, 3.5819329529473963`*^9}, 3.5819505008518505`*^9, 
   3.581997263354965*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"For", "[", 
   RowBox[{
    RowBox[{"j", "=", "1"}], ",", 
    RowBox[{"j", "\[LessEqual]", "nc"}], ",", 
    RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", 
    RowBox[{"For", "[", 
     RowBox[{
      RowBox[{"i", "=", 
       RowBox[{"bc", "+", "1"}]}], ",", 
      RowBox[{"i", "\[LessEqual]", "n"}], ",", 
      RowBox[{"i", "++"}], ",", 
      RowBox[{"If", "[", " ", 
       RowBox[{
        RowBox[{
         RowBox[{"Dij", "[", 
          RowBox[{"[", 
           RowBox[{"i", ",", "j"}], "]"}], "]"}], "\[NotEqual]", "0"}], " ", 
        ",", 
        RowBox[{
         RowBox[{"cnt", "[", 
          RowBox[{"[", "j", "]"}], "]"}], "=", 
         RowBox[{
          RowBox[{"cnt", "[", 
           RowBox[{"[", "j", "]"}], "]"}], "+", "1"}]}], ",", "0"}], "]"}]}], 
     "]"}]}], "]"}], ";"}], "\[IndentingNewLine]", "cnt"}], "Input"],

Cell[BoxData[
 RowBox[{"{", "}"}]], "Output",
 CellChangeTimes->{
  3.581932641539649*^9, 3.581932749242238*^9, 3.5819327950283184`*^9, {
   3.581932924914147*^9, 3.581932952962996*^9}, 3.5819505008674507`*^9, 
   3.5819972633705654`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"con", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"{", "i", "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";", " ", 
  RowBox[{"conj", "=", 
   RowBox[{"Table", "[", 
    RowBox[{"0", ",", 
     RowBox[{"{", "nc", "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"For", "[", 
   RowBox[{
    RowBox[{"i", "=", "1"}], ",", 
    RowBox[{"i", "\[LessEqual]", "n"}], ",", 
    RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
    RowBox[{"For", "[", 
     RowBox[{
      RowBox[{"j", "=", "1"}], ",", 
      RowBox[{"j", "\[LessEqual]", "n"}], ",", 
      RowBox[{"j", "++"}], ",", 
      RowBox[{"If", "[", " ", 
       RowBox[{
        RowBox[{
         RowBox[{"Dij", "[", 
          RowBox[{"[", 
           RowBox[{"i", ",", "j"}], "]"}], "]"}], "\[NotEqual]", "0"}], " ", 
        ",", 
        RowBox[{
         RowBox[{"con", "[", 
          RowBox[{"[", "i", "]"}], "]"}], "=", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"con", "[", 
            RowBox[{"[", "i", "]"}], "]"}], ",", "j"}], "}"}]}]}], "]"}]}], 
     "]"}]}], "]"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"conn", "=", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"Flatten", "[", 
       RowBox[{"con", "[", 
        RowBox[{"[", "i", "]"}], "]"}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "n"}], "}"}]}], "]"}]}], ";"}], 
  "\[IndentingNewLine]", 
  RowBox[{"(*", " ", 
   RowBox[{
   "identify", " ", "C", " ", "atoms", " ", "at", " ", "double", " ", 
    "bonds"}], " ", "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"j", "=", "0"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"For", "[", 
   RowBox[{
    RowBox[{"i", "=", "1"}], ",", 
    RowBox[{"i", "\[LessEqual]", "nc"}], ",", 
    RowBox[{"i", "++"}], ",", 
    RowBox[{"If", "[", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{
         RowBox[{"Length", "[", 
          RowBox[{"conn", "[", 
           RowBox[{"[", "i", "]"}], "]"}], "]"}], "\[Equal]", "4"}], " ", "&&",
         " ", 
        RowBox[{
         RowBox[{"u", "[", 
          RowBox[{"[", 
           RowBox[{"conn", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "4"}], "]"}], "]"}], "]"}], "]"}], 
         "\[NotEqual]", "\"\<O\>\""}]}], ")"}], ",", 
      RowBox[{
       RowBox[{"j", "=", 
        RowBox[{"j", "+", "1"}]}], ";", 
       RowBox[{
        RowBox[{"conj", "[", 
         RowBox[{"[", "j", "]"}], "]"}], "=", "i"}]}]}], "]"}]}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"conj", "=", 
  RowBox[{"Take", "[", 
   RowBox[{"conj", ",", 
    RowBox[{"{", 
     RowBox[{"1", ",", "j"}], "}"}]}], "]"}]}], "\[IndentingNewLine]", 
 RowBox[{"nj", "=", 
  RowBox[{"Length", "[", "conj", "]"}]}], "\[IndentingNewLine]", 
 RowBox[{"MatrixForm", "[", 
  RowBox[{"conn", ",", " ", 
   RowBox[{"TableAlignments", "\[Rule]", "Left"}]}], "]"}]}], "Input",
 CellChangeTimes->{{3.5804504793099327`*^9, 3.580450480916735*^9}}],

Cell[BoxData[
 RowBox[{"{", "}"}]], "Output",
 CellChangeTimes->{
  3.581932641555249*^9, 3.5819327492578382`*^9, 3.5819327950439186`*^9, {
   3.5819329249297466`*^9, 3.5819329529785957`*^9}, 3.581950500883051*^9, 
   3.5819972633705654`*^9}],

Cell[BoxData["0"], "Output",
 CellChangeTimes->{
  3.581932641555249*^9, 3.5819327492578382`*^9, 3.5819327950439186`*^9, {
   3.5819329249297466`*^9, 3.5819329529785957`*^9}, 3.581950500883051*^9, 
   3.581997263386165*^9}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", 
   TagBox[GridBox[{
      {
       RowBox[{"{", 
        RowBox[{"1", ",", "2", ",", "3", ",", "4"}], "}"}]},
      {
       RowBox[{"{", 
        RowBox[{"2", ",", "1"}], "}"}]},
      {
       RowBox[{"{", 
        RowBox[{"3", ",", "1"}], "}"}]},
      {
       RowBox[{"{", 
        RowBox[{"4", ",", "1"}], "}"}]}
     },
     GridBoxAlignment->{
      "Columns" -> {{Left}}, "ColumnsIndexed" -> {}, "Rows" -> {{Baseline}}, 
       "RowsIndexed" -> {}},
     GridBoxSpacings->{"Columns" -> {
         Offset[0.27999999999999997`], {
          Offset[0.5599999999999999]}, 
         Offset[0.27999999999999997`]}, "ColumnsIndexed" -> {}, "Rows" -> {
         Offset[0.2], {
          Offset[0.4]}, 
         Offset[0.2]}, "RowsIndexed" -> {}}],
    Column], "\[NoBreak]", ")"}],
  Function[BoxForm`e$, 
   MatrixForm[BoxForm`e$, TableAlignments -> Left]]]], "Output",
 CellChangeTimes->{
  3.581932641555249*^9, 3.5819327492578382`*^9, 3.5819327950439186`*^9, {
   3.5819329249297466`*^9, 3.5819329529785957`*^9}, 3.581950500883051*^9, 
   3.581997263386165*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"j", "=", "0"}], ";", " ", 
  RowBox[{"h", "=", "0"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"coco", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"0", ",", "0", ",", "0", ",", "0"}], "}"}], ",", 
     RowBox[{"{", "no", "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"coch", "=", "coco"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"For", "[", 
   RowBox[{
    RowBox[{"i", "=", 
     RowBox[{"nc", "+", "1"}]}], ",", " ", 
    RowBox[{"i", "\[LessEqual]", 
     RowBox[{"nc", "+", "no"}]}], ",", " ", 
    RowBox[{"i", "++"}], ",", 
    RowBox[{"If", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"Length", "[", 
        RowBox[{"conn", "[", 
         RowBox[{"[", "i", "]"}], "]"}], "]"}], "\[Equal]", "2"}], ",", 
      RowBox[{
       RowBox[{"j", "=", 
        RowBox[{"j", "+", "1"}]}], ";", 
       RowBox[{
        RowBox[{"coco", "[", 
         RowBox[{"[", "j", "]"}], "]"}], "=", 
        RowBox[{"Join", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{"conn", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "1"}], "]"}], "]"}], "}"}], ",", 
          RowBox[{"Take", "[", 
           RowBox[{
            RowBox[{"conn", "[", 
             RowBox[{"[", 
              RowBox[{"conn", "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}], ",", 
            "3"}], "]"}]}], "]"}]}]}], ",", 
      RowBox[{
       RowBox[{"h", "=", 
        RowBox[{"h", "+", "1"}]}], ";", 
       RowBox[{
        RowBox[{"coch", "[", 
         RowBox[{"[", "h", "]"}], "]"}], "=", 
        RowBox[{"conn", "[", 
         RowBox[{"[", "i", "]"}], "]"}]}]}]}], "]"}]}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"coco", "=", 
   RowBox[{"Take", "[", 
    RowBox[{"coco", ",", 
     RowBox[{"{", 
      RowBox[{"1", ",", "j"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"coch", "=", 
   RowBox[{"Take", "[", 
    RowBox[{"coch", ",", 
     RowBox[{"{", 
      RowBox[{"1", ",", "h"}], "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"oj", "=", 
  RowBox[{"Length", "[", "coco", "]"}]}], "\[IndentingNewLine]", 
 RowBox[{"hj", "=", 
  RowBox[{"Length", "[", "coch", "]"}]}]}], "Input",
 CellChangeTimes->{{3.5804505432700443`*^9, 3.5804506082285585`*^9}, {
   3.5804506582734466`*^9, 3.5804506742478747`*^9}, 3.5804507096443367`*^9, {
   3.5804508095313125`*^9, 3.5804508306849494`*^9}, 3.5804510247492905`*^9, 
   3.5804511224366617`*^9, {3.580460502510331*^9, 3.580460514147952*^9}, {
   3.58046055297642*^9, 3.5804605556908245`*^9}, 3.580460689991461*^9, {
   3.5804610478404894`*^9, 3.5804610672625237`*^9}, {3.5804611544354763`*^9, 
   3.5804611697547035`*^9}, {3.580464611183148*^9, 3.5804646160191565`*^9}}],

Cell[BoxData["0"], "Output",
 CellChangeTimes->{
  3.5819326415708494`*^9, 3.5819327492734385`*^9, 3.581932795059519*^9, {
   3.581932924945347*^9, 3.581932952994196*^9}, 3.581950500898651*^9, 
   3.581997263401765*^9}],

Cell[BoxData["0"], "Output",
 CellChangeTimes->{
  3.5819326415708494`*^9, 3.5819327492734385`*^9, 3.581932795059519*^9, {
   3.581932924945347*^9, 3.581932952994196*^9}, 3.581950500898651*^9, 
   3.581997263401765*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"(*", 
    RowBox[{"--", 
     RowBox[{"--", 
      RowBox[{"--", 
       RowBox[{"--", 
        RowBox[{"--", 
         RowBox[{"--", 
          RowBox[{"--", 
           RowBox[{"--", 
            RowBox[{"--", 
             RowBox[{"--", 
              RowBox[{"--", 
               RowBox[{"--", 
                RowBox[{"--", 
                 RowBox[{"--", 
                  RowBox[{"--", 
                   RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{"--", 
                    RowBox[{
                    "--", "-"}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}]}\
]}]}]}]}]}]}]}]}]}]}]}], "*)"}], "Text"], "\[IndentingNewLine]", 
  StyleBox[
   RowBox[{"(*", 
    RowBox[{"next", ":", " ", 
     RowBox[{"N", " ", "LP"}]}], "*)"}], "Text"], "\n", 
  RowBox[{
   StyleBox[
    RowBox[{
     RowBox[{"k", " ", "=", " ", "1"}], ";"}], "Input"], "\n", 
   StyleBox[
    RowBox[{
     RowBox[{"Clear", "[", 
      RowBox[{"conn", ",", "P"}], "]"}], ";"}], "Input"], "\n", 
   RowBox[{
    StyleBox[
     RowBox[{"P", " ", "=", " ", "4"}], "Input"], 
    StyleBox[";", "Input"], 
    RowBox[{"rln", "=", "1.154"}], ";", " ", 
    RowBox[{"pn", "=", "1"}], ";"}], "\n", 
   RowBox[{
    RowBox[{
     StyleBox["tpyrn", "Input"], 
     StyleBox[" ", "Input"], 
     StyleBox["=", "Input"], 
     StyleBox[" ", "Input"], 
     RowBox[{
      StyleBox["{", "Input"], 
      StyleBox[
       RowBox[{"{", 
        RowBox[{"0.0", ",", " ", "0.0", ",", " ", "0.0"}], "}"}], "Input"], 
      "}"}]}], 
    StyleBox[";", "Input"]}], "\n", 
   StyleBox[
    RowBox[{
     RowBox[{"conn", " ", "=", " ", 
      RowBox[{"{", 
       RowBox[{"{", 
        RowBox[{"2", ",", "3", ",", "4", ",", "1"}], "}"}], "}"}]}], ";"}], 
    "Input"], 
   StyleBox["\[IndentingNewLine]", "Input"], 
   StyleBox[
    RowBox[{
     RowBox[{"f", "=", 
      RowBox[{"{", 
       RowBox[{"-", "1"}], "}"}]}], ";"}], "Input"]}]}]], "Input",
 CellChangeTimes->{{3.5804505432700443`*^9, 3.5804506082285585`*^9}, {
   3.5804506582734466`*^9, 3.5804506742478747`*^9}, 3.5804507096443367`*^9, {
   3.5804508095313125`*^9, 3.5804508306849494`*^9}, 3.5804510247492905`*^9, 
   3.5804511224366617`*^9, {3.580460502510331*^9, 3.580460514147952*^9}, {
   3.58046055297642*^9, 3.5804605556908245`*^9}, 3.580460689991461*^9, {
   3.5804610478404894`*^9, 3.5804610672625237`*^9}, {3.5804611544354763`*^9, 
   3.5804611697547035`*^9}, 3.580464611183148*^9, {3.5804885900917444`*^9, 
   3.5804886187021947`*^9}, {3.5804887019907413`*^9, 3.580488708293152*^9}, {
   3.580489378080329*^9, 3.5804893866759434`*^9}, {3.580525866734337*^9, 
   3.5805259166232247`*^9}, {3.5806227077140474`*^9, 
   3.5806227086032486`*^9}, {3.581918474448312*^9, 3.5819184886287365`*^9}, {
   3.5819188422033577`*^9, 3.581918845245363*^9}, {3.5819313655042076`*^9, 
   3.5819313663622093`*^9}, {3.581931789294552*^9, 3.581931790339754*^9}, {
   3.5819322877154274`*^9, 3.581932364795163*^9}, 3.5819324025160294`*^9, {
   3.5819972346041145`*^9, 3.581997250188542*^9}}],

Cell[BoxData[
 StyleBox[
  RowBox[{"<<", " ", "\"\<C:\\\\Kimb_Mathem\\\\LpyrNtrans.m\>\""}], 
  "Code"]], "Input",
 CellChangeTimes->{{3.5804505432700443`*^9, 3.5804506082285585`*^9}, {
   3.5804506582734466`*^9, 3.5804506742478747`*^9}, 3.5804507096443367`*^9, {
   3.5804508095313125`*^9, 3.5804508306849494`*^9}, 3.5804510247492905`*^9, 
   3.5804511224366617`*^9, {3.580460502510331*^9, 3.580460514147952*^9}, {
   3.58046055297642*^9, 3.5804605556908245`*^9}, 3.580460689991461*^9, {
   3.5804610478404894`*^9, 3.5804610672625237`*^9}, {3.5804611544354763`*^9, 
   3.5804611697547035`*^9}, 3.580464611183148*^9, {3.5804885900917444`*^9, 
   3.5804886187021947`*^9}, {3.5804887019907413`*^9, 3.580488708293152*^9}, {
   3.580489378080329*^9, 3.5804893866759434`*^9}, {3.580525866734337*^9, 
   3.5805259166232247`*^9}, {3.5806227077140474`*^9, 
   3.5806227086032486`*^9}, {3.581918474448312*^9, 3.5819184886287365`*^9}, {
   3.5819188422033577`*^9, 3.581918845245363*^9}, {3.5819313655042076`*^9, 
   3.5819313663622093`*^9}, {3.581931789294552*^9, 3.581931790339754*^9}, {
   3.5819322877154274`*^9, 3.581932364795163*^9}, 3.5819324025160294`*^9}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"DeleteFile", "[", " ", "\"\<C:\\\\Kimb_Mathem\\\\tmp4\>\"", "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"For", "[", 
  RowBox[{
   RowBox[{"i", "=", "2"}], ",", 
   RowBox[{"i", "<", 
    RowBox[{"n", "+", "2"}]}], ",", 
   RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
   RowBox[{"For", "[", 
    RowBox[{
     RowBox[{"j", "=", "2"}], ",", 
     RowBox[{"j", "\[LessEqual]", "i"}], ",", 
     RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", 
     RowBox[{"If", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"Pij", "[", 
         RowBox[{"[", 
          RowBox[{"i", ",", "j"}], "]"}], "]"}], ">", "0"}], ",", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{
          RowBox[{"i", "-", "1"}], ",", 
          RowBox[{"j", "-", "1"}], ",", 
          RowBox[{"Extract", "[", 
           RowBox[{"Pij", ",", 
            RowBox[{"{", 
             RowBox[{"i", ",", "j"}], "}"}]}], "]"}]}], "}"}], ">>>", " ", 
        "\"\<C:\\\\Kimb_Mathem\\\\tmp4\>\""}]}], "]"}]}], "]"}]}], 
  "]"}]}], "Input",
 CellChangeTimes->{{3.5681154444347754`*^9, 3.568115497147268*^9}, {
   3.568127214924899*^9, 3.568127265796589*^9}, {3.568128215276656*^9, 
   3.568128238536297*^9}, {3.568130718566253*^9, 3.5681307485339055`*^9}, 
   3.5681308175952272`*^9, {3.5681308916497574`*^9, 3.5681309458598523`*^9}, {
   3.5681309983383446`*^9, 3.5681310118947687`*^9}, {3.5681310723448744`*^9, 
   3.568131074357278*^9}, {3.5681315215476637`*^9, 3.568131584509374*^9}, {
   3.56813162789305*^9, 3.5681316427754765`*^9}, {3.5681317113063965`*^9, 
   3.568131715268804*^9}, {3.5681323203938665`*^9, 3.5681323260098763`*^9}, {
   3.568132387567585*^9, 3.568132440014877*^9}, {3.5681325121026034`*^9, 
   3.568132590945142*^9}, {3.56813274298981*^9, 3.5681327548302307`*^9}, {
   3.568132877867647*^9, 3.568132888335265*^9}, {3.568133052345154*^9, 
   3.568133072344389*^9}, {3.5681331740533676`*^9, 3.5681331951446047`*^9}, {
   3.568133267793932*^9, 3.5681332684803333`*^9}, {3.5681335079111547`*^9, 
   3.568133536833606*^9}, {3.5681336272201643`*^9, 3.5681336314945717`*^9}, {
   3.568133737994359*^9, 3.568133744000369*^9}, {3.568133978515581*^9, 
   3.5681339917912045`*^9}, {3.5681340324916763`*^9, 
   3.5681340418204927`*^9}, {3.568134244262048*^9, 3.5681342454944506`*^9}, {
   3.5681343354754086`*^9, 3.568134336099409*^9}, {3.5681345828762426`*^9, 
   3.568134583952645*^9}, {3.5681348406666956`*^9, 3.5681348474683075`*^9}, 
   3.568134889151581*^9, {3.5681352914606876`*^9, 3.568135296265496*^9}, {
   3.568135518877887*^9, 3.568135519595488*^9}, 3.5681357593679094`*^9, {
   3.5681787035360527`*^9, 3.568178728449296*^9}, 3.5682779128076954`*^9, 
   3.5804489138003826`*^9}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"a", "=", 
  RowBox[{
  "ReadList", "[", "\"\<C:\\\\Kimb_Mathem\\\\tmp4\>\"", 
   "]"}]}], "\[IndentingNewLine]", 
 RowBox[{"nb", "=", 
  RowBox[{"Length", "[", "a", "]"}]}], "\[IndentingNewLine]", 
 RowBox[{"cw", "=", 
  RowBox[{"nb", "-", "bh"}]}]}], "Input",
 CellChangeTimes->{{3.5681343270045934`*^9, 3.568134331481801*^9}, 
   3.5681343823222904`*^9, {3.5681345934218616`*^9, 3.5681345986322703`*^9}, {
   3.568134635479535*^9, 3.5681346521559644`*^9}, {3.5681351451792307`*^9, 
   3.568135145428831*^9}, {3.5681351781888885`*^9, 3.5681352349417877`*^9}, {
   3.568135340163973*^9, 3.568135340335573*^9}, {3.568135523495495*^9, 
   3.5681355244314966`*^9}, {3.5681357674799232`*^9, 
   3.5681357771675406`*^9}, {3.568136284433632*^9, 3.5681362987700567`*^9}, {
   3.568136395053426*^9, 3.5681363965510283`*^9}, {3.568303192134327*^9, 
   3.568303196237134*^9}, {3.568307596598462*^9, 3.5683075983924656`*^9}, 
   3.5804489330664167`*^9, {3.580526129501199*^9, 3.580526130624401*^9}, {
   3.5819194987382727`*^9, 3.581919499346674*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"2", ",", "1", ",", "1.1526553314099903`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"3", ",", "1", ",", "1.1526552850062577`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"4", ",", "1", ",", "1.1526552844103315`"}], "}"}]}], 
  "}"}]], "Output",
 CellChangeTimes->{
  3.581932641602049*^9, 3.5819327493046384`*^9, 3.5819327950907187`*^9, {
   3.581932924976547*^9, 3.581932953009796*^9}, 3.5819505009142513`*^9, 
   3.581997263417365*^9}],

Cell[BoxData["3"], "Output",
 CellChangeTimes->{
  3.581932641602049*^9, 3.5819327493046384`*^9, 3.5819327950907187`*^9, {
   3.581932924976547*^9, 3.581932953009796*^9}, 3.5819505009142513`*^9, 
   3.581997263417365*^9}],

Cell[BoxData["0"], "Output",
 CellChangeTimes->{
  3.581932641602049*^9, 3.5819327493046384`*^9, 3.5819327950907187`*^9, {
   3.581932924976547*^9, 3.581932953009796*^9}, 3.5819505009142513`*^9, 
   3.581997263417365*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{" ", 
  RowBox[{
   RowBox[{
    RowBox[{"wx", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{
        RowBox[{
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"a", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "3"}], "]"}], "]"}], "+", 
            RowBox[{"u2", "[", 
             RowBox[{"[", 
              RowBox[{"a", "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}]}], ")"}], 
          "/", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"u2", "[", 
             RowBox[{"[", 
              RowBox[{"a", "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "+", 
            RowBox[{"2", "*", 
             RowBox[{"a", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "3"}], "]"}], "]"}]}], "+", 
            RowBox[{"u2", "[", 
             RowBox[{"[", 
              RowBox[{"a", "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], 
           ")"}]}], "*", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"x", "[", 
            RowBox[{"[", 
             RowBox[{"a", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}], "-", 
           RowBox[{"x", "[", 
            RowBox[{"[", 
             RowBox[{"a", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}]}], ")"}]}],
         "+", 
        RowBox[{"x", "[", 
         RowBox[{"[", 
          RowBox[{"a", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}]}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "cw"}], "}"}]}], "]"}]}], ";", "      ", 
    RowBox[{"wy", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{
        RowBox[{
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"a", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "3"}], "]"}], "]"}], "+", 
            RowBox[{"u2", "[", 
             RowBox[{"[", 
              RowBox[{"a", "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}]}], ")"}], 
          "/", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"u2", "[", 
             RowBox[{"[", 
              RowBox[{"a", "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "+", 
            RowBox[{"2", "*", 
             RowBox[{"a", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "3"}], "]"}], "]"}]}], "+", 
            RowBox[{"u2", "[", 
             RowBox[{"[", 
              RowBox[{"a", "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], 
           ")"}]}], "*", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"y", "[", 
            RowBox[{"[", 
             RowBox[{"a", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}], "-", 
           RowBox[{"y", "[", 
            RowBox[{"[", 
             RowBox[{"a", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}]}], ")"}]}],
         "+", 
        RowBox[{"y", "[", 
         RowBox[{"[", 
          RowBox[{"a", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}]}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "cw"}], "}"}]}], "]"}]}], ";", "  ", 
    RowBox[{"wz", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{
        RowBox[{
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"a", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "3"}], "]"}], "]"}], "+", 
            RowBox[{"u2", "[", 
             RowBox[{"[", 
              RowBox[{"a", "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}]}], ")"}], 
          "/", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"u2", "[", 
             RowBox[{"[", 
              RowBox[{"a", "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "+", 
            RowBox[{"2", "*", 
             RowBox[{"a", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "3"}], "]"}], "]"}]}], "+", 
            RowBox[{"u2", "[", 
             RowBox[{"[", 
              RowBox[{"a", "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], 
           ")"}]}], "*", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"z", "[", 
            RowBox[{"[", 
             RowBox[{"a", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}], "-", 
           RowBox[{"z", "[", 
            RowBox[{"[", 
             RowBox[{"a", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}]}], ")"}]}],
         "+", 
        RowBox[{"z", "[", 
         RowBox[{"[", 
          RowBox[{"a", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}]}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "cw"}], "}"}]}], "]"}]}], ";"}], "\n", 
   RowBox[{
    RowBox[{"Plot7", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"Graphics3D", "[", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"Opacity", "[", "0.5", "]"}], ",", 
          RowBox[{"Sphere", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{"wx", "[", 
               RowBox[{"[", "i", "]"}], "]"}], ",", 
              RowBox[{"wy", "[", 
               RowBox[{"[", "i", "]"}], "]"}], ",", 
              RowBox[{"wz", "[", 
               RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", 
            RowBox[{"a", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "3"}], "]"}], "]"}]}], "]"}]}], "}"}], "]"}], 
       ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "cw"}], "}"}]}], "]"}]}], 
    ";"}]}]}]], "Input",
 CellChangeTimes->{{3.5681365021632137`*^9, 3.5681365171392403`*^9}, {
   3.568136560351316*^9, 3.5681365863097615`*^9}, {3.5681369135895367`*^9, 
   3.5681369265999594`*^9}, {3.568137821916732*^9, 3.5681378378443604`*^9}, {
   3.5682072492149415`*^9, 3.5682072529901485`*^9}, 3.568209802830227*^9, {
   3.5682216503836985`*^9, 3.5682217583046885`*^9}, {3.5682217912831464`*^9, 
   3.568221843465238*^9}, {3.5682219760654707`*^9, 3.5682221133457117`*^9}, {
   3.568222159615393*^9, 3.5682221693498106`*^9}, {3.5682222565227633`*^9, 
   3.5682222578487663`*^9}, 3.5682224736283445`*^9, 3.5682227227607822`*^9, 
   3.5682227558640404`*^9, {3.568222818810151*^9, 3.5682228486062036`*^9}, {
   3.5682230057764797`*^9, 3.5682230079916835`*^9}, {3.568223041952943*^9, 
   3.5682230548229656`*^9}, 3.5682232619601297`*^9, 3.5682233000709963`*^9, {
   3.568223536645412*^9, 3.5682235381118145`*^9}, {3.568223591151908*^9, 
   3.568223591526308*^9}, {3.5682236491528096`*^9, 3.568223657483224*^9}, {
   3.568223713222122*^9, 3.5682237532361927`*^9}, {3.568223800379475*^9, 
   3.568223827648323*^9}, {3.5682238887380304`*^9, 3.5682238893152313`*^9}, {
   3.5682240287170763`*^9, 3.5682240612431335`*^9}, 3.5682241088076167`*^9, {
   3.568224198772975*^9, 3.5682242396450467`*^9}, {3.568224282170721*^9, 
   3.5682243185031853`*^9}, {3.568224437890195*^9, 3.568224439590598*^9}, {
   3.568224572253231*^9, 3.5682245867144566`*^9}, {3.5682247081762695`*^9, 
   3.568224733557514*^9}, {3.56822480520844*^9, 3.5682248433505073`*^9}, {
   3.568224884206979*^9, 3.5682248959069996`*^9}, {3.568224982830352*^9, 
   3.5682250047639904`*^9}, {3.56822503870965*^9, 3.5682250803305235`*^9}, {
   3.5682251223413973`*^9, 3.568225141451431*^9}, {3.568225234177994*^9, 
   3.5682252929588966`*^9}, {3.5682253775890455`*^9, 3.568225482140429*^9}, {
   3.56822551694409*^9, 3.5682255297205124`*^9}, {3.568225573197789*^9, 
   3.5682255809666023`*^9}, 3.568225660698343*^9, 3.5682257301964645`*^9, {
   3.568225765078126*^9, 3.5682257828933573`*^9}, {3.568225927026843*^9, 
   3.568225927416844*^9}, {3.568267062420238*^9, 3.568267081358671*^9}, {
   3.568267188328059*^9, 3.5682672830358257`*^9}, {3.5682676367040462`*^9, 
   3.568267747979042*^9}, 3.5682678404716043`*^9, 3.568267882030078*^9, {
   3.5682679252577534`*^9, 3.568267996612279*^9}, {3.5682680488879704`*^9, 
   3.568268061133992*^9}, {3.568268100914062*^9, 3.5682681150944867`*^9}, {
   3.5682681672609787`*^9, 3.5682681709113846`*^9}, {3.568268272904364*^9, 
   3.5682683182536435`*^9}, {3.5682685202739983`*^9, 3.568268599194537*^9}, {
   3.5682686348873997`*^9, 3.568268637773405*^9}, {3.5682686944951043`*^9, 
   3.568268872273017*^9}, {3.5682689420207396`*^9, 3.568269016947671*^9}, {
   3.5682695040741262`*^9, 3.568269601153097*^9}, {3.5682696652380095`*^9, 
   3.5682696728352227`*^9}, {3.5682697489633565`*^9, 
   3.5682698308635006`*^9}, {3.568271156491429*^9, 3.568271160063835*^9}, {
   3.568271244927984*^9, 3.568271369010602*^9}, {3.5682716250070515`*^9, 
   3.5682717312276382`*^9}, {3.5682717843301315`*^9, 
   3.5682718576814604`*^9}, {3.568271940673606*^9, 3.5682720536490045`*^9}, {
   3.5682722259669075`*^9, 3.5682723223126764`*^9}, {3.568272370813162*^9, 
   3.568272411466833*^9}, {3.5682725893539457`*^9, 3.568272622238803*^9}, {
   3.5682726936713285`*^9, 3.5682726960113325`*^9}, {3.568272741079812*^9, 
   3.56827275143823*^9}, {3.5682729218373294`*^9, 3.568272991070251*^9}, {
   3.5682734285574193`*^9, 3.5682734302890224`*^9}, {3.5682734898811274`*^9, 
   3.5682735111595645`*^9}, {3.5682736115925407`*^9, 
   3.5682736313577757`*^9}, {3.5682738601945777`*^9, 3.568273987272401*^9}, {
   3.5682763944346285`*^9, 3.5682764313910933`*^9}, 3.5682765128232365`*^9, {
   3.5682766527398825`*^9, 3.568276660976697*^9}, {3.5682774496048822`*^9, 
   3.568277461663703*^9}, {3.5682852978767633`*^9, 3.568285382460112*^9}, {
   3.568285675366227*^9, 3.568285707861084*^9}, {3.568285985822372*^9, 
   3.568285991485182*^9}, {3.568286358600627*^9, 3.568286369505046*^9}, {
   3.5682869898869357`*^9, 3.568287001134555*^9}, {3.5682977174462996`*^9, 
   3.568297735994732*^9}, {3.568305934649544*^9, 3.5683059800144234`*^9}, {
   3.5804489722380853`*^9, 3.580449074449465*^9}, {3.580449125945155*^9, 
   3.580449181777653*^9}, 3.580449245191765*^9, 3.580449330789115*^9, {
   3.5804589977940884`*^9, 3.580459000102893*^9}, {3.5804590597729974`*^9, 
   3.5804590762778263`*^9}, 3.5804591100518856`*^9, 3.580459208191658*^9, {
   3.5804592899358015`*^9, 3.580459292322606*^9}, {3.580459332664277*^9, 
   3.5804593337250786`*^9}, {3.5804599640909853`*^9, 
   3.5804599656041884`*^9}, {3.580460783217224*^9, 3.580460832840912*^9}, {
   3.580460867847373*^9, 3.5804608846018023`*^9}, {3.580465026299877*^9, 
   3.5804650871087837`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"wxh", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"a", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "3"}], "]"}], "]"}], "+", 
          RowBox[{"u2", "[", 
           RowBox[{"[", 
            RowBox[{"a", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ")"}], 
        "/", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"u2", "[", 
           RowBox[{"[", 
            RowBox[{"a", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "+", 
          RowBox[{"a", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "3"}], "]"}], "]"}], "+", 
          RowBox[{"u2", "[", 
           RowBox[{"[", 
            RowBox[{"a", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ")"}]}], 
       "*", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"x", "[", 
          RowBox[{"[", 
           RowBox[{"a", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "-", 
         RowBox[{"x", "[", 
          RowBox[{"[", 
           RowBox[{"a", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ")"}]}], 
      "+", 
      RowBox[{"x", "[", 
       RowBox[{"[", 
        RowBox[{"a", "[", 
         RowBox[{"[", 
          RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", 
       RowBox[{"cw", "+", "1"}], ",", "nb"}], "}"}]}], "]"}]}], ";", "    ", 
  RowBox[{"wyh", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"a", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "3"}], "]"}], "]"}], "+", 
          RowBox[{"u2", "[", 
           RowBox[{"[", 
            RowBox[{"a", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ")"}], 
        "/", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"u2", "[", 
           RowBox[{"[", 
            RowBox[{"a", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "+", 
          RowBox[{"a", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "3"}], "]"}], "]"}], "+", 
          RowBox[{"u2", "[", 
           RowBox[{"[", 
            RowBox[{"a", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ")"}]}], 
       "*", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"y", "[", 
          RowBox[{"[", 
           RowBox[{"a", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "-", 
         RowBox[{"y", "[", 
          RowBox[{"[", 
           RowBox[{"a", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ")"}]}], 
      "+", 
      RowBox[{"y", "[", 
       RowBox[{"[", 
        RowBox[{"a", "[", 
         RowBox[{"[", 
          RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", 
       RowBox[{"cw", "+", "1"}], ",", "nb"}], "}"}]}], "]"}]}], ";", "    ", 
  RowBox[{"wzh", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{
      RowBox[{
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"a", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "3"}], "]"}], "]"}], "+", 
          RowBox[{"u2", "[", 
           RowBox[{"[", 
            RowBox[{"a", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ")"}], 
        "/", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"u2", "[", 
           RowBox[{"[", 
            RowBox[{"a", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "+", 
          RowBox[{"a", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "3"}], "]"}], "]"}], "+", 
          RowBox[{"u2", "[", 
           RowBox[{"[", 
            RowBox[{"a", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ")"}]}], 
       "*", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"z", "[", 
          RowBox[{"[", 
           RowBox[{"a", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "-", 
         RowBox[{"z", "[", 
          RowBox[{"[", 
           RowBox[{"a", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ")"}]}], 
      "+", 
      RowBox[{"z", "[", 
       RowBox[{"[", 
        RowBox[{"a", "[", 
         RowBox[{"[", 
          RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", 
       RowBox[{"cw", "+", "1"}], ",", "nb"}], "}"}]}], "]"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.5681365021632137`*^9, 3.5681365171392403`*^9}, {
   3.568136560351316*^9, 3.5681365863097615`*^9}, {3.5681369135895367`*^9, 
   3.5681369265999594`*^9}, {3.568137821916732*^9, 3.5681378378443604`*^9}, {
   3.5682072492149415`*^9, 3.5682072529901485`*^9}, 3.568209802830227*^9, {
   3.5682216503836985`*^9, 3.5682217583046885`*^9}, {3.5682217912831464`*^9, 
   3.568221843465238*^9}, {3.5682219760654707`*^9, 3.5682221133457117`*^9}, {
   3.568222159615393*^9, 3.5682221693498106`*^9}, {3.5682222565227633`*^9, 
   3.5682222578487663`*^9}, 3.5682224736283445`*^9, 3.5682227227607822`*^9, 
   3.5682227558640404`*^9, {3.568222818810151*^9, 3.5682228486062036`*^9}, {
   3.5682230057764797`*^9, 3.5682230079916835`*^9}, {3.568223041952943*^9, 
   3.5682230548229656`*^9}, 3.5682232619601297`*^9, 3.5682233000709963`*^9, {
   3.568223536645412*^9, 3.5682235381118145`*^9}, {3.568223591151908*^9, 
   3.568223591526308*^9}, {3.5682236491528096`*^9, 3.568223657483224*^9}, {
   3.568223713222122*^9, 3.5682237532361927`*^9}, {3.568223800379475*^9, 
   3.568223827648323*^9}, {3.5682238887380304`*^9, 3.5682238893152313`*^9}, {
   3.5682240287170763`*^9, 3.5682240612431335`*^9}, 3.5682241088076167`*^9, {
   3.568224198772975*^9, 3.5682242396450467`*^9}, {3.568224282170721*^9, 
   3.5682243185031853`*^9}, {3.568224437890195*^9, 3.568224439590598*^9}, {
   3.568224572253231*^9, 3.5682245867144566`*^9}, {3.5682247081762695`*^9, 
   3.568224733557514*^9}, {3.56822480520844*^9, 3.5682248433505073`*^9}, {
   3.568224884206979*^9, 3.5682248959069996`*^9}, {3.568224982830352*^9, 
   3.5682250047639904`*^9}, {3.56822503870965*^9, 3.5682250803305235`*^9}, {
   3.5682251223413973`*^9, 3.568225141451431*^9}, {3.568225234177994*^9, 
   3.5682252929588966`*^9}, {3.5682253775890455`*^9, 3.568225482140429*^9}, {
   3.56822551694409*^9, 3.5682255297205124`*^9}, {3.568225573197789*^9, 
   3.5682255809666023`*^9}, 3.568225660698343*^9, 3.5682257301964645`*^9, {
   3.568225765078126*^9, 3.5682257828933573`*^9}, {3.568225927026843*^9, 
   3.568225927416844*^9}, {3.568267062420238*^9, 3.568267081358671*^9}, {
   3.568267188328059*^9, 3.5682672830358257`*^9}, {3.5682676367040462`*^9, 
   3.568267747979042*^9}, 3.5682678404716043`*^9, 3.568267882030078*^9, {
   3.5682679252577534`*^9, 3.568267996612279*^9}, {3.5682680488879704`*^9, 
   3.568268061133992*^9}, {3.568268100914062*^9, 3.5682681150944867`*^9}, {
   3.5682681672609787`*^9, 3.5682681709113846`*^9}, {3.568268272904364*^9, 
   3.5682683182536435`*^9}, {3.5682685202739983`*^9, 3.568268599194537*^9}, {
   3.5682686348873997`*^9, 3.568268637773405*^9}, {3.5682686944951043`*^9, 
   3.568268872273017*^9}, {3.5682689420207396`*^9, 3.568269016947671*^9}, {
   3.5682695040741262`*^9, 3.568269601153097*^9}, {3.5682696652380095`*^9, 
   3.5682696728352227`*^9}, {3.5682697489633565`*^9, 
   3.5682698308635006`*^9}, {3.568271156491429*^9, 3.568271160063835*^9}, {
   3.568271244927984*^9, 3.568271369010602*^9}, {3.5682716250070515`*^9, 
   3.5682717312276382`*^9}, {3.5682717843301315`*^9, 
   3.5682718576814604`*^9}, {3.568271940673606*^9, 3.5682720536490045`*^9}, {
   3.5682722259669075`*^9, 3.5682723223126764`*^9}, {3.568272370813162*^9, 
   3.568272411466833*^9}, {3.5682725893539457`*^9, 3.568272622238803*^9}, {
   3.5682726936713285`*^9, 3.5682726960113325`*^9}, {3.568272741079812*^9, 
   3.56827275143823*^9}, {3.5682729218373294`*^9, 3.568272991070251*^9}, {
   3.5682734285574193`*^9, 3.5682734302890224`*^9}, {3.5682734898811274`*^9, 
   3.5682735111595645`*^9}, {3.5682736115925407`*^9, 
   3.5682736313577757`*^9}, {3.5682738601945777`*^9, 3.568273987272401*^9}, {
   3.5682763944346285`*^9, 3.5682764313910933`*^9}, 3.5682765128232365`*^9, {
   3.5682766527398825`*^9, 3.568276660976697*^9}, {3.5682774496048822`*^9, 
   3.568277461663703*^9}, {3.5682852978767633`*^9, 3.568285382460112*^9}, {
   3.568285675366227*^9, 3.568285707861084*^9}, {3.568285985822372*^9, 
   3.568285991485182*^9}, {3.568286358600627*^9, 3.568286369505046*^9}, {
   3.5682869898869357`*^9, 3.568287001134555*^9}, {3.5682977174462996`*^9, 
   3.568297735994732*^9}, {3.568305934649544*^9, 3.5683059800144234`*^9}, {
   3.5804489722380853`*^9, 3.580449074449465*^9}, {3.580449125945155*^9, 
   3.580449181777653*^9}, 3.580449245191765*^9, 3.580449330789115*^9, {
   3.5804589977940884`*^9, 3.580459000102893*^9}, {3.5804590597729974`*^9, 
   3.5804590762778263`*^9}, 3.5804591100518856`*^9, 3.580459208191658*^9, {
   3.5804592899358015`*^9, 3.580459292322606*^9}, {3.580459332664277*^9, 
   3.5804593337250786`*^9}, {3.5804599640909853`*^9, 
   3.5804599656041884`*^9}, {3.580460783217224*^9, 3.580460832840912*^9}, {
   3.580460867847373*^9, 3.5804608846018023`*^9}, {3.580465026299877*^9, 
   3.5804650896203885`*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"Plot8", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"Graphics3D", "[", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"Opacity", "[", "0.5", "]"}], ",", 
        RowBox[{"Sphere", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"wxh", "[", 
             RowBox[{"[", "i", "]"}], "]"}], ",", 
            RowBox[{"wyh", "[", 
             RowBox[{"[", "i", "]"}], "]"}], ",", 
            RowBox[{"wzh", "[", 
             RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", 
          RowBox[{"a", "[", 
           RowBox[{"[", 
            RowBox[{
             RowBox[{"i", "+", "cw"}], ",", "3"}], "]"}], "]"}]}], "]"}]}], 
       "}"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "bh"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Plot8ib", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"Graphics3D", "[", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"Opacity", "[", "0.4", "]"}], ",", 
        RowBox[{"Lighter", "[", "Blue", "]"}], ",", 
        RowBox[{"Sphere", "[", 
         RowBox[{
          RowBox[{"tpyrn", "[", 
           RowBox[{"[", "i", "]"}], "]"}], ",", "rln"}], "]"}]}], "}"}], 
      "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "1"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Plot9", "=", " ", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"Graphics3D", "[", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"Lighter", "[", 
         RowBox[{"Gray", ",", "0.8"}], "]"}], ",", 
        RowBox[{"Tube", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{"x", "[", 
               RowBox[{"[", 
                RowBox[{"a", "[", 
                 RowBox[{"[", 
                  RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], ",", 
              RowBox[{"y", "[", 
               RowBox[{"[", 
                RowBox[{"a", "[", 
                 RowBox[{"[", 
                  RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], ",", 
              RowBox[{"z", "[", 
               RowBox[{"[", 
                RowBox[{"a", "[", 
                 RowBox[{"[", 
                  RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}]}], 
             "}"}], ",", 
            RowBox[{"{", 
             RowBox[{
              RowBox[{"x", "[", 
               RowBox[{"[", 
                RowBox[{"a", "[", 
                 RowBox[{"[", 
                  RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}], ",", 
              RowBox[{"y", "[", 
               RowBox[{"[", 
                RowBox[{"a", "[", 
                 RowBox[{"[", 
                  RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}], ",", 
              RowBox[{"z", "[", 
               RowBox[{"[", 
                RowBox[{"a", "[", 
                 RowBox[{"[", 
                  RowBox[{"i", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], 
             "}"}]}], "}"}], ",", "0.07"}], "]"}]}], "}"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "nb"}], "}"}]}], "]"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.5681365021632137`*^9, 3.5681365171392403`*^9}, {
   3.568136560351316*^9, 3.5681365863097615`*^9}, {3.5681369135895367`*^9, 
   3.5681369265999594`*^9}, {3.568137821916732*^9, 3.5681378378443604`*^9}, {
   3.5682072492149415`*^9, 3.5682072529901485`*^9}, 3.568209802830227*^9, {
   3.5682216503836985`*^9, 3.5682217583046885`*^9}, {3.5682217912831464`*^9, 
   3.568221843465238*^9}, {3.5682219760654707`*^9, 3.5682221133457117`*^9}, {
   3.568222159615393*^9, 3.5682221693498106`*^9}, {3.5682222565227633`*^9, 
   3.5682222578487663`*^9}, 3.5682224736283445`*^9, 3.5682227227607822`*^9, 
   3.5682227558640404`*^9, {3.568222818810151*^9, 3.5682228486062036`*^9}, {
   3.5682230057764797`*^9, 3.5682230079916835`*^9}, {3.568223041952943*^9, 
   3.5682230548229656`*^9}, 3.5682232619601297`*^9, 3.5682233000709963`*^9, {
   3.568223536645412*^9, 3.5682235381118145`*^9}, {3.568223591151908*^9, 
   3.568223591526308*^9}, {3.5682236491528096`*^9, 3.568223657483224*^9}, {
   3.568223713222122*^9, 3.5682237532361927`*^9}, {3.568223800379475*^9, 
   3.568223827648323*^9}, {3.5682238887380304`*^9, 3.5682238893152313`*^9}, {
   3.5682240287170763`*^9, 3.5682240612431335`*^9}, 3.5682241088076167`*^9, {
   3.568224198772975*^9, 3.5682242396450467`*^9}, {3.568224282170721*^9, 
   3.5682243185031853`*^9}, {3.568224437890195*^9, 3.568224439590598*^9}, {
   3.568224572253231*^9, 3.5682245867144566`*^9}, {3.5682247081762695`*^9, 
   3.568224733557514*^9}, {3.56822480520844*^9, 3.5682248433505073`*^9}, {
   3.568224884206979*^9, 3.5682248959069996`*^9}, {3.568224982830352*^9, 
   3.5682250047639904`*^9}, {3.56822503870965*^9, 3.5682250803305235`*^9}, {
   3.5682251223413973`*^9, 3.568225141451431*^9}, {3.568225234177994*^9, 
   3.5682252929588966`*^9}, {3.5682253775890455`*^9, 3.568225482140429*^9}, {
   3.56822551694409*^9, 3.5682255297205124`*^9}, {3.568225573197789*^9, 
   3.5682255809666023`*^9}, 3.568225660698343*^9, 3.5682257301964645`*^9, {
   3.568225765078126*^9, 3.5682257828933573`*^9}, {3.568225927026843*^9, 
   3.568225927416844*^9}, {3.568267062420238*^9, 3.568267081358671*^9}, {
   3.568267188328059*^9, 3.5682672830358257`*^9}, {3.5682676367040462`*^9, 
   3.568267747979042*^9}, 3.5682678404716043`*^9, 3.568267882030078*^9, {
   3.5682679252577534`*^9, 3.568267996612279*^9}, {3.5682680488879704`*^9, 
   3.568268061133992*^9}, {3.568268100914062*^9, 3.5682681150944867`*^9}, {
   3.5682681672609787`*^9, 3.5682681709113846`*^9}, {3.568268272904364*^9, 
   3.5682683182536435`*^9}, {3.5682685202739983`*^9, 3.568268599194537*^9}, {
   3.5682686348873997`*^9, 3.568268637773405*^9}, {3.5682686944951043`*^9, 
   3.568268872273017*^9}, {3.5682689420207396`*^9, 3.568269016947671*^9}, {
   3.5682695040741262`*^9, 3.568269601153097*^9}, {3.5682696652380095`*^9, 
   3.5682696728352227`*^9}, {3.5682697489633565`*^9, 
   3.5682698308635006`*^9}, {3.568271156491429*^9, 3.568271160063835*^9}, {
   3.568271244927984*^9, 3.568271369010602*^9}, {3.5682716250070515`*^9, 
   3.5682717312276382`*^9}, {3.5682717843301315`*^9, 
   3.5682718576814604`*^9}, {3.568271940673606*^9, 3.5682720536490045`*^9}, {
   3.5682722259669075`*^9, 3.5682723223126764`*^9}, {3.568272370813162*^9, 
   3.568272411466833*^9}, {3.5682725893539457`*^9, 3.568272622238803*^9}, {
   3.5682726936713285`*^9, 3.5682726960113325`*^9}, {3.568272741079812*^9, 
   3.56827275143823*^9}, {3.5682729218373294`*^9, 3.568272991070251*^9}, {
   3.5682734285574193`*^9, 3.5682734302890224`*^9}, {3.5682734898811274`*^9, 
   3.5682735111595645`*^9}, {3.5682736115925407`*^9, 
   3.5682736313577757`*^9}, {3.5682738601945777`*^9, 3.568273987272401*^9}, {
   3.5682763944346285`*^9, 3.5682764313910933`*^9}, 3.5682765128232365`*^9, {
   3.5682766527398825`*^9, 3.568276660976697*^9}, {3.5682774496048822`*^9, 
   3.568277461663703*^9}, {3.5682852978767633`*^9, 3.568285382460112*^9}, {
   3.568285675366227*^9, 3.568285707861084*^9}, {3.568285985822372*^9, 
   3.568285991485182*^9}, {3.568286358600627*^9, 3.568286369505046*^9}, {
   3.5682869898869357`*^9, 3.568287001134555*^9}, {3.5682977174462996`*^9, 
   3.568297735994732*^9}, {3.568305934649544*^9, 3.5683059800144234`*^9}, {
   3.5804489722380853`*^9, 3.580449074449465*^9}, {3.580449125945155*^9, 
   3.580449181777653*^9}, 3.580449245191765*^9, 3.580449330789115*^9, {
   3.5804589977940884`*^9, 3.580459000102893*^9}, {3.5804590597729974`*^9, 
   3.5804590762778263`*^9}, 3.5804591100518856`*^9, 3.580459208191658*^9, {
   3.5804592899358015`*^9, 3.580459292322606*^9}, {3.580459332664277*^9, 
   3.5804593337250786`*^9}, {3.5804599640909853`*^9, 
   3.5804599656041884`*^9}, {3.580460783217224*^9, 3.580460832840912*^9}, {
   3.580460867847373*^9, 3.5804608846018023`*^9}, {3.580465026299877*^9, 
   3.5804651023032103`*^9}, {3.5804653879553123`*^9, 3.580465390466917*^9}, {
   3.5804654257853785`*^9, 3.580465427017781*^9}, {3.581920158213831*^9, 
   3.5819201805842705`*^9}, {3.581920216495534*^9, 3.581920218461137*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Show", "[", 
  RowBox[{
  "Plot1o", ",", "Plot2", ",", "Plot7", ",", "Plot8", ",", "Plot8ib", ",", 
   "Plot9", ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"AspectRatio", "\[Rule]", "Automatic"}], ",", 
     RowBox[{"Boxed", "\[Rule]", "False"}], ",", 
     RowBox[{"SphericalRegion", "->", "True"}]}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.5681365021632137`*^9, 3.5681365171392403`*^9}, {
   3.568136560351316*^9, 3.5681365863097615`*^9}, {3.5681369135895367`*^9, 
   3.5681369265999594`*^9}, {3.568137821916732*^9, 3.5681378378443604`*^9}, {
   3.5682072492149415`*^9, 3.5682072529901485`*^9}, 3.568209802830227*^9, {
   3.5682216503836985`*^9, 3.5682217583046885`*^9}, {3.5682217912831464`*^9, 
   3.568221843465238*^9}, {3.5682219760654707`*^9, 3.5682221133457117`*^9}, {
   3.568222159615393*^9, 3.5682221693498106`*^9}, {3.5682222565227633`*^9, 
   3.5682222578487663`*^9}, 3.5682224736283445`*^9, 3.5682227227607822`*^9, 
   3.5682227558640404`*^9, {3.568222818810151*^9, 3.5682228486062036`*^9}, {
   3.5682230057764797`*^9, 3.5682230079916835`*^9}, {3.568223041952943*^9, 
   3.5682230548229656`*^9}, 3.5682232619601297`*^9, 3.5682233000709963`*^9, {
   3.568223536645412*^9, 3.5682235381118145`*^9}, {3.568223591151908*^9, 
   3.568223591526308*^9}, {3.5682236491528096`*^9, 3.568223657483224*^9}, {
   3.568223713222122*^9, 3.5682237532361927`*^9}, {3.568223800379475*^9, 
   3.568223827648323*^9}, {3.5682238887380304`*^9, 3.5682238893152313`*^9}, {
   3.5682240287170763`*^9, 3.5682240612431335`*^9}, 3.5682241088076167`*^9, {
   3.568224198772975*^9, 3.5682242396450467`*^9}, {3.568224282170721*^9, 
   3.5682243185031853`*^9}, {3.568224437890195*^9, 3.568224439590598*^9}, {
   3.568224572253231*^9, 3.5682245867144566`*^9}, {3.5682247081762695`*^9, 
   3.568224733557514*^9}, {3.56822480520844*^9, 3.5682248433505073`*^9}, {
   3.568224884206979*^9, 3.5682248959069996`*^9}, {3.568224982830352*^9, 
   3.5682250047639904`*^9}, {3.56822503870965*^9, 3.5682250803305235`*^9}, {
   3.5682251223413973`*^9, 3.568225141451431*^9}, {3.568225234177994*^9, 
   3.5682252929588966`*^9}, {3.5682253775890455`*^9, 3.568225482140429*^9}, {
   3.56822551694409*^9, 3.5682255297205124`*^9}, {3.568225573197789*^9, 
   3.5682255809666023`*^9}, 3.568225660698343*^9, 3.5682257301964645`*^9, {
   3.568225765078126*^9, 3.5682257828933573`*^9}, {3.568225927026843*^9, 
   3.568225927416844*^9}, {3.568267062420238*^9, 3.568267081358671*^9}, {
   3.568267188328059*^9, 3.5682672830358257`*^9}, {3.5682676367040462`*^9, 
   3.568267747979042*^9}, 3.5682678404716043`*^9, 3.568267882030078*^9, {
   3.5682679252577534`*^9, 3.568267996612279*^9}, {3.5682680488879704`*^9, 
   3.568268061133992*^9}, {3.568268100914062*^9, 3.5682681150944867`*^9}, {
   3.5682681672609787`*^9, 3.5682681709113846`*^9}, {3.568268272904364*^9, 
   3.5682683182536435`*^9}, {3.5682685202739983`*^9, 3.568268599194537*^9}, {
   3.5682686348873997`*^9, 3.568268637773405*^9}, {3.5682686944951043`*^9, 
   3.568268872273017*^9}, {3.5682689420207396`*^9, 3.568269016947671*^9}, {
   3.5682695040741262`*^9, 3.568269601153097*^9}, {3.5682696652380095`*^9, 
   3.5682696728352227`*^9}, {3.5682697489633565`*^9, 
   3.5682698308635006`*^9}, {3.568271156491429*^9, 3.568271160063835*^9}, {
   3.568271244927984*^9, 3.568271369010602*^9}, {3.5682716250070515`*^9, 
   3.5682717312276382`*^9}, {3.5682717843301315`*^9, 
   3.5682718576814604`*^9}, {3.568271940673606*^9, 3.5682720536490045`*^9}, {
   3.5682722259669075`*^9, 3.5682723223126764`*^9}, {3.568272370813162*^9, 
   3.568272411466833*^9}, {3.5682725893539457`*^9, 3.568272622238803*^9}, {
   3.5682726936713285`*^9, 3.5682726960113325`*^9}, {3.568272741079812*^9, 
   3.56827275143823*^9}, {3.5682729218373294`*^9, 3.568272991070251*^9}, {
   3.5682734285574193`*^9, 3.5682734302890224`*^9}, {3.5682734898811274`*^9, 
   3.5682735111595645`*^9}, {3.5682736115925407`*^9, 
   3.5682736313577757`*^9}, {3.5682738601945777`*^9, 3.568273987272401*^9}, {
   3.5682763944346285`*^9, 3.5682764313910933`*^9}, 3.5682765128232365`*^9, {
   3.5682766527398825`*^9, 3.568276660976697*^9}, {3.5682774496048822`*^9, 
   3.568277461663703*^9}, {3.5682852978767633`*^9, 3.568285382460112*^9}, {
   3.568285675366227*^9, 3.568285707861084*^9}, {3.568285985822372*^9, 
   3.568285991485182*^9}, {3.568286358600627*^9, 3.568286369505046*^9}, {
   3.5682869898869357`*^9, 3.568287001134555*^9}, {3.5682977174462996`*^9, 
   3.568297735994732*^9}, {3.568305934649544*^9, 3.5683059800144234`*^9}, {
   3.5804489722380853`*^9, 3.580449074449465*^9}, {3.580449125945155*^9, 
   3.580449181777653*^9}, 3.580449245191765*^9, 3.580449330789115*^9, {
   3.5804589977940884`*^9, 3.580459000102893*^9}, {3.5804590597729974`*^9, 
   3.5804590762778263`*^9}, 3.5804591100518856`*^9, 3.580459208191658*^9, {
   3.5804592899358015`*^9, 3.580459292322606*^9}, {3.580459332664277*^9, 
   3.5804593337250786`*^9}, {3.5804599640909853`*^9, 
   3.5804599656041884`*^9}, {3.580460783217224*^9, 3.580460832840912*^9}, {
   3.580460867847373*^9, 3.5804608846018023`*^9}, {3.580465026299877*^9, 
   3.580465098200403*^9}, 3.580465296975952*^9, {3.580488471874737*^9, 
   3.580488483902358*^9}, 3.5819327408494234`*^9}],

Cell[BoxData[
 Graphics3DBox[{
   {RGBColor[
     NCache[
      Rational[1, 3], 0.3333333333333333], 
     NCache[
      Rational[1, 3], 0.3333333333333333], 1], Opacity[0.8], 
    SphereBox[{{1.7007534*^-9, 3.5904794*^-9, -0.1499958505958386}}, 0.24]}, 
   {RGBColor[1, 1, 1], Opacity[0.8], 
    SphereBox[{{-2.6456164*^-9, -1.7760899227176337`, 0.5845065307586286}}, 
     0.12]}, 
   {RGBColor[1, 1, 1], Opacity[0.8], 
    SphereBox[{{-1.5381389454241607`, 0.8880449396269678, 
     0.5845065436087654}}, 0.12]}, 
   {RGBColor[1, 1, 1], Opacity[0.8], 
    SphereBox[{{1.5381389463690238`, 0.8880449432174472, 0.5845065428528751}},
      0.12]}, Text3DBox["1", {-0.4999999982992466, 
    3.5904794*^-9, -0.1499958505958386}, {1, 0}], 
   Text3DBox["2", {-0.5000000026456164, -1.7760899227176337`, 
    0.5845065307586286}, {1, 0}], 
   Text3DBox["3", {-2.0381389454241607`, 0.8880449396269678, 
    0.5845065436087654}, {1, 0}], 
   Text3DBox["4", {1.0381389463690238`, 0.8880449432174472, 
    0.5845065428528751}, {1, 0}], 
   {Opacity[0.5], 
    SphereBox[{-1.392570397368313*^-9, -1.2640482684445002`, 
     0.3727515915264284}, 1.1526553314099903`]}, 
   {Opacity[0.5], 
    SphereBox[{-1.0946978682933797`, 0.6320241141270103, 0.3727515955593445}, 
     1.1526552850062577`]}, 
   {Opacity[0.5], 
    SphereBox[{1.094697869808993, 0.6320241166029831, 0.3727515949557186}, 
     1.1526552844103315`]}, 
   {RGBColor[
     NCache[
      Rational[1, 3], 0.3333333333333333], 
     NCache[
      Rational[1, 3], 0.3333333333333333], 1], Opacity[0.4], 
    SphereBox[{1.364313995111388*^-9, 
     1.0000505534790705`*^-8, -1.5192158505958386`}, 1.154]}, 
   {RGBColor[0.9, 0.9, 0.9], 
    TubeBox[{{-2.6456164*^-9, -1.7760899227176337`, 0.5845065307586286}, {
     1.7007534*^-9, 3.5904794*^-9, -0.1499958505958386}}, 0.07]}, 
   {RGBColor[0.9, 0.9, 0.9], 
    TubeBox[{{-1.5381389454241607`, 0.8880449396269678, 0.5845065436087654}, {
     1.7007534*^-9, 3.5904794*^-9, -0.1499958505958386}}, 0.07]}, 
   {RGBColor[0.9, 0.9, 0.9], 
    TubeBox[{{1.5381389463690238`, 0.8880449432174472, 0.5845065428528751}, {
     1.7007534*^-9, 3.5904794*^-9, -0.1499958505958386}}, 0.07]}},
  AspectRatio->Automatic,
  AutomaticImageSize->True,
  Boxed->False,
  ImageSize->{360., 360.},
  SphericalRegion->True,
  ViewPoint->{-1.7925865707839068`, 1.1850661001709968`, -2.613857632785439},
  ViewVertical->{0.29882112383014037`, 
   0.7417348550373625, -0.7019865681424263}]], "Output",
 CellChangeTimes->{
  3.5819326420700502`*^9, 3.5819327493514385`*^9, 3.5819327951219187`*^9, {
   3.5819329250077467`*^9, 3.581932953040996*^9}, 3.581950500961051*^9, 
   3.5819972634485655`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"tc", "=", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        RowBox[{"x", "[", 
         RowBox[{"[", "i", "]"}], "]"}], ",", 
        RowBox[{"y", "[", 
         RowBox[{"[", "i", "]"}], "]"}], ",", 
        RowBox[{"z", "[", 
         RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "bc"}], "}"}]}], "]"}]}], ";"}], "   ", 
  RowBox[{"(*", " ", 
   RowBox[{"core", " ", "clouds"}], " ", "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"twc", "=", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        RowBox[{"wx", "[", 
         RowBox[{"[", "i", "]"}], "]"}], ",", 
        RowBox[{"wy", "[", 
         RowBox[{"[", "i", "]"}], "]"}], ",", 
        RowBox[{"wz", "[", 
         RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "cw"}], "}"}]}], "]"}]}], ";"}], "  ", 
  RowBox[{"(*", " ", 
   RowBox[{"Bond", " ", "c\[ODoubleDot]ouds", " ", "between", " ", "cores"}], 
   " ", "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"twh", "=", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        RowBox[{"wxh", "[", 
         RowBox[{"[", "i", "]"}], "]"}], ",", 
        RowBox[{"wyh", "[", 
         RowBox[{"[", "i", "]"}], "]"}], ",", 
        RowBox[{"wzh", "[", 
         RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "bh"}], "}"}]}], "]"}]}], ";"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{"H", "-", 
    RowBox[{"bonding", " ", "clouds"}]}], " ", "*)"}], 
  " "}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"tw", "=", 
   RowBox[{"Join", "[", 
    RowBox[{"tc", ",", "twc", ",", "twh", ",", "tpyrn"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"bt", "=", 
  RowBox[{"Length", "[", "tw", "]"}]}]}], "Input",
 CellChangeTimes->{{3.5804614754216404`*^9, 3.5804614810532503`*^9}, {
   3.580461527603732*^9, 3.5804615546541796`*^9}, 3.5804647573086047`*^9, 
   3.5804649196424894`*^9, {3.58048949832534*^9, 3.5804895043313503`*^9}, 
   3.580489609803136*^9, {3.5819186625846424`*^9, 3.5819186639886446`*^9}}],

Cell[BoxData["5"], "Output",
 CellChangeTimes->{
  3.58193264208565*^9, 3.5819327494294386`*^9, 3.581932795137519*^9, {
   3.5819329250545473`*^9, 3.581932953087796*^9}, 3.581950500961051*^9, 
   3.5819972634641657`*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{"cloud", "-", 
    RowBox[{"cloud", " ", "distances"}]}], " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"dw", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"Norm", "[", 
        RowBox[{
         RowBox[{"tw", "[", 
          RowBox[{"[", "i", "]"}], "]"}], "-", 
         RowBox[{"tw", "[", 
          RowBox[{"[", "j", "]"}], "]"}]}], "]"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "bt"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"j", ",", "1", ",", "bt"}], "}"}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{"MatrixForm", "[", "dw", "]"}], "*)"}]}]}]], "Input",
 CellChangeTimes->{{3.5804614754216404`*^9, 3.5804614810532503`*^9}, {
   3.580461527603732*^9, 3.5804615546541796`*^9}, {3.5804647573086047`*^9, 
   3.58046477754184*^9}, {3.580464923870097*^9, 3.5804649422313294`*^9}, 
   3.5804897528397865`*^9, 3.580489783634241*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{"nucleus", "-", 
    RowBox[{"nucleus", " ", "distances"}]}], " ", "*)"}], "\n", 
  RowBox[{
   RowBox[{
    RowBox[{"tv1", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"t", "[", 
        RowBox[{"[", "i", "]"}], "]"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "bc"}], "}"}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"tv2", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "cw"}], "}"}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"tv3", "=", " ", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"t", "[", 
        RowBox[{"[", "i", "]"}], "]"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", 
         RowBox[{"bc", "+", "1"}], ",", "n"}], "}"}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"tv6", " ", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "1"}], "}"}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"tv", "  ", "=", " ", 
     RowBox[{"Join", "[", 
      RowBox[{"tv1", ",", "tv2", ",", "tv3", ",", "tv6"}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{"MatrixForm", "[", "tv", "]"}], "*)"}], "\[IndentingNewLine]", 
   RowBox[{"Length", "[", "tv", "]"}]}]}]], "Input",
 CellChangeTimes->{{3.5804614754216404`*^9, 3.5804614810532503`*^9}, {
   3.580461527603732*^9, 3.5804615546541796`*^9}, {3.5804647573086047`*^9, 
   3.58046477754184*^9}, {3.580464923870097*^9, 3.5804649550077515`*^9}, {
   3.5804899034268513`*^9, 3.5804899043784533`*^9}, 3.58191868394108*^9}],

Cell[BoxData["5"], "Output",
 CellChangeTimes->{
  3.5819326421012497`*^9, 3.5819327494450383`*^9, 3.581932795153119*^9, {
   3.581932925085747*^9, 3.5819329531189966`*^9}, 3.581950500976651*^9, 
   3.5819972634797654`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"ow", "=", 
   RowBox[{"Join", "[", 
    RowBox[{
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"-", "2.0"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", 
         RowBox[{"bc", "+", "nb"}]}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"-", "2.0"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "1"}], "}"}]}], "]"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{" ", "ow"}], "\[IndentingNewLine]", 
 RowBox[{"cnv", "=", 
  RowBox[{"Join", "[", 
   RowBox[{
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"cn", "[", 
       RowBox[{"[", "i", "]"}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "bc"}], "}"}]}], "]"}], ",", 
    RowBox[{"Table", "[", 
     RowBox[{"0.0", ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "cw"}], "}"}]}], "]"}], ",", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"cn", "[", 
       RowBox[{"[", "i", "]"}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", 
        RowBox[{"bc", "+", "1"}], ",", "n"}], "}"}]}], "]"}], ",", 
    RowBox[{"Table", "[", 
     RowBox[{"0.0", ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "1"}], "}"}]}], "]"}]}], "]"}]}]}], "Input",
 CellChangeTimes->{{3.580461670421983*^9, 3.5804617339140944`*^9}, {
  3.5819187786488466`*^9, 3.5819187833756547`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", "2.`"}], ",", 
   RowBox[{"-", "2.`"}], ",", 
   RowBox[{"-", "2.`"}], ",", 
   RowBox[{"-", "2.`"}], ",", 
   RowBox[{"-", "2.`"}]}], "}"}]], "Output",
 CellChangeTimes->{
  3.58193264211685*^9, 3.5819327494606385`*^9, 3.581932795168719*^9, {
   3.581932925085747*^9, 3.5819329531345963`*^9}, 3.5819505009922514`*^9, 
   3.5819972634797654`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"7.`", ",", "1.`", ",", "1.`", ",", "1.`", ",", "0.`"}], 
  "}"}]], "Output",
 CellChangeTimes->{
  3.58193264211685*^9, 3.5819327494606385`*^9, 3.581932795168719*^9, {
   3.581932925085747*^9, 3.5819329531345963`*^9}, 3.5819505009922514`*^9, 
   3.581997263495365*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{"nucleus", "-", 
    RowBox[{"cloud", " ", "distances"}]}], " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"dcw", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"sq", "=", 
         RowBox[{"Norm", "[", 
          RowBox[{
           RowBox[{"tw", "[", 
            RowBox[{"[", "i", "]"}], "]"}], "-", 
           RowBox[{"tv", "[", 
            RowBox[{"[", "j", "]"}], "]"}]}], "]"}]}], ";", 
        RowBox[{"If", "[", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"sq", ">", " ", "0.1"}], " ", "&&", " ", 
            RowBox[{"i", "\[NotEqual]", "j"}]}], ")"}], ",", 
          RowBox[{"1", "/", "sq"}], ",", "0"}], "]"}]}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "bt"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"j", ",", "1", ",", "bt"}], "}"}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{"Length", "[", "dcw", "]"}]}]}]], "Input",
 CellChangeTimes->{3.5804649749445868`*^9}],

Cell[BoxData["5"], "Output",
 CellChangeTimes->{
  3.58193264213245*^9, 3.5819327494762383`*^9, 3.581932795168719*^9, {
   3.581932925116947*^9, 3.581932953150196*^9}, 3.5819505009922514`*^9, 
   3.581997263495365*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"Elcharge", "=", 
  RowBox[{"Sum", "[", 
   RowBox[{
    RowBox[{"ow", "[", 
     RowBox[{"[", "i", "]"}], "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"i", ",", "1", ",", 
      RowBox[{"Length", "[", "ow", "]"}]}], "}"}]}], 
   "]"}]}], "\[IndentingNewLine]", 
 RowBox[{"Length", "[", "ow", "]"}], "\[IndentingNewLine]", 
 RowBox[{"Nuccharge", "=", 
  RowBox[{"Sum", "[", 
   RowBox[{
    RowBox[{"cnv", "[", 
     RowBox[{"[", "i", "]"}], "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"i", ",", "1", ",", 
      RowBox[{"Length", "[", "cnv", "]"}]}], "}"}]}], 
   "]"}]}], "\[IndentingNewLine]", 
 RowBox[{"Length", "[", "cnv", "]"}]}], "Input",
 CellChangeTimes->{3.580489056610564*^9}],

Cell[BoxData[
 RowBox[{"-", "10.`"}]], "Output",
 CellChangeTimes->{
  3.5819326421480503`*^9, 3.5819327494762383`*^9, 3.5819327951843185`*^9, {
   3.581932925116947*^9, 3.581932953150196*^9}, 3.581950501007851*^9, 
   3.5819972635109653`*^9}],

Cell[BoxData["5"], "Output",
 CellChangeTimes->{
  3.5819326421480503`*^9, 3.5819327494762383`*^9, 3.5819327951843185`*^9, {
   3.581932925116947*^9, 3.581932953150196*^9}, 3.581950501007851*^9, 
   3.5819972635109653`*^9}],

Cell[BoxData["10.`"], "Output",
 CellChangeTimes->{
  3.5819326421480503`*^9, 3.5819327494762383`*^9, 3.5819327951843185`*^9, {
   3.581932925116947*^9, 3.581932953150196*^9}, 3.581950501007851*^9, 
   3.5819972635109653`*^9}],

Cell[BoxData["5"], "Output",
 CellChangeTimes->{
  3.5819326421480503`*^9, 3.5819327494762383`*^9, 3.5819327951843185`*^9, {
   3.581932925116947*^9, 3.581932953150196*^9}, 3.581950501007851*^9, 
   3.5819972635109653`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"Vee", "=", "0.0"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"For", "[", 
  RowBox[{
   RowBox[{"i", "=", "1"}], ",", 
   RowBox[{"i", "<", "bt"}], ",", 
   RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
   RowBox[{"For", "[", 
    RowBox[{
     RowBox[{"j", "=", 
      RowBox[{"i", "+", "1"}]}], ",", 
     RowBox[{"j", "\[LessEqual]", "bt"}], ",", 
     RowBox[{"j", "++"}], ",", "\[IndentingNewLine]", 
     RowBox[{"Vee", "=", 
      RowBox[{"Vee", "+", 
       RowBox[{
        RowBox[{"ow", "[", 
         RowBox[{"[", "i", "]"}], "]"}], "*", 
        RowBox[{
         RowBox[{"ow", "[", 
          RowBox[{"[", "j", "]"}], "]"}], "/", 
         RowBox[{"dw", "[", 
          RowBox[{"[", 
           RowBox[{"i", ",", "j"}], "]"}], "]"}]}]}]}]}]}], "]"}]}], 
  "]"}], "\[IndentingNewLine]", "Vee", "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"sc", "=", 
    StyleBox["0.34885",
     FontColor->RGBColor[1, 0, 0]]}], ";"}], "   ", 
  RowBox[{"(*", " ", 
   RowBox[{"core", " ", "clouds"}], " ", "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{"For", "[", 
  RowBox[{
   RowBox[{"i", "=", "1"}], " ", ",", 
   RowBox[{"i", "\[LessEqual]", "bc"}], ",", 
   RowBox[{"i", "++"}], ",", 
   RowBox[{"Vee", "=", 
    RowBox[{"Vee", "+", 
     RowBox[{"3", "*", 
      RowBox[{"sc", "/", 
       RowBox[{"u2", "[", 
        RowBox[{"[", "i", "]"}], "]"}]}]}]}]}]}], 
  "]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"sw", "=", 
    StyleBox["0.34885",
     FontColor->RGBColor[1, 0, 0]]}], ";"}], "  ", 
  RowBox[{"(*", " ", 
   RowBox[{"bonding", " ", "clouds"}], " ", "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{"For", "[", 
  RowBox[{
   RowBox[{"i", "=", "1"}], " ", ",", 
   RowBox[{"i", "\[LessEqual]", "nb"}], ",", 
   RowBox[{"i", "++"}], ",", 
   RowBox[{"Vee", "=", 
    RowBox[{"Vee", "+", 
     RowBox[{"3", "*", 
      RowBox[{"sw", "/", 
       RowBox[{"a", "[", 
        RowBox[{"[", 
         RowBox[{"i", ",", "3"}], "]"}], "]"}]}]}]}]}]}], 
  "]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Vee", "=", 
   RowBox[{"Vee", "+", 
    RowBox[{"3", "*", 
     RowBox[{"sw", "/", "rln"}]}]}]}], 
  ";"}], "\[IndentingNewLine]", "Vee"}], "Input",
 CellChangeTimes->{{3.5804618141918354`*^9, 3.580461913283209*^9}, 
   3.5804623255451336`*^9, {3.5804628735428963`*^9, 3.5804628958041353`*^9}, 
   3.580476850943878*^9, {3.580477105068324*^9, 3.5804771083443303`*^9}, {
   3.5804771721796427`*^9, 3.5804771783728533`*^9}, {3.5804772740010214`*^9, 
   3.5804772770274267`*^9}, {3.580482892684083*^9, 3.5804828962564898`*^9}, {
   3.580482954194991*^9, 3.5804829592026*^9}, {3.5804831929846106`*^9, 
   3.580483196182616*^9}, {3.5804872622642126`*^9, 3.5804872631222143`*^9}, {
   3.5804873217939167`*^9, 3.580487322527118*^9}, {3.5804873526819715`*^9, 
   3.5804873772676144`*^9}, {3.580487884424505*^9, 3.5804878847053056`*^9}, 
   3.5804890446921434`*^9, {3.580622257965257*^9, 3.580622261210063*^9}, {
   3.5806223411134033`*^9, 3.580622378725069*^9}, {3.5806224544944024`*^9, 
   3.5806224929172697`*^9}, {3.5806225414489555`*^9, 
   3.5806225448029613`*^9}, {3.5806230023205643`*^9, 3.580623005799371*^9}, {
   3.580623383101633*^9, 3.580623389170044*^9}, {3.580623584887988*^9, 
   3.580623618568447*^9}, {3.581919849114888*^9, 3.5819198617977104`*^9}, {
   3.5819213195266705`*^9, 3.581921348761122*^9}, {3.5819213915675974`*^9, 
   3.5819213966532063`*^9}, {3.581921532482645*^9, 3.5819215357118506`*^9}, {
   3.581921576911523*^9, 3.581921610919583*^9}, {3.581931513876068*^9, 
   3.5819315452789235`*^9}, {3.5819315930930076`*^9, 3.581931684431168*^9}, {
   3.5819318585118732`*^9, 3.581931861210678*^9}, {3.5819318914747314`*^9, 
   3.5819318939551363`*^9}, {3.581932919235737*^9, 3.581932948345388*^9}, 
   3.5819504342397337`*^9}],

Cell[BoxData["22.448907086329864`"], "Output",
 CellChangeTimes->{
  3.58193264216365*^9, 3.5819327494918385`*^9, 3.5819327951999187`*^9, {
   3.5819329251325474`*^9, 3.5819329531813965`*^9}, 3.581950501023451*^9, 
   3.5819972635265656`*^9}],

Cell[BoxData["30.942336004961277`"], "Output",
 CellChangeTimes->{
  3.58193264216365*^9, 3.5819327494918385`*^9, 3.5819327951999187`*^9, {
   3.5819329251325474`*^9, 3.5819329531813965`*^9}, 3.581950501023451*^9, 
   3.5819972635265656`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{"all", " ", "nuclei", " ", "with", " ", "all", " ", "clouds"}], 
    ",", " ", 
    RowBox[{
     RowBox[{"except", " ", "i"}], " ", "=", " ", "j"}]}], " ", "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"Vne", "=", 
    RowBox[{"cnv", ".", 
     RowBox[{"(", 
      RowBox[{"ow", ".", "dcw"}], ")"}]}]}], "\[IndentingNewLine]", 
   RowBox[{"For", "[", 
    RowBox[{
     RowBox[{"i", "=", "1"}], ",", 
     RowBox[{"i", "\[LessEqual]", "bc"}], ",", 
     RowBox[{"i", "++"}], ",", 
     RowBox[{"Vne", "=", 
      RowBox[{"Vne", "-", 
       RowBox[{
        StyleBox["3.0",
         FontColor->RGBColor[1, 0, 0]], "*", 
        RowBox[{
         RowBox[{"cn", "[", 
          RowBox[{"[", "i", "]"}], "]"}], "/", 
         RowBox[{"u2", "[", 
          RowBox[{"[", "i", "]"}], "]"}]}]}]}]}]}], "]"}], 
   "\[IndentingNewLine]", "Vne"}]}]], "Input",
 CellChangeTimes->{{3.5804618141918354`*^9, 3.580461913283209*^9}, {
   3.5804623255451336`*^9, 3.5804623320503445`*^9}, 3.5804623790064273`*^9, 
   3.5804648933564434`*^9, {3.580486684704398*^9, 3.5804866854531994`*^9}, 
   3.5804867335324836`*^9, {3.580486767821344*^9, 3.5804868108774195`*^9}, {
   3.5804868414066734`*^9, 3.5804868626539106`*^9}, {3.580487002024555*^9, 
   3.5804870027109566`*^9}, 3.580487037218217*^9, {3.580487070789476*^9, 
   3.5804870718814783`*^9}, {3.5804874240520964`*^9, 3.580487452054146*^9}, {
   3.580487501646633*^9, 3.5804875425967045`*^9}, {3.580487596245199*^9, 
   3.5804878336464157`*^9}, 3.5804879072785454`*^9, 3.5806228673959274`*^9, {
   3.580622903821992*^9, 3.580622904133992*^9}, {3.5806229371124496`*^9, 
   3.580622962181694*^9}, {3.580623418685296*^9, 3.580623468745784*^9}, {
   3.5806235308650928`*^9, 3.5806235316294947`*^9}, {3.5806236760545483`*^9, 
   3.580623676772149*^9}, 3.5806237132606134`*^9, {3.5806237775483265`*^9, 
   3.580623778655928*^9}, 3.5806238112911854`*^9, 3.5806238549244623`*^9, {
   3.581919900282978*^9, 3.581919900984979*^9}, {3.5819203545557756`*^9, 
   3.581920375444213*^9}, {3.581920506156842*^9, 3.581920507732445*^9}, {
   3.581920542536106*^9, 3.5819206207234435`*^9}, {3.58192136462635*^9, 
   3.5819213648135505`*^9}}],

Cell[BoxData[
 RowBox[{"-", "50.75237693226853`"}]], "Output",
 CellChangeTimes->{
  3.58193264217925*^9, 3.5819327495074387`*^9, 3.581932795215519*^9, {
   3.581932925148147*^9, 3.5819329531813965`*^9}, 3.581950501039051*^9, 
   3.5819972635421658`*^9}],

Cell[BoxData[
 RowBox[{"-", "148.3269517859067`"}]], "Output",
 CellChangeTimes->{
  3.58193264217925*^9, 3.5819327495074387`*^9, 3.581932795215519*^9, {
   3.581932925148147*^9, 3.5819329531813965`*^9}, 3.581950501039051*^9, 
   3.5819972635421658`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"For", "[", 
  RowBox[{
   RowBox[{"i", "=", 
    RowBox[{"cw", "+", "1"}]}], ",", 
   RowBox[{"i", "\[LessEqual]", "nb"}], ",", 
   RowBox[{"i", "++"}], ",", 
   RowBox[{"Vne", "=", 
    RowBox[{"Vne", "-", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"3", "-", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"u2", "[", 
            RowBox[{"[", 
             RowBox[{"a", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "1"}], "]"}], "]"}], "]"}], "]"}], "/", 
           RowBox[{"a", "[", 
            RowBox[{"[", 
             RowBox[{"i", ",", "3"}], "]"}], "]"}]}], ")"}], "^", "2"}]}], 
       ")"}], "/", 
      RowBox[{"a", "[", 
       RowBox[{"[", 
        RowBox[{"i", ",", "3"}], "]"}], "]"}]}]}]}]}], 
  "]"}], "\[IndentingNewLine]", "Vne"}], "Input",
 CellChangeTimes->{{3.580461986228937*^9, 3.5804620091453776`*^9}, {
   3.5804624463517456`*^9, 3.580462451655755*^9}, 3.5804648450587587`*^9, 
   3.580464883887227*^9}],

Cell[BoxData[
 RowBox[{"-", "155.5335607287605`"}]], "Output",
 CellChangeTimes->{
  3.58193264219485*^9, 3.581932749523039*^9, 3.581932795231119*^9, {
   3.5819329251637473`*^9, 3.5819329531969967`*^9}, 3.5819505010546513`*^9, 
   3.5819972635577655`*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"(*", " ", 
  RowBox[{"nucleus", "-", "nucleus"}], " ", "*)"}]], "Input",
 CellChangeTimes->{{3.580461986228937*^9, 3.5804620091453776`*^9}, {
  3.5804624463517456`*^9, 3.580462451655755*^9}, {3.5804648450587587`*^9, 
  3.5804648790824184`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"Vnn", "=", "0.0"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"For", "[", 
  RowBox[{
   RowBox[{"i", "=", "1"}], ",", 
   RowBox[{"i", "<", "n"}], ",", 
   RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
   RowBox[{"For", "[", 
    RowBox[{
     RowBox[{"j", "=", 
      RowBox[{"i", "+", "1"}]}], ",", 
     RowBox[{"j", "<", 
      RowBox[{"n", "+", "1"}]}], ",", 
     RowBox[{"j", "++"}], ",", 
     RowBox[{"Vnn", "=", 
      RowBox[{"Vnn", "+", 
       RowBox[{
        RowBox[{"cn", "[", 
         RowBox[{"[", "i", "]"}], "]"}], "*", 
        RowBox[{
         RowBox[{"cn", "[", 
          RowBox[{"[", "j", "]"}], "]"}], "/", 
         RowBox[{"dc", "[", 
          RowBox[{"[", 
           RowBox[{"i", ",", "j"}], "]"}], "]"}]}]}]}]}]}], "]"}]}], 
  "]"}], "\[IndentingNewLine]", "Vnn"}], "Input",
 CellChangeTimes->{{3.580461986228937*^9, 3.5804620091453776`*^9}, {
  3.5804624463517456`*^9, 3.580462451655755*^9}, {3.5804648450587587`*^9, 
  3.580464866945597*^9}}],

Cell[BoxData["11.901463515165869`"], "Output",
 CellChangeTimes->{
  3.58193264221045*^9, 3.581932749538639*^9, 3.581932795246719*^9, {
   3.581932925179347*^9, 3.5819329532125964`*^9}, 3.5819505010702515`*^9, 
   3.5819972635733657`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"ratio", "=", 
  RowBox[{
   RowBox[{"-", 
    RowBox[{"(", 
     RowBox[{"Vne", "+", 
      RowBox[{"2", "*", "Vnn"}]}], ")"}]}], "/", 
   RowBox[{"(", 
    RowBox[{"Vee", "-", "Vnn"}], ")"}]}]}]], "Input",
 CellChangeTimes->{{3.580461986228937*^9, 3.5804620091453776`*^9}, {
  3.5804624463517456`*^9, 3.580462451655755*^9}, {3.5804648450587587`*^9, 
  3.580464857866381*^9}}],

Cell[BoxData["6.9183086945741215`"], "Output",
 CellChangeTimes->{
  3.58193264221045*^9, 3.581932749538639*^9, 3.581932795246719*^9, {
   3.5819329251949472`*^9, 3.581932953228196*^9}, 3.5819505010702515`*^9, 
   3.5819972635889654`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"ratiomp2", "=", 
  RowBox[{
   RowBox[{"-", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"-", "155.6621124"}], "+", 
      RowBox[{"2", "*", "11.9014636"}]}], ")"}]}], "/", 
   RowBox[{"(", 
    RowBox[{"31.1706102", "-", "11.9014636"}], ")"}]}]}]], "Input",
 CellChangeTimes->{{3.580461986228937*^9, 3.5804620091453776`*^9}, {
  3.5804624463517456`*^9, 3.580462451655755*^9}, {3.5804648450587587`*^9, 
  3.580464857866381*^9}, {3.5804657856624107`*^9, 3.580465868716957*^9}, {
  3.58052513033449*^9, 3.580525195386604*^9}, {3.580526224146565*^9, 
  3.5805263240959406`*^9}, {3.580533563696274*^9, 3.580533628046387*^9}, {
  3.5805336590436416`*^9, 3.580533659995243*^9}, {3.5806213276576233`*^9, 
  3.5806213282504244`*^9}, {3.5819204074086685`*^9, 3.581920480401197*^9}}],

Cell[BoxData["6.843021537860945`"], "Output",
 CellChangeTimes->{
  3.5819326422260504`*^9, 3.581932749554239*^9, 3.581932795262319*^9, {
   3.5819329252105474`*^9, 3.5819329532437963`*^9}, 3.581950501085851*^9, 
   3.5819972635889654`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{"Vee", "\[IndentingNewLine]", "Vnn", "\[IndentingNewLine]", \
"Vne", "\[IndentingNewLine]", 
 RowBox[{"Tlp", " ", "=", 
  RowBox[{
   StyleBox["1.59",
    FontColor->RGBColor[1, 0, 0]], "*", 
   RowBox[{"2.25", "/", 
    RowBox[{"rln", "^", "2"}]}]}]}]}], "Input",
 CellChangeTimes->{{3.5804620346046224`*^9, 3.5804620359462247`*^9}, {
   3.580462097113932*^9, 3.5804620981435337`*^9}, 3.580462151105627*^9, {
   3.5804763243493533`*^9, 3.580476329840563*^9}, {3.580476444937565*^9, 
   3.580476458634389*^9}, {3.58047649332885*^9, 3.580476496043255*^9}, 
   3.5804767076888266`*^9, {3.58047745555414*^9, 3.5804774557413406`*^9}, 
   3.580483254167918*^9, {3.5804834637194357`*^9, 3.580483464452637*^9}, {
   3.580486296076515*^9, 3.5804864867556505`*^9}, {3.5804868827935457`*^9, 
   3.5804869648652897`*^9}, {3.580487096810322*^9, 3.5804871692724495`*^9}, {
   3.5804886599642673`*^9, 3.5804886611810694`*^9}, {3.580488733908397*^9, 
   3.5804887603348436`*^9}, {3.5804887932197013`*^9, 
   3.5804888678658323`*^9}, {3.580488905711499*^9, 3.5804889063667*^9}, 
   3.5806222929561186`*^9, {3.580622404231114*^9, 3.580622425166351*^9}, 
   3.580622513587306*^9, {3.5806225751606145`*^9, 3.5806225758626156`*^9}, {
   3.580623056265459*^9, 3.5806230569674606`*^9}, 3.580623087886715*^9, {
   3.580623139132805*^9, 3.5806231663236523`*^9}, {3.580623197944908*^9, 
   3.580623217819343*^9}, {3.5806233326979446`*^9, 3.5806233335871463`*^9}, {
   3.58191872406435*^9, 3.581918738197975*^9}, 3.5819196404645214`*^9, 
   3.581921229935714*^9, {3.5819212715253863`*^9, 3.5819212726641884`*^9}, {
   3.581921419866047*^9, 3.581921466556929*^9}, {3.5819215532618814`*^9, 
   3.5819215537766824`*^9}, 3.5819216479072475`*^9, 3.5819216983421364`*^9, {
   3.5819313972034636`*^9, 3.581931417249499*^9}, {3.5819314782144055`*^9, 
   3.581931480164409*^9}, {3.581931919648381*^9, 3.581932022015761*^9}, 
   3.581950576667984*^9}],

Cell[BoxData["30.942336004961277`"], "Output",
 CellChangeTimes->{
  3.58193264224165*^9, 3.5819327495698385`*^9, 3.581932795277919*^9, {
   3.5819329252105474`*^9, 3.5819329532437963`*^9}, 3.581950501101451*^9, 
   3.5819972636045656`*^9}],

Cell[BoxData["11.901463515165869`"], "Output",
 CellChangeTimes->{
  3.58193264224165*^9, 3.5819327495698385`*^9, 3.581932795277919*^9, {
   3.5819329252105474`*^9, 3.5819329532437963`*^9}, 3.581950501101451*^9, 
   3.5819972636045656`*^9}],

Cell[BoxData[
 RowBox[{"-", "155.5335607287605`"}]], "Output",
 CellChangeTimes->{
  3.58193264224165*^9, 3.5819327495698385`*^9, 3.581932795277919*^9, {
   3.5819329252105474`*^9, 3.5819329532437963`*^9}, 3.581950501101451*^9, 
   3.5819972636045656`*^9}],

Cell[BoxData["2.686383583286527`"], "Output",
 CellChangeTimes->{
  3.58193264224165*^9, 3.5819327495698385`*^9, 3.581932795277919*^9, {
   3.5819329252105474`*^9, 3.5819329532437963`*^9}, 3.581950501101451*^9, 
   3.5819972636045656`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"Vtot", "=", 
  RowBox[{"Vne", "+", "Vee", "+", "Vnn"}]}], "\[IndentingNewLine]", 
 RowBox[{"Ttot", " ", "=", 
  RowBox[{"Ttot", "+", "Tlp"}]}], "\[IndentingNewLine]", 
 RowBox[{"Vir", "   ", "=", 
  RowBox[{
   RowBox[{"-", "Vtot"}], "/", "Ttot"}]}], "\[IndentingNewLine]", 
 RowBox[{"Etot", "=", 
  RowBox[{"Ttot", "+", "Vtot"}]}]}], "Input",
 CellChangeTimes->{
  3.580462058613065*^9, {3.5804621657228527`*^9, 3.580462173351266*^9}}],

Cell[BoxData[
 RowBox[{"-", "112.68976120863336`"}]], "Output",
 CellChangeTimes->{
  3.5819326422572503`*^9, 3.5819327495854387`*^9, 3.581932795293519*^9, {
   3.581932925226147*^9, 3.5819329532593966`*^9}, 3.581950501117051*^9, 
   3.581997263620166*^9}],

Cell[BoxData["56.342363282528154`"], "Output",
 CellChangeTimes->{
  3.5819326422572503`*^9, 3.5819327495854387`*^9, 3.581932795293519*^9, {
   3.581932925226147*^9, 3.5819329532593966`*^9}, 3.581950501117051*^9, 
   3.581997263620166*^9}],

Cell[BoxData["2.000089358047546`"], "Output",
 CellChangeTimes->{
  3.5819326422572503`*^9, 3.5819327495854387`*^9, 3.581932795293519*^9, {
   3.581932925226147*^9, 3.5819329532593966`*^9}, 3.581950501117051*^9, 
   3.581997263620166*^9}],

Cell[BoxData[
 RowBox[{"-", "56.3473979261052`"}]], "Output",
 CellChangeTimes->{
  3.5819326422572503`*^9, 3.5819327495854387`*^9, 3.581932795293519*^9, {
   3.581932925226147*^9, 3.5819329532593966`*^9}, 3.581950501117051*^9, 
   3.581997263620166*^9}]
}, Open  ]]
},
WindowSize->{850, 881},
WindowMargins->{{Automatic, 415}, {11, Automatic}},
PrintingCopies->1,
PrintingPageRange->{32000, 32000},
PrintingOptions->{"Magnification"->1.,
"PaperOrientation"->"Portrait",
"PaperSize"->{595.2, 841.8}},
Magnification->0.8999999761581421,
FrontEndVersion->"9.0 for Microsoft Windows (64-bit) (January 25, 2013)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[557, 20, 1250, 26, 50, "Input"],
Cell[1810, 48, 1519, 31, 50, "Input"],
Cell[CellGroupData[{
Cell[3354, 83, 1077, 16, 29, "Input"],
Cell[4434, 101, 321, 7, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[4792, 113, 5774, 156, 350, "Input"],
Cell[10569, 271, 223, 4, 28, "Output"],
Cell[10795, 277, 307, 7, 28, "Output"],
Cell[11105, 286, 307, 7, 28, "Output"],
Cell[11415, 295, 296, 6, 28, "Output"],
Cell[11714, 303, 225, 4, 28, "Output"],
Cell[11942, 309, 225, 4, 28, "Output"]
}, Open  ]],
Cell[12182, 316, 5105, 104, 190, "Input"],
Cell[17290, 422, 3790, 74, 70, "Input"],
Cell[21083, 498, 1254, 35, 70, "Input"],
Cell[CellGroupData[{
Cell[22362, 537, 3042, 63, 70, "Input"],
Cell[25407, 602, 1004, 25, 97, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[26448, 632, 3060, 74, 130, "Input"],
Cell[29511, 708, 316, 7, 28, "Output"]
}, Open  ]],
Cell[29842, 718, 1708, 43, 90, "Input"],
Cell[31553, 763, 511, 11, 50, "Input"],
Cell[32067, 776, 2048, 48, 90, "Input"],
Cell[34118, 826, 2106, 50, 90, "Input"],
Cell[36227, 878, 1261, 33, 50, "Input"],
Cell[CellGroupData[{
Cell[37513, 915, 2727, 71, 130, "Input"],
Cell[40243, 988, 234, 4, 28, "Output"],
Cell[40480, 994, 234, 4, 28, "Output"],
Cell[40717, 1000, 234, 4, 28, "Output"],
Cell[40954, 1006, 234, 4, 28, "Output"],
Cell[41191, 1012, 219, 4, 28, "Output"]
}, Open  ]],
Cell[41425, 1019, 1971, 35, 90, "Input"],
Cell[43399, 1056, 2207, 39, 29, "Input"],
Cell[CellGroupData[{
Cell[45631, 1099, 1371, 43, 150, "Input"],
Cell[47005, 1144, 223, 4, 28, "Output"],
Cell[47231, 1150, 223, 4, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[47491, 1159, 883, 26, 70, "Input"],
Cell[48377, 1187, 238, 5, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[48652, 1197, 3081, 94, 210, "Input"],
Cell[51736, 1293, 242, 5, 28, "Output"],
Cell[51981, 1300, 223, 4, 28, "Output"],
Cell[52207, 1306, 1119, 33, 82, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[53363, 1344, 2861, 78, 210, "Input"],
Cell[56227, 1424, 219, 4, 28, "Output"],
Cell[56449, 1430, 219, 4, 28, "Output"]
}, Open  ]],
Cell[56683, 1437, 3682, 100, 170, "Input"],
Cell[60368, 1539, 1153, 17, 29, "Input"],
Cell[61524, 1558, 2724, 54, 90, "Input"],
Cell[CellGroupData[{
Cell[64273, 1616, 1065, 19, 70, "Input"],
Cell[65341, 1637, 510, 13, 28, "Output"],
Cell[65854, 1652, 221, 4, 28, "Output"],
Cell[66078, 1658, 221, 4, 28, "Output"]
}, Open  ]],
Cell[66314, 1665, 11150, 255, 210, "Input"],
Cell[77467, 1922, 9944, 229, 190, "Input"],
Cell[87414, 2153, 8221, 161, 150, "Input"],
Cell[CellGroupData[{
Cell[95660, 2318, 5239, 75, 50, "Input"],
Cell[100902, 2395, 2674, 60, 337, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[103613, 2460, 2266, 65, 110, "Input"],
Cell[105882, 2527, 220, 4, 28, "Output"]
}, Open  ]],
Cell[106117, 2534, 1001, 26, 70, "Input"],
Cell[CellGroupData[{
Cell[107143, 2564, 1898, 54, 170, "Input"],
Cell[109044, 2620, 223, 4, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[109304, 2629, 1422, 43, 90, "Input"],
Cell[110729, 2674, 402, 11, 28, "Output"],
Cell[111134, 2687, 308, 7, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[111479, 2699, 1082, 31, 70, "Input"],
Cell[112564, 2732, 218, 4, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[112819, 2741, 711, 21, 90, "Input"],
Cell[113533, 2764, 243, 5, 28, "Output"],
Cell[113779, 2771, 223, 4, 28, "Output"],
Cell[114005, 2777, 226, 4, 28, "Output"],
Cell[114234, 2783, 223, 4, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[114494, 2792, 3815, 92, 230, "Input"],
Cell[118312, 2886, 242, 4, 28, "Output"],
Cell[118557, 2892, 242, 4, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[118836, 2901, 2238, 47, 90, "Input"],
Cell[121077, 2950, 254, 5, 28, "Output"],
Cell[121334, 2957, 254, 5, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[121625, 2967, 1016, 30, 50, "Input"],
Cell[122644, 2999, 256, 5, 28, "Output"]
}, Open  ]],
Cell[122915, 3007, 268, 5, 29, "Input"],
Cell[CellGroupData[{
Cell[123208, 3016, 1016, 29, 90, "Input"],
Cell[124227, 3047, 238, 4, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[124502, 3056, 400, 11, 29, "Input"],
Cell[124905, 3069, 238, 4, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[125180, 3078, 799, 16, 29, "Input"],
Cell[125982, 3096, 240, 4, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[126259, 3105, 1929, 30, 90, "Input"],
Cell[128191, 3137, 240, 4, 28, "Output"],
Cell[128434, 3143, 240, 4, 28, "Output"],
Cell[128677, 3149, 256, 5, 28, "Output"],
Cell[128936, 3156, 239, 4, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[129212, 3165, 461, 11, 90, "Input"],
Cell[129676, 3178, 256, 5, 28, "Output"],
Cell[129935, 3185, 239, 4, 28, "Output"],
Cell[130177, 3191, 238, 4, 28, "Output"],
Cell[130418, 3197, 254, 5, 60, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)