function z=f(x)
   z=sin(x);
%  z=1./(x.^2+1);