function y=f(x)
%y=-2+p(x)*(2*x-1)+q(x)*(x^2-x);
y=sin(2*pi.*x);
%y=0;