Problem

Experiment with different boundary conditions and different source terms $f(x,t)$. Plot the solution in each time step to see the time evolution of the solution. Either use the command pause each time you plot the solution (otherwise you will not see the solution) or use the commands movie and getframe to create a movie of your solution.