[SciPy-user] Solving Two boundary value problems with Python/Scipy

John Hassler hasslerjc at adelphia.net
Fri Oct 6 08:32:36 CDT 2006

I don't think that there's a function specifically for boundary value 
ode problems in SciPy, but it is pretty simple to combine an ode solver 
with a zero finder to use the "shooting method."  I'm not familiar with 
how rlab or Matlab does it, but the Scilab "bvode" uses a finite 
difference method.  If the problem isn't too complex, the shooting 
method is easy to set up, and works as well as the finite difference 
method.  (In my vast experience ... well, ok, two applications ... the 
shooting method was more forgiving and less brittle than the finite 
difference method.)

What sort of bvp do you have?


Volker Lorrmann wrote:
> Hi list,
> i´m going to write a simulation, in which it would be neccessary to 
> solve an boundary value ode.
> My question is, would it be possible to solve such ode´s with any solver 
> in scipy? Otherwise i have to choose another language for this 
> simulation (scilab,rlab).
> Many thanks
> Volker

More information about the SciPy-user mailing list