pbnorm function

Compute the probability for the bivariate normal function