Check ode45 for a certain value and start a new solver when these limitations are reached

I use the ode45 solver for my rigid multibody system. I wanted to limit the movement of the bodies to a certain max and min. I know it is possible to stop the solver completely when a condition is reached but I want it to compute the remaining values and if possible mirror the second output of the body so that a bouncing is simulated.

In my first attempt I solved the system and swapped the unwanted values with my max/min but that doesn't make sense. Could I start another ode45 with initial conditions that are my last values of the first ode45 outputs and just mirror the second output with a minus?

Answer 1

[Torsten]

That's exactly what the Event facility of the ODE integrators does.

Look at the "ballode" example:

https://uk.mathworks.com/help/matlab/math/ode-event-location.html