In this paper we investigate the structure of the global discretization error of the implicit Euler scheme applied to systems to stiff differential equations, extending earlier work on this subject (cf. [1], [9]). We restrain our considerations to the linear, self-adjoint, constant coefficient case but—in contrast to [1], [9]—we make no assumptions about the nature of the stiff spectrum; the stiff eigenvalues may be arbitrarily distributed on the negative real axis. Our main result says that the global error of the implicit Euler scheme admits an asymptotic expansion in powers of the stepsize τ which is not asymptotically correct in the conventional sense: Near the initial pointt=0 the expansion is spoiled at theO(τ2) by ‘irregular’ error components which are, however, (algebraically) damped, such that away fromt=0 the ‘pure’ asymptotic expansion reappears. We present numerical experiments confirming this result. Our considerations should be particularly helpful for a rigorous, quantitative analysis of the structure of the full (space & time) discretization error in the PDE (method of lines) context, and thus for a sound theoretical justification of extrapolation techniques for this important class of stiff problems.