Division of two polynomials
qsprayDivision(qsprayA, qsprayB)
qsprayA
: a qspray
object, the dividendqsprayB
: a qspray
object, the divisorA list with two qspray
objects, the quotient and the remainder.
library(qspray) x <- qlone(1) y <- qlone(2) z <- qlone(3) B <- x*y^2 + z*x^2 + 1 A <- B * (x^2*y^2*z^2 - 3) + x*y divis <- qsprayDivision(A, B) B * divis[["Q"]] + divis[["R"]] == A # should be TRUE