Set up diffusion MMS test case

Added Python script to run for several resolutions
and produce plot of convergence

Currently showing O(1), rather than O(2). I must
have broken something...