Para x e y pequenos, o algoritmo usa a expressão em função da função gama, de outro modo, ele aplica a função exponencial no resutado da função betaln function fornecido no DCDFLIB: Biblioteca de Rotinas FORTRAN para Funções, Inversas e Outros Parâmetros de Distribuição Cumulativa (ver cdfbet para maiores informações sobre DCDFLIB).