摘要

This article deals with the numerical solution to some models described by the system of strongly coupled reaction-diffusion equations with the Neumann boundary value conditions. A linearized three-level scheme is derived by the method of reduction of order. The uniquely solvability and second-order convergence in L-2-norm are proved by the energy method. A numerical example is presented to demonstrate the accuracy and efficiency of the proposed method.