In reply to Jun Luo.

You can use the PROBBNRM function to compute the bivariate normal CDF. For higher dimensions, there have been some discussions on the SAS Support Communities. One discusses the Monte Carlo approach for t distributions, but also applies to MVN distributions.