ptnorm function

Compute the probability for the trivariate normal function