(* 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[      8092,        233]
NotebookOptionsPosition[      7263,        205]
NotebookOutlinePosition[      7960,        228]
CellTagsIndexPosition[      7917,        225]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["\<\
(* Ar atom nach Kimball, Aussen-sp-Kugeln 10.12.2011 *)
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.4, sig2 -> 0.4, sig3 -> 0.4};

T = 2.25*k1/R1^2+9.0*k2/R2^2+9.0*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+4.0*3.0*sig3/R3+16/(R1+R2)+24/((R1+R2)*Sqrt[\
8/3])+16/a+48/(R2+R3)+16/(R1+R2+a)+24*Sqrt[3/8]/a /. c;
Vne=-3.0*z/R1-8*z/(R1+R2)-8*z/a;

func = T + Vee + Vne;

t = Table[FindMinimum[func, {R1,0.07}, {R2,0.2}, {R3,0.5}, {MaxIterations -> \
500}], {z,18,26}]

Table[N[Vne /. c /. t[[{z-17},2]],10], {z,18,26}]
Table[N[Vee /. c /. t[[{z-17},2]],10], {z,18,26}]
Table[N[-(Vee+Vne)/T /. c /. t[[{z-17},2]],10], {z,18,26}]
Table[N[Sqrt[a/(R1+R2+R3)] /. c /. t[[{z-17},2]],10], {z,18,26}]
\
\>", "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.582948845164712*^9, 3.5829488512019224`*^9}, {
  3.5829488914031925`*^9, 3.5829489011532097`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "515.6667109150261`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"R1", "\[Rule]", "0.07711800971360289`"}], ",", 
       RowBox[{"R2", "\[Rule]", "0.27033658921390286`"}], ",", 
       RowBox[{"R3", "\[Rule]", "0.7964803884161074`"}]}], "}"}]}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "587.1556187901322`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"R1", "\[Rule]", "0.07272328087812782`"}], ",", 
       RowBox[{"R2", "\[Rule]", "0.25102394606285783`"}], ",", 
       RowBox[{"R3", "\[Rule]", "0.6902359278529366`"}]}], "}"}]}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "663.6153929440742`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"R1", "\[Rule]", "0.06879809865620665`"}], ",", 
       RowBox[{"R2", "\[Rule]", "0.23421567370790466`"}], ",", 
       RowBox[{"R3", "\[Rule]", "0.6100291387719222`"}]}], "}"}]}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "745.0356278756125`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"R1", "\[Rule]", "0.06527178047448152`"}], ",", 
       RowBox[{"R2", "\[Rule]", "0.21946867000491133`"}], ",", 
       RowBox[{"R3", "\[Rule]", "0.5471302757652722`"}]}], "}"}]}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "831.4092141965693`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"R1", "\[Rule]", "0.06208698871503191`"}], ",", 
       RowBox[{"R2", "\[Rule]", "0.20643483118425401`"}], ",", 
       RowBox[{"R3", "\[Rule]", "0.4963670614995652`"}]}], "}"}]}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "922.7311296316776`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"R1", "\[Rule]", "0.05919676764998288`"}], ",", 
       RowBox[{"R2", "\[Rule]", "0.1948380401644524`"}], ",", 
       RowBox[{"R3", "\[Rule]", "0.45446751488564696`"}]}], "}"}]}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "1018.9977260360312`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"R1", "\[Rule]", "0.05656231886324758`"}], ",", 
       RowBox[{"R2", "\[Rule]", "0.1844571197830332`"}], ",", 
       RowBox[{"R3", "\[Rule]", "0.4192539086755075`"}]}], "}"}]}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "1120.2062900564856`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"R1", "\[Rule]", "0.05415132168808806`"}], ",", 
       RowBox[{"R2", "\[Rule]", "0.17511319118268998`"}], ",", 
       RowBox[{"R3", "\[Rule]", "0.3892168684421618`"}]}], "}"}]}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "1226.3547622126478`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"R1", "\[Rule]", "0.0519366459720168`"}], ",", 
       RowBox[{"R2", "\[Rule]", "0.16666025131427997`"}], ",", 
       RowBox[{"R3", "\[Rule]", "0.36327508081166954`"}]}], "}"}]}], "}"}]}], 
  "}"}]], "Output",
 CellChangeTimes->{3.5829486896168385`*^9, 3.582948860483938*^9, 
  3.5829489064884195`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"-", "1241.978402283204`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1405.5522761780305`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1578.8282482008017`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1761.8557177044797`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1954.6660142331925`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "2157.279815306124`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "2369.711292173279`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "2591.9704329849733`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "2824.064465695024`"}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.5829486896168385`*^9, 3.582948860483938*^9, 
  3.5829489064884195`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", "210.64497999445484`", "}"}], ",", 
   RowBox[{"{", "231.24103805980877`", "}"}], ",", 
   RowBox[{"{", "251.59746297120603`", "}"}], ",", 
   RowBox[{"{", "271.78446925083864`", "}"}], ",", 
   RowBox[{"{", "291.84758752936654`", "}"}], ",", 
   RowBox[{"{", "311.817556917185`", "}"}], ",", 
   RowBox[{"{", "331.7158400970783`", "}"}], ",", 
   RowBox[{"{", "351.55785267751486`", "}"}], ",", 
   RowBox[{"{", "371.3549361811258`", "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.5829486896168385`*^9, 3.582948860483938*^9, 
  3.5829489064884195`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", "1.9999999991104775`", "}"}], ",", 
   RowBox[{"{", "1.9999999990837902`", "}"}], ",", 
   RowBox[{"{", "2.000000000992372`", "}"}], ",", 
   RowBox[{"{", "2.000000009794947`", "}"}], ",", 
   RowBox[{"{", "2.0000000020318667`", "}"}], ",", 
   RowBox[{"{", "2.00000000094764`", "}"}], ",", 
   RowBox[{"{", "1.9999999999959392`", "}"}], ",", 
   RowBox[{"{", "1.9999999998263835`", "}"}], ",", 
   RowBox[{"{", "1.9999999958506274`", "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.5829486896168385`*^9, 3.582948860483938*^9, 
  3.5829489064884195`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", "0.9943720994995239`", "}"}], ",", 
   RowBox[{"{", "0.9922418837590123`", "}"}], ",", 
   RowBox[{"{", "0.9903155698099795`", "}"}], ",", 
   RowBox[{"{", "0.9885744491405891`", "}"}], ",", 
   RowBox[{"{", "0.986997760130254`", "}"}], ",", 
   RowBox[{"{", "0.9855658263310554`", "}"}], ",", 
   RowBox[{"{", "0.9842610536795957`", "}"}], ",", 
   RowBox[{"{", "0.9830680937729245`", "}"}], ",", 
   RowBox[{"{", "0.9819736984191736`", "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.5829486896168385`*^9, 3.582948860483938*^9, 
  3.5829489064884195`*^9}]
}, Open  ]]
},
PrintingStyleEnvironment->"Printout",
WindowSize->{857, 776},
WindowMargins->{{Automatic, 281}, {Automatic, 48}},
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, 958, 24, 371, "Input"],
Cell[1540, 48, 3072, 85, 192, "Output"],
Cell[4615, 135, 790, 22, 52, "Output"],
Cell[5408, 159, 610, 13, 52, "Output"],
Cell[6021, 174, 610, 13, 31, "Output"],
Cell[6634, 189, 613, 13, 52, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)