(* 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[     39801,        905]
NotebookOptionsPosition[     38365,        857]
NotebookOutlinePosition[     38969,        879]
CellTagsIndexPosition[     38926,        876]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell["\<\
(* CH4   28.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, {
   3.578720542295065*^9, 3.5787205428566666`*^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", ",", "\[Lambda]"}], "]"}], 
  ";"}]], "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}, {3.5787175511146116`*^9, 
   3.5787175547494183`*^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[{
   RowBox[{"c", " ", "=", " ", 
    RowBox[{"{", 
     RowBox[{
     "R1", ",", " ", "R2", ",", " ", "R2", ",", " ", "R2", ",", " ", "R2"}], 
     "}"}]}], ";"}], " "}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"ec", " ", "=", " ", 
    RowBox[{"{", 
     RowBox[{
     "0", ",", " ", "p", ",", " ", "p", ",", " ", "p", ",", " ", "p"}], 
     "}"}]}], ";"}], " "}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"d1", " ", "=", " ", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"R1", " ", "+", " ", "R2", " ", "+", " ", "p"}], ")"}], "/", 
     "sq"}]}], ";", " ", 
   RowBox[{"d2", " ", "=", " ", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"R1", " ", "+", " ", "R2"}], ")"}], "/", "sq"}]}], ";"}], 
  " "}], "\n", 
 RowBox[{
  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[{
   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"}]}], "}"}]}], "}"}]}], ";"}], " "}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"For", "[", 
    RowBox[{
     RowBox[{"i", " ", "=", " ", "1"}], ",", " ", 
     RowBox[{"i", " ", "<=", " ", "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", " ", "<=", " ", "nc"}], ",", " ", 
        RowBox[{"j", "++"}], ",", " ", "\n", "    ", 
        RowBox[{
         RowBox[{
          RowBox[{"dw", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "j"}], "]"}], "]"}], " ", "=", " ", 
          RowBox[{"Norm", "[", 
           RowBox[{
            RowBox[{"w", "[", 
             RowBox[{"[", "i", "]"}], "]"}], " ", "-", " ", 
            RowBox[{"w", "[", 
             RowBox[{"[", "j", "]"}], "]"}]}], "]"}]}], ";", " ", 
         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", "]"}], "]"}]}], "]"}]}], ";", " ", 
         RowBox[{
          RowBox[{"dwr", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "j"}], "]"}], "]"}], " ", "=", " ", 
          RowBox[{"If", "[", 
           RowBox[{
            RowBox[{"i", " ", "!=", " ", "j"}], ",", " ", 
            RowBox[{"1", "/", 
             RowBox[{"dw", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "j"}], "]"}], "]"}]}], ",", " ", "\n", 
            "        ", 
            RowBox[{"1.5", "*", 
             RowBox[{"(", 
              RowBox[{
               RowBox[{"sig", "[", 
                RowBox[{"[", "i", "]"}], "]"}], "/", 
               RowBox[{"c", "[", 
                RowBox[{"[", "i", "]"}], "]"}]}], ")"}]}]}], "]"}]}], ";", 
         " ", 
         RowBox[{
          RowBox[{"dnr", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "j"}], "]"}], "]"}], " ", "=", " ", 
          RowBox[{"If", "[", 
           RowBox[{
            RowBox[{"i", " ", "!=", " ", "j"}], ",", " ", 
            RowBox[{"1", "/", 
             RowBox[{"dn", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "j"}], "]"}], "]"}]}], ",", " ", "0"}], 
           "]"}]}], ";", " ", "\n", "    ", 
         RowBox[{
          RowBox[{"dwnr", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", "j"}], "]"}], "]"}], " ", "=", " ", 
          RowBox[{"If", "[", 
           RowBox[{
            RowBox[{"i", " ", "!=", " ", "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", "]"}], "]"}]}]}], "]"}]}]}]}], "]"}]}]}], 
    "]"}], ";"}], " "}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"Ekin", " ", "=", " ", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"-", 
       RowBox[{"(", 
        RowBox[{"9", "/", "8"}], ")"}]}], ")"}], "*", 
     RowBox[{"ow", " ", ".", " ", "c2"}]}]}], ";"}], " "}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"vee", " ", "=", " ", 
    RowBox[{"0.5", "*", 
     RowBox[{"ow", " ", ".", " ", 
      RowBox[{"(", 
       RowBox[{"ow", " ", ".", " ", "dwr"}], ")"}]}]}]}], ";"}], " "}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"vnn", " ", "=", " ", 
    RowBox[{"0.5", "*", 
     RowBox[{"cn", " ", ".", " ", 
      RowBox[{"(", 
       RowBox[{"cn", " ", ".", " ", "dnr"}], ")"}]}]}]}], ";"}], " "}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"vne", " ", "=", " ", 
    RowBox[{"cn", " ", ".", " ", 
     RowBox[{"(", 
      RowBox[{"ow", " ", ".", " ", "dwnr"}], ")"}]}]}], ";"}], " "}], "\n", 
 RowBox[{
  RowBox[{
   RowBox[{"Epot", " ", "=", " ", 
    RowBox[{"vne", " ", "+", " ", "vee", " ", "+", " ", "vnn"}]}], ";"}], 
  " "}], "\n", 
 RowBox[{
  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}, {
   3.5786677972082167`*^9, 3.5786678267702684`*^9}, 3.578667911260017*^9, {
   3.5786679871229506`*^9, 3.578668026169819*^9}, {3.578668087041126*^9, 
   3.5786680886791286`*^9}, {3.578668159472053*^9, 3.578668179845689*^9}, {
   3.578668332086356*^9, 3.578668361976009*^9}, {3.578668412691698*^9, 
   3.578668413346899*^9}, {3.5786685853060007`*^9, 3.5786686488449125`*^9}, {
   3.5786687638483143`*^9, 3.5786687758447356`*^9}, {3.578715368421178*^9, 
   3.5787153700123806`*^9}, {3.5787154222100725`*^9, 3.5787154258136787`*^9}, 
   3.5787155896451664`*^9, {3.578715692932948*^9, 3.5787157128073826`*^9}, 
   3.5787157851447096`*^9, {3.5787159581490135`*^9, 3.57871599001987*^9}, {
   3.5787161271441107`*^9, 3.5787161562693615`*^9}, 3.578716187438216*^9, {
   3.5787174390284147`*^9, 3.578717456391245*^9}, 3.5787175063113327`*^9, {
   3.5787175762774553`*^9, 3.5787175778374586`*^9}, {3.578717639629167*^9, 
   3.578717644340375*^9}, 3.5787176797368374`*^9, {3.5787177641485853`*^9, 
   3.578717771792599*^9}, {3.5787179056720343`*^9, 3.5787179109604435`*^9}, {
   3.5787179462945056`*^9, 3.5787179462945056`*^9}, 3.578718000379801*^9, {
   3.5787181347896366`*^9, 3.578718154133671*^9}, {3.578718222680191*^9, 
   3.578718237827818*^9}, {3.578718362955638*^9, 3.578718418881736*^9}, {
   3.578718450612191*^9, 3.5787184740434327`*^9}},
 FormatType->"InputForm"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"t", " ", "=", 
  RowBox[{"FindMinimum", "[", 
   RowBox[{"func", ",", 
    RowBox[{"{", 
     RowBox[{"R1", ",", "0.26"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"R2", ",", "1.24"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"p", ",", "0.53"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"Method", "\[Rule]", "Automatic"}], ",", 
      RowBox[{"MaxIterations", "\[Rule]", "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, {3.5786568600064273`*^9, 
   3.578656911704918*^9}, {3.578657115706476*^9, 3.5786571617265573`*^9}, {
   3.5786674217311573`*^9, 3.5786675482473793`*^9}, 3.5786675866390467`*^9, {
   3.5786676279167194`*^9, 3.578667652237162*^9}, {3.5786676890688267`*^9, 
   3.5786676933588343`*^9}, {3.5787174121339674`*^9, 3.5787174314000015`*^9}, 
   3.578717511833742*^9, {3.5787175687894425`*^9, 3.5787175967290916`*^9}, 
   3.578717656804797*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", "40.56089226557409`"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"R1", "\[Rule]", "0.26537991269064937`"}], ",", 
     RowBox[{"R2", "\[Rule]", "1.2343054127201176`"}], ",", 
     RowBox[{"p", "\[Rule]", "0.5334356115481017`"}]}], "}"}]}], 
  "}"}]], "Output",
 CellChangeTimes->{3.5787205765839252`*^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[{
   RowBox[{"d1", "*", "sq", "*", "bohr"}], "     ", "/.", " ", 
   RowBox[{"t", "[", 
    RowBox[{"[", "2", "]"}], "]"}]}], "  "}], "\[IndentingNewLine]", 
 RowBox[{"d1", "/.", " ", 
  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}, {3.5786686933673906`*^9, 3.5786687128518248`*^9}}],

Cell[BoxData["40.56089386077506`"], "Output",
 CellChangeTimes->{3.5787205765995255`*^9}],

Cell[BoxData[
 RowBox[{"-", "121.15886817891526`"}]], "Output",
 CellChangeTimes->{3.5787205765995255`*^9}],

Cell[BoxData["26.425380893475342`"], "Output",
 CellChangeTimes->{3.5787205765995255`*^9}],

Cell[BoxData["13.611701159090778`"], "Output",
 CellChangeTimes->{3.5787205765995255`*^9}],

Cell[BoxData["1.999999960671454`"], "Output",
 CellChangeTimes->{3.5787205765995255`*^9}],

Cell[BoxData["1.0758808380570837`"], "Output",
 CellChangeTimes->{3.5787205765995255`*^9}],

Cell[BoxData["1.1738229202482673`"], "Output",
 CellChangeTimes->{3.5787205765995255`*^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.26537991269064937`], 
    {RGBColor[0., 0., 0.], SphereBox[{{0, 0, 0}}, 0.08]}}, 
   {Opacity[0.6], 
    SphereBox[{{0.8658437263256379, 0.8658437263256379, 0.8658437263256379}}, 
     1.2343054127201176`], 
    {RGBColor[0., 0., 0.], 
     SphereBox[{{1.1738229202482673`, 1.1738229202482673`, 
      1.1738229202482673`}}, 0.08]}}, 
   {Opacity[0.6], 
    SphereBox[{{-0.8658437263256379, -0.8658437263256379, 
     0.8658437263256379}}, 1.2343054127201176`], 
    {RGBColor[0., 0., 0.], 
     SphereBox[{{-1.1738229202482673`, -1.1738229202482673`, 
      1.1738229202482673`}}, 0.08]}}, 
   {Opacity[0.6], 
    SphereBox[{{
     0.8658437263256379, -0.8658437263256379, -0.8658437263256379}}, 
     1.2343054127201176`], 
    {RGBColor[0., 0., 0.], 
     SphereBox[{{
      1.1738229202482673`, -1.1738229202482673`, -1.1738229202482673`}}, 
      0.08]}}, 
   {Opacity[0.6], 
    SphereBox[{{-0.8658437263256379, 
     0.8658437263256379, -0.8658437263256379}}, 1.2343054127201176`], 
    {RGBColor[0., 0., 0.], 
     SphereBox[{{-1.1738229202482673`, 
      1.1738229202482673`, -1.1738229202482673`}}, 0.08]}}},
  AspectRatio->Automatic,
  Axes->True]], "Output",
 CellChangeTimes->{3.5787205766151257`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"wd", "=", 
   RowBox[{"w", ".", 
    RowBox[{"RotationMatrix", "[", 
     RowBox[{
      RowBox[{"\[Pi]", "/", "4"}], ",", 
      RowBox[{"{", 
       RowBox[{"0", ",", "1", ",", "0"}], "}"}]}], "]"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"nd", "=", 
   RowBox[{"n", ".", 
    RowBox[{"RotationMatrix", "[", 
     RowBox[{
      RowBox[{"\[Pi]", "/", "4"}], ",", 
      RowBox[{"{", 
       RowBox[{"0", ",", "1", ",", "0"}], "}"}]}], "]"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"plot2", "=", 
    RowBox[{
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"Graphics", "[", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"Circle", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{"wd", "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", "1"}], "]"}], "]"}], ",", 
              RowBox[{"wd", "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", "2"}], "]"}], "]"}]}], "}"}], ",", 
            RowBox[{"c", "[", 
             RowBox[{"[", "i", "]"}], "]"}]}], "]"}], ",", 
          RowBox[{"Disk", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{
              RowBox[{"nd", "[", 
               RowBox[{"[", 
                RowBox[{"i", ",", "1"}], "]"}], "]"}], ",", 
              RowBox[{"nd", "[", 
               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, {
   3.578719135172594*^9, 3.578719270019231*^9}, {3.578719304526491*^9, 
   3.578719384336231*^9}, {3.57871945197795*^9, 3.5787195085280495`*^9}, {
   3.5787195871833878`*^9, 3.578719605965821*^9}, {3.578719665869926*^9, 
   3.578719681298353*^9}, {3.5787197271936336`*^9, 3.578719882491906*^9}, {
   3.5787199375132027`*^9, 3.5787199417720103`*^9}, {3.578720091875474*^9, 
   3.578720131187543*^9}, {3.5787202075184774`*^9, 3.5787202448649426`*^9}, {
   3.578720313005862*^9, 3.5787203458127203`*^9}, {3.578720440161686*^9, 
   3.578720453359309*^9}, {3.5787204882409706`*^9, 3.5787204989425893`*^9}}],

Cell[BoxData[
 GraphicsBox[{{CircleBox[{0, 0}, 0.26537991269064937`], 
    DiskBox[{0, 0}, 0.08]}, {
    CircleBox[{0, 0.8658437263256379}, 1.2343054127201176`], 
    DiskBox[{0, 1.1738229202482673`}, 0.08]}, {
    CircleBox[{-1.2244879406653755`, -0.8658437263256379}, 
     1.2343054127201176`], 
    DiskBox[{-1.6600362936394912`, -1.1738229202482673`}, 0.08]}, {
    CircleBox[{1.2244879406653755`, -0.8658437263256379}, 
     1.2343054127201176`], 
    DiskBox[{1.6600362936394912`, -1.1738229202482673`}, 0.08]}, {
    CircleBox[{0, 0.8658437263256379}, 1.2343054127201176`], 
    DiskBox[{0, 1.1738229202482673`}, 0.08]}},
  AspectRatio->Automatic,
  Axes->True,
  Frame->True,
  GridLines->Automatic,
  PlotRange->{{-3, 3}, {-3, 3}}]], "Output",
 CellChangeTimes->{3.5787205766931257`*^9}]
}, Open  ]]
},
WindowToolbars->"EditBar",
WindowSize->{837, 831},
WindowMargins->{{Automatic, 495}, {Automatic, 2}},
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, 749, 13, 49, "Text"],
Cell[1309, 35, 972, 17, 52, "Input"],
Cell[2284, 54, 869, 15, 31, "Input"],
Cell[3156, 71, 1452, 32, 92, "Input"],
Cell[4611, 105, 1166, 28, 72, "Input"],
Cell[5780, 135, 3540, 57, 52, "Input"],
Cell[9323, 194, 12468, 288, 304, "Input"],
Cell[CellGroupData[{
Cell[21816, 486, 4179, 65, 52, "Input"],
Cell[25998, 553, 368, 10, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[26403, 568, 2178, 47, 152, "Input"],
Cell[28584, 617, 89, 1, 31, "Output"],
Cell[28676, 620, 107, 2, 31, "Output"],
Cell[28786, 624, 90, 1, 31, "Output"],
Cell[28879, 627, 90, 1, 31, "Output"],
Cell[28972, 630, 89, 1, 31, "Output"],
Cell[29064, 633, 90, 1, 31, "Output"],
Cell[29157, 636, 90, 1, 31, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[29284, 642, 2595, 55, 92, "Input"],
Cell[31882, 699, 1288, 32, 422, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[33207, 736, 4342, 98, 152, "Input"],
Cell[37552, 836, 797, 18, 411, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)