this runs fine but it is not the program I want; instead want to use DMDA and MatSetValuesStencil(); i.e. see ksp example ex46.c
this is [close to: add actual exact solution] a good first PDE code for the book; note extensive comments about runtime at top; restructuring to have a linear PDE chapter 2 is appropriate