function y = f_poly(x)
y = sin(pi/2*x);
