function x=cauchyrnd(a,b,M,N)
F = rand(M,N);
x = tan((F-.5)*pi)*b+a;