function y=fp3(x)
y=-sin(3*pi*x)*27*pi^2-27*x*pi^3*cos(3*pi*x);
