Revision history [back]

Many thanks! finally i decided to use the scipy odeint, since supports the backwards integration and is the fastest ode solver that i found

fa= lambda y,t: [y[1],y[0]^2-y[1]]

t=scipy.linspace(10, 0, 1/0.01)
len(t)
time l=integrate.odeint(fa,[0,.1],-t)