R/fun.fmklb.R

In GLDEX: Fitting Single and Mixture of Generalised Lambda Distributions

"fun.fmklb" <-
function(L3, L4, k)
{
j <- 0:k
result <- sum(((-1)^j * (L3^(k - j) * L4^(j))^-1 * (gamma(k + 1)/(gamma(k - j + 1) * gamma(j + 1))) * fun.beta(L3 * (k - j) + 1, L4 *
j + 1)))
}