% Script to test our crude zeta summation function. Barnett 4/6/14

help zeta      % check the documentation adequate

try
  zeta(1)
catch me
  'zeta(1) failed as it should'
end

fprintf('error in zeta(4) = %.3g\n', zeta(4) - pi^4/90)

disp('please wait 10 seconds...')
fprintf('error in zeta(2) = %.3g\n', zeta(2) - pi^2/6)
