Salve, ragazzi, ho un dubbio.
Sto analizzando i diversi metodi di integrazione proposti da Py.

Scipy ci viene in soccorso con quad, dblquad e in generale con nquad.

Poi ci sono altri metodi di integrazione, quello della somma di riemann, quello di simpson o quello di romberg.

Tuttavia questi ultimi metodi sono relativi ad una sola variabile, almeno negli esempi che fino ad ora ho sempre visto.




Ora, a me serve costruire uno script che data una funzione, qualsiasi essa sia, mi risolva l'integrale doppio di questa.

E sono partito da un esempio basic, tipo: \int_{0}^{\frac{1}{2}}\int_{0}^{1-2y} xydxdy.

Posto che l'integrale ad una variabile ormai è chiaro, come posso fare una doppia integrazione qui usando i due cicli for??

Il problema nasce soprattutto nell'estremo superiore dell'integrale in x (cioè 1-2y)