Abstract We propose a simple and efficient way to calculate bivariate normal probabilities. The algorithm is based on a formula for the partial derivative of the bivariate probability with respect to the correlation coefficient. Keywords: Bivariate normal, computation, approximation