https://github.com/cran/pracma
Raw File
Tip revision: c79a04b5074656b36e591191eb8137b70a349932 authored by Hans W. Borchers on 30 June 2014, 00:00:00 UTC
version 1.7.0
Tip revision: c79a04b
taylor.R
##
##  t a y l o r . R  Taylor Series Approximation
##


taylor <- function(f, x0, n = 4, ...) {
    if (!is.numeric(x0))
        stop("Argument 'x0' must be a numeric value or vector.")
    if (!is.numeric(n) || floor(n) != ceil(n) || length(n) != 1)
        stop("Order 'n' must be an scalar integer between 1 and 8.")

    fun <- match.fun(f)
    f <- function(x) fun(x, ...)

    if (n > 8) {
        # n <- 8
        warning("Order 'n' is too high; should not be greater than 'n=8'.")
    } else if (n <= 0) {
        n <- 4
        warning("Order 'n' is too low; will be reset to the default 'n=4'.")
    }

    T <- f(x0)
    for (i in 1:n) {
        T <- polyadd(T, fderiv(f, x0, i)/fact(i) * polypow(c(1, -x0), i))
    }
    return(T)
}
back to top