I do not know if something is wrong is my computer. I use the function pmvnorm from the package TruncatedNormal to compute the cdf of multivariate
pmvnorm
TruncatedNormal