% FILE d012hlp.m begins.
s=0;
disp('n     p(n)     Error')
for k=0:10
%    fprintf('%3d  %12.8f  %12.4e\n',k,s,s-pi);
disp(  [k,s,s-pi]);
    s=s+(16^(-k))*(4/(8*k+1)- 2/(8*k+4)-1/(8*k+5)-1/(8*k+6));
end
% FILE d012hlp.m ends.
