(*Be - 
        Atom and homologues 1  s (2) 2  
          s (2) concentric spherical, 
      with exchange correction page 40;   4.12  .2011 *) 
Clear[z, sig];  sig = 0.29;
T = 9/(4*P^2) + 9*4/(4*Q^2);
Vne = -3*z/P - 3*z/Q;
Vee = +(2/Q)*(3 - (6/10)*(P/Q)^2) - 
          6*sig*P^2/Q^3 + 3*sig/P;
func = T + Vne + Vee;  
    Vir  =  (Vne + Vee)/T;
  
t  =  
      Table[FindMinimum[func, {P, 0.4}, {Q, 3}], {z, 4, 
          18}]
  
Table[
      N[Vir  /.  t[[{z - 3 }, 2]], 10], {z , 4, 
        18}]
Table[
      N[Vne  /.  t[[{z - 3 }, 2]], 10], {z , 4, 
        18}]
Table[
      N[Vee /.  t[[{z - 3 }, 2]], 10], {z , 4, 
        18}]
  
Table[
        N[t[[{z - 3}, 1]]*27.21138, 10], {z, 4, 
          18}];
PA = 
        Table[t[[{z - 3}, 2, 1]], {z, 4, 
            18}];
QA = 
        Table[t[[{z - 3}, 2, 2]], {z, 4, 
            18}];

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