(* 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[     12442,        325]
NotebookOptionsPosition[     11529,        295]
NotebookOutlinePosition[     12225,        318]
CellTagsIndexPosition[     12182,        315]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

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

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

func = T + Vee + Vne;

t = Table[FindMinimum[func, {R1,0.04}, {R2,0.25}, {R3,0.6}, {MaxIterations -> \
500, Method -> \"Newton\"}], {z,15,26}]

Table[N[Vne /. c /. t[[{z-14},2]],10], {z,15,26}]
Table[N[Vee /. c /. t[[{z-14},2]],10], {z,15,26}]
Table[N[-(Vee+Vne)/T /. c /. t[[{z-14},2]],10], {z,15,26}]
Table[N[a /. c /. t[[{z-14},2]],10], {z,15,26}]
\
\>", "Input",
 PageWidth->WindowWidth,
 CellChangeTimes->{{3.5746123826363306`*^9, 3.5746125087781525`*^9}, {
  3.574612573175065*^9, 3.5746126286487627`*^9}, {3.5746126911736727`*^9, 
  3.5746127207357244`*^9}, {3.5746127858970385`*^9, 3.5746128390775323`*^9}, {
  3.57461297443897*^9, 3.5746133186379747`*^9}, {3.5746133562652407`*^9, 
  3.574613363316453*^9}, {3.5746134449981966`*^9, 3.5746134456533976`*^9}}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"FindMinimum", "::", "nrnum"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"The function value \
\[NoBreak]\\!\\(\\(\\(-88.87980613256218`\\)\\) + \
\\(\\(507.5077041500493`\\\\ \[ImaginaryI]\\)\\)\\)\[NoBreak] is not a real \
number at \[NoBreak]\\!\\({R1, R2, R3}\\)\[NoBreak] = \
\[NoBreak]\\!\\({0.04031896713285391`, 0.16912981486857842`, \
\\(\\(-0.09309045543977634`\\)\\)}\\)\[NoBreak]. \\!\\(\\*ButtonBox[\\\"\
\[RightSkeleton]\\\", ButtonStyle->\\\"Link\\\", ButtonFrame->None, \
ButtonData:>\\\"paclet:ref/FindMinimum\\\", ButtonNote -> \
\\\"FindMinimum::nrnum\\\"]\\)\"\>"}]], "Message", "MSG",
 CellChangeTimes->{{3.574613081361558*^9, 3.5746133238015833`*^9}, 
   3.5746134536094117`*^9}],

Cell[BoxData[
 RowBox[{
  StyleBox[
   RowBox[{"FindMinimum", "::", "nrnum"}], "MessageName"], 
  RowBox[{
  ":", " "}], "\<\"The function value \
\[NoBreak]\\!\\(\\(\\(-22460.994600724775`\\)\\) + \\(\\(655.0905616155553`\\\
\\ \[ImaginaryI]\\)\\)\\)\[NoBreak] is not a real number at \
\[NoBreak]\\!\\({R1, R2, R3}\\)\[NoBreak] = \
\[NoBreak]\\!\\({0.04562684890741754`, 0.11375419332683194`, \
\\(\\(-0.11513947986484785`\\)\\)}\\)\[NoBreak]. \\!\\(\\*ButtonBox[\\\"\
\[RightSkeleton]\\\", ButtonStyle->\\\"Link\\\", ButtonFrame->None, \
ButtonData:>\\\"paclet:ref/FindMinimum\\\", ButtonNote -> \
\\\"FindMinimum::nrnum\\\"]\\)\"\>"}]], "Message", "MSG",
 CellChangeTimes->{{3.574613081361558*^9, 3.5746133238015833`*^9}, 
   3.5746134536874123`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "341.932491956393`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"R1", "\[Rule]", "0.09254471858595847`"}], ",", 
       RowBox[{"R2", "\[Rule]", "0.34061626425294206`"}], ",", 
       RowBox[{"R3", "\[Rule]", "1.0704822686217355`"}]}], "}"}]}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "398.4698878791991`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"R1", "\[Rule]", "0.08628398088003397`"}], ",", 
       RowBox[{"R2", "\[Rule]", "0.31145330204809685`"}], ",", 
       RowBox[{"R3", "\[Rule]", "0.8796843089702161`"}]}], "}"}]}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "459.6087484299541`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"R1", "\[Rule]", "0.08080853039640279`"}], ",", 
       RowBox[{"R2", "\[Rule]", "0.2867851270813821`"}], ",", 
       RowBox[{"R3", "\[Rule]", "0.7500679863879944`"}]}], "}"}]}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "525.3352786928641`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"R1", "\[Rule]", "0.07598087392250386`"}], ",", 
       RowBox[{"R2", "\[Rule]", "0.2656729099596509`"}], ",", 
       RowBox[{"R3", "\[Rule]", "0.6555219771085457`"}]}], "}"}]}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "595.6409583760214`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"R1", "\[Rule]", "0.07169348189650378`"}], ",", 
       RowBox[{"R2", "\[Rule]", "0.2474142885608197`"}], ",", 
       RowBox[{"R3", "\[Rule]", "0.5831350763718222`"}]}], "}"}]}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "670.520219482222`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"R1", "\[Rule]", "0.0678611496498792`"}], ",", 
       RowBox[{"R2", "\[Rule]", "0.2314762040542397`"}], ",", 
       RowBox[{"R3", "\[Rule]", "0.5257341435805841`"}]}], "}"}]}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "749.969258993716`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"R1", "\[Rule]", "0.0644155491888764`"}], ",", 
       RowBox[{"R2", "\[Rule]", "0.21744821615704818`"}], ",", 
       RowBox[{"R3", "\[Rule]", "0.4789885925649421`"}]}], "}"}]}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "833.9853853537335`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"R1", "\[Rule]", "0.06130128469791216`"}], ",", 
       RowBox[{"R2", "\[Rule]", "0.20500991056187373`"}], ",", 
       RowBox[{"R3", "\[Rule]", "0.4401160047839307`"}]}], "}"}]}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "922.5666370337392`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"R1", "\[Rule]", "0.05847299422891402`"}], ",", 
       RowBox[{"R2", "\[Rule]", "0.19390790835137317`"}], ",", 
       RowBox[{"R3", "\[Rule]", "0.40724009721171894`"}]}], "}"}]}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "1015.711549005993`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"R1", "\[Rule]", "0.055893191030741425`"}], ",", 
       RowBox[{"R2", "\[Rule]", "0.1839394398351013`"}], ",", 
       RowBox[{"R3", "\[Rule]", "0.3790457988838496`"}]}], "}"}]}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "946.0723163464886`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"R1", "\[Rule]", "0.04`"}], ",", 
       RowBox[{"R2", "\[Rule]", "0.25`"}], ",", 
       RowBox[{"R3", "\[Rule]", "0.6`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "1113.577545400156`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"R1", "\[Rule]", "0.03982729683809949`"}], ",", 
       RowBox[{"R2", "\[Rule]", "0.19184155424793503`"}], ",", 
       RowBox[{"R3", "\[Rule]", "0.47189263229782463`"}]}], "}"}]}], "}"}]}], 
  "}"}]], "Output",
 CellChangeTimes->{{3.5746123888919415`*^9, 3.5746125149401627`*^9}, {
   3.574612578525875*^9, 3.574612636651577*^9}, 3.5746127303141413`*^9, 
   3.574612805865074*^9, 3.574612846347145*^9, {3.5746129812561817`*^9, 
   3.574613323941984*^9}, 3.5746134536874123`*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"-", "813.448633356039`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "941.4567560139374`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1078.385642175064`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1224.3100802076776`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1379.2712414526688`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1543.2935095283942`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1716.3921818718757`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "1898.5773072290817`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "2089.8557408897314`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "2290.2323103381987`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "2703.3141355167472`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"-", "3040.9509584275293`"}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{{3.5746123888919415`*^9, 3.5746125149401627`*^9}, {
   3.574612578525875*^9, 3.574612636651577*^9}, 3.5746127303141413`*^9, 
   3.574612805865074*^9, 3.574612846347145*^9, {3.5746129812561817`*^9, 
   3.574613323941984*^9}, 3.574613453703012*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", "129.58364944167562`", "}"}], ",", 
   RowBox[{"{", "144.5169802555331`", "}"}], ",", 
   RowBox[{"{", "159.16814368075003`", "}"}], ",", 
   RowBox[{"{", "173.6395227787889`", "}"}], ",", 
   RowBox[{"{", "187.9893246999328`", "}"}], ",", 
   RowBox[{"{", "202.25307056367706`", "}"}], ",", 
   RowBox[{"{", "216.4536638777743`", "}"}], ",", 
   RowBox[{"{", "230.60653651858973`", "}"}], ",", 
   RowBox[{"{", "244.722466821536`", "}"}], ",", 
   RowBox[{"{", "258.8092116927815`", "}"}], ",", 
   RowBox[{"{", "206.20476917025863`", "}"}], ",", 
   RowBox[{"{", "255.07195542655276`", "}"}]}], "}"}]], "Output",
 CellChangeTimes->{{3.5746123888919415`*^9, 3.5746125149401627`*^9}, {
   3.574612578525875*^9, 3.574612636651577*^9}, 3.5746127303141413`*^9, 
   3.574612805865074*^9, 3.574612846347145*^9, {3.5746129812561817`*^9, 
   3.574613323941984*^9}, 3.574613453703012*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", "1.9999999999953875`", "}"}], ",", 
   RowBox[{"{", "1.9999999999999847`", "}"}], ",", 
   RowBox[{"{", "1.9999999964439188`", "}"}], ",", 
   RowBox[{"{", "1.9999999999178415`", "}"}], ",", 
   RowBox[{"{", "1.9999999999988358`", "}"}], ",", 
   RowBox[{"{", "1.9999999999995928`", "}"}], ",", 
   RowBox[{"{", "1.9999999999911064`", "}"}], ",", 
   RowBox[{"{", "1.9999999999963731`", "}"}], ",", 
   RowBox[{"{", "1.9999999999992237`", "}"}], ",", 
   RowBox[{"{", "1.999999999376367`", "}"}], ",", 
   RowBox[{"{", "1.6099611330022636`", "}"}], ",", 
   RowBox[{"{", "1.6658952190341079`", "}"}]}], "}"}]], "Output",
 CellChangeTimes->{{3.5746123888919415`*^9, 3.5746125149401627`*^9}, {
   3.574612578525875*^9, 3.574612636651577*^9}, 3.5746127303141413`*^9, 
   3.574612805865074*^9, 3.574612846347145*^9, {3.5746129812561817`*^9, 
   3.574613323941984*^9}, 3.574613453703012*^9}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", "1.4950972654126444`", "}"}], ",", 
   RowBox[{"{", "1.2631502236060952`", "}"}], ",", 
   RowBox[{"{", "1.099748278753444`", "}"}], ",", 
   RowBox[{"{", "0.976927127074678`", "}"}], ",", 
   RowBox[{"{", "0.8805119421297811`", "}"}], ",", 
   RowBox[{"{", "0.8024316149522565`", "}"}], ",", 
   RowBox[{"{", "0.7376971697794413`", "}"}], ",", 
   RowBox[{"{", "0.6830312650144732`", "}"}], ",", 
   RowBox[{"{", "0.6361778268418187`", "}"}], ",", 
   RowBox[{"{", "0.5955252077527108`", "}"}], ",", 
   RowBox[{"{", "0.9014923896927417`", "}"}], ",", 
   RowBox[{"{", "0.7039892528134529`", "}"}]}], "}"}]], "Output",
 CellChangeTimes->{{3.5746123888919415`*^9, 3.5746125149401627`*^9}, {
   3.574612578525875*^9, 3.574612636651577*^9}, 3.5746127303141413`*^9, 
   3.574612805865074*^9, 3.574612846347145*^9, {3.5746129812561817`*^9, 
   3.574613323941984*^9}, 3.574613453703012*^9}]
}, Open  ]]
},
PrintingStyleEnvironment->"Printout",
WindowSize->{967, 853},
WindowMargins->{{209, Automatic}, {Automatic, 1}},
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, 1244, 29, 354, "Input"],
Cell[1826, 53, 753, 15, 21, "Message"],
Cell[2582, 70, 754, 15, 21, "Message"],
Cell[3339, 87, 4164, 113, 252, "Output"],
Cell[7506, 202, 1179, 30, 52, "Output"],
Cell[8688, 234, 936, 18, 52, "Output"],
Cell[9627, 254, 942, 18, 31, "Output"],
Cell[10572, 274, 941, 18, 108, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)