Integral dupla com limites variáveis em python Scipy + sympy (?)
O problema matemático completo éaqui.
Resumidamente, quero integrar uma função com uma integral dupla. A integral interna tem limites20 ex-2, enquanto o exterior tem limites22 e30.
Eu sei que com Scipy eu posso calcular a integral dupla comscipy.integrate.nquad. Eu gostaria de fazer algo assim:
def f(x, y):
return (x ** 2 + y ** 2)
res = sp.integrate.nquad(f, [[22, 30], [20, x-2]])
É possível? Talvez usando tambémsympy?