This paper develops an algorithm to calculate the Brownian multivariate normal probability subject to any preset error tolerance criteria. The algorithm is founded upon the computational simplicity of the tridiagonal structure of the inverse of the Brownian correlation matrix. Compared with...