

function y=f(x,a)
y=a./((a.^2)+(x.^2));


