Integral doble con límites variables en python Scipy + sympy (?)

El problema matemático completo esaquí.

Brevemente quiero integrar una función con una integral doble. La integral interna tiene límites20 yx-2, mientras que el exterior tiene límites22 y30.

Sé que con Scipy puedo calcular la integral doble conscipy.integrate.nquad. Me gustaría hacer algo como esto:

def f(x, y):
    return (x ** 2 + y ** 2)
res = sp.integrate.nquad(f, [[22, 30], [20, x-2]])

¿Es posible? Tal vez usando tambiénsympy?