How to use `nquad` with specific weight function?

Question

Right now I need to calculate the integral of f(x)(x-a)(y-c), with x from a to b and y from c to d numerically. I\'ve checked the docs of scipy.integrate.nquad, and