function y=fp4(x)
y=-cos(3*pi*x)*108*pi^3+81*x*pi^4*sin(3*pi*x);