function y=f(x)
y=1./(1+x.^2);
% y=1./(2+sin(x+1));
%y=(x.*(x+2)+x.^4)./x;
return