A numerical integrator for a model with a discontinuous sink term: the dynamics of the sexual phase of monogonont rotifera

We present an explicit numerical method especially adapted to approximate the solution of a nonlinear model with a discontinuous sink term: the dynamics of the sexual phase of monogonont rotifera. We show the effectiveness of this numerical technique in the simulation of the dynamics of the solutions. In particular, we show that this method provides a good approximation to the equilibrium solution of the problem. On the other hand, the numerical simulation presented validates the asymptotic behaviour of the numerical solution with regard to the theoretical one.