  We propose an approach to solve the stochastic neutron point kinetics
equations using an adaptation of the diagonalization-decomposition method
(DDM). This new approach (Double-DDM) yields a nonstiff solution for the
stochastic formulation, allowing the calculation of the neutron and precursor
densities at any time of interest without the need of using progressive time
steps. We use Double-DDM to compute results for stochastic problems with
constant, linear, and sinusoidal reactivities. We show that these results
strongly agree with those obtained by other approaches established in the
literature. We also compute and analyze the first four statistical moments of
the solutions.
