Comparing Numerical Methods for the Solution of Stiff Systems of ODEs Arising in Chemistry