taylor(sin, 0, 4)  #=> -0.1666666  0.0000000  1.0000000  0.0000000
taylor(exp, 1, 4)  #=>  0.04166657 0.16666673 0.50000000 1.00000000 1.00000000
f <- function(x) log(1+x)
p <- taylor(f, 0, 4)
p                     # log(1+x) = 0 + x - 1/2 x^2 + 1/3 x^3 - 1/4 x^4 +- ...
                      # [1] -0.250004  0.333334 -0.500000  1.000000  0.000000
## Not run: 
# x <- seq(-1.0, 1.0, length.out=100)
# yf <- f(x)
# yp <- polyval(p, x)
# plot(x, yf, type = "l", col = "gray", lwd = 3)
# lines(x, yp, col = "red")
# grid()## End(Not run)
Run the code above in your browser using DataLab