FINAL u: sin(pi*x)*exp(-k*pi**2*t)