Hi This bites me on Python 2.5 (never noticed it before) scipy.integrate.odeint changes the argument y0. The fix would be to make a copy in line 118, as is also done with the t argument. I opened ticket #370 for this bug. Best, tnoo -- tnoo answer at tnoo.net