(* 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[     80774,       1872]
NotebookOptionsPosition[     77170,       1756]
NotebookOutlinePosition[     77681,       1775]
CellTagsIndexPosition[     77638,       1772]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell["\<\
(* CH4   24.05.2013; all with matrix operations; e.g. -4*(3-(p/R2)^2)/R2 for \
H atoms is included in cn.(ow.dwnr), see below *)\
\>", "Text",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5783693559648447`*^9, 3.5783693662920628`*^9}, 
   3.578369554833994*^9, {3.5783729062258806`*^9, 3.5783729062258806`*^9}, {
   3.578414671505531*^9, 3.5784147130484037`*^9}, {3.578415068885029*^9, 
   3.5784152004400597`*^9}, 3.578415236117323*^9, 3.5784152716541853`*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"Clear", "[", 
   RowBox[{
   "z", ",", "sig1", ",", "sig2", ",", "sig4", ",", "k1", ",", "k2", ",", 
    "k4", ",", "nc", ",", "R1", ",", "R2", ",", "R3", ",", "R4", ",", "w", 
    ",", "p", ",", "vee", ",", "vne", ",", "vnn", ",", "n", ",", "c", ",", 
    "dw", ",", "dn", ",", "dwn", ",", "dwnr", ",", "dwr", ",", "dnr", ",", 
    "a1", ",", "a2", ",", "a3", ",", "b1", ",", "b2", ",", "b3", ",", "d1", 
    ",", "d2", ",", "d3", ",", "d4", ",", "d5", ",", "pi", ",", "i", ",", "j",
     ",", "d", ",", "t"}], "]"}], ";"}]], "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5783693559648447`*^9, 3.5783693662920628`*^9}, 
   3.578369554833994*^9, 3.578370271764453*^9, {3.5783729062258806`*^9, 
   3.5783729062258806`*^9}, {3.5783896086269007`*^9, 3.5783896345697465`*^9}, 
   3.5783944253753614`*^9, 3.57840293976919*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"R1", "=", "0.26537992"}], ";", 
  RowBox[{"R2", "=", "1.2343054"}], ";", 
  RowBox[{"p", "=", "0.5334356688274884"}], ";"}]], "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5783693559648447`*^9, 3.5783693662920628`*^9}, 
   3.578369554833994*^9, 3.578370271764453*^9, {3.578370376315837*^9, 
   3.578370376315837*^9}, {3.5783729062258806`*^9, 3.5783729062258806`*^9}, {
   3.5783780090233684`*^9, 3.578378019537787*^9}, {3.5783786664709234`*^9, 
   3.5783786985289793`*^9}, {3.57839148096379*^9, 3.5783914901678057`*^9}, {
   3.578391539635493*^9, 3.5783915617251315`*^9}, 3.578408000583571*^9, {
   3.578408075073702*^9, 3.5784080792389092`*^9}, {3.578409085674677*^9, 
   3.578409089059883*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"z", "=", "6.0"}], ";", "  ", 
  RowBox[{"vne", "=", "0.0"}], ";"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"nc", "=", "5"}], ";", " ", 
   RowBox[{"pi", "=", "0.0"}], ";", " ", 
   RowBox[{"bohr", "=", "0.529177"}], ";", 
   RowBox[{"sq", "=", 
    RowBox[{"Sqrt", "[", "3", "]"}]}], ";"}], "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{"k", "=", 
    RowBox[{"Table", "[", 
     RowBox[{"1.0", ",", 
      RowBox[{"{", "nc", "}"}]}], "]"}]}], "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"k", "=", 
    RowBox[{"{", 
     RowBox[{
     "1.04605075", ",", "1.20889711", ",", "1.20889711", ",", "1.20889711", 
      ",", "1.20889711"}], "}"}]}], ";"}], " ", "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{"sig", "=", 
    RowBox[{"Table", "[", 
     RowBox[{"0.4", ",", 
      RowBox[{"{", "nc", "}"}]}], "]"}]}], "*)"}]}], "\n", 
 RowBox[{
  RowBox[{"sig", "=", 
   RowBox[{"{", 
    RowBox[{
    "0.2650207", ",", "0.30451569", ",", "0.30451569", ",", "0.30451569", ",",
      "0.30451569"}], "}"}]}], ";"}]}], "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5783693559648447`*^9, 3.5783693662920628`*^9}, 
   3.578369554833994*^9, 3.578370271764453*^9, {3.5783703674238214`*^9, 
   3.5783703954726706`*^9}, {3.578371324750703*^9, 3.5783713343291197`*^9}, 
   3.5783717590710654`*^9, {3.5783729062258806`*^9, 3.5783729062258806`*^9}, 
   3.5783735444386015`*^9, {3.5783737801394157`*^9, 3.578373791355835*^9}, {
   3.578401369954033*^9, 3.578401375242442*^9}, {3.578402976710055*^9, 
   3.578403050903785*^9}, {3.57840781804805*^9, 3.5784079245802374`*^9}, {
   3.5784080374308357`*^9, 3.5784080457144504`*^9}, {3.5784081607802525`*^9, 
   3.5784081643058586`*^9}, {3.578409120493938*^9, 3.578409156998002*^9}}],

Cell[BoxData[
 RowBox[{"(*", " ", 
  RowBox[{
   RowBox[{"C", " ", "He", " ", "shell"}], ",", " ", 
   RowBox[{"i", "=", "j"}]}], " ", "*)"}]], "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5783693559648447`*^9, 3.5783693662920628`*^9}, {
   3.578369554833994*^9, 3.5783695604344034`*^9}, {3.5783729062258806`*^9, 
   3.5783729062258806`*^9}, 3.578373808328665*^9, {3.5783944107425356`*^9, 
   3.57839441348814*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"ow", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "2"}], ",", 
     RowBox[{"-", "2"}], ",", 
     RowBox[{"-", "2"}], ",", 
     RowBox[{"-", "2"}], ",", 
     RowBox[{"-", "2"}]}], "}"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"cn", "=", 
   RowBox[{"{", 
    RowBox[{"z", ",", "1", ",", "1", ",", "1", ",", "1"}], "}"}]}], 
  ";"}], "\n", 
 RowBox[{
  RowBox[{"c", "=", 
   RowBox[{"{", 
    RowBox[{"R1", ",", "R2", ",", "R2", ",", "R2", ",", "R2"}], "}"}]}], 
  ";"}], "\n", 
 RowBox[{"Ekin", " ", "=", " ", 
  RowBox[{
   RowBox[{"-", 
    RowBox[{"(", 
     RowBox[{"9", "/", "8"}], ")"}]}], "*", 
   RowBox[{"Sum", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"k", "[", 
       RowBox[{"[", "i", "]"}], "]"}], "*", 
      RowBox[{
       RowBox[{"ow", "[", 
        RowBox[{"[", "i", "]"}], "]"}], "/", 
       RowBox[{
        RowBox[{"c", "[", 
         RowBox[{"[", "i", "]"}], "]"}], "^", "2"}]}]}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "nc"}], "}"}]}], "]"}]}]}], "\n", 
 RowBox[{"vee", "   ", "=", " ", 
  RowBox[{
   RowBox[{"-", 
    RowBox[{"(", 
     RowBox[{"3", "/", "2"}], ")"}]}], "*", 
   RowBox[{"Sum", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"sig", "[", 
       RowBox[{"[", "i", "]"}], "]"}], "*", 
      RowBox[{
       RowBox[{"ow", "[", 
        RowBox[{"[", "i", "]"}], "]"}], "/", 
       RowBox[{"c", "[", 
        RowBox[{"[", "i", "]"}], "]"}]}]}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "nc"}], "}"}]}], "]"}]}]}]}], "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5783693559648447`*^9, 3.5783693662920628`*^9}, {
   3.578369554833994*^9, 3.5783695604344034`*^9}, 3.5783701337978106`*^9, {
   3.578371289853441*^9, 3.578371302442663*^9}, {3.57837134613834*^9, 
   3.57837135733916*^9}, {3.578371392314421*^9, 3.5783716741757164`*^9}, {
   3.578371712130583*^9, 3.578371724501405*^9}, {3.5783726153945694`*^9, 
   3.5783726153945694`*^9}, 3.5783726553930397`*^9, 3.578372749586005*^9, {
   3.5783728157925215`*^9, 3.5783728573041945`*^9}, {3.578373248116081*^9, 
   3.5783732528896894`*^9}, {3.5783732892065535`*^9, 
   3.5783733600150776`*^9}, {3.578373605637509*^9, 3.5783736465095806`*^9}, 
   3.578373878918789*^9, {3.5783778323218584`*^9, 3.5783779683072968`*^9}, {
   3.5783782452545834`*^9, 3.5783782551606007`*^9}, {3.5783783031618853`*^9, 
   3.5783783174047103`*^9}, {3.5783792814240036`*^9, 
   3.5783793528877287`*^9}, {3.578379408611027*^9, 3.5783794493738985`*^9}, {
   3.578379483085558*^9, 3.578379496922782*^9}, {3.5783795638312993`*^9, 
   3.5783797837136855`*^9}, {3.578380015670493*^9, 3.578380570345067*^9}, {
   3.5783806138535433`*^9, 3.5783806149455457`*^9}, {3.578380675645252*^9, 
   3.5783807269069424`*^9}, 3.5783807789174337`*^9, 3.5783894188369675`*^9, {
   3.5783894739674644`*^9, 3.5783894894270916`*^9}, {3.578389884825386*^9, 
   3.5783899482706976`*^9}, {3.5783916231892395`*^9, 3.5783916258256445`*^9}, 
   3.578402886229896*^9, {3.578403430655252*^9, 3.5784034604825044`*^9}, 
   3.5784035095133905`*^9, {3.578403576437508*^9, 3.578403610585968*^9}, {
   3.578403655342447*^9, 3.5784036911445093`*^9}, {3.578403738755793*^9, 
   3.5784037505962143`*^9}, {3.578403806834313*^9, 3.578403809704718*^9}, {
   3.5784039070800886`*^9, 3.578403917859708*^9}, 3.5784081931347094`*^9, {
   3.5784082628512316`*^9, 3.578408283412068*^9}, {3.578408320524533*^9, 
   3.5784083687130175`*^9}, {3.578408401691476*^9, 3.5784084583975754`*^9}, {
   3.578408505088457*^9, 3.578408509612465*^9}, {3.578408598517021*^9, 
   3.578408635988287*^9}, {3.578408829693827*^9, 3.5784088678670945`*^9}, {
   3.578408960640457*^9, 3.5784090167225556`*^9}, {3.578409331016308*^9, 
   3.5784093564755526`*^9}, {3.578409544440283*^9, 3.578409547076687*^9}, {
   3.5784143031576843`*^9, 3.5784143198653135`*^9}}],

