This paper presents a numerical solution to the equations of point kinetics for nuclear power reactors, a set of seven coupled differential equations that describe the temporal variation of neutron density and the concentration of delayed neutron precursors. Due to the nature of the system, we propose to numerically solve the point kinetics equations by implementing the Adams-Bashforth and Adams-Moulton methods, which are predictor-corrector schemes with their respective modifiers to increase precision. The proposed method was tested computationally for different forms of reactivity with up to six groups of delayed neutron precursors. This method was used in a recent publication to solve the inverse problem of finding the reactivity. In this work, it is shown that it can also be used for the calculation of nuclear power, that it is simple and easy to implement, and that it produces good results when compared with those in the literature for neutron population density and concentration of delayed neutron precursors.
nuclear reactor power, nuclear density, point kinetics equations, numerical methods, Adams-Bashforth method, Adams-Moulton method
 Sánchez J. On the numerical solution of the point reactor kinetics equations by generalized Runge-Kutta methods, Nuclear Science and Engineering, 103: 94-99, 1989.
 Aboanber AE, Nahla AA. Solution of the point kinetics equations in the presence of Newtonian temperature feedback by Padé approximation via the analytical inversion method, Journal of Physics A: Mathematical and General, 35(45):9609-9627, 2002b.
 Aboanber AE, Nahla AA. Generalization of the analytical inverse method for the solution of point kinetics equations, Journal of Physics A: Mathematical and General, 35(14): 3245-3263, 2002a.
 Aboanber AE. Analytical solution of the point kinetics equations by exponential mode analysis, Progress in Nuclear Energy, 42(2): 179-197, 2003.
 Kinard, M.; Allen, E. J.: Efficient numerical solution of the point kinetics equations in nuclear reactor dynamics, Annals of Nuclear Energy, 31(9): 1039-1051, 2004.
 Quintero LB. CORE: a numerical algorithm to solve the point kinetics equations, Annals of Nuclear Energy, 35(11): 2136-2138, 2008.
 Li H, Chen W, Luo L, Zhu Q. A new integral method for solving the point reactor neutron kinetics equations, Annals of Nuclear Energy, 36(4): 427-432, 2009.
 Nahla, A. A.: Taylor series method for solving the nonlinear point kinetics equations, Nuclear Engineering and Design, 241(5): 1592-1595, 2011.
 Hamada, Y. M.: Generalized power series method with step size control for neutron kinetics equations, Nuclear Engineering and Design, 241(8): 3032-3041, 2011.
 Hamada YM. Confirmation of accuracy of generalized power series method for the solution of point kinetics equations with feedback, Annals of Nuclear Energy, 55: 184-193, 2013.
 Ganapol BD. A highly accurate algorithm for the solution of the point kinetics equations, Annals of Nuclear Energy, 62: 564- 571, 2013.
 Picca P, Furfaro R, Ganapol B. A highly accurate technique for the solution of the non-linear point kinetics equations, Annals of Nuclear Energy, 58: 43-53, 2013.
 Salah A. Hassan SA. Samia.: The Analytical Algorithm for the Differential Transform Method to Solution of the Reactor Point kinetics Equations, World Applied Sciences Journal, 27(3):367-370, 2013.
 Kim HT, Park Y, Kazantzis N, Parlos A, Vista IV F, Chong KT. A numerical solution to the point kinetic equations using Taylor-Lie series combined with a scaling and squaring technique, Nuclear Engineering and Design, 272: 1-10, 2014.
 Patra A, Ray SS. A numerical approach based on Haar wavelet operational method to solve neutron point kinetics equation involving imposed reactivity insertions, Annals of Nuclear Energy, 68: 112-117, 2014.
 Leite QB, Palma AP, Vilhena MT, Bodmann EJ. Analytical representation of the solution of the point reactor kinetics equations with adaptive time step, Progress in Nuclear Energy, 70: 112-118, 2014.
 Hamada YM. Trigonometric Fourier-series solutions of the point reactor kinetics equations. Nuclear Engineering and Design, 281: 142-153, 2015.
 Razak MA, Devan K, Sathiyasheela T. The modified exponential time differencing (ETD) method for solving the reactor point kinetics equations, Annals of Nuclear Energy, 76: 193-199, 2015.
 Nahla AA. Numerical treatment for the point reactor kinetics equations using theta method, eigenvalues and eigenvectors, Progress in Nuclear Energy, 85: 756-763, 2015.
 Suescún DD, Narváez PM, Lozano PH. Calculation of Nuclear Reactivity Using the Generalised Adams Bashforth-Moulton Predictor-Corrector Method, Kerntechnik, 81(1): 86-93, 2016.
 Yun C, Xingjie P, Qing L, Kan W. A numerical solution to the nonlinear point kinetics equations using Magnus expansion, Annals of Nuclear Energy, 89: 84-89, 2016.
 Duderstadt JJ, Hamilton LJ. Nuclear Reactor Analysis, second ed. John Wiley & Sons Inc., New York, 1976