% generatore per la distribuzione di cauchy
clear;
clf;
gamma=1;
n=10000;
rd=rand(1,n);
xx=-gamma*tan(pi*(rd-0.5));
a=-7;
b=-a;
dx=(b-a)/20;
xgr=a:dx:b;
nd=histc(xx,xgr);
nd=nd/(n*dx);
y=(gamma/pi)*1./(gamma^2+(xgr+dx/2).^2);
bar(xgr,nd);
hold on
plot(xgr,y,'rs--','LineWidth',3,'MarkerSize',2);
hold off