def integral_bounded_mp(fn, lb, ub):
    (integral, _) = mp.quad(fn, [lb, ub], error=True)
    return integral
def integral_bounded_mp(fn, lb, ub):
    (integral, _) = mp.quad(fn, [lb, ub], error=True)
    return integral