The Complete Flux scheme (CFS) (J.H.M. ten Thije Boonkkamp et al., J. Sci. Comput. 46 (2011) 47–70) is an extension of the widely used exponential di¿erence scheme for advection-di¿usion-reaction equations. In the present paper we provide a rigorous proof that the convergence order of this scheme is 2 for all grid Péclet numbers, whereas that of the exponential scheme reduces to 1 for high grid Péclet numbers in the presence of source terms. The performance of both schemes is compared in two case studies: a model system and a real-world model of a parallel-plate glow discharge. The results indicate that the usage of CFS allows a considerable reduction of the number of grid points that is required to obtain the same accuracy. The MATLAB/Octave source code that has been used in these studies has been made available.