read MarkovWZ:

r:=5:
print(`The acceleration formula for Zeta(`, r, `) is `):
ACCpaper(((z+x)!/(z+2*x+1)!)^r,z,x);

print(`The error in using 30 terms is`):
print(evalf(Appxf(((z+x)!/(z+2*x+1)!)^r,z,x,30)[2]-Zeta(r))):


quit;