compute cell areas
add function to evaluate source
add function to evaluate residual, boundary conditions, etc.; finish first working version