% §3.5 Summation Formulas: The Real Case
% MATLAB Session p. 64

zeta4 = pi^4/90;
alpha = 4; f = inline('1./x.^alpha','x','alpha');
error = [];
for n = [2,20,200]
    [w,c] = SummationFormula(n,alpha);
    error = [error; abs(zeta4 - w*f(c,alpha)')];
end
error