The standard bivariate normal distribution function is given by \Phi( \bf b}, \rho)= \frac{1}{2 \pi \sqrt{1- \rho^{2}}} \int_{- \infty}^{b_1} \int_{- \infty}^{b_2} e^{-(x^{2}-2 \rho xy +y^{2})/(2(1- \rho^{2}))}dydx }
This implementation mainly differs from the original code in two regards;
|
|
|
|