y[i]=func_BaryRat(x[i]);
}
ymax=maxel(y);
Ran myran = new Ran(17);
for (j=0;j<M;j++) {
BaryRat_interp z = new BaryRat_interp(x,y,j);
for (i=0;i<N;i++) {
xx[i]=2.0*myran.doub();
yy[i]=z.interp(xx[i]); // interpolated values
zz[i]=func_BaryRat(xx[i]); // Actual values
}
dy=maxel(vecsub(zz,yy));
System.out.printf(" BaryRat_interp: Order: %d Max. error: %f\n", j, dy);
sbeps=pow(10.0,-(j+4.0)/2)*ymax;