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