read MarkovWZ:

for r from 2 to 8 do
print(`The best kernel,as far as appx. after 20 terms, for Zeta(`,r,`) is`):
lu:=BestKernelAppx(r,z,x,20):
H:=lu[1]:
print(H);
print(`the error is`, evalf(Appx(H,z,x,20)[2]-Zeta(r))):
print(`The error with 60 terms is`, evalf(Appx(H,z,x,60)[2]-Zeta(r))):
print(`The correponding MWZ pair is`):
print(Markov(H,z,x)):
print(`The corresponding Acceleration scheme is`):
print( Acc(H,z,x)):
od:

print(`The whole thing too`, time() , `seconds of CPU time `):

quit: