Abstract: In this study, a numerical study of the stochastic dynamics of a nuclear reactor is carried out using the Euler-Maruyama method to resolve the stochastic point kinetic equations. The method was tested using different uncertainties, initial conditions and different forms of reactivity with one and six groups of delayed neutron precursors using up to 5000 Brownian trajectories for each numerical results. The average expected values, standard deviations of neutron density and delayed neutron precursors were calculated in each computational experiment. It was compared with other methods reported in literature and on average with the deterministic model of point kinetics, demonstrating that the Euler-Maruyama method for this proposal, apart from being easy to implement is extremely precise in obtaining expected values.