


function w = u1(x,t)

     w = 40/(pi.^2)*sin(pi*x/10)*exp(-(pi/10)^2*t);
