famfun.R
famfun <- function(fam, deriv, theta) {
setfam(list(fam))
result <- .C(C_aster_family,
index = as.integer(1),
deriv = as.integer(deriv),
theta = as.double(theta),
value = double(1))$value
clearfam()
return(result)
}