I want to solve a boundary value ODE problem, where I know some boundary conditions at each end of the interval. As far as I know, odeint is just for initial value problems where you know boundary conditions at one end of the interval. Is there a relaxation solver in SciPy?