## E6test.txt    check of P-Q pairs
#  if following X** =0, then the relation EPQE holds

with(DEtools):
read "equations.txt":
read "E6PQ.txt":
sube9:={e9=3*r+6-e1-e2-e3-e4-e5-e6-e7-e8}:

e6pq:=proc(A,B,EE) 
        return(subs(sube9,factor(mult(EE,A,[dx,x])-mult(B,E6,[dx,x])))):
      end proc:

Xn00:=e6pq(Pn00,Qn00,subs({e1=e1-1,e2=e2-1,e3=e3-1,r=r-1},E6));
X0n0:=e6pq(P0n0,Q0n0,subs({e4=e4-1,e5=e5-1,e6=e6-1,r=r-1},E6));
Xp00:=e6pq(Pp00,Qp00,subs({e1=e1+1,e2=e2+1,e3=e3+1,r=r+1},E6));
X0p0:=e6pq(P0p0,Q0p0,subs({e4=e4+1,e5=e5+1,e6=e6+1,r=r+1},E6));
Xnnp:=e6pq(Pnnp,Qnnp,
            subs({e1=e1-1,e2=e2-1,e3=e3-1,e4=e4-1,e5=e5-1,e6=e6-1,
            e7=e7+1,e8=e8+1,e9=e9+1,r=r-1},E6));
Xppn:=e6pq(Pppn,Qppn,
            subs({e1=e1+1,e2=e2+1,e3=e3+1,e4=e4+1,e5=e5+1,e6=e6+1,
            e7=e7-1,e8=e8-1,e9=e9-1,r=r+1},E6));
Xn0p:=e6pq(Pn0p,Qn0p,
            subs({e1=e1-1,e2=e2-1,e3=e3-1,e7=e7+1,e8=e8+1,e9=e9+1},E6));
Xp0n:=e6pq(Pp0n,Qp0n,
            subs({e1=e1+1,e2=e2+1,e3=e3+1,e7=e7-1,e8=e8-1,e9=e9-1},E6));
X0np:=e6pq(P0np,Q0np,
            subs({e4=e4-1,e5=e5-1,e6=e6-1,e7=e7+1,e8=e8+1,e9=e9+1},E6));
X0pn:=e6pq(P0pn,Q0pn,
            subs({e4=e4+1,e5=e5+1,e6=e6+1,e7=e7-1,e8=e8-1,e9=e9-1},E6));
X00n:=e6pq(P00n,Q00n,subs({e7=e7-1,e8=e8-1,e9=e9-1,r=r-1},E6));
X00p:=e6pq(P00p,Q00p,subs({e7=e7+1,e8=e8+1,e9=e9+1,r=r+1},E6));