function y=f(x)
     y=(pi*pi+q(x))*sin(pi*x)+p(x)*pi*cos(pi*x);