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