quad = quadfun(f, g, dg, a, b, precision=10, endpoints=False,