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

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

Xap:=se5pq(Pap,Qap,subs({a=a+1},SE5));
Xan:=se5pq(Pan,Qan,subs({a=a-1},SE5));
Xbp:=se5pq(Pbp,Qbp,subs({b=b+1},SE5));
Xbn:=se5pq(Pbn,Qbn,subs({b=b-1},SE5));
Xpp:=se5pq(Ppp,Qpp,subs({p=p+1},SE5));
Xpn:=se5pq(Ppn,Qpn,subs({p=p-1},SE5));
Xqp:=se5pq(Pqp,Qqp,subs({q=q+1},SE5));
Xqn:=se5pq(Pqn,Qqn,subs({q=q-1},SE5));
Xcp:=se5pq(Pcp,Qcp,subs({c=c+1},SE5));
Xgp:=se5pq(Pgp,Qgp,subs({g=g+2},SE5));
Xag:=se5pq(Pag,Qag,subs({a=a+1,g=g-2},SE5));
Xrag:=se5pq(Prag,Qrag,subs({a=a-1,g=g+2},SE5));
Xcpq:=se5pq(Pcpq,Qcpq,subs({c=c-1,p=p+1,q=q+1},SE5));  
Xrcpq:=se5pq(Prcpq,Qrcpq,subs({c=c+1,p=p-1,q=q-1},SE5));

## remark: cn and gn not checked because Qcn and Qgn are not computed