On the numerical approximation of the bivariate normal (tetrachoric) correlation coefficient