You should go to the source...this dude Professor Alan Genz has been researching how to do this and other multivariate integrals numerically since the 1980's. All of the code that others have implemented are derived from his algorithm and papers. His code can compute the CDF and expectation of the multivariate normal and T distributions for dimensions up to 1000.
http://www.math.wsu.edu/faculty/genz/software/software.html
I've also written code to call those subroutines from Java: Compute the multivariate normal CDF in Java