The package pbivnorm will certainly solve your problem. It allows a vectorized call to the cdf of the bivariate normal distribution. -- View this message in context: http://r.789695.n4.nabble.com/Bivariate-normal-integral-tp4571018p4639041.html Sent from the R help mailing list archive at Nabble.com.