function y=uexact(x,t);
y=exp(-4*t*pi^2)*sin(2*pi*x);