A method of calculating space-dependent reactor kinetics by direct numerical integration of the time-dependent multigroup equations is described. In particular, this method is applied to study the buildup and decay of neutron flux in a subcritical reactor with a step-function point external source.