# NOT RUN {
p <- poly.from.zeros(-2:5)
ip <- integral(p)
ipv <- integral(p, limits = c(-2, 5))
plot(polylist(p, deriv(p)))
x <- polynom()
H <- polylist(1, x)
for(n in 2:10)
H[[n+1]] <- x * H[[n]] - (n-1)*H[[n-1]]
solve(deriv(H))
# }
Run the code above in your browser using DataLab