from scipy.stats import mvn
import numpy as np
low = np.array([-10, -10])
upp = np.array([0.1, -0.2])
mu = np.array([-0.3, 0.17])
S = np.array([[1.2, 0.35], [0.35, 2.1]])
p, i = mvn.mvnun(low, upp, mu, S)
print(p)
0.2881578675080012
Multivariate Normal CDF in Python using scipy