Chebyshev cheb1 = new Chebyshev(bj1,aa,bb,NN);
Chebyshev cint=cheb1.integral();
for (i=0;i<MM;i++) {
x=-5.0+i;
y[i]=cint.eval(x,m);
yy[i]=b.j0(aa)-b.j0(x);
// System.out.printf(y[i] << " %f\n", yy[i] << endl;
}
System.out.printf("Chebyshev (integral): Maximum discrepancy = %f\n", maxel(vecsub(y,yy)));
sbeps=2.0e-6;
localflag = maxel(vecsub(y,yy)) > sbeps;