robrighter

import pycxsimulator
from pylab import *

n = 200     # size of grid: n * n
Dh = 1. / n # spatial resolution, assuming space is [0,1] * [0,1]
Dt = 0.02   # temporal resolution

a, b = 12.0, 16.0 # parameter values

Du = 0.0001 # diffusion constant of u