The RightEuclidean(Poly1, Poly2, A) calling sequence returns a list [m, S] where m is a positive integer and S is an array with m elements storing the subresultant sequence of the first kind of Poly1 and Poly2.

If Poly1 and Poly2 are fraction-free, and the commutation rule of the Ore algebra A are also fraction-free, then the output of RightEuclidean is fraction-free.

If the optional fourth argument to the RightEuclidean command c1 is specified, the first co-sequence of Poly1 and Poly2 is assigned to c1, so that:

and c1[m+1] Poly2 is a least common left multiple (LCLM) of Poly1 and Poly2.

If the optional fifth argument to the RightEuclidean command c2 is specified, the second co-sequence of Poly1 and Poly2 is assigned to c2, so that:

and c1[m+1] Poly2 = - c2[m+1] Poly1 is an LCLM of Poly1 and Poly2.

The RightPseudoRemainder(Poly1, Poly2, A) calling sequence returns the right pseudo-remainder R of Poly1 and Poly2 such that:

where the degree of R is less than that of Poly2, pquo is the right pseudo-quotient, and mult (the multiplier) is the ${\left(\mathrm{degree}\left(\mathrm{Poly1}\right)-\mathrm{degree}\left(\mathrm{Poly2}\right)+1\right)}^{\mathrm{th}}$ sigma power of the leading coefficient of Poly2.

If the fourth argument 'mult' is specified, it is assigned the multiplier defined above. If the fifth argument 'pquo' is specified, it is assigned the right pseudo-quotient defined above.

The RightPseudoQuotient(Poly1, Poly2, A) calling sequence returns the right pseudo-quotient Q of Poly1 and Poly2 such that:

where the degree of the right pseudo-remainder prem is less than that of Poly2 and mult (the multiplier) is the ${\left(\mathrm{degree}\left(\mathrm{Poly1}\right)-\mathrm{degree}\left(\mathrm{Poly2}\right)+1\right)}^{\mathrm{th}}$ sigma power of the leading coefficient of Poly2.

If the fourth argument 'mult' is specified, it is assigned the multiplier defined above. If the fifth argument 'prem' is specified, it is assigned the right pseudo-remainder defined above.

$\mathrm{with}\left(\mathrm{OreTools}\right)\:$

$\mathrm{with}\left(\mathrm{OreTools}\left[\mathrm{FractionFree}\right]\right)\:$

$A≔\mathrm{SetOreRing}\left(x\,'\mathrm{differential}'\right)$

$A≔\mathrm{UnivariateOreRing}\left(x\,\mathrm{differential}\right)$

$\mathrm{Ore1}≔\mathrm{OrePoly}\left(1\,3x\,x^2-\left(1+x\right)\,23x+1\right)$

$\mathrm{Ore1}≔\mathrm{OrePoly}\left(1\,3x\,x^2-x-1\,23x+1\right)$

$\mathrm{Ore2}≔\mathrm{OrePoly}\left(x\,x\,x^2+x+1\right)$

$\mathrm{Ore2}≔\mathrm{OrePoly}\left(x\,x\,x^2+x+1\right)$

$V≔\mathrm{RightEuclidean}\left(\mathrm{Ore2}\,\mathrm{Ore1}\,A\,'\mathrm{c1}'\,'\mathrm{c2}'\right)$

$V≔\left[4\,S\right]$

$m≔V\left[1\right]\;$$S≔V\left[2\right]$

$m≔4$

$S≔S$

$\mathrm{print}\left(S\right)$

$\left[\begin{array}{cccc}\mathrm{OrePoly}\left(1\,3x\,x^2-x-1\,23x+1\right)& \mathrm{OrePoly}\left(x\,x\,x^2+x+1\right)& \mathrm{OrePoly}\left(-x^5+x^4+49x^3+7x^2-20x\,2x^5-17x^4+32x^3-14x^2-20x-1\right)& \mathrm{OrePoly}\left(x^8+3x^7-209x^6+1005x^5+61x^4-587x^3+971x^2+35x-20\right)\end{array}\right]$

