In my study, I often need to verify that two operator is symmetric i.e. [P,Q]=PQ-QP=0, where A and D are operator polynomial such like  D²+4u+u_xD^-1 multiply with D³+uD+u_x,where D is differential operator.

I tried to use the Ore_package which can easily deal with the operator polynomial without integral(i.e. D^-1 term), so in my case , how to deal with operator with both differential and integral?