Bessj1_Adapt bj1 = new Bessj1_Adapt();
Adapt adp = new Adapt(0.0);
for (i=0;i<N;i++) {
x1=1.0*i;
x2=x1+2.0;
y[i]=adp.integrate(bj1,x1,x2);
yy[i]=bj0.funk(x1)-bj0.funk(x2);
}
System.out.printf("Adapt: Maximum discrepancy = %f\n", maxel(vecsub(y,yy)));
localflag = maxel(vecsub(y,yy)) > sbeps;
globalflag = globalflag || localflag;