function y = uexact(x,t)

    y = exp(-t)*sin(pi*x);