Doing repeated Euler-steps in a simulation, Differentialequations.jl

https://docs.sciml.ai/DiffEqDocs/stable/basics/integrator/

Or use discrete callbacks