function y=myfunc(x)
global alpha
y=sin(alpha*x.^2).*cos(x);

