def romberg(fun, a, b, releps=1e-3, abseps=1e-3):