The continuous indication of the neutron density and its rate of change are important for the safe startup and operation of reactors. The best way to achieve this is to obtain analytical solutions of the neutron kinetics equations because none of the developed numerical methods can well satisfy the need for real-time or even super-time computation for the safe startup and operation of reactors in practice. In this paper, an accurate analytical solution of point kinetics equations with one-group delayed neutrons and an extraneous neutron source is proposed to calculate the change in neutron density, where the whole process from the subcritical stage to critical and supercritical stages is considered for step reactivity insertions. The accurate analytical solution can also be used as a benchmark of all numerical methods employed to solve stiff neutron kinetics equations.
