Numerical simulation of stiff systems with a diagonal splitting method