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

read "equations.txt":
read "E5PQ.txt":
subE5r:={r=(e1+e2+e3+e4+e5+e6+e7+e8-6)/3}:

e5pq:=proc(A,B,EE) 
        return(subs(subE5r,factor(mult(EE,A,[dx,x])-mult(B,E5,[dx,x])))):
      end proc:

Xn0:=e5pq(Pn0,Qn0,subs({e1=e1-1,e2=e2-1,e3=e3-1,r=r-1},E5));
X0n:=e5pq(P0n,Q0n,subs({e4=e4-1,e5=e5-1,e6=e6-1,r=r-1},E5));
Xp0:=e5pq(Pp0,Qp0,subs({e1=e1+1,e2=e2+1,e3=e3+1,r=r+1},E5));
X0p:=e5pq(P0p,Q0p,subs({e4=e4+1,e5=e5+1,e6=e6+1,r=r+1},E5));