I have the following code
def f(x, y, ro): return x*y*ro def bounds_y(): return [0, 0.5] def bounds_x(y): return [0, 1-2*y] scipy.integrate.nqu