(* 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[      5525,        173]
NotebookOptionsPosition[      4792,        147]
NotebookOutlinePosition[      5474,        170]
CellTagsIndexPosition[      5431,        167]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["\<\
(* Ne like ions F- to Ar(+8) nach Kimball 30.04.2011 *)

Clear[k0,k2,sig0,sig2,c]; n = 1.;
 c = {k0 -> 0.98, k2 -> 1.10,
     sig0 -> 0.3, sig2 -> 0.3};

T = 2.25*n*k0/P^2+4.5*(n+1.0)*k2/R^2;

Vee=3.0*n*sig0/P+(n+1.0)*6.0*sig2/R+16/(P+R)+24/Sqrt[8/3]/(P+R);

Vne=-3.0*n*z/P-8*z/(P+R);
      
func = T + Vee + Vne /. c;

t=Table[FindMinimum[func, {P, 0.09},{R, 0.5}], {z,9,18}]
                        

u = t[[{1,2,3,4,5,6,7,8,9,10},2]];
v = u[[{1,2,3,4,5,6,7,8,9,10},1]];
w = u[[{1,2,3,4,5,6,7,8,9,10},2]];
x = Table[N[P+R /. v[[i]] /. w[[i]]],{i,1,10}]
(* Kriterium f\[UDoubleDot]r Overlap der Aussenwolken: muss <= 0 sein *)
y = Table[N[R/(P+R)-Sqrt[2/3] /. v[[i]] /. w[[i]]],{i,1,10}]
\
\>", "Input"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "100.05896511275769`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.16165430228051653`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.7945844298616241`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "128.8896014714696`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.14385224998030405`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.6657800718264498`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "161.5284292507662`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.12953151862351328`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.5734098795474808`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "197.97064938870105`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.11777446574264976`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.503806123520395`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "238.21330316387844`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.10795604648631822`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.4494156060652847`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "282.25446843807487`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.09963716107832683`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.4057101660416223`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "330.09284282975267`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.09250097066732882`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.36980529721361505`"}]}], "}"}]}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "381.727512843433`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.08631344497946385`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.3397734506169151`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "437.15781895067903`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.0808980489184948`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.3142761740722062`"}]}], "}"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "496.3832731471916`"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"P", "\[Rule]", "0.07611945652756581`"}], ",", 
       RowBox[{"R", "\[Rule]", "0.2923546151786939`"}]}], "}"}]}], "}"}]}], 
  "}"}]], "Output"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "0.9562387321421406`", ",", "0.8096323218067538`", ",", 
   "0.7029413981709941`", ",", "0.6215805892630447`", ",", 
   "0.5573716525516029`", ",", "0.5053473271199491`", ",", 
   "0.46230626788094387`", ",", "0.42608689559637897`", ",", 
   "0.395174222990701`", ",", "0.3684740717062597`"}], "}"}]], "Output"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"0.014451176314461023`", ",", "0.0058273973698705905`", ",", 
   RowBox[{"-", "0.0007673024424148522`"}], ",", 
   RowBox[{"-", "0.0059723589314839565`"}], ",", 
   RowBox[{"-", "0.010184304356257168`"}], ",", 
   RowBox[{"-", "0.013662284061398577`"}], ",", 
   RowBox[{"-", "0.016582491705709024`"}], ",", 
   RowBox[{"-", "0.01906898076333685`"}], ",", 
   RowBox[{"-", "0.021211474288567622`"}], ",", 
   RowBox[{"-", "0.023076805623282737`"}]}], "}"}]], "Output"]
}, Open  ]]
},
PrintingStyleEnvironment->"Printout",
WindowSize->{623, 717},
WindowMargins->{{240, Automatic}, {Automatic, 92}},
PrintingCopies->1,
PrintingPageRange->{2, 6},
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) (November 20, 2012)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[579, 22, 714, 25, 422, "Input"],
Cell[1296, 49, 2614, 74, 212, "Output"],
Cell[3913, 125, 353, 7, 52, "Output"],
Cell[4269, 134, 507, 10, 52, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)