R/getJac.R

"getJac" <- function(k, x) {
  jac <- array(0, dim = c(length(x), length(k), length(k)))

  for (i in 1:length(k)) {
    jac[, i, i] <- -x * exp(-k[i] * x)
  }

  jac
}
glotaran/TIMP documentation built on Feb. 11, 2023, 4:11 a.m.