(* 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[     35713,        821]
NotebookOptionsPosition[     34310,        774]
NotebookOutlinePosition[     34915,        796]
CellTagsIndexPosition[     34872,        793]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell["\<\
(* CH4   25.05.2013; all interactions with matrix operations; e.g. \
-4*(3-(p/R2)^2)/R2 for H atoms
 is included in dwnr  *)\
\>", "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, {
   3.5784553214077992`*^9, 3.578455336883026*^9}, {3.578461864956892*^9, 
   3.5784618814461207`*^9}, {3.5784625589397106`*^9, 3.578462637454649*^9}, {
   3.5784819580459146`*^9, 3.5784819597931175`*^9}, 3.578482692495205*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"Clear", "[", 
   RowBox[{
   "z", ",", "sig", ",", "k", ",", "nc", ",", "R1", ",", "R2", ",", "p", ",", 
    "Ekin", ",", "Epot", ",", "vne", ",", "vnn", ",", "vee", ",", "ow", ",", 
    "cn", ",", "n", ",", "w", ",", "c", ",", "c2", ",", "ec", ",", "dw", ",", 
    "dn", ",", "dwn", ",", "dwnr", ",", "dwr", ",", "dnr", ",", "d1", ",", 
    "d2", ",", "pi", ",", "i", ",", "j", ",", "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, {3.5784594637094746`*^9, 
   3.578459599913314*^9}, {3.5784609547444935`*^9, 3.578461001996976*^9}, {
   3.5784612573826246`*^9, 3.578461273450653*^9}}],

