def expected_value_integral(func, bounds, parameters): return quad(func, *bounds, parameters)[0] / (bounds[1] - bounds[0])