i2p/haiku-beta6
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Use 5 factor Stirling's Series

Browse Source
This commit is contained in:
John Scipione
2016-03-03 10:07:25 -08:00
parent fb4dcb4965
commit 261bdbab6d
1 changed files with 6 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files
src/kits/shared/ExpressionParser.cpp
+6 -4
View File
@@ -712,12 +712,14 @@ ExpressionParser::_ParseFactorial(MAPM value)
if (value < 1000) if (value < 1000)
return value.factorial(); return value.factorial();
else { else {
// Use Stirling's approximation (with extra term) // Use Stirling's approximation (5 term expansion)
// http://hyperphysics.phy-astr.gsu.edu/hbase/math/stirling.html // http://en.wikipedia.org/wiki/Stirling%27s_approximation
// n! ≈ (n/e)^n * sqrt(2πn) * (1 + (1/12n))
return value.pow(value) / value.exp() return value.pow(value) / value.exp()
* (MAPM(2) * MAPM(MM_PI) * value).sqrt() * (MAPM(2) * MAPM(MM_PI) * value).sqrt()
* (MAPM(1) + (MAPM(1) / (MAPM(12) * value))); * (MAPM(1) + (MAPM(1) / (MAPM(12) * value))
+ (MAPM(1) / (MAPM(288) * value.pow(2)))
- (MAPM(139) / (MAPM(51840) * value.pow(3)))
- (MAPM(571) / (MAPM(2488320) * value.pow(4))));
} }
} }