For a multi-variate Gaussian, it is straight-forward to compute the CDF or PDF given a point. rv = scipy.stats.multivariate_normal(mean, cov) and then rv.pdf(point) or rv.cdf(upper)