A discretization scheme is introduced for a set of convection-diffusion equations with a non-linear reaction term, where the convection velocity is constant for each reactant. This constancy allows a transformation to new spatial variables, which ensures the global stability of discretization. Convection-diffusion equations are notorious for their lack of stability, arising from the algebraic interaction of the convection and diffusion terms. Unexpectedly, our implemented numerical algorithm proves to be faster than computing exact solutions derived for a special case, while remaining reasonably accurate, as demonstrated in our runtime and error analysis.MSC class: 65M08 (primary); 65M12, 35C05 (secondary).