tidy up
write node locations more accurately, to eliminate this as possible source of error
This version *seems* to work.  It uses manufactured Dirichlet solution for test.  Convergence is not evident, though perhaps this is because quadrature is not accurate enough for right-hand side?