Cell[BoxData[
 RowBox[{"(*", 
  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}, 3.5784586223399963`*^9, {3.578459802323669*^9, 
   3.5784598064576764`*^9}, 3.5784623423177304`*^9}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"z", "=", "6.0"}], ";"}], "\n", 
 RowBox[{
  RowBox[{"nc", "=", "5"}], ";", " ", 
  RowBox[{"bohr", "=", "0.529177"}], ";", 
  RowBox[{"sq", "=", 
   RowBox[{"Sqrt", "[", "3", "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"k", "=", 
    RowBox[{"{", 
     RowBox[{
     "1.04605075", ",", "1.20889711", ",", "1.20889711", ",", "1.20889711", 
      ",", "1.20889711"}], "}"}]}], ";"}], " "}], "\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}, 
   3.5784587030521383`*^9, {3.578481973630342*^9, 3.5784819824287577`*^9}, 
   3.578482310544134*^9}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"dw", "=", 
   RowBox[{"Table", "[", 
    RowBox[{"0", ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "nc"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"j", ",", "1", ",", "nc"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"dn", "=", "dw"}], ";", " ", 
  RowBox[{"dwn", "=", "dw"}], ";", " ", 
  RowBox[{"dwr", "=", "dw"}], ";", " ", 
  RowBox[{"dnr", "=", "dw"}], ";", 
  RowBox[{"dwnr", "=", "dw"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"c2", "=", 
   RowBox[{"Table", "[", 
    RowBox[{"0", ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", "1", ",", "nc"}], "}"}]}], "]"}]}], ";"}]}], "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}, {3.5784554299839897`*^9, 3.578455463695649*^9}, {
   3.578461345319979*^9, 3.5784613583616023`*^9}, {3.578480998503829*^9, 
   3.5784811442236853`*^9}, {3.5784814608262415`*^9, 3.5784814881262894`*^9}}],

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"}], "}"}]}], 
  ";"}]}], "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}, {3.578454093592043*^9, 
   3.5784541407977257`*^9}, {3.578455373418291*^9, 3.5784553903911204`*^9}, {
   3.5784554730556655`*^9, 3.578455482384482*^9}, {3.578456998330344*^9, 
   3.5784570025111513`*^9}, {3.578457293950864*^9, 3.5784573307357283`*^9}, {
   3.578458063921416*^9, 3.5784580686794243`*^9}, 3.5784590454259396`*^9, {
   3.5784592151386375`*^9, 3.5784592662443275`*^9}, {3.5784596927334766`*^9, 
   3.578459697023484*^9}, {3.578459839498534*^9, 3.57845989952744*^9}, {
   3.5784599603987465`*^9, 3.578459968651161*^9}, 3.578460029319668*^9, 
   3.578460147583475*^9, {3.578460262025276*^9, 3.5784604020355225`*^9}, 
   3.578460573261423*^9, {3.5784619324270105`*^9, 3.578462010099547*^9}, {
   3.5784622442403584`*^9, 3.5784622854868307`*^9}, {3.578462453093525*^9, 
   3.578462453670726*^9}, {3.578479608042187*^9, 3.5784797011587505`*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"c", "=", 
   RowBox[{"{", 
    RowBox[{"R1", ",", "R2", ",", "R2", ",", "R2", ",", "R2"}], "}"}]}], 
  ";"}], "\n", 
 RowBox[{
  RowBox[{"ec", "=", 
   RowBox[{"{", 
    RowBox[{"0", ",", "p", ",", "p", ",", "p", ",", "p"}], "}"}]}], 
  ";"}], "\n", 
 RowBox[{
  RowBox[{"d1", "=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"R1", "+", "R2", "+", "p"}], ")"}], "/", "sq"}]}], ";", "  ", 
  RowBox[{"d2", "=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"R1", "+", "R2"}], ")"}], "/", "sq"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"n", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"0", ",", "0", ",", "0"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"d1", ",", "d1", ",", "d1"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"-", "d1"}], ",", 
       RowBox[{"-", "d1"}], ",", "d1"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"d1", ",", 
       RowBox[{"-", "d1"}], ",", 
       RowBox[{"-", "d1"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"-", "d1"}], ",", "d1", ",", 
       RowBox[{"-", "d1"}]}], "}"}]}], "}"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"w", "=", 
   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[{
  RowBox[{
   RowBox[{"For", "[", 
    RowBox[{
     RowBox[{"i", "=", "1"}], ",", 
     RowBox[{"i", "\[LessEqual]", "nc"}], ",", 
     RowBox[{"i", "++"}], ",", 
     RowBox[{
      RowBox[{
       RowBox[{"c2", "[", 
        RowBox[{"[", "i", "]"}], "]"}], "=", 
       RowBox[{
        RowBox[{"k", "[", 
         RowBox[{"[", "i", "]"}], "]"}], "/", 
        RowBox[{
         RowBox[{"c", "[", 
          RowBox[{"[", "i", "]"}], "]"}], "^", "2"}]}]}], ";", 
      "\[IndentingNewLine]", 
      RowBox[{"For", "[", 
       RowBox[{
        RowBox[{"j", "=", "1"}], ",", 
        RowBox[{"j", "\[LessEqual]", "nc"}], ",", 
        RowBox[{"j", "++"}], ",", "\n", 
        RowBox[{
         RowBox[{
          RowBox[{"dw", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "j"}], "]"}], "]"}], "=", 
          RowBox[{"Norm", "[", 
           RowBox[{
            RowBox[{"w", "[", 
             RowBox[{"[", "i", "]"}], "]"}], "-", 
            RowBox[{"w", "[", 
             RowBox[{"[", "j", "]"}], "]"}]}], "]"}]}], ";", "\n", 
         RowBox[{
          RowBox[{"dn", " ", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "j"}], "]"}], "]"}], "=", 
          RowBox[{"Norm", "[", 
           RowBox[{
            RowBox[{"n", "[", 
             RowBox[{"[", "i", "]"}], "]"}], "-", 
            RowBox[{"n", "[", 
             RowBox[{"[", "j", "]"}], "]"}]}], "]"}]}], ";", "\n", 
         RowBox[{
          RowBox[{"dwn", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "j"}], "]"}], "]"}], "=", 
          RowBox[{"Norm", "[", 
           RowBox[{
            RowBox[{"w", "[", 
             RowBox[{"[", "i", "]"}], "]"}], "-", 
            RowBox[{"n", "[", 
             RowBox[{"[", "j", "]"}], "]"}]}], "]"}]}], ";", "\n", 
         RowBox[{
          RowBox[{"dwr", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "j"}], "]"}], "]"}], "=", 
          RowBox[{"If", "[", 
           RowBox[{
            RowBox[{"i", "\[NotEqual]", "j"}], ",", 
            RowBox[{"1", "/", 
             RowBox[{"dw", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "j"}], "]"}], "]"}]}], ",", 
            RowBox[{"1.5", "*", 
             RowBox[{
              RowBox[{"sig", "[", 
               RowBox[{"[", "i", "]"}], "]"}], "/", 
              RowBox[{"c", "[", 
               RowBox[{"[", "i", "]"}], "]"}]}]}]}], "]"}]}], ";", "\n", 
         RowBox[{
          RowBox[{"dnr", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "j"}], "]"}], "]"}], "=", 
          RowBox[{"If", "[", 
           RowBox[{
            RowBox[{"i", "\[NotEqual]", "j"}], ",", 
            RowBox[{"1", "/", 
             RowBox[{"dn", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "j"}], "]"}], "]"}]}], ",", "0"}], "]"}]}], 
         ";", "\n", 
         RowBox[{
          RowBox[{"dwnr", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "j"}], "]"}], "]"}], "=", 
          RowBox[{"If", "[", 
           RowBox[{
            RowBox[{"i", "\[NotEqual]", "j"}], ",", 
            RowBox[{"1", "/", 
             RowBox[{"dwn", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "j"}], "]"}], "]"}]}], ",", 
            RowBox[{
             RowBox[{"(", 
              RowBox[{"1.5", "-", 
               RowBox[{"0.5", "*", 
                RowBox[{
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{"ec", "[", 
                    RowBox[{"[", "i", "]"}], "]"}], "/", 
                   RowBox[{"c", "[", 
                    RowBox[{"[", "i", "]"}], "]"}]}], ")"}], "^", "2"}]}]}], 
              ")"}], "/", 
             RowBox[{"c", "[", 
              RowBox[{"[", "i", "]"}], "]"}]}]}], "]"}]}]}]}], "]"}]}]}], 
    "]"}], ";"}], "\[IndentingNewLine]"}], "\n", 
 RowBox[{
  RowBox[{"Ekin", "=", 
   RowBox[{
    RowBox[{"-", 
     RowBox[{"(", 
      RowBox[{"9", "/", "8"}], ")"}]}], "*", 
    RowBox[{"(", 
     RowBox[{"ow", ".", "c2"}], ")"}]}]}], ";"}], "\n", 
 RowBox[{
  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[{"vne", "=", 
   RowBox[{"cn", ".", 
    RowBox[{"(", 
     RowBox[{"ow", ".", "dwnr"}], ")"}]}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"Epot", "=", 
   RowBox[{"vne", "+", "vee", "+", "vnn"}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"func", "=", 
   RowBox[{"Ekin", "+", "Epot"}]}], ";"}]}], "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}, {3.578454093592043*^9, 
   3.5784541407977257`*^9}, {3.578455373418291*^9, 3.5784553903911204`*^9}, {
   3.5784554730556655`*^9, 3.578455482384482*^9}, {3.578456998330344*^9, 
   3.5784570025111513`*^9}, {3.578457293950864*^9, 3.5784573307357283`*^9}, {
   3.578458063921416*^9, 3.5784580686794243`*^9}, 3.5784590454259396`*^9, {
   3.5784592151386375`*^9, 3.5784592662443275`*^9}, {3.5784596927334766`*^9, 
   3.578459697023484*^9}, {3.578459839498534*^9, 3.57845989952744*^9}, {
   3.5784599603987465`*^9, 3.578459968651161*^9}, 3.578460029319668*^9, 
   3.578460147583475*^9, {3.578460262025276*^9, 3.5784604020355225`*^9}, 
   3.578460573261423*^9, {3.5784619324270105`*^9, 3.578462010099547*^9}, {
   3.5784622442403584`*^9, 3.5784622854868307`*^9}, {3.578462453093525*^9, 
   3.578462453670726*^9}, {3.578479608042187*^9, 3.578479714013173*^9}, 
   3.5784797501896367`*^9, {3.578479970758424*^9, 3.5784799822244444`*^9}, 
   3.578480246878909*^9, {3.5784802810897694`*^9, 3.5784803143178277`*^9}, {
   3.578480533950613*^9, 3.578480636520794*^9}, {3.5784807191541386`*^9, 
   3.578480924388099*^9}, {3.578480956180955*^9, 3.578480956539756*^9}, {
   3.5784812251254272`*^9, 3.57848128353193*^9}, {3.578481359925264*^9, 
   3.5784813929661217`*^9}, {3.5784817072754745`*^9, 
   3.5784817101146793`*^9}, {3.57848235969982*^9, 3.5784824488851767`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"t", " ", "=", " ", 
  RowBox[{"FindMinimum", "[", 
   RowBox[{"func", ",", 
    RowBox[{"{", 
     RowBox[{"R1", ",", "0.265"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"R2", ",", "1.23"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"p", ",", "0.53"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"Method", " ", "\[Rule]", "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.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}, {3.578454093592043*^9, 
   3.5784541407977257`*^9}, {3.578455373418291*^9, 3.5784553903911204`*^9}, {
   3.5784554730556655`*^9, 3.578455482384482*^9}, {3.578456998330344*^9, 
   3.5784570025111513`*^9}, {3.578457293950864*^9, 3.5784573307357283`*^9}, {
   3.578458063921416*^9, 3.5784580686794243`*^9}, 3.5784590454259396`*^9, {
   3.5784592151386375`*^9, 3.5784592662443275`*^9}, {3.5784596927334766`*^9, 
   3.578459697023484*^9}, {3.578459839498534*^9, 3.57845989952744*^9}, {
   3.5784599603987465`*^9, 3.578459968651161*^9}, 3.578460029319668*^9, 
   3.578460147583475*^9, {3.578460262025276*^9, 3.5784604020355225`*^9}, 
   3.578460573261423*^9, {3.5784619324270105`*^9, 3.578462010099547*^9}, {
   3.5784622442403584`*^9, 3.5784622854868307`*^9}, {3.578462453093525*^9, 
   3.578462453670726*^9}, {3.578479608042187*^9, 3.5784796373858385`*^9}, {
   3.578481907111825*^9, 3.5784819149742393`*^9}, {3.5784821295682163`*^9, 
   3.5784821328754215`*^9}, 3.5784822034499454`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", "40.56089226557407`"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"R1", "\[Rule]", "0.26537991410835604`"}], ",", 
     RowBox[{"R2", "\[Rule]", "1.2343054251474905`"}], ",", 
     RowBox[{"p", "\[Rule]", "0.5334356944092935`"}]}], "}"}]}], 
  "}"}]], "Output",
 CellChangeTimes->{3.5784825261209126`*^9, 3.578482584387015*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"Ekin", "    ", "/.", " ", 
   RowBox[{"t", "[", 
    RowBox[{"[", "2", "]"}], "]"}]}], "  "}], "\n", 
 RowBox[{
  RowBox[{"vne", "     ", "/.", " ", 
   RowBox[{"t", "[", 
    RowBox[{"[", "2", "]"}], "]"}]}], "  "}], "\n", 
 RowBox[{
  RowBox[{"vee", "     ", "/.", " ", 
   RowBox[{"t", "[", 
    RowBox[{"[", "2", "]"}], "]"}]}], "  "}], "\n", 
 RowBox[{
  RowBox[{"vnn", "     ", "/.", " ", 
   RowBox[{"t", "[", 
    RowBox[{"[", "2", "]"}], "]"}]}], " ", " "}], "\n", 
 RowBox[{
  RowBox[{"vir", "=", 
   RowBox[{
    RowBox[{
     RowBox[{"-", "Epot"}], "/", "Ekin"}], "     ", "/.", " ", 
    RowBox[{"t", "[", 
     RowBox[{"[", "2", "]"}], "]"}]}]}], "  "}], "\n", 
 RowBox[{
  RowBox[{"d1", "*", "sq", "*", "bohr"}], "     ", "/.", " ", 
  RowBox[{"t", "[", 
   RowBox[{"[", "2", "]"}], "]"}], " ", " "}]}], "Input",
 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}, {
   3.5784584627517166`*^9, 3.5784584719089327`*^9}, 3.5784585415318546`*^9, {
   3.5784585720767083`*^9, 3.578458572092308*^9}, {3.578460080503358*^9, 
   3.5784600805189576`*^9}, {3.5784797917169094`*^9, 
   3.5784797917325096`*^9}, {3.5784820752021203`*^9, 3.578482083704135*^9}, {
   3.5784824699608135`*^9, 3.578482470148014*^9}, {3.578482574324997*^9, 
   3.578482601063444*^9}}],

