Finite Element Method for Schnakenberg Model

In the study, solution of an initial boundary value problem for the Schnakenberg reaction-diffusion model is considered in numerical meaning. The approximate solution is assumed to be a finite series, some of exponential form of the cubic B-spline basis. After adapting the boundary data, the system is integrated in time variable by using Crank-Nicolson implicit method. The resultant iteration algorithm is initiated by the aid of the initial data adaption. The numerical results are compared with the analytical solution by using the relative error calculation that measures the ratio of the difference in the two successive time levels to the previous time level to validate the accuracy and the reliability.