Horrible limit with factorials. Need to use Stirling formula?

Click For Summary
SUMMARY

The discussion centers on the application of Stirling's formula for approximating factorials in limit calculations. Users debate the necessity of using Stirling's formula, with one participant asserting its utility in probability theory and asymptotic analysis. The formula, St(n) = √(2πn) n^n e^(-n), is confirmed to be a reliable approximation for large n, with bounds established for n!. The conversation highlights the importance of understanding Stirling's formula for accurate factorial approximations in mathematical problems.

PREREQUISITES
  • Understanding of factorials and their properties
  • Familiarity with Stirling's approximation and its derivation
  • Basic knowledge of asymptotic analysis
  • Experience with limit calculations in calculus
NEXT STEPS
  • Study the derivation of Stirling's formula from integral representations
  • Learn about asymptotic analysis techniques in mathematical proofs
  • Explore the applications of Stirling's formula in probability theory
  • Investigate numerical methods for approximating factorials
USEFUL FOR

Mathematicians, students studying calculus and probability, and anyone interested in advanced techniques for approximating factorials in mathematical analysis.

tsuwal
Messages
105
Reaction score
0

Homework Statement



2013-03-16 17.17.22.jpg



Homework Equations





The Attempt at a Solution



Do I need to used the boring Stirling formula?
 
Physics news on Phys.org
tsuwal said:

Homework Statement



View attachment 56798


Homework Equations





The Attempt at a Solution



Do I need to used the boring Stirling formula?

What is stopping you from trying the use of Stirling's formula to see if it works?

BTW: Stirling is far from boring; it is extremely useful and plays a crucial role in many areas such as probability theory and practice.
 
Actually, i think this problem can be done without stirrings. For the top factorial, you have (2n-2) terms, giving a leading term of (2n)^(2n-2) for the expansion. Then the bottom will have (n^(n-2)+...)(n^n+...), giving a leading coefficient n^(2n-2)... This your limit becomes (2n-2)!/(2^(2n-2)*n!*n-2!)~O(2^(2n-2)*n^(2n-2)/(2^(2n-2)*n^(2n-2))~1
 
tt2348 said:
Actually, i think this problem can be done without stirrings. For the top factorial, you have (2n-2) terms, giving a leading term of (2n)^(2n-2) for the expansion. Then the bottom will have (n^(n-2)+...)(n^n+...), giving a leading coefficient n^(2n-2)... This your limit becomes (2n-2)!/(2^(2n-2)*n!*n-2!)~O(2^(2n-2)*n^(2n-2)/(2^(2n-2)*n^(2n-2))~1

Too bad this answer is wrong.
 
I don't use the stirling formula because it is just a aproximation and I'm not sure if I can substitute n! by n^n/e^n
In fact, wolfram alpha says the limit of n!*e^n/n^n as n goes to infinity is not 1, is infinity, so it is not a good aproximation, right?
 
tsuwal said:
I don't use the stirling formula because it is just a aproximation and I'm not sure if I can substitute n! by n^n/e^n
In fact, wolfram alpha says the limit of n!*e^n/n^n as n goes to infinity is not 1, is infinity, so it is not a good aproximation, right?

Your approximation is not good, but Stirling's formula IS good. The Stirling formula is
St(n) = \sqrt{2 \pi n}\: n^{n} e^{-n}. This is an asymptotic formula for n!, in the sense that
\frac{n!}{St(n)} \to 1 \; \text{ as } \ \to \infty.
St(n) is a useably good approximation to n! for large, finite n, in the sense that the % errors are small.

In fact, St(n) is a lower bound on n! for any positive integer n. An upper bound St1(n) can be obtained simply as
St1(n) = St(n) e^{1/(12n)}.
We have ## St(n) < n! < St1(n)## for all integers ##n \geq 1##. The fancier form St1 improves the accuracy of Stirling remarkably, even for small n. For example, here are some numerical values:
\begin{array}{cccc}<br /> n &amp; St(n) &amp; n\! &amp; St1(n) \\ <br /> 1 &amp; 0.922137 &amp; 1 &amp; 1.00227 \\<br /> 2 &amp; 1.919 &amp; 2 &amp; 2.00065 \\<br /> 5 &amp; 118.019 &amp; 120 &amp; 120.003 \\<br /> 10 &amp; 3.59870e+06 &amp; 3.62880e+06 &amp; 3.62881e+06 \\<br /> 15 &amp; 1.30043e+12 &amp; 1.30767e+12 &amp; 1.30768e+12 \\<br /> 20 &amp; 2.42279e+18 &amp; 2.43290e+18 &amp; 2.43290e+18 \\<br /> 30 &amp; 2.64517e+32 &amp; 2.65253e+32 &amp; 2.65253e+32 \\<br /> 40 &amp; 8.14217e+47 &amp; 8.15915e+47 &amp; 8.15915e+47 \\<br /> 50 &amp; 3.03634e+64 &amp; 3.04141e+64 &amp; 3.04141e+64 \end{array}<br />
Anyhow, for limit results, ordinary Stirling is good enough.
 
Last edited:
Now it makes sense. I didn't knew "that" stirling formula, mine is deduced from the numerical integration of the logarithmic function by ln(n!). Would like to know how do you get to your formula though.. Thanks!
 
tsuwal said:
Now it makes sense. I didn't knew "that" stirling formula, mine is deduced from the numerical integration of the logarithmic function by ln(n!). Would like to know how do you get to your formula though.. Thanks!

One way is to do an asymptotic analysis of the integral representation
n\! = \int_0^{\infty} e^{-x} x^n \, dx,
as is done in http://mathworld.wolfram.com/StirlingsApproximation.html . Another way (as done in Feller, Introduction to Probability Theory and its Applications, Volume I) is to look directly at the series for ln(n!). Feller gives a nice, elementary argument that shows
C + \left( n + \frac{1}{2}\right) \ln(n) - n &lt; \ln (n!)<br /> &lt; C + \left( n + \frac{1}{2}\right) \ln(n) - n +\frac{1}{12n}.
The argument leaves the constant C unevaluated, but other methods (such as the above integral method, or the limiting theorems of Probability theory itself) serve to give
C = \frac{1}{2} \ln(2 \pi). Noter that Feller's method gives the lower and upper bounds I mentioned.
 

Similar threads

  • · Replies 2 ·
Replies
2
Views
1K
Replies
10
Views
4K
  • · Replies 9 ·
Replies
9
Views
8K
  • · Replies 1 ·
Replies
1
Views
3K
Replies
3
Views
1K
  • · Replies 2 ·
Replies
2
Views
2K
Replies
1
Views
1K
Replies
6
Views
3K
  • · Replies 3 ·
Replies
3
Views
2K
Replies
6
Views
4K