I think it works! Next step, compare to an analytic solution! Do some
convergence testing, etc.