# NOT RUN {
(x <- polynomial())
(p <- polynomial(c(1, 5, 2, 2)/10))
plot(p, xlim = 0:1, ylim = 0:1)
P <- p
for(j in 1:7) {
lines(P, col = j)
P <- p(P)
}
(r <- solve(p-x))
segments(r, 0, r, p(r), lty = "dashed")
# }
Run the code above in your browser using DataLab