%
%  loop.m -- djm, 29 june 2005
%

inp = 1;  Ntest = 14;

max_data = zeros(Ntest,1);
for N = 1:Ntest
    L2poly
    max_data(N) = max(error);
    pause(1)
end

figure(2);  clf
plot(1:Ntest,log(abs(max_data)),'*')
ylabel('\bf log( max|f(x)-P_N(x)|)')
xlabel('\bf N')

clear inp

nh = 1:14;  detH = zeros(size(nh));
for j = nh
    detH(j) = det(hilb(j));
end

figure(3);  clf
plot(nh,log(abs(detH)),'*')
title('/bf det(hilb(N))')
figure(2)