Cell[BoxData["40.56089216702621`"], "Output",
 CellChangeTimes->{3.5784143359489417`*^9, 3.5784155274478345`*^9, 
  3.578415560863093*^9, 3.5784156294876137`*^9, 3.5784156751332936`*^9, 
  3.5784157213405747`*^9, 3.57841583231917*^9, 3.5784171204942327`*^9, 
  3.5784172813773155`*^9, 3.57841733550941*^9, 3.578417735353712*^9, 
  3.578418066947895*^9, 3.5784181925125155`*^9, 3.578418290839488*^9, 
  3.5784183398391743`*^9}],

Cell[BoxData["5.956461171075636`"], "Output",
 CellChangeTimes->{3.5784143359489417`*^9, 3.5784155274478345`*^9, 
  3.578415560863093*^9, 3.5784156294876137`*^9, 3.5784156751332936`*^9, 
  3.5784157213405747`*^9, 3.57841583231917*^9, 3.5784171204942327`*^9, 
  3.5784172813773155`*^9, 3.57841733550941*^9, 3.578417735353712*^9, 
  3.578418066947895*^9, 3.5784181925125155`*^9, 3.578418290839488*^9, 
  3.5784183398391743`*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"d2", "=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"R1", "+", "R2", "+", "p"}], ")"}], "/", "sq"}]}], ";", "  ", 
  RowBox[{"d4", "=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"R1", "+", "R2"}], ")"}], "/", "sq"}]}], ";"}]], "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5783693559648447`*^9, 3.5783693662920628`*^9}, {
   3.578369554833994*^9, 3.5783695604344034`*^9}, 3.5783701337978106`*^9, {
   3.578371289853441*^9, 3.578371302442663*^9}, {3.57837134613834*^9, 
   3.57837135733916*^9}, {3.578371392314421*^9, 3.5783716741757164`*^9}, {
   3.578371712130583*^9, 3.578371724501405*^9}, {3.5783726153945694`*^9, 
   3.5783726153945694`*^9}, 3.5783726553930397`*^9, 3.578372749586005*^9, {
   3.5783728157925215`*^9, 3.5783728573041945`*^9}, {3.578373248116081*^9, 
   3.5783732528896894`*^9}, {3.5783732892065535`*^9, 
   3.5783733600150776`*^9}, {3.578373605637509*^9, 3.5783736465095806`*^9}, 
   3.578373878918789*^9, {3.5783778323218584`*^9, 3.5783779683072968`*^9}, {
   3.5783782452545834`*^9, 3.5783782551606007`*^9}, {3.5783783031618853`*^9, 
   3.5783783174047103`*^9}, {3.5783792814240036`*^9, 
   3.5783793528877287`*^9}, {3.578379408611027*^9, 3.5783794493738985`*^9}, {
   3.578379483085558*^9, 3.578379496922782*^9}, {3.5783795638312993`*^9, 
   3.5783797837136855`*^9}, {3.578380015670493*^9, 3.578380570345067*^9}, {
   3.5783806138535433`*^9, 3.5783806149455457`*^9}, {3.578380675645252*^9, 
   3.5783807269069424`*^9}, 3.5783807789174337`*^9, 3.5783894188369675`*^9, {
   3.5783894739674644`*^9, 3.5783894894270916`*^9}, {3.578389884825386*^9, 
   3.5783899482706976`*^9}, {3.5783916231892395`*^9, 3.5783916258256445`*^9}, 
   3.578402886229896*^9, {3.578403430655252*^9, 3.5784034604825044`*^9}, 
   3.5784035095133905`*^9, {3.578403576437508*^9, 3.578403610585968*^9}, {
   3.578403655342447*^9, 3.5784036911445093`*^9}, {3.578403738755793*^9, 
   3.5784037505962143`*^9}, {3.578403806834313*^9, 3.578403809704718*^9}, {
   3.5784039070800886`*^9, 3.578403917859708*^9}, 3.5784081931347094`*^9, {
   3.5784082628512316`*^9, 3.578408283412068*^9}, {3.578408320524533*^9, 
   3.5784083687130175`*^9}, {3.578408401691476*^9, 3.5784084583975754`*^9}, {
   3.578408505088457*^9, 3.578408509612465*^9}, {3.578408598517021*^9, 
   3.578408635988287*^9}, {3.578408829693827*^9, 3.5784088678670945`*^9}, {
   3.578408960640457*^9, 3.5784090167225556`*^9}, {3.578409331016308*^9, 
   3.5784093564755526`*^9}, {3.578409544440283*^9, 3.578409547076687*^9}, {
   3.5784143031576843`*^9, 3.578414314998105*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"n", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"d2", ",", "d2", ",", "d2"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"-", "d2"}], ",", 
       RowBox[{"-", "d2"}], ",", "d2"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"d2", ",", 
       RowBox[{"-", "d2"}], ",", 
       RowBox[{"-", "d2"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"-", "d2"}], ",", "d2", ",", 
       RowBox[{"-", "d2"}]}], "}"}]}], "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"MatrixForm", "[", "n", "]"}]}], "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5783693559648447`*^9, 3.5783693662920628`*^9}, {
   3.578369554833994*^9, 3.5783695604344034`*^9}, 3.5783701337978106`*^9, {
   3.578371289853441*^9, 3.578371302442663*^9}, {3.57837134613834*^9, 
   3.57837135733916*^9}, {3.578371392314421*^9, 3.5783716741757164`*^9}, {
   3.578371712130583*^9, 3.578371724501405*^9}, {3.5783726153945694`*^9, 
   3.5783726153945694`*^9}, 3.5783726553930397`*^9, 3.578372749586005*^9, {
   3.5783728157925215`*^9, 3.5783728573041945`*^9}, {3.578373248116081*^9, 
   3.5783732528896894`*^9}, {3.5783732892065535`*^9, 
   3.5783733600150776`*^9}, {3.578373605637509*^9, 3.5783736465095806`*^9}, 
   3.578373878918789*^9, {3.5783778323218584`*^9, 3.5783779683072968`*^9}, {
   3.5783782452545834`*^9, 3.5783782551606007`*^9}, {3.5783783031618853`*^9, 
   3.5783783174047103`*^9}, {3.5783792814240036`*^9, 
   3.5783793528877287`*^9}, {3.578379408611027*^9, 3.5783794493738985`*^9}, {
   3.578379483085558*^9, 3.578379496922782*^9}, {3.5783795638312993`*^9, 
   3.5783797837136855`*^9}, {3.578380015670493*^9, 3.578380570345067*^9}, {
   3.5783806138535433`*^9, 3.5783806149455457`*^9}, {3.578380675645252*^9, 
   3.5783807269069424`*^9}, 3.5783807789174337`*^9, 3.5783894188369675`*^9, {
   3.5783894739674644`*^9, 3.5783894894270916`*^9}, {3.578389884825386*^9, 
   3.5783899482706976`*^9}, 3.578390559479771*^9}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"0", "0", "0"},
     {"1.173822950194629`", "1.173822950194629`", "1.173822950194629`"},
     {
      RowBox[{"-", "1.173822950194629`"}], 
      RowBox[{"-", "1.173822950194629`"}], "1.173822950194629`"},
     {"1.173822950194629`", 
      RowBox[{"-", "1.173822950194629`"}], 
      RowBox[{"-", "1.173822950194629`"}]},
     {
      RowBox[{"-", "1.173822950194629`"}], "1.173822950194629`", 
      RowBox[{"-", "1.173822950194629`"}]}
    },
    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.5784143359801416`*^9, 3.5784155274790344`*^9, 
  3.578415560894293*^9, 3.578415629503214*^9, 3.5784156751644936`*^9, 
  3.5784157213717747`*^9, 3.5784158323347697`*^9, 3.5784171205254326`*^9, 
  3.578417281408515*^9, 3.57841733552501*^9, 3.5784177353693123`*^9, 
  3.578418066963495*^9, 3.5784181925437155`*^9, 3.5784182908706884`*^9, 
  3.5784183398703747`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"w", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"d4", ",", "d4", ",", "d4"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"-", "d4"}], ",", 
       RowBox[{"-", "d4"}], ",", "d4"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"d4", ",", 
       RowBox[{"-", "d4"}], ",", 
       RowBox[{"-", "d4"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"-", "d4"}], ",", "d4", ",", 
       RowBox[{"-", "d4"}]}], "}"}]}], "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"MatrixForm", "[", "w", "]"}]}], "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5783693559648447`*^9, 3.5783693662920628`*^9}, {
   3.578369554833994*^9, 3.5783695604344034`*^9}, 3.5783701337978106`*^9, {
   3.578371289853441*^9, 3.578371302442663*^9}, {3.57837134613834*^9, 
   3.57837135733916*^9}, {3.578371392314421*^9, 3.5783716741757164`*^9}, {
   3.578371712130583*^9, 3.578371724501405*^9}, {3.5783726153945694`*^9, 
   3.5783726153945694`*^9}, 3.5783726553930397`*^9, 3.578372749586005*^9, {
   3.5783728157925215`*^9, 3.5783728573041945`*^9}, {3.578373248116081*^9, 
   3.5783732528896894`*^9}, {3.5783732892065535`*^9, 
   3.5783733600150776`*^9}, {3.578373605637509*^9, 3.5783736465095806`*^9}, 
   3.578373878918789*^9, {3.5783778323218584`*^9, 3.5783779683072968`*^9}, {
   3.5783782452545834`*^9, 3.5783782551606007`*^9}, {3.5783783031618853`*^9, 
   3.5783783174047103`*^9}, {3.5783792814240036`*^9, 
   3.5783793528877287`*^9}, {3.578379408611027*^9, 3.5783794493738985`*^9}, {
   3.578379483085558*^9, 3.578379496922782*^9}, {3.5783795638312993`*^9, 
   3.5783797837136855`*^9}, {3.578380015670493*^9, 3.578380570345067*^9}, {
   3.5783806138535433`*^9, 3.5783806149455457`*^9}, {3.578380675645252*^9, 
   3.5783807269069424`*^9}, 3.5783807789174337`*^9, 3.5783894188369675`*^9, {
   3.5783894739674644`*^9, 3.5783894894270916`*^9}, {3.578389884825386*^9, 
   3.578389953637107*^9}, 3.578390556656166*^9}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"0", "0", "0"},
     {"0.8658437232017302`", "0.8658437232017302`", "0.8658437232017302`"},
     {
      RowBox[{"-", "0.8658437232017302`"}], 
      RowBox[{"-", "0.8658437232017302`"}], "0.8658437232017302`"},
     {"0.8658437232017302`", 
      RowBox[{"-", "0.8658437232017302`"}], 
      RowBox[{"-", "0.8658437232017302`"}]},
     {
      RowBox[{"-", "0.8658437232017302`"}], "0.8658437232017302`", 
      RowBox[{"-", "0.8658437232017302`"}]}
    },
    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.578414335995742*^9, 3.5784155274946346`*^9, 
  3.578415560894293*^9, 3.5784156295188136`*^9, 3.578415675180094*^9, 
  3.578415721387375*^9, 3.5784158323503695`*^9, 3.578417120541033*^9, 
  3.578417281424115*^9, 3.5784173355406103`*^9, 3.5784177353849125`*^9, 
  3.578418066979095*^9, 3.5784181925593157`*^9, 3.578418290886288*^9, 
  3.5784183398859744`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"dw", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"Norm", "[", 
      RowBox[{
       RowBox[{"w", "[", 
        RowBox[{"[", "i", "]"}], "]"}], "-", 
       RowBox[{"w", "[", 
        RowBox[{"[", "j", "]"}], "]"}]}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "nc"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"j", ",", "1", ",", "nc"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"MatrixForm", "[", "dw", "]"}]}], "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5783693559648447`*^9, 3.5783693662920628`*^9}, {
   3.578369554833994*^9, 3.5783695604344034`*^9}, 3.5783701337978106`*^9, {
   3.578371289853441*^9, 3.578371302442663*^9}, {3.57837134613834*^9, 
   3.57837135733916*^9}, {3.578371392314421*^9, 3.5783716741757164`*^9}, {
   3.578371712130583*^9, 3.578371724501405*^9}, {3.5783726153945694`*^9, 
   3.5783726153945694`*^9}, 3.5783726553930397`*^9, 3.578372749586005*^9, {
   3.5783728157925215`*^9, 3.5783728573041945`*^9}, {3.578373248116081*^9, 
   3.5783732528896894`*^9}, {3.5783732892065535`*^9, 
   3.5783733600150776`*^9}, {3.578373605637509*^9, 3.5783736465095806`*^9}, 
   3.578373878918789*^9, {3.5783778323218584`*^9, 3.5783779683072968`*^9}, {
   3.5783782452545834`*^9, 3.5783782551606007`*^9}, {3.5783783031618853`*^9, 
   3.5783783174047103`*^9}, {3.5783792814240036`*^9, 
   3.5783793528877287`*^9}, {3.578379408611027*^9, 3.5783794493738985`*^9}, {
   3.578379483085558*^9, 3.578379496922782*^9}, {3.5783795638312993`*^9, 
   3.5783797837136855`*^9}, {3.578380015670493*^9, 3.578380570345067*^9}, {
   3.5783806138535433`*^9, 3.5783806149455457`*^9}, {3.578380675645252*^9, 
   3.5783807269069424`*^9}, 3.5783807789174337`*^9, 3.5783894188369675`*^9, {
   3.5783894739674644`*^9, 3.5783894894270916`*^9}, {3.578389884825386*^9, 
   3.5783899611875205`*^9}, {3.5783900128548107`*^9, 
   3.5783900278308372`*^9}, {3.5783900678917074`*^9, 
   3.5783900722753153`*^9}, {3.578390104099371*^9, 3.5783901302606173`*^9}, 
   3.578390540541338*^9}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"0", "1.4996853200000002`", "1.4996853200000002`", 
      "1.4996853200000002`", "1.4996853200000002`"},
     {"1.4996853200000002`", "0.`", "2.448975872495006`", 
      "2.448975872495006`", "2.448975872495006`"},
     {"1.4996853200000002`", "2.448975872495006`", "0.`", 
      "2.448975872495006`", "2.448975872495006`"},
     {"1.4996853200000002`", "2.448975872495006`", "2.448975872495006`", 
      "0.`", "2.448975872495006`"},
     {"1.4996853200000002`", "2.448975872495006`", "2.448975872495006`", 
      "2.448975872495006`", "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.5784143360113416`*^9, 3.5784155274946346`*^9, 
  3.578415560909893*^9, 3.5784156295344133`*^9, 3.578415675195694*^9, 
  3.578415721402975*^9, 3.5784158323659697`*^9, 3.578417120556633*^9, 
  3.5784172814397154`*^9, 3.5784173355562105`*^9, 3.5784177354005127`*^9, 
  3.5784180669946947`*^9, 3.578418192574916*^9, 3.5784182909018884`*^9, 
  3.578418339901574*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"dn", " ", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"Norm", "[", 
      RowBox[{
       RowBox[{"n", "[", 
        RowBox[{"[", "i", "]"}], "]"}], "-", 
       RowBox[{"n", "[", 
        RowBox[{"[", "j", "]"}], "]"}]}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "nc"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"j", ",", "1", ",", "nc"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"MatrixForm", "[", "dn", "]"}]}], "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5783693559648447`*^9, 3.5783693662920628`*^9}, {
   3.578369554833994*^9, 3.5783695604344034`*^9}, 3.5783701337978106`*^9, {
   3.578371289853441*^9, 3.578371302442663*^9}, {3.57837134613834*^9, 
   3.57837135733916*^9}, {3.578371392314421*^9, 3.5783716741757164`*^9}, {
   3.578371712130583*^9, 3.578371724501405*^9}, {3.5783726153945694`*^9, 
   3.5783726153945694`*^9}, 3.5783726553930397`*^9, 3.578372749586005*^9, {
   3.5783728157925215`*^9, 3.5783728573041945`*^9}, {3.578373248116081*^9, 
   3.5783732528896894`*^9}, {3.5783732892065535`*^9, 
   3.5783733600150776`*^9}, {3.578373605637509*^9, 3.5783736465095806`*^9}, 
   3.578373878918789*^9, {3.5783778323218584`*^9, 3.5783779683072968`*^9}, {
   3.5783782452545834`*^9, 3.5783782551606007`*^9}, {3.5783783031618853`*^9, 
   3.5783783174047103`*^9}, {3.5783792814240036`*^9, 
   3.5783793528877287`*^9}, {3.578379408611027*^9, 3.5783794493738985`*^9}, {
   3.578379483085558*^9, 3.578379496922782*^9}, {3.5783795638312993`*^9, 
   3.5783797837136855`*^9}, {3.578380015670493*^9, 3.578380570345067*^9}, {
   3.5783806138535433`*^9, 3.5783806149455457`*^9}, {3.578380675645252*^9, 
   3.5783807269069424`*^9}, 3.5783807789174337`*^9, 3.5783894188369675`*^9, {
   3.5783894739674644`*^9, 3.5783894894270916`*^9}, {3.578389884825386*^9, 
   3.5783899611875205`*^9}, {3.5783900128548107`*^9, 
   3.5783900278308372`*^9}, {3.5783901592766685`*^9, 3.5783901772323*^9}, {
   3.578390300753317*^9, 3.5783903380529823`*^9}, 3.5783905350189285`*^9}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"0", "2.033120988827489`", "2.033120988827489`", "2.033120988827489`", 
      "2.033120988827489`"},
     {"2.033120988827489`", "0.`", "3.3200726719800846`", 
      "3.3200726719800846`", "3.3200726719800846`"},
     {"2.033120988827489`", "3.3200726719800846`", "0.`", 
      "3.3200726719800846`", "3.3200726719800846`"},
     {"2.033120988827489`", "3.3200726719800846`", "3.3200726719800846`", 
      "0.`", "3.3200726719800846`"},
     {"2.033120988827489`", "3.3200726719800846`", "3.3200726719800846`", 
      "3.3200726719800846`", "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.5784143360269413`*^9, 3.5784155275102344`*^9, 
  3.5784155609254932`*^9, 3.5784156295500135`*^9, 3.5784156752112937`*^9, 
  3.5784157214185753`*^9, 3.57841583238157*^9, 3.5784171205722327`*^9, 
  3.578417281455315*^9, 3.5784173355718107`*^9, 3.578417735416113*^9, 
  3.578418067010295*^9, 3.5784181925905156`*^9, 3.578418290917488*^9, 
  3.5784183399171743`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"dwn", "=", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"Norm", "[", 
       RowBox[{
        RowBox[{"w", "[", 
         RowBox[{"[", "i", "]"}], "]"}], "-", 
        RowBox[{"n", "[", 
         RowBox[{"[", "j", "]"}], "]"}]}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "nc"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"j", ",", "1", ",", "nc"}], "}"}]}], "]"}]}], ";"}], 
  "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
     RowBox[{"dwn", "[", 
      RowBox[{"[", 
       RowBox[{"1", ",", "1"}], "]"}], "]"}], "=", 
     RowBox[{
      RowBox[{"1", "/", "1.5"}], "*", "R1"}]}], ";"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{"MatrixForm", "[", "dwn", "]"}], "\[IndentingNewLine]"}], "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5783693559648447`*^9, 3.5783693662920628`*^9}, {
   3.578369554833994*^9, 3.5783695604344034`*^9}, 3.5783701337978106`*^9, {
   3.578371289853441*^9, 3.578371302442663*^9}, {3.57837134613834*^9, 
   3.57837135733916*^9}, {3.578371392314421*^9, 3.5783716741757164`*^9}, {
   3.578371712130583*^9, 3.578371724501405*^9}, {3.5783726153945694`*^9, 
   3.5783726153945694`*^9}, 3.5783726553930397`*^9, 3.578372749586005*^9, {
   3.5783728157925215`*^9, 3.5783728573041945`*^9}, {3.578373248116081*^9, 
   3.5783732528896894`*^9}, {3.5783732892065535`*^9, 
   3.5783733600150776`*^9}, {3.578373605637509*^9, 3.5783736465095806`*^9}, 
   3.578373878918789*^9, {3.5783778323218584`*^9, 3.5783779683072968`*^9}, {
   3.5783782452545834`*^9, 3.5783782551606007`*^9}, {3.5783783031618853`*^9, 
   3.5783783174047103`*^9}, {3.5783792814240036`*^9, 
   3.5783793528877287`*^9}, {3.578379408611027*^9, 3.5783794493738985`*^9}, {
   3.578379483085558*^9, 3.578379496922782*^9}, {3.5783795638312993`*^9, 
   3.5783797837136855`*^9}, {3.578380015670493*^9, 3.578380570345067*^9}, {
   3.5783806138535433`*^9, 3.5783806149455457`*^9}, {3.578380675645252*^9, 
   3.5783807269069424`*^9}, 3.5783807789174337`*^9, 3.5783894188369675`*^9, {
   3.5783894739674644`*^9, 3.5783894894270916`*^9}, {3.578389884825386*^9, 
   3.5783899611875205`*^9}, {3.5783900128548107`*^9, 
   3.5783900278308372`*^9}, {3.5783901592766685`*^9, 3.5783901772323*^9}, {
   3.578390300753317*^9, 3.57839034279539*^9}, 3.578390533505726*^9, {
   3.578395020188806*^9, 3.578395024541214*^9}, {3.5783952278875713`*^9, 
   3.5783952342679825`*^9}, {3.578399945411931*^9, 3.5784000197616615`*^9}, {
   3.5784000519601183`*^9, 3.5784000616165347`*^9}, {3.5784002276164265`*^9, 
   3.578400228146827*^9}, {3.5784002591752815`*^9, 3.5784002634184895`*^9}, 
   3.578400302434158*^9, {3.578400542112979*^9, 3.578400543220581*^9}, {
   3.578400693074444*^9, 3.578400700125656*^9}, {3.5784007327921133`*^9, 
   3.578400734008916*^9}, 3.5784007912922163`*^9, {3.578400901038409*^9, 
   3.5784009063736186`*^9}, 3.5784010782771206`*^9, 3.578401329222361*^9, 
   3.578402693397957*^9}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"0", "2.033120988827489`", "2.033120988827489`", "2.033120988827489`", 
      "2.033120988827489`"},
     {"1.4996853200000002`", "0.5334356688274887`", "2.9009190752909126`", 
      "2.9009190752909126`", "2.9009190752909126`"},
     {"1.4996853200000002`", "2.9009190752909126`", "0.5334356688274887`", 
      "2.9009190752909126`", "2.9009190752909126`"},
     {"1.4996853200000002`", "2.9009190752909126`", "2.9009190752909126`", 
      "0.5334356688274887`", "2.9009190752909126`"},
     {"1.4996853200000002`", "2.9009190752909126`", "2.9009190752909126`", 
      "2.9009190752909126`", "0.5334356688274887`"}
    },
    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.5784143360581417`*^9, 3.578415527525834*^9, 
  3.5784155609410934`*^9, 3.5784156295656137`*^9, 3.578415675226894*^9, 
  3.578415721434175*^9, 3.57841583239717*^9, 3.5784171205878325`*^9, 
  3.5784172814709153`*^9, 3.5784173355874104`*^9, 3.5784177354317126`*^9, 
  3.5784180670414953`*^9, 3.5784181926061153`*^9, 3.5784182909330883`*^9, 
  3.5784183399327745`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"dwr", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"If", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"dw", "[", 
         RowBox[{"[", 
          RowBox[{"i", ",", "j"}], "]"}], "]"}], ">", " ", "0"}], ",", 
       RowBox[{"1", "/", 
        RowBox[{"dw", "[", 
         RowBox[{"[", 
          RowBox[{"i", ",", "j"}], "]"}], "]"}]}], ",", "0"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "nc"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"j", ",", "1", ",", "nc"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"MatrixForm", "[", "dwr", "]"}]}], "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5783693559648447`*^9, 3.5783693662920628`*^9}, {
   3.578369554833994*^9, 3.5783695604344034`*^9}, 3.5783701337978106`*^9, {
   3.578371289853441*^9, 3.578371302442663*^9}, {3.57837134613834*^9, 
   3.57837135733916*^9}, {3.578371392314421*^9, 3.5783716741757164`*^9}, {
   3.578371712130583*^9, 3.578371724501405*^9}, {3.5783726153945694`*^9, 
   3.5783726153945694`*^9}, 3.5783726553930397`*^9, 3.578372749586005*^9, {
   3.5783728157925215`*^9, 3.5783728573041945`*^9}, {3.578373248116081*^9, 
   3.5783732528896894`*^9}, {3.5783732892065535`*^9, 
   3.5783733600150776`*^9}, {3.578373605637509*^9, 3.5783736465095806`*^9}, 
   3.578373878918789*^9, {3.5783778323218584`*^9, 3.5783779683072968`*^9}, {
   3.5783782452545834`*^9, 3.5783782551606007`*^9}, {3.5783783031618853`*^9, 
   3.5783783174047103`*^9}, {3.5783792814240036`*^9, 
   3.5783793528877287`*^9}, {3.578379408611027*^9, 3.5783794493738985`*^9}, {
   3.578379483085558*^9, 3.578379496922782*^9}, {3.5783795638312993`*^9, 
   3.5783797837136855`*^9}, {3.578380015670493*^9, 3.578380570345067*^9}, {
   3.5783806138535433`*^9, 3.5783806149455457`*^9}, {3.578380675645252*^9, 
   3.5783807269069424`*^9}, 3.5783807789174337`*^9, 3.5783894188369675`*^9, {
   3.5783894739674644`*^9, 3.5783894894270916`*^9}, {3.578389884825386*^9, 
   3.5783899611875205`*^9}, {3.5783900128548107`*^9, 
   3.5783900278308372`*^9}, {3.5783901592766685`*^9, 3.5783901772323*^9}, {
   3.578390300753317*^9, 3.5783903881290703`*^9}, {3.578390483616838*^9, 
   3.578390531493322*^9}, {3.5783906225974817`*^9, 3.5783906229094825`*^9}}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"0", "0.6668065537908978`", "0.6668065537908978`", 
      "0.6668065537908978`", "0.6668065537908978`"},
     {"0.6668065537908978`", "0", "0.4083339534828509`", 
      "0.4083339534828509`", "0.4083339534828509`"},
     {"0.6668065537908978`", "0.4083339534828509`", "0", 
      "0.4083339534828509`", "0.4083339534828509`"},
     {"0.6668065537908978`", "0.4083339534828509`", "0.4083339534828509`", 
      "0", "0.4083339534828509`"},
     {"0.6668065537908978`", "0.4083339534828509`", "0.4083339534828509`", 
      "0.4083339534828509`", "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.578414336073742*^9, 3.5784155275414343`*^9, 
  3.578415560956693*^9, 3.578415629581214*^9, 3.578415675242494*^9, 
  3.5784157214497747`*^9, 3.5784158324127703`*^9, 3.5784171206034327`*^9, 
  3.578417281486515*^9, 3.57841733560301*^9, 3.578417735447313*^9, 
  3.5784180670570955`*^9, 3.5784181926217155`*^9, 3.5784182909642887`*^9, 
  3.5784183399483747`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"dnr", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"If", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"dn", "[", 
         RowBox[{"[", 
          RowBox[{"i", ",", "j"}], "]"}], "]"}], ">", " ", "0"}], ",", 
       RowBox[{"1", "/", 
        RowBox[{"dn", "[", 
         RowBox[{"[", 
          RowBox[{"i", ",", "j"}], "]"}], "]"}]}], ",", "0"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "nc"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"j", ",", "1", ",", "nc"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"MatrixForm", "[", "dnr", "]"}]}], "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5783693559648447`*^9, 3.5783693662920628`*^9}, {
   3.578369554833994*^9, 3.5783695604344034`*^9}, 3.5783701337978106`*^9, {
   3.578371289853441*^9, 3.578371302442663*^9}, {3.57837134613834*^9, 
   3.57837135733916*^9}, {3.578371392314421*^9, 3.5783716741757164`*^9}, {
   3.578371712130583*^9, 3.578371724501405*^9}, {3.5783726153945694`*^9, 
   3.5783726153945694`*^9}, 3.5783726553930397`*^9, 3.578372749586005*^9, {
   3.5783728157925215`*^9, 3.5783728573041945`*^9}, {3.578373248116081*^9, 
   3.5783732528896894`*^9}, {3.5783732892065535`*^9, 
   3.5783733600150776`*^9}, {3.578373605637509*^9, 3.5783736465095806`*^9}, 
   3.578373878918789*^9, {3.5783778323218584`*^9, 3.5783779683072968`*^9}, {
   3.5783782452545834`*^9, 3.5783782551606007`*^9}, {3.5783783031618853`*^9, 
   3.5783783174047103`*^9}, {3.5783792814240036`*^9, 
   3.5783793528877287`*^9}, {3.578379408611027*^9, 3.5783794493738985`*^9}, {
   3.578379483085558*^9, 3.578379496922782*^9}, {3.5783795638312993`*^9, 
   3.5783797837136855`*^9}, {3.578380015670493*^9, 3.578380570345067*^9}, {
   3.5783806138535433`*^9, 3.5783806149455457`*^9}, {3.578380675645252*^9, 
   3.5783807269069424`*^9}, 3.5783807789174337`*^9, 3.5783894188369675`*^9, {
   3.5783894739674644`*^9, 3.5783894894270916`*^9}, {3.578389884825386*^9, 
   3.5783899611875205`*^9}, {3.5783900128548107`*^9, 
   3.5783900278308372`*^9}, {3.5783901592766685`*^9, 3.5783901772323*^9}, {
   3.578390300753317*^9, 3.5783903881290703`*^9}, {3.578390483616838*^9, 
   3.578390514941693*^9}}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"0", "0.4918546439170376`", "0.4918546439170376`", 
      "0.4918546439170376`", "0.4918546439170376`"},
     {"0.4918546439170376`", "0", "0.301198226303764`", "0.301198226303764`", 
      "0.301198226303764`"},
     {"0.4918546439170376`", "0.301198226303764`", "0", "0.301198226303764`", 
      "0.301198226303764`"},
     {"0.4918546439170376`", "0.301198226303764`", "0.301198226303764`", "0", 
      "0.301198226303764`"},
     {"0.4918546439170376`", "0.301198226303764`", "0.301198226303764`", 
      "0.301198226303764`", "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.578414336089342*^9, 3.5784155275570345`*^9, 
  3.578415560972293*^9, 3.5784156295968137`*^9, 3.578415675258094*^9, 
  3.578415721465375*^9, 3.57841583242837*^9, 3.578417120619033*^9, 
  3.5784172815021152`*^9, 3.5784173356186104`*^9, 3.5784177354629126`*^9, 
  3.5784180670726953`*^9, 3.5784181926373158`*^9, 3.5784182909798884`*^9, 
  3.5784183399639745`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"dwnr", "=", 
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"If", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"dwn", "[", 
         RowBox[{"[", 
          RowBox[{"i", ",", "j"}], "]"}], "]"}], ">", " ", "0"}], ",", 
       RowBox[{"1", "/", 
        RowBox[{"dwn", "[", 
         RowBox[{"[", 
          RowBox[{"i", ",", "j"}], "]"}], "]"}]}], ",", 
       RowBox[{"1.5", "/", 
        RowBox[{"c", "[", 
         RowBox[{"[", "i", "]"}], "]"}]}]}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "nc"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"j", ",", "1", ",", "nc"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"MatrixForm", "[", "dwnr", "]"}]}], "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5783693559648447`*^9, 3.5783693662920628`*^9}, {
   3.578369554833994*^9, 3.5783695604344034`*^9}, 3.5783701337978106`*^9, {
   3.578371289853441*^9, 3.578371302442663*^9}, {3.57837134613834*^9, 
   3.57837135733916*^9}, {3.578371392314421*^9, 3.5783716741757164`*^9}, {
   3.578371712130583*^9, 3.578371724501405*^9}, {3.5783726153945694`*^9, 
   3.5783726153945694`*^9}, 3.5783726553930397`*^9, 3.578372749586005*^9, {
   3.5783728157925215`*^9, 3.5783728573041945`*^9}, {3.578373248116081*^9, 
   3.5783732528896894`*^9}, {3.5783732892065535`*^9, 
   3.5783733600150776`*^9}, {3.578373605637509*^9, 3.5783736465095806`*^9}, 
   3.578373878918789*^9, {3.5783778323218584`*^9, 3.5783779683072968`*^9}, {
   3.5783782452545834`*^9, 3.5783782551606007`*^9}, {3.5783783031618853`*^9, 
   3.5783783174047103`*^9}, {3.5783792814240036`*^9, 
   3.5783793528877287`*^9}, {3.578379408611027*^9, 3.5783794493738985`*^9}, {
   3.578379483085558*^9, 3.578379496922782*^9}, {3.5783795638312993`*^9, 
   3.5783797837136855`*^9}, {3.578380015670493*^9, 3.578380570345067*^9}, {
   3.5783806138535433`*^9, 3.5783806149455457`*^9}, {3.578380675645252*^9, 
   3.5783807269069424`*^9}, 3.5783807789174337`*^9, 3.5783894188369675`*^9, {
   3.5783894739674644`*^9, 3.5783894894270916`*^9}, {3.578389884825386*^9, 
   3.5783899611875205`*^9}, 3.5783905374681325`*^9, {3.5783952781976595`*^9, 
   3.578395330286151*^9}, {3.5784026468786755`*^9, 3.578402661386701*^9}, 
   3.5784028964167137`*^9, {3.57841816081326*^9, 3.578418176038886*^9}, {
   3.578418233431387*^9, 3.5784182455370083`*^9}}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", GridBox[{
     {"5.652273917333309`", "0.4918546439170376`", "0.4918546439170376`", 
      "0.4918546439170376`", "0.4918546439170376`"},
     {"0.6668065537908978`", "1.874640295798061`", "0.34471833720481054`", 
      "0.34471833720481054`", "0.34471833720481054`"},
     {"0.6668065537908978`", "0.34471833720481054`", "1.874640295798061`", 
      "0.34471833720481054`", "0.34471833720481054`"},
     {"0.6668065537908978`", "0.34471833720481054`", "0.34471833720481054`", 
      "1.874640295798061`", "0.34471833720481054`"},
     {"0.6668065537908978`", "0.34471833720481054`", "0.34471833720481054`", 
      "0.34471833720481054`", "1.874640295798061`"}
    },
    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.578414336104942*^9, 3.5784155275726347`*^9, 
  3.578415560987893*^9, 3.5784156296124134`*^9, 3.578415675273694*^9, 
  3.578415721480975*^9, 3.5784158324439697`*^9, 3.578417120650233*^9, 
  3.5784172815177155`*^9, 3.5784173356342106`*^9, 3.578417735478513*^9, 
  3.578418067088295*^9, 3.578418192652916*^9, 3.578418290995488*^9, 
  3.5784183399951744`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"vee", "=", 
   RowBox[{"vee", "+", 
    RowBox[{"0.5", "*", 
     RowBox[{"ow", ".", 
      RowBox[{"(", 
       RowBox[{"ow", ".", "dwr"}], ")"}]}]}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"vnn", "=", 
   RowBox[{"0.5", "*", 
    RowBox[{"cn", ".", 
     RowBox[{"(", 
      RowBox[{"cn", ".", "dnr"}], ")"}]}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"vne", "=", 
    RowBox[{"(*", 
     RowBox[{"vne", "+"}], "*)"}], 
    RowBox[{"cn", ".", 
     RowBox[{"(", 
      RowBox[{"ow", ".", "dwnr"}], ")"}]}]}], ";"}], "  ", 
  RowBox[{"(*", "0.90069", "**)"}]}], "\[IndentingNewLine]", 
 RowBox[{"MatrixForm", "[", 
  RowBox[{"ow", ".", "dwnr"}], "]"}], "\[IndentingNewLine]", 
 RowBox[{"cn", ".", "%"}]}], "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5783693559648447`*^9, 3.5783693662920628`*^9}, {
   3.578369554833994*^9, 3.5783695604344034`*^9}, 3.5783701337978106`*^9, {
   3.578371289853441*^9, 3.578371302442663*^9}, {3.57837134613834*^9, 
   3.57837135733916*^9}, {3.578371392314421*^9, 3.5783716741757164`*^9}, {
   3.578371712130583*^9, 3.578371724501405*^9}, {3.5783726153945694`*^9, 
   3.5783726153945694`*^9}, 3.5783726553930397`*^9, 3.578372749586005*^9, {
   3.5783728157925215`*^9, 3.5783728573041945`*^9}, {3.578373248116081*^9, 
   3.5783732528896894`*^9}, {3.5783732892065535`*^9, 
   3.5783733600150776`*^9}, {3.578373605637509*^9, 3.5783736465095806`*^9}, 
   3.578373878918789*^9, {3.5783778323218584`*^9, 3.5783779683072968`*^9}, {
   3.5783782452545834`*^9, 3.5783782551606007`*^9}, {3.5783783031618853`*^9, 
   3.5783783174047103`*^9}, {3.5783792814240036`*^9, 
   3.5783793528877287`*^9}, {3.578379408611027*^9, 3.5783794493738985`*^9}, {
   3.578379483085558*^9, 3.578379496922782*^9}, {3.5783795638312993`*^9, 
   3.5783797837136855`*^9}, {3.578380015670493*^9, 3.578380570345067*^9}, {
   3.5783806138535433`*^9, 3.5783806149455457`*^9}, {3.578380675645252*^9, 
   3.5783807269069424`*^9}, 3.5783807789174337`*^9, 3.5783894188369675`*^9, {
   3.5783894739674644`*^9, 3.5783894894270916`*^9}, {3.578389884825386*^9, 
   3.5783899611875205`*^9}, {3.578390760304924*^9, 3.578390770117341*^9}, 
   3.5783925300500326`*^9, {3.578392846309388*^9, 3.5783928501469946`*^9}, {
   3.578392901798685*^9, 3.5783929100823*^9}, {3.578394026841461*^9, 
   3.5783940274966626`*^9}, {3.578394057979116*^9, 3.5783940627215242`*^9}, {
   3.578395120184982*^9, 3.578395126596593*^9}, {3.5783954002366734`*^9, 
   3.5783955994802237`*^9}, {3.5783956359218874`*^9, 3.578395682285169*^9}, {
   3.578396078042264*^9, 3.578396117759934*^9}, {3.5783961994260774`*^9, 
   3.5783962032168837`*^9}, 3.578400136059866*^9, {3.5784008039750385`*^9, 
   3.578400817313062*^9}, {3.578401035065045*^9, 3.5784010405718546`*^9}, 
   3.5784094487341146`*^9, {3.578418030584231*^9, 3.578418055949876*^9}, {
   3.578418260684635*^9, 3.5784182806058702`*^9}}],

