Numerical simulation of reaction–diffusion equations on spherical domains

The two-variable reaction diffusion equations on the spherical domain is considered and simulated, using the semi-implicit Euler finite difference method. It is shown that the method keeps the kinetics from overshooting the stable branches when a large time step is used in the simulation.