(* 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[     11113,        284]
NotebookOptionsPosition[     10242,        255]
NotebookOutlinePosition[     10939,        278]
CellTagsIndexPosition[     10896,        275]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["\<\
(* Cl atom nach Kimball, Aussen-sp-Kugeln 18.12.2011/24.03.2013 *)
Clear[k1,k2,k3,sig1,sig2,sig3,c,z,R1,R2,R3,a];
c = {k1 -> 1.0, k2 -> 1.0, k3->1.0, sig1 -> 0.32, sig2 -> 0.32, sig3 -> 0.32};

T = 2.25*k1/R1^2+9.0*k2/R2^2+7.875*k3/R3^2 /. c;
a = (R1+R2)/3+Sqrt[(R2+R3)^2-8*(R1+R2)^2/9];
Vee = 3.0*sig1/R1+4.0*3.0*sig2/R2+3.0*3.0*sig3/R3+16/(R1+R2)+24/((R1+R2)*Sqrt[\
8/3])+
    14/a+42/(R2+R3)+14/(R1+R2+a)+18*Sqrt[3/8]/a /. c;
Vne=-3.0*z/R1-8*z/(R1+R2)-7*z/a;

func = T + Vee + Vne;

t = Table[FindMinimum[func, {R1,0.06}, {R2,0.2}, {R3,0.43}, {Method->\"Newton\
\",MaxIterations -> 500}], {z,17,26}]

Table[N[Vne /. c /. t[[{z-16},2]],10], {z,17,26}]
Table[N[Vee /. c /. t[[{z-16},2]],10], {z,17,26}]
Table[N[-(Vee+Vne)/T /. c /. t[[{z-16},2]],10], {z,17,26}]
Table[N[a /. c /. t[[{z-16},2]],10], {z,17,26}]
\
\>", "Input",
 CellChangeTimes->{{3.5731434931493855`*^9, 3.573143509935015*^9}, {
   3.5731435457994785`*^9, 3.573143660210079*^9}, {3.573143706635761*^9, 
   3.573143733327408*^9}, {3.5731438491108108`*^9, 3.5731439199505353`*^9}, {
   3.573143965284215*^9, 3.573143979932641*^9}, 3.57314404192715*^9, {
   3.573144123343693*^9, 3.573144163950564*^9}, {3.573144202732232*^9, 
   3.573144203434233*^9}, 3.5731442369742923`*^9, {3.573144317688834*^9, 
   3.5731443254420476`*^9}}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"FindMinimum", "::", "lstol"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"The line search decreased the step size to within the \
tolerance specified by AccuracyGoal and PrecisionGoal but was unable to find \
a sufficient decrease in the function. You may need more than \
\[NoBreak]\\!\\(MachinePrecision\\)\[NoBreak] digits of working precision to \
meet these tolerances. \\!\\(\\*ButtonBox[\\\"\[RightSkeleton]\\\", \
ButtonStyle->\\\"Link\\\", ButtonFrame->None, \
ButtonData:>\\\"paclet:ref/message/FindMinimum/lstol\\\", ButtonNote -> \
\\\"FindMinimum::lstol\\\"]\\)\"\>"}]], "Message", "MSG",
 CellChangeTimes->{
  3.57314404771476*^9, {3.5731441289129024`*^9, 3.5731441705961757`*^9}, 
   3.5731442086914425`*^9, 3.5731442409054995`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "456.2099110511343`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"R1", "\[Rule]", "0.0816329703655353`"}], ",", 
       RowBox[{"R2", "\[Rule]", "0.288942546139099`"}], ",", 
       RowBox[{"R3", "\[Rule]", "0.8550326474433549`"}]}], "}"}]}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "522.7403636427763`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"R1", "\[Rule]", "0.0767380212319812`"}], ",", 
       RowBox[{"R2", "\[Rule]", "0.26727214004346267`"}], ",", 
       RowBox[{"R3", "\[Rule]", "0.7328113705855436`"}]}], "}"}]}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "594.1025400618049`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"R1", "\[Rule]", "0.07239167875866229`"}], ",", 
       RowBox[{"R2", "\[Rule]", "0.2485465010103451`"}], ",", 
       RowBox[{"R3", "\[Rule]", "0.6426209370636539`"}]}], "}"}]}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "670.2852787838276`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"R1", "\[Rule]", "0.06850753294268336`"}], ",", 
       RowBox[{"R2", "\[Rule]", "0.23222053259227926`"}], ",", 
       RowBox[{"R3", "\[Rule]", "0.5730318593023928`"}]}], "}"}]}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "751.2811461661835`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"R1", "\[Rule]", "0.06501621069968305`"}], ",", 
       RowBox[{"R2", "\[Rule]", "0.21787123926268312`"}], ",", 
       RowBox[{"R3", "\[Rule]", "0.5175457460785067`"}]}], "}"}]}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "837.0849953102896`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"R1", "\[Rule]", "0.06186142781357465`"}], ",", 
       RowBox[{"R2", "\[Rule]", "0.205166851981688`"}], ",", 
       RowBox[{"R3", "\[Rule]", "0.4721752539684697`"}]}], "}"}]}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "927.6931472342006`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"R1", "\[Rule]", "0.05899707424938501`"}], ",", 
       RowBox[{"R2", "\[Rule]", "0.1938443063864721`"}], ",", 
       RowBox[{"R3", "\[Rule]", "0.43432794344239`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "1023.1029012249353`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"R1", "\[Rule]", "0.05638503574038974`"}], ",", 
       RowBox[{"R2", "\[Rule]", "0.18369279485453743`"}], ",", 
       RowBox[{"R3", "\[Rule]", "0.40223988419368145`"}]}], "}"}]}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "1123.3122293593829`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"R1", "\[Rule]", "0.053993546178689796`"}], ",", 
       RowBox[{"R2", "\[Rule]", "0.17454167184196454`"}], ",", 
       RowBox[{"R3", "\[Rule]", "0.3746660882157774`"}]}], "}"}]}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "1228.3195789889048`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"R1", "\[Rule]", "0.05179592807034305`"}], ",", 
       RowBox[{"R2", "\[Rule]", "0.1662514908256886`"}], ",", 
       RowBox[{"R3", "\[Rule]", "0.3507009538725807`"}]}], "}"}]}], "}"}]}], 
  "}"}]], "Output",
 CellChangeTimes->{{3.573143529201049*^9, 3.5731436644532866`*^9}, {
   3.5731437157461767`*^9, 3.5731437378046155`*^9}, {3.5731438579404263`*^9, 
   3.5731439258629456`*^9}, {3.573143970806625*^9, 3.5731439856734505`*^9}, 
   3.5731440477771597`*^9, {3.5731441289597025`*^9, 3.573144170642976*^9}, 
   3.5731442087382426`*^9, 3.573144240952299*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"-", "1089.8609361823233`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1241.0746205798764`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1401.704946344715`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1571.8063328408016`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1751.412039946955`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1940.543432007893`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "2139.2147988217216`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "2347.436020487187`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "2565.214138082292`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "2792.554278453123`"}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{{3.573143529201049*^9, 3.5731436644532866`*^9}, {
   3.5731437157461767`*^9, 3.5731437378046155`*^9}, {3.5731438579404263`*^9, 
   3.5731439258629456`*^9}, {3.573143970806625*^9, 3.5731439856734505`*^9}, 
   3.5731440477771597`*^9, {3.5731441289597025`*^9, 3.573144170642976*^9}, 
   3.5731442087382426`*^9, 3.573144240952299*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", "177.44111407522837`", "}"}], ",", 
   RowBox[{"{", "195.5938927152515`", "}"}], ",", 
   RowBox[{"{", "213.49986622099584`", "}"}], ",", 
   RowBox[{"{", "231.23577518517115`", "}"}], ",", 
   RowBox[{"{", "248.84974761458295`", "}"}], ",", 
   RowBox[{"{", "266.373441385487`", "}"}], ",", 
   RowBox[{"{", "283.8285031912724`", "}"}], ",", 
   RowBox[{"{", "301.23021785313426`", "}"}], ",", 
   RowBox[{"{", "318.58967788044896`", "}"}], ",", 
   RowBox[{"{", "335.91511991457253`", "}"}]}], "}"}]], "Output",
 CellChangeTimes->{{3.573143529201049*^9, 3.5731436644532866`*^9}, {
   3.5731437157461767`*^9, 3.5731437378046155`*^9}, {3.5731438579404263`*^9, 
   3.5731439258629456`*^9}, {3.573143970806625*^9, 3.5731439856734505`*^9}, 
   3.5731440477771597`*^9, {3.5731441289597025`*^9, 3.573144170642976*^9}, 
   3.5731442087382426`*^9, 3.573144240952299*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", "1.9999999999894211`", "}"}], ",", 
   RowBox[{"{", "1.9999999988922368`", "}"}], ",", 
   RowBox[{"{", "1.9999999999998164`", "}"}], ",", 
   RowBox[{"{", "1.9999999998687499`", "}"}], ",", 
   RowBox[{"{", "1.9999999999999933`", "}"}], ",", 
   RowBox[{"{", "1.9999999999978175`", "}"}], ",", 
   RowBox[{"{", "1.999999998747379`", "}"}], ",", 
   RowBox[{"{", "1.9999999998199764`", "}"}], ",", 
   RowBox[{"{", "1.9999999986797277`", "}"}], ",", 
   RowBox[{"{", "1.9999999995434896`", "}"}]}], "}"}]], "Output",
 CellChangeTimes->{{3.573143529201049*^9, 3.5731436644532866`*^9}, {
   3.5731437157461767`*^9, 3.5731437378046155`*^9}, {3.5731438579404263`*^9, 
   3.5731439258629456`*^9}, {3.573143970806625*^9, 3.5731439856734505`*^9}, 
   3.5731440477771597`*^9, {3.5731441289597025`*^9, 3.573144170642976*^9}, 
   3.5731442087382426`*^9, 3.5731442409678993`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", "1.2128421592103316`", "}"}], ",", 
   RowBox[{"{", "1.0607003166995004`", "}"}], ",", 
   RowBox[{"{", "0.9452053043515635`", "}"}], ",", 
   RowBox[{"{", "0.8539287251686525`", "}"}], ",", 
   RowBox[{"{", "0.7796457521035838`", "}"}], ",", 
   RowBox[{"{", "0.7178264940401031`", "}"}], ",", 
   RowBox[{"{", "0.6654644307985011`", "}"}], ",", 
   RowBox[{"{", "0.6204738052829738`", "}"}], ",", 
   RowBox[{"{", "0.5813555703498943`", "}"}], ",", 
   RowBox[{"{", "0.5470008207747593`", "}"}]}], "}"}]], "Output",
 CellChangeTimes->{{3.573143529201049*^9, 3.5731436644532866`*^9}, {
   3.5731437157461767`*^9, 3.5731437378046155`*^9}, {3.5731438579404263`*^9, 
   3.5731439258629456`*^9}, {3.573143970806625*^9, 3.5731439856734505`*^9}, 
   3.5731440477771597`*^9, {3.5731441289597025`*^9, 3.573144170642976*^9}, 
   3.5731442087382426`*^9, 3.5731442409678993`*^9}]
}, Open  ]]
},
PrintingStyleEnvironment->"Printout",
WindowSize->{973, 776},
WindowMargins->{{Automatic, 283}, {Automatic, 14}},
PrintingCopies->1,
PrintingPageRange->{Automatic, Automatic},
PrintingOptions->{"PrintCellBrackets"->False,
"PrintMultipleHorizontalPages"->False,
"PrintRegistrationMarks"->False,
"PrintingMargins"->{{34, 14.125}, {56.6875, 56.6875}}},
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, 1303, 29, 337, "Input"],
Cell[1885, 53, 799, 15, 56, "Message"],
Cell[2687, 70, 3642, 96, 212, "Output"],
Cell[6332, 168, 1117, 27, 52, "Output"],
Cell[7452, 197, 919, 17, 52, "Output"],
Cell[8374, 216, 924, 17, 31, "Output"],
Cell[9301, 235, 925, 17, 52, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)