Analysis of Turing Instability for Biological Models

Reaction-diffusion equations are often used to model biological phenomena. This type of system can produce stable spatial patterns from an uniform initial distribution. This phenomenon is known as Turing instability. This paper presents an analysis of the Turing instability for three biological models: a) Schnakenberg model, b) glycolysis model and c) blood coagulation model. The method of lines is used in the numerical solution, and the spatial discretization is done using a finite difference scheme. The resulting system of ordinary differential equations is then solved by an adaptive integration scheme with the use of SciPy, a Python package for scientific computing.