read `MultiAlmkvistZeilberger.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 a3, annihilatating the (a1,a2,2*a3) mixed moment of the Ekhad-Zeilberger tri-variate distribution`): print(`where the shift operator is a3 is denoted by A3 is`): ope:=MAZ(1,F*x^a1*y^a2,z^2*(1-c)*(1+2*c)/(c+1),[x,y,z],a3,A3,{a1,a2},c)[1]; lprint(ope): print(`-----------------------------------------`): print(`This took`, time(), `seconds. `): quit: