}
Ran myran = new Ran(17);
Poly_interp z = new Poly_interp(x,y,3);
for (i=0;i<N;i++) {
xx[i]=myran.doub();
yy[i]=z.interp(xx[i]); // interpolated values
dyy[i]=z.dy; // Estimated errors
zz[i]=1.0+xx[i]*(1.0+xx[i]*(1.0+xx[i]*(1.0+xx[i]))); // Actual Values
}
System.out.printf(" Poly_interp: Max. estimated error: %f\n", maxel(dyy));
System.out.printf(" Poly_interp: Max. actual error: %f\n", maxel(vecsub(zz,yy)));