function y=fp(x)
y=sin(3*pi*x)+x*3*pi*cos(3*pi*x);