Cell[BoxData["40.560893359905116`"], "Output",
 CellChangeTimes->{3.5784825261365128`*^9, 3.5784825844026146`*^9}],

Cell[BoxData[
 RowBox[{"-", "121.15886679581777`"}]], "Output",
 CellChangeTimes->{3.5784825261365128`*^9, 3.5784825844026146`*^9}],

Cell[BoxData["26.42538065869418`"], "Output",
 CellChangeTimes->{3.5784825261365128`*^9, 3.5784825844026146`*^9}],

Cell[BoxData["13.611700511644393`"], "Output",
 CellChangeTimes->{3.5784825261365128`*^9, 3.5784825844026146`*^9}],

Cell[BoxData["1.9999999730200457`"], "Output",
 CellChangeTimes->{3.5784825261365128`*^9, 3.5784825844026146`*^9}],

Cell[BoxData["1.075880889231818`"], "Output",
 CellChangeTimes->{3.5784825261365128`*^9, 3.578482584418215*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"plot1", "=", 
   RowBox[{
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"Graphics3D", "[", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"Opacity", "[", "0.6", "]"}], ",", 
         RowBox[{"Sphere", "[", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{"w", "[", 
             RowBox[{"[", "i", "]"}], "]"}], "}"}], ",", 
           RowBox[{"c", "[", 
            RowBox[{"[", "i", "]"}], "]"}]}], "]"}], ",", 
         RowBox[{"{", 
          RowBox[{
           RowBox[{"Darker", "[", 
            RowBox[{"Red", ",", "1.0"}], "]"}], ",", 
           RowBox[{"Sphere", "[", 
            RowBox[{
             RowBox[{"{", 
              RowBox[{"n", "[", 
               RowBox[{"[", "i", "]"}], "]"}], "}"}], ",", "0.08"}], "]"}]}], 
          "}"}]}], "}"}], " ", "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "1", ",", "nc"}], "}"}]}], "]"}], "   ", "/.", " ", 
    RowBox[{"t", "[", 
     RowBox[{"[", "2", "]"}], "]"}]}]}], " ", " ", ";"}], "\n", 
 RowBox[{"Show", "[", 
  RowBox[{"plot1", ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"AspectRatio", " ", "\[Rule]", " ", "Automatic"}], ",", 
     RowBox[{"Axes", " ", "->", " ", "True"}]}], "}"}]}], "]"}]}], "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}, {3.578458572092308*^9, 3.5784585723263087`*^9}, {
   3.578458946662166*^9, 3.578458947800968*^9}, {3.5784600805345573`*^9, 
   3.5784600807529583`*^9}, {3.5784605921686563`*^9, 3.578460617471901*^9}, {
   3.578460728622096*^9, 3.5784608330798798`*^9}, 3.5784610501074605`*^9, {
   3.57847953008885*^9, 3.5784795603997035`*^9}, 3.57847979174811*^9}],

