The neutron point kinetics equations model the time-dependent behavior of nuclear reactors and are used to understand the dynamics of nuclear reactor operations (e.g. power fluctuations caused by control rod motions during start-up and shut-down procedures). These equations are deterministic in nature, and therefore can provide only average values of the modeled populations. However, the actual dynamical process is stochastic: the neutron density and the delayed neutron precursor concentrations vary randomly with time.

A system of stochastic differential equations can accurately model the random behavior of the neutron density and the precursor concentrations in a point reactor. Due to the issue of stiffness, numerical implementation may be challenging. This research focused on studying the influence of stochastic fluctuations upon the results of the neutron point kinetics equations, and to find an approach to obtain a nonstiff solution for the stochastic neutron point kinetics equations.