phi(t+1,n) = phi(t,n) + I(n) - <>. I(n) = f(phi(t,n)) - phi(t,n) + D(<> - phi(t,n)). where <<*>> denotes average of neighbour cells. f(x) = 1.3 * tanh(x)