function y=fn(x)
  y=pi*pi*sin(pi*x)a+exp(sin(pi*x));