Stirling's factorial approximation

diff --git a/extra/math/analysis/analysis.factor b/extra/math/analysis/analysis.factor
index bde5cad7a019015d2bc2fde67ef8ec29d5838c5d..2973d3d3f2fb6572cebdb084d0af724e79b3927b 100644 (file)
--- a/extra/math/analysis/analysis.factor
+++ b/extra/math/analysis/analysis.factor
@@ -108,3 +108,12 @@ PRIVATE>
         swap -1.0 * exp
         *
     ] if ;
+
+! James Stirling's approximation for N!:
+! http://www.csse.monash.edu.au/~lloyd/tildeAlgDS/Numerical/Stirling/
+
+: stirling-fact ( n -- fact )
+  [ pi 2 * * sqrt ]
+  [ dup e / swap ^ ]
+  [ 12 * recip 1 + ]
+  tri * * ;