Schnakenberg Turing Patterns

ut = Du∇²u + a − u + u²v  |  vt = Dv∇²v + b − u²v

Step: 0
Turing instability:
Dv/Du ≫ 1 needed
Homogeneous state:
u* = a+b, v* = b/(a+b)²
Short-range activation
(u), long-range inhibition (v)
→ spots when b/a² > 1
→ stripes near boundary
Display