An implicit-explicit (IMEX) method is developed for the numerical solution of reaction-diffusion equations with pure Neumann boundary conditions. The corresponding method of lines scheme with finite differences is analyzed: ...