Cell[BoxData[
 Graphics3DBox[{
   {Opacity[0.6], SphereBox[{{0, 0, 0}}, 0.26537991410835604`], 
    {RGBColor[0., 0., 0.], SphereBox[{{0, 0, 0}}, 0.08]}}, 
   {Opacity[0.6], 
    SphereBox[{{0.8658437343190983, 0.8658437343190983, 0.8658437343190983}}, 
     1.2343054251474905`], 
    {RGBColor[0., 0., 0.], 
     SphereBox[{{1.173822976081659, 1.173822976081659, 1.173822976081659}}, 
      0.08]}}, 
   {Opacity[0.6], 
    SphereBox[{{-0.8658437343190983, -0.8658437343190983, 
     0.8658437343190983}}, 1.2343054251474905`], 
    {RGBColor[0., 0., 0.], 
     SphereBox[{{-1.173822976081659, -1.173822976081659, 1.173822976081659}}, 
      0.08]}}, 
   {Opacity[0.6], 
    SphereBox[{{
     0.8658437343190983, -0.8658437343190983, -0.8658437343190983}}, 
     1.2343054251474905`], 
    {RGBColor[0., 0., 0.], 
     SphereBox[{{1.173822976081659, -1.173822976081659, -1.173822976081659}}, 
      0.08]}}, 
   {Opacity[0.6], 
    SphereBox[{{-0.8658437343190983, 
     0.8658437343190983, -0.8658437343190983}}, 1.2343054251474905`], 
    {RGBColor[0., 0., 0.], 
     SphereBox[{{-1.173822976081659, 1.173822976081659, -1.173822976081659}}, 
      0.08]}}},
  AspectRatio->Automatic,
  Axes->True]], "Output",
 CellChangeTimes->{3.5784825261521125`*^9, 3.578482584433815*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"plot2", "=", 
    RowBox[{
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"Graphics", "[", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"Circle", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{"w", "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", "1"}], "]"}], "]"}], ",", 
              RowBox[{"w", "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", "2"}], "]"}], "]"}]}], "}"}], ",", 
            RowBox[{"c", "[", 
             RowBox[{"[", "i", "]"}], "]"}]}], "]"}], ",", 
          RowBox[{"Disk", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{"n", "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", "1"}], "]"}], "]"}], ",", 
              RowBox[{"n", "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", "2"}], "]"}], "]"}]}], "}"}], ",", "0.08"}],
            "]"}]}], "}"}], " ", "]"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "1", ",", "nc"}], "}"}]}], "]"}], "   ", "/.", " ", 
     RowBox[{"t", "[", 
      RowBox[{"[", "2", "]"}], "]"}]}]}], " ", " ", ";"}], "  "}], "\n", 
 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"}]}], "}"}]}], "]"}]}], "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}, {3.578458572092308*^9, 3.5784585723263087`*^9}, {
   3.578458946662166*^9, 3.578458947800968*^9}, {3.5784600805345573`*^9, 
   3.5784600807529583`*^9}, {3.5784605921686563`*^9, 3.578460617471901*^9}, {
   3.578460728622096*^9, 3.5784608330798798`*^9}, 3.5784610501074605`*^9, {
   3.57847953008885*^9, 3.578479548138082*^9}, 3.57847979174811*^9}],

