x[i]=(i)*pi/(N-1);
y[i]=sin(x[i]);
}
yp1=1.0;
ypn=-1.0;
Spline_interp z = new Spline_interp(x,y,yp1,ypn);
Ran myran = new Ran(17);
for(i=0;i<N;i++) {
xx[i]=pi*myran.doub();
yy[i]=z.interp(xx[i]);
zz[i]=sin(xx[i]);
}
System.out.printf(" Spline_interp: Max. actual error: %f\n", maxel(vecsub(zz,yy)));
sbeps=1.e-5*maxel(y);
localflag = maxel(vecsub(zz,yy)) > sbeps;