Looking into using `dense_output=True` in `scipy.integrate.solve_ivp`
Looking into using
dense_output=Trueinscipy.integrate.solve_ivp