Cell[BoxData[
 TagBox[
  RowBox[{"(", "\[NoBreak]", 
   TagBox[GridBox[{
      {
       RowBox[{"-", "16.6390002649938`"}]},
      {
       RowBox[{"-", "6.801299902659061`"}]},
      {
       RowBox[{"-", "6.801299902659061`"}]},
      {
       RowBox[{"-", "6.801299902659061`"}]},
      {
       RowBox[{"-", "6.801299902659061`"}]}
     },
     GridBoxAlignment->{
      "Columns" -> {{Center}}, "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$]]]], "Output",
 CellChangeTimes->{3.5784143361205416`*^9, 3.5784155276038346`*^9, 
  3.5784155610034933`*^9, 3.5784156296280136`*^9, 3.5784156752892942`*^9, 
  3.5784157214965754`*^9, 3.57841583245957*^9, 3.5784171206658325`*^9, 
  3.5784172815333157`*^9, 3.5784173356498103`*^9, 3.578417735494113*^9, 
  3.578418067103895*^9, 3.578418192684116*^9, 3.5784182910110884`*^9, 
  3.5784183400107746`*^9}],

Cell[BoxData[
 RowBox[{"-", "127.03920120059902`"}]], "Output",
 CellChangeTimes->{3.5784143361205416`*^9, 3.5784155276038346`*^9, 
  3.5784155610034933`*^9, 3.5784156296280136`*^9, 3.5784156752892942`*^9, 
  3.5784157214965754`*^9, 3.57841583245957*^9, 3.5784171206658325`*^9, 
  3.5784172815333157`*^9, 3.5784173356498103`*^9, 3.578417735494113*^9, 
  3.578418067103895*^9, 3.578418192684116*^9, 3.5784182910110884`*^9, 
  3.5784183400107746`*^9}]
}, Open  ]],

Cell[BoxData[{
 RowBox[{
  RowBox[{"b1", "=", 
   RowBox[{"{", 
    RowBox[{"0", ",", "d2", ",", 
     RowBox[{"-", "d2"}], ",", "d2", ",", 
     RowBox[{"-", "d2"}]}], "}"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"b2", "=", 
   RowBox[{"{", 
    RowBox[{"0", ",", "d2", ",", 
     RowBox[{"-", "d2"}], ",", 
     RowBox[{"-", "d2"}], ",", "d2"}], "}"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"b3", "=", 
    RowBox[{"{", 
     RowBox[{"0", ",", "d2", ",", "d2", ",", 
      RowBox[{"-", "d2"}], ",", 
      RowBox[{"-", "d2"}]}], "}"}]}], ";"}], "        "}], "\n", 
 RowBox[{
  RowBox[{"a1", "=", 
   RowBox[{"{", 
    RowBox[{"0", ",", "d4", ",", 
     RowBox[{"-", "d4"}], ",", "d4", ",", 
     RowBox[{"-", "d4"}]}], "}"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"a2", "=", 
   RowBox[{"{", 
    RowBox[{"0", ",", "d4", ",", 
     RowBox[{"-", "d4"}], ",", 
     RowBox[{"-", "d4"}], ",", "d4"}], "}"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"a3", "=", 
   RowBox[{"{", 
    RowBox[{"0", ",", "d4", ",", "d4", ",", 
     RowBox[{"-", "d4"}], ",", 
     RowBox[{"-", "d4"}]}], "}"}]}], ";"}]}], "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5783693559648447`*^9, 3.5783693662920628`*^9}, {
   3.578369554833994*^9, 3.5783695604344034`*^9}, 3.5783701337978106`*^9, {
   3.578371289853441*^9, 3.578371302442663*^9}, {3.57837134613834*^9, 
   3.57837135733916*^9}, {3.578371392314421*^9, 3.5783716741757164`*^9}, {
   3.578371712130583*^9, 3.578371724501405*^9}, {3.5783726153945694`*^9, 
   3.5783726153945694`*^9}, 3.5783726553930397`*^9, 3.578372749586005*^9, {
   3.5783728157925215`*^9, 3.5783728573041945`*^9}, {3.578373248116081*^9, 
   3.5783732528896894`*^9}, {3.5783732892065535`*^9, 
   3.5783733600150776`*^9}, {3.578373605637509*^9, 3.5783736465095806`*^9}, 
   3.578373878918789*^9, {3.5783778323218584`*^9, 3.5783779683072968`*^9}, {
   3.5783782452545834`*^9, 3.5783782551606007`*^9}, {3.5783783031618853`*^9, 
   3.5783783174047103`*^9}, {3.5783792814240036`*^9, 
   3.5783793528877287`*^9}, {3.578379408611027*^9, 3.5783794493738985`*^9}, {
   3.578379483085558*^9, 3.578379496922782*^9}, {3.5783795638312993`*^9, 
   3.5783797837136855`*^9}, {3.578380015670493*^9, 3.578380570345067*^9}, {
   3.5783806138535433`*^9, 3.5783806149455457`*^9}, {3.578380675645252*^9, 
   3.5783807269069424`*^9}, 3.5783807789174337`*^9, {3.578391576435958*^9, 
   3.5783915916927843`*^9}, {3.5783918945985165`*^9, 
   3.5783919055497355`*^9}, {3.578392550954069*^9, 3.5783925648224936`*^9}, {
   3.5783926112949753`*^9, 3.578392616895385*^9}, {3.578392688624311*^9, 
   3.5783927172503614`*^9}, {3.578392777310467*^9, 3.5783927862492824`*^9}, 
   3.57839286456142*^9, 3.5783945876156464`*^9, {3.5783947732247725`*^9, 
   3.578394784456792*^9}, {3.578409205545287*^9, 3.578409212144099*^9}, {
   3.5784144291591053`*^9, 3.578414434197914*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"Epot", "=", 
   RowBox[{"vne", "+", "vee", "+", "vnn"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"func", "=", 
   RowBox[{"Ekin", "+", "Epot"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"t", " ", "=", " ", 
    RowBox[{"FindMinimum", "[", 
     RowBox[{"func", ",", 
      RowBox[{"{", 
       RowBox[{"R1", ",", "0.267"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"R2", ",", "1.244"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"p", ",", "0.65"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"Method", " ", "->", " ", "Automatic"}], "}"}], ",", " ", 
      RowBox[{"{", 
       RowBox[{"MaxIterations", " ", "->", " ", "500"}], "}"}]}], "]"}]}], 
   "*)"}]}]}], "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5783693559648447`*^9, 3.5783693662920628`*^9}, {
   3.578369554833994*^9, 3.5783695604344034`*^9}, 3.5783701337978106`*^9, {
   3.578371289853441*^9, 3.578371302442663*^9}, {3.57837134613834*^9, 
   3.57837135733916*^9}, {3.578371392314421*^9, 3.5783716741757164`*^9}, {
   3.578371712130583*^9, 3.578371724501405*^9}, {3.5783726153945694`*^9, 
   3.5783726153945694`*^9}, 3.5783726553930397`*^9, 3.578372749586005*^9, {
   3.5783728157925215`*^9, 3.5783728573041945`*^9}, {3.578373248116081*^9, 
   3.5783732528896894`*^9}, {3.5783732892065535`*^9, 
   3.5783733600150776`*^9}, {3.578373605637509*^9, 3.5783736465095806`*^9}, 
   3.578373878918789*^9, {3.5783778323218584`*^9, 3.5783779683072968`*^9}, {
   3.5783782452545834`*^9, 3.5783782551606007`*^9}, {3.5783783031618853`*^9, 
   3.5783783174047103`*^9}, {3.5783792814240036`*^9, 
   3.5783793528877287`*^9}, {3.578379408611027*^9, 3.5783794493738985`*^9}, {
   3.578379483085558*^9, 3.578379496922782*^9}, {3.5783795638312993`*^9, 
   3.5783797837136855`*^9}, {3.578380015670493*^9, 3.578380570345067*^9}, {
   3.5783806138535433`*^9, 3.5783806149455457`*^9}, {3.578380675645252*^9, 
   3.5783807269069424`*^9}, 3.5783807789174337`*^9, {3.578391576435958*^9, 
   3.5783915916927843`*^9}, {3.5783918945985165`*^9, 
   3.5783919055497355`*^9}, {3.578392550954069*^9, 3.5783925648224936`*^9}, {
   3.5783926112949753`*^9, 3.578392616895385*^9}, {3.578392688624311*^9, 
   3.5783927172503614`*^9}, {3.578392777310467*^9, 3.5783927862492824`*^9}, 
   3.57839286456142*^9, 3.5783945876156464`*^9, {3.5783947732247725`*^9, 
   3.578394784456792*^9}, {3.578409205545287*^9, 3.578409212144099*^9}, 
   3.5784144291591053`*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{
  "func", "\[IndentingNewLine]", "Ekin", "\[IndentingNewLine]", "vne", "  ", 
   RowBox[{"(*", 
    RowBox[{"/.", " ", 
     RowBox[{"t", "[", 
      RowBox[{"[", "2", "]"}], "]"}]}], "*)"}], "\n", "vee", "  ", 
   RowBox[{"(*", 
    RowBox[{"/.", " ", 
     RowBox[{"t", "[", 
      RowBox[{"[", "2", "]"}], "]"}]}], "*)"}], "\n", "vnn", "  ", 
   RowBox[{"(*", 
    RowBox[{"/.", " ", 
     RowBox[{"t", "[", 
      RowBox[{"[", "2", "]"}], "]"}]}], "*)"}], "\n", 
   RowBox[{"vir", "=", 
    RowBox[{
     RowBox[{"-", "Epot"}], "/", "Ekin"}]}], "  ", 
   RowBox[{"(*", 
    RowBox[{"/.", " ", 
     RowBox[{"t", "[", 
      RowBox[{"[", "2", "]"}], "]"}]}], "*)"}], "\n", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"R1", "+", "R2", "+", "p"}], ")"}], "*", "bohr", "  ", 
    RowBox[{"(*", 
     RowBox[{"/.", " ", 
      RowBox[{"t", "[", 
       RowBox[{"[", "2", "]"}], "]"}]}], "*)"}]}]}]}]], "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5783693559648447`*^9, 3.5783693662920628`*^9}, {
   3.578369554833994*^9, 3.5783695604344034`*^9}, 3.5783696029912786`*^9, {
   3.5783698599861298`*^9, 3.5783699151166267`*^9}, {3.5783699654735155`*^9, 
   3.5783699759255333`*^9}, 3.578370096825746*^9, {3.5783723049696245`*^9, 
   3.578372305827626*^9}, {3.5783729062258806`*^9, 3.5783729062258806`*^9}, {
   3.578373477623684*^9, 3.578373497186118*^9}, {3.5783736524999914`*^9, 
   3.5783736582720013`*^9}, 3.578373849497137*^9, {3.57837838017922*^9, 
   3.5783783892272367`*^9}, {3.578378465870171*^9, 3.578378465885771*^9}, {
   3.578378532607088*^9, 3.5783785430279064`*^9}, {3.5783916033460045`*^9, 
   3.5783916076984124`*^9}, {3.578391690815358*^9, 3.578391693436163*^9}}],

Cell[BoxData[
 RowBox[{"-", "46.4412273064229`"}]], "Output",
 CellChangeTimes->{3.578414336151742*^9, 3.5784155276194344`*^9, 
  3.5784155610346932`*^9, 3.578415629659214*^9, 3.5784156753204937`*^9, 
  3.5784157215277753`*^9, 3.5784158324907703`*^9, 3.5784171206814327`*^9, 
  3.5784172815645156`*^9, 3.5784173356810102`*^9, 3.5784177355097127`*^9, 
  3.578418067135095*^9, 3.578418192715316*^9, 3.5784182910422883`*^9, 
  3.5784183400263743`*^9}],

Cell[BoxData["40.56089216702621`"], "Output",
 CellChangeTimes->{3.578414336151742*^9, 3.5784155276194344`*^9, 
  3.5784155610346932`*^9, 3.578415629659214*^9, 3.5784156753204937`*^9, 
  3.5784157215277753`*^9, 3.5784158324907703`*^9, 3.5784171206814327`*^9, 
  3.5784172815645156`*^9, 3.5784173356810102`*^9, 3.5784177355097127`*^9, 
  3.578418067135095*^9, 3.578418192715316*^9, 3.5784182910422883`*^9, 
  3.5784183400263743`*^9}],

Cell[BoxData[
 RowBox[{"-", "127.03920120059902`"}]], "Output",
 CellChangeTimes->{3.578414336151742*^9, 3.5784155276194344`*^9, 
  3.5784155610346932`*^9, 3.578415629659214*^9, 3.5784156753204937`*^9, 
  3.5784157215277753`*^9, 3.5784158324907703`*^9, 3.5784171206814327`*^9, 
  3.5784172815645156`*^9, 3.5784173356810102`*^9, 3.5784177355097127`*^9, 
  3.578418067135095*^9, 3.578418192715316*^9, 3.5784182910422883`*^9, 
  3.5784183400419745`*^9}],

Cell[BoxData["26.425380915318428`"], "Output",
 CellChangeTimes->{3.578414336151742*^9, 3.5784155276194344`*^9, 
  3.5784155610346932`*^9, 3.578415629659214*^9, 3.5784156753204937`*^9, 
  3.5784157215277753`*^9, 3.5784158324907703`*^9, 3.5784171206814327`*^9, 
  3.5784172815645156`*^9, 3.5784173356810102`*^9, 3.5784177355097127`*^9, 
  3.578418067135095*^9, 3.578418192715316*^9, 3.5784182910422883`*^9, 
  3.5784183400419745`*^9}],

Cell[BoxData["13.611700811831488`"], "Output",
 CellChangeTimes->{3.578414336151742*^9, 3.5784155276194344`*^9, 
  3.5784155610346932`*^9, 3.578415629659214*^9, 3.5784156753204937`*^9, 
  3.5784157215277753`*^9, 3.5784158324907703`*^9, 3.5784171206814327`*^9, 
  3.5784172815645156`*^9, 3.5784173356810102`*^9, 3.5784177355097127`*^9, 
  3.578418067135095*^9, 3.578418192715316*^9, 3.5784182910422883`*^9, 
  3.5784183400419745`*^9}],

Cell[BoxData["2.144975488092865`"], "Output",
 CellChangeTimes->{3.578414336151742*^9, 3.5784155276194344`*^9, 
  3.5784155610346932`*^9, 3.578415629659214*^9, 3.5784156753204937`*^9, 
  3.5784157215277753`*^9, 3.5784158324907703`*^9, 3.5784171206814327`*^9, 
  3.5784172815645156`*^9, 3.5784173356810102`*^9, 3.5784177355097127`*^9, 
  3.578418067135095*^9, 3.578418192715316*^9, 3.5784182910422883`*^9, 
  3.5784183400419745`*^9}],

Cell[BoxData["1.075880865504764`"], "Output",
 CellChangeTimes->{3.578414336151742*^9, 3.5784155276194344`*^9, 
  3.5784155610346932`*^9, 3.578415629659214*^9, 3.5784156753204937`*^9, 
  3.5784157215277753`*^9, 3.5784158324907703`*^9, 3.5784171206814327`*^9, 
  3.5784172815645156`*^9, 3.5784173356810102`*^9, 3.5784177355097127`*^9, 
  3.578418067135095*^9, 3.578418192715316*^9, 3.5784182910422883`*^9, 
  3.5784183400419745`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"plot1", "=", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"Graphics3D", "[", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"Opacity", "[", "0.6", "]"}], ",", 
         RowBox[{"Sphere", "[", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{
             RowBox[{"a1", "[", 
              RowBox[{"[", "i", "]"}], "]"}], ",", 
             RowBox[{"a2", "[", 
              RowBox[{"[", "i", "]"}], "]"}], ",", 
             RowBox[{"a3", "[", 
              RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", 
           RowBox[{"c", "[", 
            RowBox[{"[", "i", "]"}], "]"}]}], "]"}], ",", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"Darker", "[", 
            RowBox[{"Red", ",", "1.0"}], "]"}], ",", 
           RowBox[{"Sphere", "[", 
            RowBox[{
             RowBox[{"{", 
              RowBox[{
               RowBox[{"b1", "[", 
                RowBox[{"[", "i", "]"}], "]"}], ",", 
               RowBox[{"b2", "[", 
                RowBox[{"[", "i", "]"}], "]"}], ",", 
               RowBox[{"b3", "[", 
                RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", "0.08"}], 
            "]"}]}], "}"}]}], "}"}], " ", "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "nc"}], "}"}]}], "]"}]}], 
   RowBox[{"(*", 
    RowBox[{"/.", " ", 
     RowBox[{"t", "[", 
      RowBox[{"[", "2", "]"}], "]"}]}], "*)"}], ";"}], "\n"}], "\n", 
 RowBox[{
  RowBox[{"Show", "[", 
   RowBox[{"plot1", ",", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"AspectRatio", " ", "\[Rule]", " ", "Automatic"}], ",", 
      RowBox[{"Axes", " ", "->", " ", "True"}]}], "}"}]}], "]"}], 
  "\n"}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"plot2", "=", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"Graphics", "[", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"Circle", "[", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{
             RowBox[{"a1", "[", 
              RowBox[{"[", "i", "]"}], "]"}], ",", 
             RowBox[{"a3", "[", 
              RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", 
           RowBox[{"c", "[", 
            RowBox[{"[", "i", "]"}], "]"}]}], "]"}], ",", 
         RowBox[{"Disk", "[", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{
             RowBox[{"b1", "[", 
              RowBox[{"[", "i", "]"}], "]"}], ",", 
             RowBox[{"b3", "[", 
              RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", "0.08"}], "]"}]}],
         "}"}], " ", "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "nc"}], "}"}]}], "]"}]}], 
   RowBox[{"(*", 
    RowBox[{"/.", " ", 
     RowBox[{"t", "[", 
      RowBox[{"[", "2", "]"}], "]"}]}], "*)"}], ";", 
   RowBox[{"(*", 
    RowBox[{"/.", " ", 
     RowBox[{"t", "[", 
      RowBox[{"[", "2", "]"}], "]"}]}], "*)"}], ";"}], "\n"}], "\n", 
 RowBox[{
  RowBox[{"Show", "[", 
   RowBox[{"plot2", ",", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"AspectRatio", " ", "\[Rule]", " ", "Automatic"}], ",", 
      RowBox[{"Axes", " ", "->", " ", "True"}], ",", 
      RowBox[{"GridLines", " ", "->", " ", "Automatic"}], ",", " ", 
      RowBox[{"PlotRange", " ", "\[Rule]", " ", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"{", 
          RowBox[{
           RowBox[{"-", "3"}], ",", "3"}], "}"}], ",", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"-", "3"}], ",", "3"}], "}"}]}], "}"}]}], ",", " ", 
      RowBox[{"Frame", " ", "->", " ", "True"}]}], "}"}]}], "]"}], 
  "\n"}], "\n"}], "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5783729062414804`*^9, 3.5783729062414804`*^9}, {
   3.578378465901371*^9, 3.578378466135371*^9}, {3.5783919335673847`*^9, 
   3.5783919401349964`*^9}, {3.578391990975486*^9, 3.578392058929205*^9}, {
   3.5783920945908675`*^9, 3.578392104684085*^9}, {3.5783921429509525`*^9, 
   3.5783921446201553`*^9}, 3.5783921794862165`*^9, 3.578392214445878*^9, 
   3.578392277126788*^9, {3.5783923342228885`*^9, 3.578392362256138*^9}, {
   3.578393245046088*^9, 3.578393299536984*^9}, {3.5783933461186657`*^9, 
   3.5783933685827055`*^9}, {3.578393404899569*^9, 3.578393415179987*^9}, {
   3.578393477252496*^9, 3.578393498936534*^9}, {3.578402910207138*^9, 
   3.5784029167435493`*^9}, {3.5784154216016483`*^9, 
   3.5784154620525193`*^9}, {3.5784155549506826`*^9, 3.578415555871084*^9}, {
   3.578415608786377*^9, 3.5784156211415987`*^9}, {3.5784156592056656`*^9, 
   3.5784157153033648`*^9}, {3.5784158197455473`*^9, 3.5784158281227627`*^9}}],

Cell[BoxData[
 Graphics3DBox[{
   {Opacity[0.6], SphereBox[{0, 0, 0}, 0.26537992], 
    {RGBColor[0., 0., 0.], SphereBox[{0, 0, 0}, 0.08]}}, 
   {Opacity[0.6], 
    SphereBox[{0.8658437232017302, 0.8658437232017302, 0.8658437232017302}, 
     1.2343054], 
    {RGBColor[0., 0., 0.], 
     SphereBox[{1.173822950194629, 1.173822950194629, 1.173822950194629}, 
      0.08]}}, 
   {Opacity[0.6], 
    SphereBox[{-0.8658437232017302, -0.8658437232017302, 0.8658437232017302}, 
     1.2343054], 
    {RGBColor[0., 0., 0.], 
     SphereBox[{-1.173822950194629, -1.173822950194629, 1.173822950194629}, 
      0.08]}}, 
   {Opacity[0.6], 
    SphereBox[{0.8658437232017302, -0.8658437232017302, -0.8658437232017302}, 
     1.2343054], 
    {RGBColor[0., 0., 0.], 
     SphereBox[{1.173822950194629, -1.173822950194629, -1.173822950194629}, 
      0.08]}}, 
   {Opacity[0.6], 
    SphereBox[{-0.8658437232017302, 0.8658437232017302, -0.8658437232017302}, 
     1.2343054], 
    {RGBColor[0., 0., 0.], 
     SphereBox[{-1.173822950194629, 1.173822950194629, -1.173822950194629}, 
      0.08]}}},
  AspectRatio->Automatic,
  AutomaticImageSize->True,
  Axes->True,
  ImageSize->{360.2605029682672, 386.28357292862876`},
  ViewPoint->{-0.0010229225502464412`, -1.2121588193650408`, \
-3.159219832532237},
  ViewVertical->{
   0.026601423781041356`, -0.9644183107393388, -0.2630393243671024}]], \
"Output",
 CellChangeTimes->{3.578414336167342*^9, 3.5784155282902355`*^9, 
  3.578415561065893*^9, 3.578415629690414*^9, 3.578415675351694*^9, 
  3.5784157215745754`*^9, 3.57841583252197*^9, 3.5784171207282333`*^9, 
  3.5784172815957155`*^9, 3.5784173357122107`*^9, 3.5784177355409126`*^9, 
  3.578418067166295*^9, 3.578418192746516*^9, 3.5784182910734887`*^9, 
  3.5784183400575743`*^9}],

Cell[BoxData[
 GraphicsBox[{{CircleBox[{0, 0}, 0.26537992], DiskBox[{0, 0}, 0.08]}, {
    CircleBox[{0.8658437232017302, 0.8658437232017302}, 1.2343054], 
    DiskBox[{1.173822950194629, 1.173822950194629}, 0.08]}, {
    CircleBox[{-0.8658437232017302, 0.8658437232017302}, 1.2343054], 
    DiskBox[{-1.173822950194629, 1.173822950194629}, 0.08]}, {
    CircleBox[{0.8658437232017302, -0.8658437232017302}, 1.2343054], 
    DiskBox[{1.173822950194629, -1.173822950194629}, 0.08]}, {
    CircleBox[{-0.8658437232017302, -0.8658437232017302}, 1.2343054], 
    DiskBox[{-1.173822950194629, -1.173822950194629}, 0.08]}},
  AspectRatio->Automatic,
  Axes->True,
  Frame->True,
  GridLines->Automatic,
  PlotRange->{{-3, 3}, {-3, 3}}]], "Output",
 CellChangeTimes->{3.578414336167342*^9, 3.5784155282902355`*^9, 
  3.578415561065893*^9, 3.578415629690414*^9, 3.578415675351694*^9, 
  3.5784157215745754`*^9, 3.57841583252197*^9, 3.5784171207282333`*^9, 
  3.5784172815957155`*^9, 3.5784173357122107`*^9, 3.5784177355409126`*^9, 
  3.578418067166295*^9, 3.578418192746516*^9, 3.5784182910734887`*^9, 
  3.5784183400731745`*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{"from", " ", "old", " ", "program"}], " ", "*)"}], "       ", 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"Ekin", "=", 
     RowBox[{
      RowBox[{"2.25", "*", 
       RowBox[{
        RowBox[{"k", "[", 
         RowBox[{"[", "1", "]"}], "]"}], "/", 
        RowBox[{"R1", "^", "2"}]}]}], "+", 
      RowBox[{"9.0", "*", 
       RowBox[{
        RowBox[{"k", "[", 
         RowBox[{"[", "2", "]"}], "]"}], "/", 
        RowBox[{"R2", "^", "2"}]}]}]}]}], ";"}], "\n", " "}]}]], "Input",
 CellChangeTimes->{{3.5784163520300827`*^9, 3.578416352045683*^9}, {
  3.5784167316567497`*^9, 3.578416866097786*^9}, {3.5784168996066446`*^9, 
  3.5784169997968206`*^9}, {3.5784170630237317`*^9, 3.5784170940833864`*^9}, {
  3.5784172056703825`*^9, 3.578417273873702*^9}, {3.578417315619375*^9, 
  3.5784173253381925`*^9}, {3.5784175560625973`*^9, 3.5784176402559457`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"      ", 
  RowBox[{"vne", "=", 
   RowBox[{
    RowBox[{"-", "3"}], "*", 
    RowBox[{"z", "/", "R1", "   "}]}]}]}]], "Input",
 CellChangeTimes->{{3.5784163520300827`*^9, 3.578416352045683*^9}, {
   3.5784167316567497`*^9, 3.578416866097786*^9}, {3.5784168996066446`*^9, 
   3.5784169997968206`*^9}, {3.5784170630237317`*^9, 
   3.5784170940833864`*^9}, {3.5784172056703825`*^9, 3.578417273873702*^9}, {
   3.578417315619375*^9, 3.5784173253381925`*^9}, {3.5784175560625973`*^9, 
   3.5784176402559457`*^9}, 3.5784176844508233`*^9}],

Cell[BoxData[
 RowBox[{"-", "67.8272870079997`"}]], "Output",
 CellChangeTimes->{3.5784177356189127`*^9, 3.578418067244295*^9, 
  3.578418192808916*^9, 3.5784182911046886`*^9, 3.578418340135575*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"       ", 
  RowBox[{
   RowBox[{
    RowBox[{"vee", "=", 
     RowBox[{"3", "*", 
      RowBox[{
       RowBox[{"sig", "[", 
        RowBox[{"[", "1", "]"}], "]"}], "/", "R1"}]}]}], ";"}], "\n", 
   "       ", 
   RowBox[{
    RowBox[{"p", "=", 
     RowBox[{
      RowBox[{"p", "/", 
       RowBox[{"(", 
        RowBox[{"R1", "+", "R2"}], ")"}]}], "+", "1"}]}], ";", " ", 
    RowBox[{"wd", "=", "sq"}], ";", " ", 
    RowBox[{"f", "=", 
     RowBox[{"Sqrt", "[", 
      RowBox[{"2", "/", "3"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
   "       ", 
   RowBox[{
    RowBox[{"p2", "=", 
     RowBox[{"p", "*", "p"}]}], ";", " ", 
    RowBox[{"r", "=", 
     RowBox[{"R1", "+", "R2"}]}], ";", " ", 
    RowBox[{"r2", "=", 
     RowBox[{"r", "*", "r"}]}], ";"}]}]}]], "Input",
 CellChangeTimes->{{3.5784163520300827`*^9, 3.578416352045683*^9}, {
  3.5784167316567497`*^9, 3.578416866097786*^9}, {3.5784168996066446`*^9, 
  3.5784169997968206`*^9}, {3.5784170630237317`*^9, 3.5784170940833864`*^9}, {
  3.5784172056703825`*^9, 3.578417273873702*^9}, {3.578417315619375*^9, 
  3.5784173253381925`*^9}, {3.5784175560625973`*^9, 3.578417650863964*^9}, {
  3.578417864350339*^9, 3.578417880886368*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"       ", 
  RowBox[{"vne", "=", 
   RowBox[{"vne", "-", 
    RowBox[{"4", "*", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"3", "-", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{"p", "-", "1"}], ")"}], "*", 
           RowBox[{"(", 
            RowBox[{"1", "+", 
             RowBox[{"R1", "/", "R2"}]}], ")"}]}], ")"}], "^", "2"}]}], ")"}],
       "/", "R2"}]}]}]}]}]], "Input",
 CellChangeTimes->{{3.5784163520300827`*^9, 3.578416352045683*^9}, {
   3.5784167316567497`*^9, 3.578416866097786*^9}, {3.5784168996066446`*^9, 
   3.5784169997968206`*^9}, {3.5784170630237317`*^9, 
   3.5784170940833864`*^9}, {3.5784172056703825`*^9, 3.578417273873702*^9}, {
   3.578417315619375*^9, 3.5784173253381925`*^9}, {3.5784175560625973`*^9, 
   3.578417650863964*^9}, 3.578417692516037*^9}],

Cell[BoxData[
 RowBox[{"-", "76.94407433353535`"}]], "Output",
 CellChangeTimes->{3.578417735634513*^9, 3.5784180672598953`*^9, 
  3.5784181928245163`*^9, 3.578418291135889*^9, 3.5784183401511745`*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"       ", 
  RowBox[{
   RowBox[{
    RowBox[{"vee", "=", 
     RowBox[{"vee", "+", 
      RowBox[{"4", "*", "3", "*", 
       RowBox[{
        RowBox[{"sig", "[", 
         RowBox[{"[", "2", "]"}], "]"}], "/", "R2"}]}]}]}], ";"}], "\n", 
   "       ", 
   RowBox[{
    RowBox[{"vee", "=", 
     RowBox[{"vee", "+", 
      RowBox[{"8", "*", 
       RowBox[{"2", "/", "r"}]}]}]}], ";"}]}]}]], "Input",
 CellChangeTimes->{{3.5784163520300827`*^9, 3.578416352045683*^9}, {
  3.5784167316567497`*^9, 3.578416866097786*^9}, {3.5784168996066446`*^9, 
  3.5784169997968206`*^9}, {3.5784170630237317`*^9, 3.5784170940833864`*^9}, {
  3.5784172056703825`*^9, 3.578417273873702*^9}, {3.578417315619375*^9, 
  3.5784173253381925`*^9}, {3.5784175560625973`*^9, 3.578417662049184*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"       ", 
  RowBox[{"vne", "=", 
   RowBox[{"vne", "-", 
    RowBox[{"8", "*", 
     RowBox[{"z", "/", "r"}]}]}]}]}]], "Input",
 CellChangeTimes->{{3.5784163520300827`*^9, 3.578416352045683*^9}, {
   3.5784167316567497`*^9, 3.578416866097786*^9}, {3.5784168996066446`*^9, 
   3.5784169997968206`*^9}, {3.5784170630237317`*^9, 
   3.5784170940833864`*^9}, {3.5784172056703825`*^9, 3.578417273873702*^9}, {
   3.578417315619375*^9, 3.5784173253381925`*^9}, {3.5784175560625973`*^9, 
   3.578417662049184*^9}, 3.5784177018916535`*^9}],

Cell[BoxData[
 RowBox[{"-", "108.95078891549846`"}]], "Output",
 CellChangeTimes->{3.578417735650113*^9, 3.5784180672754955`*^9, 
  3.578418192840116*^9, 3.578418291151489*^9, 3.5784183401667747`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"       ", 
  RowBox[{"vne", "=", 
   RowBox[{"vne", "-", 
    RowBox[{"8", "/", 
     RowBox[{"(", 
      RowBox[{"r", "*", "p"}], ")"}]}]}]}]}]], "Input",
 CellChangeTimes->{{3.5784163520300827`*^9, 3.578416352045683*^9}, {
   3.5784167316567497`*^9, 3.578416866097786*^9}, {3.5784168996066446`*^9, 
   3.5784169997968206`*^9}, {3.5784170630237317`*^9, 
   3.5784170940833864`*^9}, {3.5784172056703825`*^9, 3.578417273873702*^9}, {
   3.578417315619375*^9, 3.5784173253381925`*^9}, {3.5784175560625973`*^9, 
   3.578417672844403*^9}, 3.578417704262858*^9}],

Cell[BoxData[
 RowBox[{"-", "112.88562606683476`"}]], "Output",
 CellChangeTimes->{3.578417735650113*^9, 3.5784180672910957`*^9, 
  3.5784181928557158`*^9, 3.5784182911670885`*^9, 3.578418340182375*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"       ", 
  RowBox[{
   RowBox[{
    RowBox[{"vnn", "=", 
     RowBox[{"4", "*", 
      RowBox[{"z", "/", 
       RowBox[{"(", 
        RowBox[{"r", "*", "p"}], ")"}]}]}]}], ";"}], "\n", "       ", 
   RowBox[{
    RowBox[{"vee", "=", 
     RowBox[{"vee", "+", 
      RowBox[{"12", "/", 
       RowBox[{"(", 
        RowBox[{"f", "*", "r"}], ")"}]}]}]}], ";"}]}]}]], "Input",
 CellChangeTimes->{{3.5784163520300827`*^9, 3.578416352045683*^9}, {
  3.5784167316567497`*^9, 3.578416866097786*^9}, {3.5784168996066446`*^9, 
  3.5784169997968206`*^9}, {3.5784170630237317`*^9, 3.5784170940833864`*^9}, {
  3.5784172056703825`*^9, 3.578417273873702*^9}, {3.578417315619375*^9, 
  3.5784173253381925`*^9}, {3.5784175560625973`*^9, 3.578417678086012*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"       ", 
  RowBox[{"vne", "=", 
   RowBox[{"vne", "-", 
    RowBox[{"24", "*", 
     RowBox[{"wd", "/", 
      RowBox[{"Sqrt", "[", 
       RowBox[{
        RowBox[{"2", "*", "r2", "*", "p"}], " ", "+", " ", 
        RowBox[{"3", "*", "r2", "*", "p2"}], " ", "+", " ", 
        RowBox[{"3", "*", "r2"}]}], "]"}]}]}]}]}]}]], "Input",
 CellChangeTimes->{{3.5784163520300827`*^9, 3.578416352045683*^9}, {
  3.5784167316567497`*^9, 3.578416866097786*^9}, {3.5784168996066446`*^9, 
  3.5784169997968206`*^9}, {3.5784170630237317`*^9, 3.5784170940833864`*^9}, {
  3.5784172056703825`*^9, 3.578417273873702*^9}, {3.578417315619375*^9, 
  3.5784173253381925`*^9}, {3.5784175560625973`*^9, 3.578417707164463*^9}}],

Cell[BoxData[
 RowBox[{"-", "121.15886615975022`"}]], "Output",
 CellChangeTimes->{3.5784177356813135`*^9, 3.578418067306696*^9, 
  3.578418192871316*^9, 3.578418291198289*^9, 3.5784183401979747`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"       ", 
  RowBox[{
   RowBox[{
    RowBox[{"vnn", "=", 
     RowBox[{"vnn", "+", 
      RowBox[{"3", "/", 
       RowBox[{"(", 
        RowBox[{"f", "*", "r", "*", "p"}], ")"}]}]}]}], ";"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", "\[IndentingNewLine]", 
   "\[IndentingNewLine]", "Ekin", "\[IndentingNewLine]", 
   RowBox[{"Epot", "=", 
    RowBox[{"vne", "+", "vnn", "+", "vee"}]}], "\[IndentingNewLine]", 
   RowBox[{"Epot", "+", "Ekin"}], "\[IndentingNewLine]", "vne", "  ", 
   RowBox[{"(*", 
    RowBox[{"/.", " ", 
     RowBox[{"t", "[", 
      RowBox[{"[", "2", "]"}], "]"}]}], "*)"}], "\n", "vee", "  ", 
   RowBox[{"(*", 
    RowBox[{"/.", " ", 
     RowBox[{"t", "[", 
      RowBox[{"[", "2", "]"}], "]"}]}], "*)"}], "\n", "vnn", "  ", 
   RowBox[{"(*", 
    RowBox[{"/.", " ", 
     RowBox[{"t", "[", 
      RowBox[{"[", "2", "]"}], "]"}]}], "*)"}], "\n", 
   RowBox[{"vir", "=", 
    RowBox[{
     RowBox[{"-", "Epot"}], "/", "Ekin"}]}], "  ", 
   RowBox[{"(*", 
    RowBox[{"/.", " ", 
     RowBox[{"t", "[", 
      RowBox[{"[", "2", "]"}], "]"}]}], "*)"}], "\n", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"R1", "+", "R2"}], ")"}], "*", "p", "*", "bohr", Cell[""]}], 
   "  ", 
   RowBox[{"(*", 
    RowBox[{"/.", " ", 
     RowBox[{"t", "[", 
      RowBox[{"[", "2", "]"}], "]"}]}], "*)"}], "\n"}]}]], "Input",
 CellChangeTimes->{{3.5784163520300827`*^9, 3.578416352045683*^9}, {
  3.5784167316567497`*^9, 3.578416866097786*^9}, {3.5784168996066446`*^9, 
  3.5784169997968206`*^9}, {3.5784170630237317`*^9, 3.5784170940833864`*^9}, {
  3.5784172056703825`*^9, 3.578417273873702*^9}, {3.578417315619375*^9, 
  3.5784173253381925`*^9}, {3.5784175560625973`*^9, 3.578417678086012*^9}, {
  3.578417720112486*^9, 3.5784177210328875`*^9}, {3.5784183246447477`*^9, 
  3.5784183318675604`*^9}}],

Cell[BoxData["40.56089216702621`"], "Output",
 CellChangeTimes->{3.578417120743833*^9, 3.5784172816581154`*^9, 
  3.5784173357902107`*^9, 3.5784177356813135`*^9, 3.578418067306696*^9, 
  3.578418192886916*^9, 3.578418291213889*^9, 3.578418340213575*^9}],

Cell[BoxData[
 RowBox[{"-", "81.12178443260031`"}]], "Output",
 CellChangeTimes->{3.578417120743833*^9, 3.5784172816581154`*^9, 
  3.5784173357902107`*^9, 3.5784177356813135`*^9, 3.578418067306696*^9, 
  3.578418192886916*^9, 3.578418291213889*^9, 3.578418340213575*^9}],

Cell[BoxData[
 RowBox[{"-", "40.5608922655741`"}]], "Output",
 CellChangeTimes->{3.578417120743833*^9, 3.5784172816581154`*^9, 
  3.5784173357902107`*^9, 3.5784177356813135`*^9, 3.578418067306696*^9, 
  3.578418192886916*^9, 3.578418291213889*^9, 3.578418340213575*^9}],

Cell[BoxData[
 RowBox[{"-", "121.15886615975022`"}]], "Output",
 CellChangeTimes->{3.578417120743833*^9, 3.5784172816581154`*^9, 
  3.5784173357902107`*^9, 3.5784177356813135`*^9, 3.578418067306696*^9, 
  3.578418192886916*^9, 3.578418291213889*^9, 3.5784183402291746`*^9}],

Cell[BoxData["26.425380915318428`"], "Output",
 CellChangeTimes->{3.578417120743833*^9, 3.5784172816581154`*^9, 
  3.5784173357902107`*^9, 3.5784177356813135`*^9, 3.578418067306696*^9, 
  3.578418192886916*^9, 3.578418291213889*^9, 3.5784183402291746`*^9}],

Cell[BoxData["13.61170081183149`"], "Output",
 CellChangeTimes->{3.578417120743833*^9, 3.5784172816581154`*^9, 
  3.5784173357902107`*^9, 3.5784177356813135`*^9, 3.578418067306696*^9, 
  3.578418192886916*^9, 3.578418291213889*^9, 3.5784183402291746`*^9}],

Cell[BoxData["2.000000002429628`"], "Output",
 CellChangeTimes->{3.578417120743833*^9, 3.5784172816581154`*^9, 
  3.5784173357902107`*^9, 3.5784177356813135`*^9, 3.578418067306696*^9, 
  3.578418192886916*^9, 3.578418291213889*^9, 3.5784183402291746`*^9}],

Cell[BoxData[
 RowBox[{"1.075880865504764`", " ", 
  InterpretationBox[Cell[""],
   TextCell[""]]}]], "Output",
 CellChangeTimes->{3.578417120743833*^9, 3.5784172816581154`*^9, 
  3.5784173357902107`*^9, 3.5784177356813135`*^9, 3.578418067306696*^9, 
  3.578418192886916*^9, 3.578418291213889*^9, 3.5784183402291746`*^9}]
}, Open  ]]
},
WindowToolbars->"EditBar",
WindowSize->{875, 964},
WindowMargins->{{Automatic, 451}, {Automatic, 4}},
PrintingCopies->1,
PrintingPageRange->{Automatic, Automatic},
PrivateNotebookOptions->{"VersionedStylesheet"->{"Default.nb"[8.] -> False}},
FrontEndVersion->"9.0 for Microsoft Windows (64-bit) (January 25, 2013)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[557, 20, 473, 8, 30, "Text"],
Cell[1033, 30, 866, 15, 52, "Input"],
Cell[1902, 47, 741, 13, 31, "Input"],
Cell[2646, 62, 1788, 43, 132, "Input"],
Cell[4437, 107, 426, 9, 31, "Input"],
Cell[CellGroupData[{
Cell[4888, 120, 3852, 86, 112, "Input"],
Cell[8743, 208, 426, 6, 31, "Output"],
Cell[9172, 216, 426, 6, 31, "Output"]
}, Open  ]],
Cell[9613, 225, 2571, 41, 31, "Input"],
Cell[CellGroupData[{
Cell[12209, 270, 2070, 41, 52, "Input"],
Cell[14282, 313, 1381, 32, 101, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[15700, 350, 2068, 41, 52, "Input"],
Cell[17771, 393, 1388, 32, 101, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[19196, 430, 2072, 39, 52, "Input"],
Cell[21271, 471, 1485, 31, 101, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[22793, 507, 2071, 38, 52, "Input"],
Cell[24867, 547, 1490, 31, 101, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[26394, 583, 2987, 57, 92, "Input"],
Cell[29384, 642, 1560, 31, 101, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[30981, 678, 2281, 43, 52, "Input"],
Cell[33265, 723, 1486, 31, 101, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[34788, 759, 2231, 43, 52, "Input"],
Cell[37022, 804, 1476, 31, 101, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[38535, 840, 2376, 46, 52, "Input"],
Cell[40914, 888, 1580, 31, 101, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[42531, 924, 2942, 56, 112, "Input"],
Cell[45476, 982, 1251, 33, 101, "Output"],
Cell[46730, 1017, 449, 7, 31, "Output"]
}, Open  ]],
Cell[47194, 1027, 2856, 62, 132, "Input"],
Cell[50053, 1091, 2455, 47, 92, "Input"],
Cell[CellGroupData[{
Cell[52533, 1142, 1745, 40, 172, "Input"],
Cell[54281, 1184, 448, 7, 31, "Output"],
Cell[54732, 1193, 432, 6, 31, "Output"],
Cell[55167, 1201, 450, 7, 31, "Output"],
Cell[55620, 1210, 433, 6, 31, "Output"],
Cell[56056, 1218, 433, 6, 31, "Output"],
Cell[56492, 1226, 432, 6, 31, "Output"],
Cell[56927, 1234, 432, 6, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[57396, 1245, 4571, 120, 292, "Input"],
Cell[61970, 1367, 1773, 42, 401, "Output"],
Cell[63746, 1411, 1120, 20, 379, "Output"]
}, Open  ]],
Cell[64881, 1434, 935, 23, 72, "Input"],
Cell[CellGroupData[{
Cell[65841, 1461, 557, 11, 31, "Input"],
Cell[66401, 1474, 198, 3, 31, "Output"]
}, Open  ]],
Cell[66614, 1480, 1222, 33, 92, "Input"],
Cell[CellGroupData[{
Cell[67861, 1517, 873, 22, 31, "Input"],
Cell[68737, 1541, 201, 3, 31, "Output"]
}, Open  ]],
Cell[68953, 1547, 796, 20, 52, "Input"],
Cell[CellGroupData[{
Cell[69774, 1571, 556, 11, 31, "Input"],
Cell[70333, 1584, 200, 3, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[70570, 1592, 581, 12, 31, "Input"],
Cell[71154, 1606, 202, 3, 31, "Output"]
}, Open  ]],
Cell[71371, 1612, 772, 19, 52, "Input"],
Cell[CellGroupData[{
Cell[72168, 1635, 730, 15, 31, "Input"],
Cell[72901, 1652, 200, 3, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[73138, 1660, 1840, 47, 272, "Input"],
Cell[74981, 1709, 253, 3, 31, "Output"],
Cell[75237, 1714, 270, 4, 31, "Output"],
Cell[75510, 1720, 269, 4, 31, "Output"],
Cell[75782, 1726, 273, 4, 31, "Output"],
Cell[76058, 1732, 256, 3, 31, "Output"],
Cell[76317, 1737, 255, 3, 31, "Output"],
Cell[76575, 1742, 255, 3, 31, "Output"],
Cell[76833, 1747, 321, 6, 63, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)