`gammaln`

Logaritmo natural da função gama.

Categoria: Mathematics > Elementary Math

Sintaxe

Y = gammaln(X)

Retorna log(gamma(X)) para cada elemento, calculado de modo a não estourar.

Descrição

gamma(172) já estoura, mas seu logaritmo é de apenas algumas centenas. gammaln trabalha com o logaritmo o tempo todo, e por isso lida com argumentos de qualquer tamanho prático. É assim que se calculam coeficientes binomiais grandes e razões de fatoriais sem estouro: exp(gammaln(n+1) - gammaln(k+1) - gammaln(n-k+1)).

Argumentos de entrada

X: Array real.

Argumentos de saída

Y: Array do mesmo tamanho de X.

Exemplo

Além do alcance de gamma

gammaln([10 100 1000])
exp(gammaln(51) - gammaln(26) - gammaln(26))

Saída

ans = 1×3

   1.0e+03 *

    0.0128    0.3591    5.9052
ans = 1.2641e+14

Veja também

gamma factorial nchoosek psi