Furthermore, an approximate algorithm is pro-posed to compute the kinetic spectrum of each species, thus making the method effective as well for reactions with fractional orders.