Abstract: A new numerical method of solving the point reactor neutron kinetics equations by using the technique of first order Taylor polynomials for the approximation of neutron density in integral of one step is presented. The full implicit formulation is used to overcome the stiffness of the equations and to guarantee the necessary exactness of solution. The numerical results show that for reasonable reactivity inputs, a sufficient accuracy can be achieved in first order polynomial approach presented here. When the reactivity is given, this method can provide a straighforward computation procedure for solving reactor dynamics problems.

Key words: point reactor neutron kinetics, stiffness, reactivity, implicit