stehfest_inversion

Source: R/laplace_inversion_utilities.R

Function to calculate the Laplace inverse transform using the stehfest algorithm

stehfest_inversion(t, coeffs, fun, arg1 = NULL, arg2 = NULL,
  arg3 = NULL)

Arguments

t

Numeric vector with the time
coeffs

Numeric vector with the coefficients
fun

function to be used
arg1, arg2, arg3

Optional arguments to be used by fun

Value

Numeric vector

See also

Other laplace functions: stehfest_coeffs

Examples

test <- function(p,...){1/p}
coeffs8 <- stehfest_coeffs(8)
td <- logseq(-2,4,30)
stehfest_inversion(td, coeffs8,test)
#>  [1] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1