System.out.println("Testing Chebyshev (derivative)");
Chebyshev cder=cheb.derivative();
for (i=0;i<MM;i++) {
x=-5.0+i;
y[i]=cder.eval(x,m);
yy[i]=0.5*(b.jn(N-1,x)-b.jn(N+1,x));
// System.out.printf(y[i] << " %f\n", yy[i] << endl;
}
System.out.printf("Chebyshev (derivative): Maximum discrepancy = %f\n", maxel(vecsub(y,yy)));
localflag = maxel(vecsub(y,yy)) > sbeps;
globalflag = globalflag || localflag;