I'd like the use the cumulative distribution function of the multivariate normal as a moment condition in a GMM estimation. It seems like there is a mata command for this:
Code:
mvnormal(L, U, M, V)
, but I don't know how to use this function in vanilla stata. For many reasons, I'd like to use the GMM estimation command rather than write it myself in mata. Is there a way to do this?

If more detail is necessary, the parameters I'd like to estimate are parameters of the multivariate normal, M and V