T  =  
        2.25/R^2;                          
eefact  =  1.0406758;              
Vee  =  eefact*1.2/R;           
Vnn  =  Z*Z/(2*r);                 
Vne  =  -2*Z*(3 - r^2/R^2)/R;   

func  =  T  +  Vee  +  Vnn  +  Vne;     
 
    t  =  Table[
        FindMinimum[func,  {R, {0.4, 13.0}},  {r, {0.4, 6.0}}], {Z, 0.3, 
          2.2, 0.1}] ;
v  =  Table[Z, {Z, 0.3, 2.2, 0.1}]
  
u  =  
      t[[{1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 
              19, 20}, 1]] + v^2  
w =  
      t[[{1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 
            19, 20}, 2]]

ListPlot[{{v[[1]], 
          u[[1]]}, {v[[2]], u[[2]]}, {v[[3]], 
          u[[3]]}, {v[[4]], u[[4]]}, {v[[5]], 
          u[[5]]}, {v[[6]], u[[6]]}, {v[[7]], 
          u[[7]]}, {v[[8]], u[[8]]}, {v[[9]], 
          u[[9]]}, {v[[10]], u[[10]]}, {v[[11]], 
          u[[11]]}, {v[[12]], u[[12]]}, {v[[13]], 
          u[[13]]}, {v[[14]], u[[14]]}, {v[[15]], 
          u[[15]]}, {v[[16]], u[[16]]}, {v[[17]], 
          u[[17]]}, {v[[18]], u[[18]]}, {v[[19]], 
          u[[19]]}, {v[[20]], u[[20]]}}]