if (e > 1.0e-6) {
fail("norm of 1st derivative of cum minus Gaussian: " + e + "\ntolerance exceeded");
}
// check that the second derivative of cum is normal.derivative
temp = D2.applyTo(yi);
for (i=0; i < yd.size(); i++) {
diff.set(i, yd.get(i) - temp.get(i));
}
e = Utilities.norm(diff, h);
if (e > 1.0e-4) {