A method of solving reactor kinetics equations by using third-order Hermite polynomials for the approximation of neutron density (or power density) and delayed group precursors in finite time intervals is presented. When the reactivity is given, including the case of the feedback reactivity being a function of neutron density, this method can provide a straightforward procedure for computing reactor dynamics problems.