%discauchy.m
%distribuzione di Cauchy

figure(1);
clear;
gamma=1;
n=6;
x=[-n:0.1:n];
y=gamma./(pi*(gamma^2+x.^2));
plot(x,y)
axis([-n n 0 .4]);
xlabel('x')
ylabel('\rho(x)')
text(3.,.35,'\Gamma=1.0')
