I’d like to simulate a complex mecanical motion that goes from one phase to another in a strict interval of values for $deltaX$, where I’m finding a position X as a fucntion of time, which is calculated by different differential equations in two periodic domains, each for a relative change of X of length a and b. Is there a way to adapt ode or odeint from scipy, or should I program an alogrithm by hand?

Thanks!

Searched for ways to ude odeint, but it seems that I must program an Euler/Runge-Kutta method by hand.