(* Ne atom nach Kimball *)

Clear[k0,k2,sig0,sig2,u,c]; n = 1.;  z = 10.0; z2 = z*z;
 c = {k0 -> 0.963, k2 -> 1.2,
     sig0 -> 0.30, sig2 -> 0.35};

T = 2.25*n*k0/P^2+4.5*(n+1.0)*k2/R^2 /. c ;
 r = P+R; 
r2 = r*r;  
ad = 8./3.;

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

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

N[t = FindMinimum[func, {P, 0.26112},{R, 1.2769}], 10]
                        

u = t[[2]];

N[T /. u,10]

N[Vne /. u,10]

N[Vee /. u /. c, 10]

N[-(Vne+Vee)/T /. u /. c, 8] (* Virial Theorem *)

N[r /.u /. c, 8]