R/second.order.num.deriv.R

In rodd: Optimal Discriminating Designs

Defines functions second.order.num.deriv

second.order.num.deriv <- function(f, x, epsilon, ...)
{
	(f(x + epsilon, ...) - 2 * f(x, ...) + f(x - epsilon, ...)) / (epsilon * epsilon)
}

Try the rodd package in your browser