(* Li-Atom and homologues 1s(2)2s concentric spherical,with exchange correction page 40; 4.12.2011 *)
Clear[z,sig]; sig=0.29;
T=9/(4*P^2)+9*4/(8*Q^2);
Vne=-3*z/P-3*z/(2*Q);
Vee=+(2/Q)*(3/2-(3/10)*(P/Q)^2)-3*sig*P^2/Q^3 +3*sig/P;
func=T+Vne+Vee; Vir = (Vne+Vee)/T;

(* Etot, P, Q *)
t = Table[FindMinimum[func,{P,0.6},{Q,4}],{z,3,18}]

(* Virial ratio *)
Table[N[Vir /. t[[{z-2 },2]],10],{z ,3,18}]
Table[N[Vne /. t[[{z-2 },2]],10],{z ,3,18}]
Table[N[Vee /. t[[{z-2 },2]],10],{z ,3,18}]

(* Etot in eV *)
Table[N[t[[{z-2},1]]*27.21138,10],{z,3,18}];

PA=Table[t[[{z-2},2,1]],{z,3,18}];
QA=Table[t[[{z-2},2,2]],{z,3,18}];

Do[Show[Graphics[{Circle[{0,0},P/. PA[[i]]],Circle[{0,0},Q /. QA[[i]]]}],{AspectRatio -> Automatic,Axes -> True,GridLines -> Automatic, PlotRange -> {{-6,6},{-6,6}}, Frame -> True}], {i,1,16}]