The modelling of chemical reactions relies upon expressions being supplied which give the rate of reaction as a function of local temperature and local concentrations. When power law rate expressions are used with concentrations to integer powers, the analysis and numerical solution of these equations is made easier. In this note the numerical difficulties that arise when the order is non-integer are described and one method is proposed which overcomes them in a differential-algebraic equation model of a catalytic monolith reactor.