Meinhardt-Gierer (activator-inhibitor) model solved in a nonregular domain with constant boundaries.

(rho/i)*((a^2)/(1 + kappa*a^2)) - mu_a * a + rho_a rho*((a^2)/(1+kappa*a^2)) - mu_i *i 1