System.out.println("Testing bandec and banmul");
double[][] c = new double[50][17];
ranmat(c);
Bandec clu = new Bandec(c,9,7);
ranvec(r);
clu.solve(r,y);
Bandec.banmul(c,9,7,y,rr);
sbeps = 5.e15;
localflag = maxel(vecsub(r,rr)) > sbeps;
globalflag = globalflag || localflag;
if (localflag) {