I have an ODE that I can solve, but I need to also be able to solve its “inverse”. That is, my normal ode is defined as
du_dx = f(x, u, p) , and I can solve that, but I also need to be able to solve
dx_du = g(x, u, p) (because I’m trying to solve a boundary value problem w/ a shooting algorithm).
I’m hitting a weird boundary error when I try to do the inverse problem though, and I’m stuck on why. Any advice on what to try/check or what might help resolve this would be really helpful (and greatly appreciated).