function main_summation

s    = 0;
nmax = round(1e3);
for j = 1:nmax
  s = s + 1/(j*j);
  %fprintf(1,'%10d   %20.16f   %20.16f\n',j,s,pi*pi/6-s)
  %keyboard
end
fprintf(1,'Exact    answer = %20.16f\n',pi*pi/6)
fprintf(1,'Computed answer = %20.16f\n',s)
fprintf(1,'Error           = %20.16f\n',s - pi*pi/6)
fprintf(1,'Estimated error = %20.16f\n',1/(nmax+1))

keyboard

return