a=(tan(pi*(rand(1000)-0.5)));
s=reshape(a,1000000,1);
z=mean(a);
y=1./(pi*(1+t.^2));