R/numHess.R

Defines functions numHess

numHess <- function(fcn, x, ...) {
  f1 <- fcn
  n <- length(x)
  h <- matrix(0, n, n)
  f2 <- function(z, ...) { numgrad(fcn=f1, z, ...)$g}
  h <- numgrad(fcn=f2, x=x, ...)$g
  return(h)
}
    
thiyangt/fformpp documentation built on Jan. 5, 2024, 5:44 a.m.