function z=lcrcl(y)
%usage:  z=lcrcl(y);

z=-sqrt(2.^2-y.^2);