import numpy as np
from scipy.stats import norm
from scipy.stats import multivariate_normal as mvn

mean = np.array([0, 0])
covariance = np.array([[1, 0.4], [0.4, 1]])
mvna = mvn(mean, covariance)
print(mvna.cdf(np.array([-2, -4])))

multivariate_normal_frozen' object has no attribute 'cdf'