function y=fcn(x,a,b)
%    y=3*(a*x+b*x.^3)./(1+abs(x.^3));
%y=a*(1.5*x+0.5*sin(b*x));
y=a*(x+b*x.^3);
