Gray-Scott

du/dt = Du∇²u − uv² + F(1−u)
dv/dt = Dv∇²v + uv² − (F+k)v