function [Y] = rcauchy(n)

X=rand(n,1);
Y=-tan(%pi*(X-0.5));
endfunction

