In this paper it is shown how a general two-sided orthant probability for a quadrivariate normal distribution can be evaluated by a one-dimensional numerical integral calculation. The quadrivariate normal distribution can have any covariance matrix and any mean vector. This affords a practical...