def create_solver(dydt): solver = ode(dydt).set_integrator('lsoda', method='bdf', rtol=1e-2) return solver