[SciPy-User] Trajectory Integration via scipy.integrate or PyDSTool?

[Sara Gallian]

Hi all,
I'm trying to integrate particle trajectories in a finite domain, and I need to make them "reflect" at the boundaries (i.e. reverse the velocity perpendicular to the plane they collided with). I started by using  scipy.integrate.odeint , but since the integration steps are variable, simply checking the position and reversing the velocity won't work!
Can anyone suggest the quickest way to obtain this? would Vode be able to handle this, or should I try to learn to use  PyDSTool?
I'm running pretty late on a deadline, so any suggestion is more than appreciated :)
Thanks!
Sara