Gierer-Meinhardt

Activator-inhibitor: ∂a/∂t = Da∇²a + ρa²/h - μa
∂h/∂t = Di∇²h + ρa² - νh