A numerical simulation for the dynamics of a model that describes the sexual phase of Monogonont Rotifera reproduction is presented. The simulation is carried out by means of a numerical method based on the integration along the characteristic curves. The numerical experiments cover two basic situations: the existence of an asymptotic stable equilibrium state and the existence of an stable periodic solution. Our results are in agreement with the theoretical analysis made by Calsina and Ripoll (J. Math. Biol. 45 (2002) 22).