read MarkovWZ:

for r from 1 to 8 do
H:=binomial(n,k)^r:
print(`If the kernel is`,H, `then the MWZ pair is the following`):
lu:=Markov(H,k,n):
print(`consisting of a matrix of dimension`, nops(lu[2])):
print(lu):
od:

quit: