\\ generates Table 3.2 dec=28; default(realprecision,dec); tol=10^(2.0-dec); L=16; vals=vector(L); res=[1.0,1.0]; vals[1]=res[1]; for(nu=2,L, x0=concat(res[2],0*res[2]); res=OperatorNorm(x0,tol); vals[nu]=res[1]); L2=2*L-1; vv=vector(L2,j, if(j%2==1, vals[(j+1)/2], 0)); for(j=2,L, forstep(k=j, L2+1-j, 2, vv[k]=vv[k]+1/(vv[k+1]-vv[k-1]))); result=vector(L,j, vv[2*j-1]); for(j=1,L, print( vals[j]," ",result[j]));