$$\begin{gathered} \mathbf{for}i\mathbf{to}m\mathbf{do} \\ \mathrm{W1}≔\mathrm{Multiply}\left(\mathrm{c1}\left[i\right]\,\mathrm{Ore2}\,A\right); \\ \mathrm{W2}≔\mathrm{Multiply}\left(\mathrm{c2}\left[i\right]\,\mathrm{Ore1}\,A\right); \\ W≔\mathrm{Add}\left(\mathrm{W1}\,\mathrm{W2}\right); \\ C≔\mathrm{Minus}\left(W\,S\left[i\right]\right); \\ \mathrm{print}\left(C\right) \\ \mathbf{end}\mathbf{do}\: \end{gathered}$$

$\mathrm{OrePoly}\left(0\right)$

$\mathrm{OrePoly}\left(0\right)$

$\mathrm{OrePoly}\left(0\right)$

$\mathrm{OrePoly}\left(0\right)$

$\mathrm{W3}≔\mathrm{Multiply}\left(\mathrm{c1}\left[m+1\right]\,\mathrm{Ore2}\,A\right)\:$

$\mathrm{W4}≔\mathrm{Multiply}\left(\mathrm{c2}\left[m+1\right]\,\mathrm{Ore1}\,A\right)\:$

$\mathrm{Add}\left(\mathrm{W3}\,\mathrm{W4}\right)$

$\mathrm{OrePoly}\left(0\right)$

$\mathrm{Poly1}≔\mathrm{OrePoly}\left(x\,0\,x+1\,x+2\right)$

$\mathrm{Poly1}≔\mathrm{OrePoly}\left(x\,0\,x+1\,x+2\right)$

$\mathrm{Poly2}≔\mathrm{OrePoly}\left(x\,x+1\,x^2-1\right)$

$\mathrm{Poly2}≔\mathrm{OrePoly}\left(x\,x+1\,x^2-1\right)$

$B≔\mathrm{SetOreRing}\left(x\,'\mathrm{shift}'\right)$

$B≔\mathrm{UnivariateOreRing}\left(x\,\mathrm{shift}\right)$

$R≔\mathrm{RightPseudoRemainder}\left(\mathrm{Poly1}\,\mathrm{Poly2}\,B\right)$

$R≔\mathrm{OrePoly}\left(x^5+x^4-3x^3+4x\,-2x^4-6x^3-x^2+9x+6\right)$

$R≔\mathrm{RightPseudoRemainder}\left(\mathrm{Poly1}\,\mathrm{Poly2}\,B\,'m'\,'Q'\right)$

$R≔\mathrm{OrePoly}\left(x^5+x^4-3x^3+4x\,-2x^4-6x^3-x^2+9x+6\right)$

$m$

$x^4+2x^3-x^2-2x$

$Q$

$\mathrm{OrePoly}\left(x^3+2x^2-2x-4\,\left(x^2-1\right)\left(x+2\right)\right)$

$\mathrm{evalb}\left(\mathrm{factor}\left(\mathrm{Add}\left(\mathrm{Multiply}\left(Q\,\mathrm{Poly2}\,B\right)\,R\right)\right)=\mathrm{factor}\left(\mathrm{ScalarMultiply}\left(m\,\mathrm{Poly1}\right)\right)\right)$

$\mathrm{true}$

$Q≔\mathrm{RightPseudoQuotient}\left(\mathrm{Poly1}\,\mathrm{Poly2}\,A\right)$

$Q≔\mathrm{OrePoly}\left(x^3-2x^2-8x-3\,\left(x^2-1\right)\left(x+2\right)\right)$

$Q≔\mathrm{RightPseudoQuotient}\left(\mathrm{Poly1}\,\mathrm{Poly2}\,B\,'m'\,'R'\right)$

$Q≔\mathrm{OrePoly}\left(x^3+2x^2-2x-4\,\left(x^2-1\right)\left(x+2\right)\right)$

$m$

$x^4+2x^3-x^2-2x$

$R$

$\mathrm{OrePoly}\left(x^5+x^4-3x^3+4x\,-2x^4-6x^3-x^2+9x+6\right)$

$\mathrm{evalb}\left(\mathrm{factor}\left(\mathrm{Add}\left(\mathrm{Multiply}\left(Q\,\mathrm{Poly2}\,B\right)\,R\right)\right)=\mathrm{factor}\left(\mathrm{ScalarMultiply}\left(m\,\mathrm{Poly1}\right)\right)\right)$

$\mathrm{true}$

Li, Z. "A subresultant theory for Ore polynomials with applications." Proc. of ISSAC'98. pp. 132-139. Edited by O. Gloor. ACM Press, 1998.