On the numerical solution of the point reactor kinetics equations