Die Umgebung ist Python 3.7.6, Windows 10.
Die Mehrfachintegration integr.nquad scheint nicht gut zu funktionieren und das Ergebnis sollte 1.000 sein (da es sich um eine simultane Wahrscheinlichkeitsdichtefunktion handelt), aber ich erhalte 1.326.
Wie kann ich die Mehrfachintegration korrekt durchführen?
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