Given a system of time dependent ordinary differential equations, dot y i = f i(c 1, c 2, …, y 1, y 2, …, t) , where c k are rate parameters, we simultaneously solve for both y i and a set of sensitivity functions, ϱy i ϱc k , over all times t. These partial derivatives measure the sensitivity of the solution with respect to changes in the parameters c k . Often these parameters are not accurately known. An example is given from atmospheric chemical kinetics using constant as well as time varying (diurnal) rate parameters. For the purposes of this paper, our calculations considered both first- and second-order contributions to Δ y with respect to Δc. It is found that second-order sensitivity terms can be highly significant, but tend to be too costly for present widespread application.