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