This routine approximates factorials and their logarithms by Stirling's approximation.
The calls are
CALL STIRLING(FACT,I,IER)
and
CALL STIRLOG(FACTL,I,IER)
where
FACT = the computed factorial
FACTL = computed logarithm of FACT
I = the number for which the factorial is computed
IER = error return; =0 normal return, =−1 I was negative
Return to