function y=cauchy(x)
y = 1./(1+25.*x.^2);
