add the plotting infrastructure
add a test of the variable coefficient solver that uses Dirichlet
BCs on phi