double[] yy3=buildVector(y3),yy4=buildVector(y4),zz3=new double[N],zz4=new double[N];
// Test Bessjy (y0)
System.out.println("Testing Bessjy (y0)");
sbeps=1.e-14;
for (i=0;i<N;i++) zz3[i]=bess.y0(x3[i]);
System.out.printf("Bessjy: Maximum discrepancy = %f\n", maxel(vecsub(zz3,yy3)));
localflag = maxel(vecsub(zz3,yy3)) > sbeps;
globalflag = globalflag || localflag;
if (localflag) {
fail("*** Bessjy (y0): Incorrect function values");