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

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

X1p:=e3cpq(P1p,Q1p,subs(e1=e1+1,E3c));
X1n:=e3cpq(P1n,Q1n,subs(e1=e1-1,E3c));
X3pp:=e3cpq(P3pp,Q3pp,subs(e3=e3+2,E3c));
X3nn:=e3cpq(P3nn,Q3nn,subs(e3=e3-2,E3c));
X6pp:=e3cpq(P6pp,Q6pp,subs(e6=e6+2,E3c));
X6np:=e3cpq(P6nn,Q6nn,subs(e6=e6-2,E3c));