On 04/04/2008, Doreen Mbabazi <doreen at aims.ac.za> wrote: > y = odeint(f,initial_y,t,args=(p,)) > > def S(t,p): > v = y[:,2] > return v Here you are running odeint only once. You must run it every time S is evaluated. Anne