final double tolerance = 1.0e-15;
for (final double element : rho) {
for (Integer sgn=-1; sgn < 2; sgn+=2) {
final BivariateNormalDistribution bvn= new BivariateNormalDistribution(sgn*element);
final double expected = 0.25 + Math.asin(sgn*element) / (2*Math.PI) ;
final double realised = bvn.op(x,y);
if (Math.abs(realised-expected)>=tolerance)
fail(" bivariate cumulative distribution\n"
+ " rho: " + sgn*element + "\n"
+ " expected: " + expected + "\n"