Numerical solution of ordinary differential equations arising from chemical kinetics