function y = u1exact(x,t)

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