try time dependent ode solver s in scipy