## Z3test.txt     check of P-Q pairs
#  if following X** =0, then the relation EPQE holds
with(DEtools):
read "equations.txt":
read "Z3PQ.txt":

z3pq:=proc(A,B,EE)  
        return(factor(mult(EE,A,[dx,x])-mult(B,Z3,[dx,x]))): end proc:

X0p1p:=z3pq(P0p1p,Q0p1p,subs({A0=A0+1,A1=A1+1},Z3));
X0p3p:=z3pq(P0p3p,Q0p3p,subs({A0=A0+1,A3=A3+1},Z3));
X0p2p:=z3pq(P0p2p,Q0p2p,subs({A0=A0+1,A2=A2+1},Z3)); 
X0p2n:=z3pq(P0p2n,Q0p2n,subs({A0=A0+1,A2=A2-1},Z3));
X2p3p:=z3pq(P2p3p,Q2p3p,subs({A2=A2+1,A3=A3+1},Z3));
X2n3p:=z3pq(P2n3p,Q2n3p,subs({A2=A2-1,A3=A3+1},Z3));
X0n1p:=z3pq(P0n1p,Q0n1p,subs({A0=A0-1,A1=A1+1},Z3));
X0p1n:=z3pq(P0p1n,Q0p1n,subs({A0=A0+1,A1=A1-1},Z3));
X0n1n:=z3pq(P0n1n,Q0n1n,subs({A0=A0-1,A1=A1-1},Z3));
X0pp:=z3pq(P0pp,Q0pp,subs({A0=A0+2},Z3)); 
X0nn:=z3pq(P0nn,Q0nn,subs({A0=A0-2},Z3)); 
X1pp:=z3pq(P1pp,Q1pp,subs({A1=A1+2},Z3));
X1nn:=z3pq(P1nn,Q1nn,subs({A1=A1-2},Z3));
X2pp:=z3pq(P2pp,Q2pp,subs({A2=A2+2},Z3));
X2nn:=z3pq(P2nn,Q2nn,subs({A2=A2-2},Z3));
X3pp:=z3pq(P3pp,Q3pp,subs({A3=A3+2},Z3));
X3nn:=z3pq(P3nn,Q3nn,subs({A3=A3-2},Z3));