% Generates Table 8.2

format long e;

t = 0.424; table = [];
for h = 1./[4 8 16 32 64]
    u = heat([0,0],[1,1],t,0,0,[5,0,0,0],h);
    table = [table; 2/h^2 2/h-1 h u];
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]