function z=lcrcu(y) % %usage: z=lcrcu(y); % z=sqrt(2.^2-y.^2); endfunction