Double integration of copula entropy in Python using scipy.integrate.dblquad

Question

The code below uses double integration with scipy.integrate.dblquad to calculate the differential entropy, c*np.log(c) of a copula density function