How to change support of multivariate integral to [0,1]^k using scipy.integrate.quad?

Question

The following k-dimension integral has dependent limits (a support dependent on gamma):

from  scipy.integrate import nquad
import numpy as np

de