I was wondering if anybody has written R code to compute the cdf of a multivariate (or at least a bivariate) normal distribution with given covariance structure.