subroutine bivariate_normal_cdf_values ( n_data, x, y, r, fxy )
!*****************************************************************************80
!
!! BIVARIATE_NORMAL_CDF_VALUES returns some values of the bivariate normal CDF.
!
! Discussion:
!
! FXY is the probability that two variables A and B, which are
! related by a bivariate normal distribution with correlation R,
! respectively satisfy A <= X and B <= Y.
!
! Mathematica can evaluate the bivariate normal CDF via the commands:
!
! <