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).

https://github.com/shoshievass/PublicFiles/blob/master/Julia/Inv_diffeq_problem_example.ipynb