% Script to test our zeta function

disp 'zeta(1) should fail:'
zeta(1)

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