this first version of c2heat compiles and runs, but blows up; plan = back up to c2poisson.c and make it work with an exact solution and convergence plot