from scipy.stats import mvn
import numpy as np
np.random.seed(1)
inf = 1e5
mvn.mvnun([10.0, 800.0, 800.0],
 [1e5, 1e5, 1e5],
 [ 10, 1000, 1000],
 [[   5, -25, -10],
        [-25 , 1250,  500],
        [ -10,  500, 1100]])