%
% compute sigma values

% L1 point 
rho = 0.5 
rho = 0.01215
X = -0.83692;
sigmavalue = (1-rho)/abs( (X - rho)^3) + rho/abs( (X+1 - rho)^3)


% L2 point


X = -1.15568;
sigmavalue = (1-rho)/abs( (X - rho)^3) + rho/abs( (X +1 - rho)^3)

% L3 point

X = 1.00506;
sigmavalue = (1-rho)/abs( (X - rho)^3) + rho/abs( (X +1 - rho)^3)