}
// Test Bessjy (y1)
System.out.println("Testing Bessjy (y1)");
sbeps=1.5e-13;
for (i=0;i<N;i++) zz4[i]=bess.y1(x3[i]);
System.out.printf("Bessjy: Maximum discrepancy = %f\n", maxel(vecsub(zz4,yy4)));
localflag = maxel(vecsub(zz4,yy4)) > sbeps;
globalflag = globalflag || localflag;
if (localflag) {
fail("*** Bessjy (y1): Incorrect function values");