read `SMAZ.txt`:
F:=exp(-x^2/2-y^2/2-z^2/2 -(x*y+x*z+y*z)*c)/sqrt(2*Pi)^3/(c+1)*(1-c^2)*sqrt(1+2*c):


print(`The pure recurrence operator, in n,  annihilatating the (2n,2n,2n) mixed moment of the Ekhad-Zeilberger tri-variate distribution`):
print(`where the shift operator is n is denoted by N is`):

ope:=MAZ(1,F,x^2*y^2*z^2*((1-c)*(1+2*c)/(c+1))^3,[x,y,z],n,N,{c},[y1,y2],z2)[1];

lprint(ope):


print(`-----------------------------------------`):
print(`This took`, time(), `seconds. `):

quit: