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