L = 10; table = zeros(L-1,3);
alpha = {4,'exp'};

for n=2.^(2:L)
    table(log2(n)-1,1) = n;
    for mu=1:2
        [w,c]=SummationFormula(n,alpha{mu});
        m = length(w); w = sqrt(w);
        A = zeros(m,m);
        for j=1:m
            for k=1:m
                A(j,k) = w(j)*w(k)/...
                    ((c(j)+c(k)-1)*(c(j)+c(k))/2-(c(k)-1));
            end
        end
        table(log2(n)-1,mu+1) = normest(A,eps);
    end
end
table