T  =  
        1.125/R^2;
r1sq  =  x1^2 + y1^2 + z1^2;
Vne  =  -(1.5  -  0.5*r1sq/R^2)/
          R;

func  =  T  +  Vne;
   
N[t  =  
        FindMinimum[
          func,  {R, {1.0, 1.5}}, 
     {x1, {0.4, 
              2}}, {y1, {0.1, 2}}, {z1, {0.2, 2}}, {MaxIterations  ->  
              100}], 30]