x=[.01:.01:100];

for i=1:length(x),
    y2(i)=erf(x(i)/50)/x(i)^2;
end


plot(x,y);