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 negativeReturn to