% §3.2 Extrapolation
% MATLAB Session p. 51 (cont. of p. 50)

L = 11;
vals = [];
for nu = 1:L,
   n=2^(nu-1); An=A(1:n,1:n); vals=[vals;normest(An,5e-12)];
end;

L2 = 2*L-1; vv = zeros(L2,1); vv(1:2:L2) = vals;
for j=2:L
  k=j:2:L2+1-j; vv(k)=vv(k)+1./(vv(k+1)-vv(k-1));
end;
result = vv(1:2:L2);
