If you have external impulses, then none of the ode* functions can be used for the entire timespan, unless the forces can be interpolated with degree two or more, such as using interp1 'spline'.
If they cannot be made smooth enough then you need to stop the ode* function at the each location that fails C2 continuity.
This is because the mathematics of optimal point selection to calculate new points from is dependent upon intepolating polynomials that require continuity.
