Internal Function
It approximates the trivariate normal integral.
Maintainer: Giampiero Marra giampiero.marra@ucl.ac.uk