read MarkovWZ:
H:=(z+x)!/(2*x+z+1)!;

t0:=time():

for r from 2 to 6 do
 print(`Using`, H^r, `as a kernel for Zeta(`,r,`),`):
 lu:=BestDelt(H^r,Zeta(r),z,x,40);
 print(`The best delta amongst the Potential Function of the corresponding`):
 print(`MWZ pair , followed by the pair [n,k], for n+k=40 is`):
 print(lu):

t1:=time():
print(`This took `, t1-t0, `seconds of CPU time`):
t0:=t1:
od:

quit: