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

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