(* 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[      6909,        178]
NotebookOptionsPosition[      5534,        130]
NotebookOutlinePosition[      6525,        166]
CellTagsIndexPosition[      6482,        163]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell[TextData[{
 StyleBox["(* SF6 als S+6  + 6 F- in Oh 14.02.12 *)\n\
Clear[k1,k2,k3,k4,k5,k6,sig1,sig2,sig3,sig4,sig5,c,z,R1,R2,R3,R4,R5,S2,S3,S4];\
\nc = {k1 -> 0.88, k2 -> 2.435, k5 -> 0.88, k6 -> 2.43,sig1 -> 0.3, \nsig2 -> \
0.30, sig3 -> 0.30, sig4 -> 0.30, sig5 -> 0.302}; z=16.; z1=9.0;\nad = \
Sqrt[3./8.]; sq2=Sqrt[2];\n\n(* S+6 He+Ne shell *)\nT = \
2.25*k1/R1^2+9.*k2/R2^2 /. c;\n\
Vee=3.0*sig1/R1+12.*sig2/R2+16/(R1+R2)+24*ad/(R1+R2) /. c;\n\
Vne=-3.0*z/R1-8.7*z/(R1+R2);\nS2 = R2*4^(1/3);\n\n(* F- He+Ne shell *)\nTf = \
2.25*k5/R5^2+9.*k6/R6^2 /. c;\n\
Veef=3.0*sig1/R5+12*sig2/R6+16/(R5+R6)+24*ad/(R5+R6) /. c;\n\
Vnef=-3.0*z1/R5-8.3*z1/(R5+R6) /. c;\nS5 = R6*4^(1/3);\n\nT=T+6*Tf;\n\
Vee=Vee+6*Veef;\nVne=Vne+6*Vnef;\n\n(* S+6/6 F- *)\nVnn = (6*z*z1 + \
12*z1^2/sq2+1.5*z1^2)/(S2+S5);\nVee = Vee+(6*10*10 + 12*10*10/sq2 + \
1.5*10*10)/(S2+S5);\nVne = Vne \
-(6*z*10+6*z1*10+24*z1*10/sq2+3*z1*10)/(S2+S5);\n\nfunc = T + Vnn + Vne + \
Vee;\n\nt = FindMinimum[func, {R1,0.0809742}, {R2,0.546933},\n    \
{R5,0.1491164},{R6,1.3151171},{Method -> Automatic}, {MaxIterations -> 500}]\n\
    \nN[Vne /. c /. t[[2]],10]\nN[Vnn /. c /. t[[2]],10]\nN[Vee /. c /. \
t[[2]],10]\nN[-(Vee+Vne+Vnn)/T /. c /. t[[2]],10]\nN[(S2+S5) /. c /. \
t[[2]],10]\nN[0.529177*(S2+S5) /. c /. t[[2]],10]\nN[",
  FontSize->16],
 "ratio=-(Vne+2*Vnn)/(Vee-Vnn)",
 StyleBox[" /. c /. t[[2]],8]\nd=S2+S5 /.c /.t[[2]]\n",
  FontSize->16],
 "(* xy plane *)\nplot2=Graphics[{\n      Circle[{0,0},R1],Circle[{0,0},S2],\n\
      Circle[{d,0},R5],\n      Circle[{d,0},S5],\n      Circle[{-d,0},S5],\n  \
    Circle[{-d,0},R5],\n      Circle[{0,d},R5],\n      \
Circle[{0,d},S5],Circle[{0,-d},R5],\n      Circle[{0,-d},S5]} ] /. t[[2]] /. \
c\n\nShow[plot2,{AspectRatio -> Automatic,Axes -> True,\n    GridLines -> \
Automatic, PlotRange \\[Rule] {{-6,6},{-6,6}}, \n    Frame -> True}]\n"
}], "Input",
 PageWidth->Infinity,
 InitializationCell->True,
 CellChangeTimes->{{3.580216356436178*^9, 3.5802164023158584`*^9}},
 ShowSpecialCharacters->False,
 FontSize->16],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", "985.377647721967`"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"R1", "\[Rule]", "0.08097428759759646`"}], ",", 
     RowBox[{"R2", "\[Rule]", "0.5469397006476997`"}], ",", 
     RowBox[{"R5", "\[Rule]", "0.14911677762821401`"}], ",", 
     RowBox[{"R6", "\[Rule]", "1.315117047285238`"}]}], "}"}]}], 
  "}"}]], "Output",
 CellChangeTimes->{{3.580216424265097*^9, 3.5802164243118973`*^9}}],

Cell[BoxData[
 RowBox[{"-", "3322.5037480694336`"}]], "Output",
 CellChangeTimes->{{3.580216424265097*^9, 3.5802164243118973`*^9}}],

Cell[BoxData["565.934886528956`"], "Output",
 CellChangeTimes->{{3.580216424265097*^9, 3.5802164243118973`*^9}}],

Cell[BoxData["785.8135591173312`"], "Output",
 CellChangeTimes->{{3.580216424265097*^9, 3.5802164243118973`*^9}}],

Cell[BoxData["1.9999999929172212`"], "Output",
 CellChangeTimes->{{3.580216424265097*^9, 3.5802164243118973`*^9}}],

Cell[BoxData["2.9558308404932294`"], "Output",
 CellChangeTimes->{{3.580216424265097*^9, 3.5802164243118973`*^9}}],

Cell[BoxData["1.5641576966796857`"], "Output",
 CellChangeTimes->{{3.580216424265097*^9, 3.5802164243118973`*^9}}],

Cell[BoxData["9.962921593184744`"], "Output",
 CellChangeTimes->{{3.580216424265097*^9, 3.5802164243118973`*^9}}],

Cell[BoxData["2.9558308404932294`"], "Output",
 CellChangeTimes->{{3.580216424265097*^9, 3.5802164243118973`*^9}}],

Cell[BoxData[
 GraphicsBox[{CircleBox[{0, 0}, 0.08097428759759646], 
   CircleBox[{0, 0}, 0.8682126561713306], 
   CircleBox[{2.9558308404932294`, 0}, 0.14911677762821401`], 
   CircleBox[{2.9558308404932294`, 0}, 2.087618184321899], 
   CircleBox[{-2.9558308404932294`, 0}, 2.087618184321899], 
   CircleBox[{-2.9558308404932294`, 0}, 0.14911677762821401`], 
   CircleBox[{0, 2.9558308404932294`}, 0.14911677762821401`], 
   CircleBox[{0, 2.9558308404932294`}, 2.087618184321899], 
   CircleBox[{0, -2.9558308404932294`}, 0.14911677762821401`], 
   CircleBox[{0, -2.9558308404932294`}, 2.087618184321899]}]], "Output",
 CellChangeTimes->{{3.580216424265097*^9, 3.5802164243118973`*^9}}],

Cell[BoxData[
 GraphicsBox[{CircleBox[{0, 0}, 0.08097428759759646], 
   CircleBox[{0, 0}, 0.8682126561713306], 
   CircleBox[{2.9558308404932294`, 0}, 0.14911677762821401`], 
   CircleBox[{2.9558308404932294`, 0}, 2.087618184321899], 
   CircleBox[{-2.9558308404932294`, 0}, 2.087618184321899], 
   CircleBox[{-2.9558308404932294`, 0}, 0.14911677762821401`], 
   CircleBox[{0, 2.9558308404932294`}, 0.14911677762821401`], 
   CircleBox[{0, 2.9558308404932294`}, 2.087618184321899], 
   CircleBox[{0, -2.9558308404932294`}, 0.14911677762821401`], 
   CircleBox[{0, -2.9558308404932294`}, 2.087618184321899]},
  AspectRatio->Automatic,
  Axes->True,
  Frame->True,
  GridLines->Automatic,
  PlotRange->{{-6, 6}, {-6, 6}}]], "Output",
 CellChangeTimes->{{3.580216424265097*^9, 3.5802164243118973`*^9}}]
}, Open  ]]
},
AutoGeneratedPackage->Automatic,
WindowSize->{988, 740},
WindowMargins->{{Automatic, 413}, {Automatic, 37}},
PrintingCopies->1,
PrintingPageRange->{Automatic, Automatic},
PageHeaders->{{
   Cell[
    TextData[{
      CounterBox["Page"]}], "PageNumber"], Inherited, 
   Cell[
    TextData[{
      ValueBox["FileName"]}], "Header"]}, {
   Cell[
    TextData[{
      ValueBox["FileName"]}], "Header"], Inherited, 
   Cell[
    TextData[{
      CounterBox["Page"]}], "PageNumber"]}},
PrintingOptions->{"PrintCellBrackets"->False,
"PrintMultipleHorizontalPages"->False,
"PrintRegistrationMarks"->True,
"PrintingMargins"->{{36, 36}, {57.5625, 36}}},
PrivateNotebookOptions->{"VersionedStylesheet"->{"Default.nb"[8.] -> False}},
FrontEndVersion->"9.0 for Microsoft Windows (64-bit) (January 25, 2013)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[579, 22, 2042, 36, 1059, "Input",
 InitializationCell->True],
Cell[2624, 60, 453, 11, 31, "Output"],
Cell[3080, 73, 131, 2, 31, "Output"],
Cell[3214, 77, 112, 1, 31, "Output"],
Cell[3329, 80, 113, 1, 31, "Output"],
Cell[3445, 83, 114, 1, 31, "Output"],
Cell[3562, 86, 114, 1, 31, "Output"],
Cell[3679, 89, 114, 1, 31, "Output"],
Cell[3796, 92, 113, 1, 31, "Output"],
Cell[3912, 95, 114, 1, 31, "Output"],
Cell[4029, 98, 687, 11, 374, "Output"],
Cell[4719, 111, 799, 16, 411, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)