an example with inhomogeneous Dirichlet BCs -- this seems to work
(converges second order).  The notebook derives the RHS and the
boundary conditions.