Cell[BoxData[
 GraphicsBox[{{CircleBox[{0, 0}, 0.26537991410835604`], 
    DiskBox[{0, 0}, 0.08]}, {
    CircleBox[{0.8658437343190983, 0.8658437343190983}, 1.2343054251474905`], 
    DiskBox[{1.173822976081659, 1.173822976081659}, 0.08]}, {
    CircleBox[{-0.8658437343190983, -0.8658437343190983}, 
     1.2343054251474905`], 
    DiskBox[{-1.173822976081659, -1.173822976081659}, 0.08]}, {
    CircleBox[{0.8658437343190983, -0.8658437343190983}, 1.2343054251474905`],
     DiskBox[{1.173822976081659, -1.173822976081659}, 0.08]}, {
    CircleBox[{-0.8658437343190983, 0.8658437343190983}, 1.2343054251474905`],
     DiskBox[{-1.173822976081659, 1.173822976081659}, 0.08]}},
  AspectRatio->Automatic,
  Axes->True,
  Frame->True,
  GridLines->Automatic,
  PlotRange->{{-3, 3}, {-3, 3}}]], "Output",
 CellChangeTimes->{3.578482526214513*^9, 3.578482584496215*^9}]
}, Open  ]]
},
WindowToolbars->"EditBar",
WindowSize->{837, 831},
WindowMargins->{{Automatic, 269}, {Automatic, 13}},
PrintingCopies->1,
PrintingPageRange->{32000, 32000},
PrintingOptions->{"Magnification"->1.,
"PaperOrientation"->"Portrait",
"PaperSize"->{595.2, 841.8}},
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, 697, 12, 49, "Text"],
Cell[1257, 34, 897, 15, 52, "Input"],
Cell[2157, 51, 869, 15, 31, "Input"],
Cell[3029, 68, 1452, 32, 92, "Input"],
Cell[4484, 102, 1166, 28, 72, "Input"],
Cell[5653, 132, 3540, 57, 52, "Input"],
Cell[9196, 191, 10123, 249, 412, "Input"],
Cell[CellGroupData[{
Cell[19344, 444, 3771, 59, 52, "Input"],
Cell[23118, 505, 390, 10, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[23545, 520, 2002, 43, 132, "Input"],
Cell[25550, 565, 114, 1, 31, "Output"],
Cell[25667, 568, 131, 2, 31, "Output"],
Cell[25801, 572, 113, 1, 31, "Output"],
Cell[25917, 575, 114, 1, 31, "Output"],
Cell[26034, 578, 114, 1, 31, "Output"],
Cell[26151, 581, 111, 1, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[26299, 587, 2595, 55, 92, "Input"],
Cell[28897, 644, 1279, 31, 422, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[30213, 680, 3213, 72, 112, "Input"],
Cell[33429, 754, 865, 17, 379, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)