def dwdt(b, w, p, m, d, dx2): """ dw/dt of Modified Klausmeier """ return p - w - w * b ** 2 + d * laplacian(b, dx2)
def dbdt(b, w, p, m, d, dx2): """ db/dt of Modified Klausmeier """ return w * b ** 2 - m * b + laplacian(b, dx2)