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