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

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

X1pp:=e3bpq(P1pp,Q1pp,subs(e1=e1+2,E3b));
X1nn:=e3bpq(P1nn,Q1nn,subs(e1=e1-2,E3b));
X3pp:=e3bpq(P3pp,Q3pp,subs(e3=e3+2,E3b));
X3nn:=e3bpq(P3nn,Q3nn,subs(e3=e3-2,E3b));
X5pp:=e3bpq(P5pp,Q5pp,subs(e5=e5+2,E3b));
X5nn:=e3bpq(P5nn,Q5nn,subs(e5=e5-2,E3b));
X6p:=e3bpq(P6p,Q6p,subs(e6=e6+1,E3b));
X6n:=e3bpq(P6n,Q6n,subs(e6=e6-1,E3b));
X3p5p:=e3bpq(P3p5p,Q3p5p,subs({e3=e3+1,e5=e5+1},E3b));
X3n5n:=e3bpq(P3n5n,Q3n5n,subs({e3=e3-1,e5=e5-1},E3b));
X3p5n:=e3bpq(P3p5n,Q3p5n,subs({e3=e3+1,e5=e5-1},E3b));
X3n5p:=e3bpq(P3n5p,Q3n5p,subs({e3=e3-1,e5=e5+1},E3b));
X1pp6n:=e3bpq(P1pp6n,Q1pp6n,subs({e1=e1+2,e6=e6-1},E3b));
X1nn6p:=e3bpq(P1nn6p,Q1nn6p,subs({e1=e1-2,e6=e6+1},E3b));