Poly_interp polint = new Poly_interp(h,s,K);
h[0]=1.0;
for (int j=1;j<=JMAX;j++) {
s[j-1]=q.next();
if (j >= K) {
double ss=polint.rawinterp(j-K,0.0);
if (abs(polint.dy) <= eps*abs(ss)) return ss;
}
h[j]=h[j-1]/9.0;
}
throw new IllegalArgumentException("Too many steps in routine qromo");