L'environnement est python 3.7.6, windows10.
L'intégration multiple integr.nquad ne semble pas bien fonctionner et le résultat devrait être de 1.000 (car c'est une fonction de densité de probabilité simultanée), mais elle renvoie 1.326.
Comment puis-je obtenir une intégration multiple correctement?
code
import numpy as np
from scipy import integrate
def f_xy(x, y): if 0 <= y <= 1 and 0 <= x-y <= 1: return 4 * y * (x - y) else: return 0
integrate.nquad(f_xy, [[-np.inf, np.inf],[-np.inf, np.inf]])[0]
Out: 1.326
Recommended Posts