Abstract: | In this article, we discuss a scheme for dealing with Neumann and mixed boundary conditions using a compact stencil. The resulting compact algorithm for solving systems of nonlinear reaction‐diffusion equations is fourth‐order accurate in both the temporal and spatial dimensions. We also prove that the standard second‐order approximation to zero Neumann boundary conditions provides fourth‐order accuracy when the nonlinear reaction term is independent of the spatial variables. Numerical examples, including an application of this algorithm to a mathematical model describing frontal polymerization process, are presented in the article to demonstrate the accuracy and efficiency of the scheme. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005 |