function y=myfun(x)
y=x./(x.^2+1);