read MarkovWZ:

Digits:=200:

print(`For Zeta(2), using`, H2, `the 40^th term of the weighted-average with`):
print(binomial(z+2*x,x)*binomial(x+z,z)^2):
print(`is `):
lu:=WtAve(Markov(H2,z,x),z,x,binomial(z+2*x,x)*binomial(x+z,z)^2,40):

print(lu):
print(`delta is`, delt(lu,Zeta(2))):

print(`For Zeta(3), using`, H3, `the 40^th term of the weighted-average with`):
print(binomial(z+2*x,x)^2*binomial(x+z,z)^2):
print(`is `):
lu:=WtAve(Markov(H3,z,x),z,x,binomial(z+2*x,x)^2*binomial(x+z,z)^2,40):

print(lu):
print(`delta is`, delt(lu,Zeta(3)));

quit: