=> 1 + x + x^2 + x^3 + ….
fn_1 <- function(x) {1/(1-x)}
taylor(fn_1, x0 = 0, n = 4)
## [1] 1.000029 1.000003 1.000000 1.000000 1.000000
=> 1 + x + x^2/2 + x^3/6 + …
fn_2 <- function(x) {exp(x)}
taylor(fn_2, x0 = 0, n = 4)
## [1] 0.04166657 0.16666673 0.50000000 1.00000000 1.00000000
=> x - x^2/2 + x^3/3 - x^4/4 + …
fn_3 <- function(x) {log(1 + x)}
taylor(fn_3, x0 = 0, n = 4)
## [1] -0.2500044 0.3333339 -0.5000000 1.0000000 0.0000000