Двойной интеграл с переменными границами в Python Scipy + sympy (?)

Полная математическая проблемаВот.

Вкратце я хочу интегрировать функцию с двойным интегралом. Внутренний интеграл имеет границы20 а такжеx-2в то время как внешнее имеет границы22 а также30.

Я знаю, что со Сципи я могу вычислить двойной интеграл сscipy.integrate.nquad, Я хотел бы сделать что-то вроде этого:

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

Является ли это возможным? Может быть, используя такжеsympy?