% Generates Table 10.3

format long e;

table = [];
for h = 1./[8 16 32 64 128]
    p = poisson([0,0],[10,1],0,[1,1,0,0],'Cholesky',h);
    table = [table; 20/h-1 2/h-1 h p];
end
err = (table(1:end-1,4)-table(2:end,4))./(3*table(2:end,4));
errq = err(1:end-1)./err(2:end);

err = [NaN;err]; errq = [NaN;NaN;errq];
table = [table err errq]