function y=u(x)
%y=exp(-100.*x.*(1-x));
y=sin(2*pi.*x)./(1+4*pi^2);