library(resultant)
x <- qlone(1)
qspray1 <- 2*x^4 + x^3 - 3*x^2 - x + 1
qspray2 <- x^2 - 5*x + 10
division <- univariateDivision(qspray1, qspray2)
Q <- division[["Q"]]; R <- division[["R"]]
qspray1 == Q*qspray2 + R # should be TRUE
Run the code above in your browser using DataLab