Generation of random bivariate normal deviates and computation of related integrals

A method which transforms two random variables having rectangular distributions into a pair of bivariate normal deviates with prescribed covariance matrix is described. The same transformation is used for integrating the bivariate normal distribution over areas which are the intersection of the domain outside an equiprobability ellipse and a sector determined by two lines through the point of gravity of the normal distribution.