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

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

 Xap:=se3pq(Pap,Qap,subs(a=a+1,SE3));
 Xan:=se3pq(Pan,Qan,subs(a=a-1,SE3));
 Xcp:=se3pq(Pcp,Qcp,subs(c=c+1,SE3));
 Xcn:=se3pq(Pcn,Qcn,subs(c=c-1,SE3));
 Xgp:=se3pq(Pgp,Qgp,subs(g=g+2,SE3));
 Xgn:=se3pq(Pgn,Qgn,subs(g=g-2,SE3));