read WZtriple:
t0:=time():
print(`Selberg's Integral`):
F:=(x*y*z)^n*((1-x)*(1-y)*(1-z))^m*((x-y)*(x-z)*(y-z))^(2*k);
print(`WZ3(F,n,x,y,z,[n,m,k],N)`);
print(WZ3(F,n,x,y,z,[n,m,k],N));
t1:=time():
print(`It tool`, t1-t0, ` seconds `):
t0:=time():
print(`Mehta's Integral `):
F:=exp(-x^2-y^2-z^2)*((x-y)*(x-z)*(y-z))^(2*n);
print(`WZ3(F,n,x,y,z,[n],N)`);
print(WZ3(F,n,x,y,z,[n],N));
t1:=time():
print(`It tool`, t1-t0, ` seconds `):
t0:=time():
print(`Dyson's Constant-Term Identity `):
F:=((x-y)*(x-z)*(y-z)/x/y/z)^(2*n)/x/y/z;
print(`WZ3(F,n,x,y,z,[n],N)`);
print(WZ3(F,n,x,y,z,[n],N));
t1:=time():
print(`It tool`, t1-t0, ` seconds `):