Abstract
The processes of diffusion and reaction play essential roles in numerous system dynamics. Consequently, the solutions of reaction–diffusion equations have gained much attention because of not only their occurrence in many fields of science but also the existence of important properties and information in the solutions. However, despite the wide range of numerical methods explored for approximating solutions, the adoption of block methods is yet to be investigated. Hence, this article introduces a new two-step third–fourth-derivative block method as a numerical approach to solve the reaction–diffusion equation. In order to ensure improved accuracy, the method introduces the concept of nonlinearity in the solution of the linear model through the presence of higher derivatives. The method obtained accurate solutions for the model at varying values of the dimensionless diffusion parameter and saturation parameter. Furthermore, the solutions are also in good agreement with previous solutions by existing authors.