integral = quad(f,a,b,tol)