Computes the Owen's T function or the bivariate
normal integral using one of the following:
modified Euler's arctangent series, tetrachoric
series, or Vasicek's series. For the methods,
see Komelj, J. (2023)
(or reprint with better typography)
and Vasicek, O. A. (1998) .