% generates Table 3.3

n0 = 64;
table = []; v = zeros(n0,1); v(1)=1;
for n = n0*2.^(1:5)
    m = ceil(n^(4/3));
    [s0,v,u] = OperatorNorm([v;0*v],1e-14);
    s1 = ThirdOrderCorrection(s0,v,u,m);
    table = [table;s0 s1];
end

table 