Skip to main content

New upper bounds of n!

Abstract

In this article, we deduce a new family of upper bounds of n! of the form

n ! < 2 π n ( n / e ) n e M n [ m ] n , M n [ m ] = 1 2 m + 3 1 4 n + k = 1 m 2 m - 2 k + 2 2 k + 1 2 - 2 k ζ ( 2 k , n + 1 / 2 ) m = 1 , 2 , 3 , . . . .

We also proved that the approximation formula 2 π n ( n / e ) n e M n [ m ] for big factorials has a speed of convergence equal to n-2m- 3for m = 1,2,3,..., which give us a superiority over other known formulas by a suitable choice of m.

Mathematics Subject Classification (2000): 41A60; 41A25; 57Q55; 33B15; 26D07.

1 Introduction

Stirling' formula

n ! ~ 2 n π ( n / e ) n
(1)

is one of the most widely known and used in asymptotics. In other words, we have

lim n n ! e n 2 π n n n = 1 .
(2)

This formula provides an extremely accurate approximation of n! for large values of n. The first proofs of Stirling's formula was given by De Moivre (1730) [1] and Stirling (1730) [2]. Both used what is now called the Euler-MacLaurin formula to approximate log 2 + log 3 + ... + log n. The first derivation of De Moivre did not explicity determine the constant 2 π . In 1731, Stirling determine this constant using Wallis' formula

lim n 2 2 n ( n ! ) 2 ( 2 n ) ! 1 n = π .

Over the years, there have been many different upper and lower bounds for n! by various authors [310]. Artin [11] show that μ ( n ) =ln n ! e n n n 2 π n lies between any two successive partial sums of the Stirling's series

B 2 1 . 2 . n + B 4 3 . 4 . n 3 + B 6 5 . 6 . n 5 + ,
(3)

where the numbers B i are called the Bernoulli numbers and are defined by

B 0 = 1 , k = 0 n 1 ( k n ) B k = 0 , n 2.
(4)

We can't take infinite sum of the series (3) because the series diverges. Also, Impens [12] deduce Artin result with different proof and show that the Bernoulli numbers in this series cannot be improved by any method whatsoever.

The organization of this article is as follows. In Section 2, we deduce a general double inequality of n!, which already obtained in [9] with different proof. Section 3 is devoted to getting a new family of upper bounds of n! different from the partial sums of the Stirling's series. In Section 4, we measure the speed of convergence of our approximation formula 2 π n ( n / e ) n e M n [ m ] for big factorials. Also, we offer some numerical computations to prove the superiority of our formula over other known formulas.

2 A double inequality of n!

In view of the relation (2), we begin with the two sequences K n and f n defined by

f n = n ! e n 2 π n n n e - K n n 1 ,
(5)

where

lim n K n = 0 .
(6)

Then we have

lim n f n = 1 .
(7)

Now define the sequence g n by

g n = f n + 1 f n = e K n - K n + 1 + 1 1 + 1 n - n - 1 / 2 ,
(8)

which satisfies

lim n g n = 1
(9)

and

g n = g n d d n K n - K n + 1 - ( n + 1 / 2 ) ln 1 + 1 n .
(10)

There are two cases. The first case if K n = a n such that

d d n a n - a n + 1 - ( n + 1 / 2 ) ln 1 + 1 n > 0
(11)

then we get g n > 0 and hence g n is strictly increasing function. But g n → 1 as n → ∞, then g n < 1. Hence f n + 1 f n < 1 which give us that fn+1< f n . Then f n is strictly decreasing function. Also, f n → 1 as n → ∞, then we obtain f n > 1. Then

n ! e n 2 π n n n e - a n > 1 n 1
(12)

or

2 π n ( n / e ) n e a n < n ! n 1 .
(13)

The condition (11) means that the function a n - an+1-(n+ 1/2) ln 1 + 1 n is strictly increasing function also it tends to -1 as n → ∞. Then

a n - a n + 1 - ( n + 1 / 2 ) ln 1 + 1 n < - 1

or

a n - a n + 1 < ( n + 1 / 2 ) ln 1 + 1 n - 1 .
(14)

The second case if K n = b n such that

d d n b n - b n + 1 - ( n + 1 / 2 ) ln 1 + 1 n < 0 .
(15)

Similarly, we can prove that f n < 1. Then

n ! e n 2 π n n n e - b n < 1 n 1
(16)

or

n ! < 2 π n ( n / e ) n e b n n 1 .
(17)

The condition (15) means that the function b n - bn+1- (n + 1/2) ln 1 + 1 n is strictly decreasing function also it tends to -1 as n → ∞. Then

b n - b n + 1 - ( n + 1 / 2 ) ln 1 + 1 n > - 1

or

( n + 1 / 2 ) ln 1 + 1 n - 1 < b n - b n + 1 .
(18)

From the well-known expansion

1 2 ln 1 + y 1 - y = k = 1 y 2 k - 1 2 k - 1 y < 1

in which we substitute y= 1 2 n + 1 , so 1 + y 1 - y =1+ 1 n . Then

( n + 1 / 2 ) ln 1 + 1 n - 1 = k = 1 1 2 k + 1 1 ( 2 n + 1 ) 2 k .
(19)

Now we obtain the following result

Theorem 1.

2 π n ( n / e ) n e a n < n ! < 2 π n ( n / e ) n e b n n 1 ,
(20)

where the two sequences a n , b n → 0 as n → ∞ and satisfy

a n - a n + 1 < k = 1 1 2 k + 1 1 ( 2 n + 1 ) 2 k < b n - b n + 1 .
(21)

A q-analog of the inequality (20) was introduced in [13].

3 A new family of upper bounds of n!

By manipulating the series k = 1 1 2 k + 1 1 ( 2 n + 1 ) 2 k to find upper bounds, we get

k = 1 1 2 k + 1 1 ( 2 n + 1 ) 2 k < k = 1 m 1 ( 2 k + 1 ) ( 2 n + 1 ) 2 k + 1 2 m + 3 k = m + 1 1 ( 2 n + 1 ) 2 k < k = 1 1 ( 2 k + 1 ) ( 2 n + 1 ) 2 k + 1 4 ( 2 m + 3 ) n ( n + 1 ) ( 2 n + 1 ) 2 m , m = 1 , 2 , 3 , . . . .

Let's solve the following recurrence relation w.r.t n

M n [ m ] - M n + 1 [ m ] = k = 1 m 1 ( 2 k + 1 ) ( 2 n + 1 ) 2 k + 1 4 ( 2 m + 3 ) n ( n + 1 ) ( 2 n + 1 ) 2 m , m = 1 , 2 , 3 , . . .
(22)

which give us

M n [ m ] = M 0 [ m ] - i = 1 n - 1 k = 1 m 1 ( 2 k + 1 ) ( 2 i + 1 ) 2 k + 1 4 ( 2 m + 3 ) i ( i + 1 ) ( 2 i + 1 ) 2 m = M 0 [ m ] - k = 1 m 1 2 k + 1 i = 1 n - 1 1 ( 2 i + 1 ) 2 k - 1 4 ( 2 m + 3 ) i = 1 n - 1 1 i ( i + 1 ) ( 2 i + 1 ) 2 m .

But

i = 1 n - 1 1 i ( i + 1 ) ( 2 i + 1 ) 2 m = i = 1 n - 1 1 i - 1 i + 1 - 4 ( 2 i + 1 ) 2 - 4 ( 2 i + 1 ) 4 - - 4 ( 2 i + 1 ) 2 m = n - 1 n - 4 k = 1 m i = 1 n - 1 1 ( 2 i + 1 ) 2 k .

Then

M n [ m ] = M 0 [ m ] - n - 1 4 n ( 2 m + 3 ) - k = 1 m 1 2 k + 1 - 1 2 m + 3 i = 1 n - 1 1 ( 2 i + 1 ) 2 k .

The series

i = 1 1 ( 2 i + 1 ) 2 k = ζ ( 2 k ) ( 1 - 2 - 2 k ) - 1 ,

where ζ(x) is the Riemann Zeta function. By using the relation [14]

ζ ( 2 k ) = ( - 1 ) k - 1 2 2 k - 1 ( 2 k ) ! B 2 k π 2 k ,

where B r s are Bernoulli's numbers. Then

i = 1 1 ( 2 i + 1 ) 2 k = ( - 1 ) k - 1 ( 2 2 k - 1 ) 2 ( 2 k ) ! B 2 k π 2 k - 1 .

Hence, we can choose

M 0 [ m ] = 1 4 ( 2 m + 3 ) + k = 1 m 1 2 k + 1 - 1 2 m + 3 ( - 1 ) k - 1 ( 2 2 k - 1 ) 2 ( 2 k ) ! B 2 k π 2 k - 1 ,
(23)

which satisfies

lim n M n [ m ] = 0 , m = 1 , 2 , 3 , . . . .

Then

M n [ m ] = 1 2 m + 3 1 4 n + k = 1 m 2 m - 2 k + 2 2 k + 1 ( - 1 ) k - 1 ( 2 2 k - 1 ) 2 ( 2 k ) ! B 2 k π 2 k - 1 - i = 1 n - 1 1 ( 2 i + 1 ) 2 k .
(24)

By using the relation

i = 1 n - 1 1 ( 2 i + 1 ) 2 k = - 1 - ( 2 - 2 k - 1 ) ζ ( 2 k ) - 2 - 2 k ζ ( 2 k , n + 1 / 2 ) = - 1 - ( - 1 ) k - 1 ( 1 - 2 2 k ) 2 ( 2 k ) ! B 2 k π 2 k - 2 - 2 k ζ ( 2 k , n + 1 / 2 ) ,

we get

M n [ m ] = 1 2 m + 3 1 4 n + k = 1 m 2 m - 2 k + 2 2 k + 1 2 - 2 k ζ ( 2 k , n + 1 / 2 ) .
(25)

In the following Lemma, we will see that the upper bound M n m will improved with increasing the value of m.

Lemma 3.1.

M n [ m + 1 ] < M n [ m ] , m , n = 1 , 2 , 3 , . . . .
(26)

Proof. From (25), we get

M n [ m + 1 ] = 1 2 m + 5 1 4 n + k = 1 m 2 m - 2 k + 2 2 k + 1 2 - 2 k ζ ( 2 k , n + 1 / 2 ) + 1 2 m + 5 k = 1 m + 1 1 2 k + 1 2 1 - 2 k ζ ( 2 k , n + 1 / 2 ) .

Then

M n [ m + 1 ] - M n [ m ] = 2 ( 2 m + 3 ) ( 2 m + 5 ) - 1 4 n + k = 1 m + 1 2 - 2 k ζ ( 2 k , n + 1 / 2 ) < 2 ( 2 m + 3 ) ( 2 m + 5 ) - 1 4 n + k = 1 2 - 2 k ζ ( 2 k , n + 1 / 2 ) .

Now

k = 1 2 - 2 k ζ ( 2 k , n + 1 / 2 ) = k = 1 r = 0 2 - 2 k ( n + r + 1 / 2 ) 2 k = r = 0 k = 1 2 - 2 k ( n + r + 1 / 2 ) 2 k = r = 0 1 1 - 1 ( 2 n + 2 r + 1 ) 2 - 1 = r = 0 1 4 ( n + r ) ( n + r + 1 ) = 1 4 r = 0 1 n + r - 1 n + r + 1 = 1 4 n .

Then

k = 1 2 - 2 k ζ ( 2 k , n + 1 / 2 ) = 1 4 n .
(27)

and hence

M n [ m + 1 ] - M n [ m ] < 2 ( 2 m + 3 ) ( 2 m + 5 ) - 1 4 n + k = 1 2 - 2 k ζ ( 2 k , n + 1 / 2 ) < 2 ( 2 m + 3 ) ( 2 m + 5 ) - 1 4 n + 1 4 n < 0 .

The following Lemma show that M n [ 1 ] is an improvement of the upper bound 1 12 n .

Lemma 3.2.

M n [ 1 ] < 1 12 n
(28)

Proof. From Eq. (25), we have M n [ 1 ] = 1 20 n + 1 30 ζ ( 2 , n + 1 / 2 ) . By using (27), we obtain

1 n = k = 1 2 2 - 2 k ζ ( 2 k , n + 1 / 2 ) = ζ ( 2 , n + 1 / 2 ) + k = 2 2 2 - 2 k ζ ( 2 k , n + 1 / 2 ) .

Then

ζ ( 2 , n + 1 / 2 ) < 1 n

which give us that

1 30 ζ ( 2 , n + 1 / 2 ) + 1 20 n < 1 30 n + 1 20 n = 1 12 n .

4 The speed of convergence of the approximation formula 2 π n ( n / e ) n e M n [ m ] for big factorials

In what follows, we need the following result, which represents a powerful tool to measure the rate of convergence.

Lemma 4.1. If (w n )n≥1is convergent to zero and there exists the limit

lim n n k ( w n - w n + 1 ) = l
(29)

with k > 1, then there exists the limit:

lim n n k - 1 w n = l k - 1

This Lemma was first used by Mortici for constructing asymptotic expansions, or to accelerate some convergences [1521]. By using Lemma (4.1), clearly the sequence (w n )n≥1converges more quickly when the value of k satisfying (29) is larger.

To measure the accuracy of approximation formula 2 π n ( n / e ) n e M n [ m ] , define the sequence (w n )n≥1by the relation

n ! = 2 π n ( n / e ) n e M n [ m ] e w n ; n = 1 , 2 , 3 , . . . .
(30)

This approximation formula will be better as (w n )n≥1converges faster to zero. Using the relation (30), we get

w n = ln n ! - ln 2 π - ( n + 1 / 2 ) ln n + n - M n [ m ]

Then

w n - w n + 1 = ( n + 1 / 2 ) ln ( 1 + 1 / n ) - 1 + M n + 1 [ m ] - M n [ m ]

By using the relations (19) and (22), we have

w n - w n + 1 = k = m + 1 1 2 k + 1 1 ( 2 n + 1 ) 2 k - 1 4 ( 2 m + 3 ) n ( n + 1 ) ( 2 n + 1 ) 2 m = k = m + 3 1 2 k + 1 1 ( 2 n + 1 ) 2 k - ( 5 + 2 m + 8 n + 8 n 2 ) ( 2 n + 1 ) - 4 - 2 m 4 n ( n + 1 ) ( 15 + 16 m + 4 m 2 )

Then

lim n n 2 ( m + 2 ) ( w n - w n + 1 ) = - 1 ( 2 m + 3 ) ( 2 m + 5 ) 2 2 m + 3 ; n , m = 1 , 2 , 3 , . . . .
(31)

Theorem 2. The rate of convergence of the sequence w n is equal to n-2m-3, since

lim n n 2 m + 3 w n = 1 ( 2 m + 3 ) 2 ( 2 m + 5 ) 2 2 m + 3 .

In 2011, Mortici [22] shows by numerical computations that his formula

n ! ~ 2 π n ( n / e ) n e 1 12 n + 2 5 n = μ n
(32)

is much stronger than other known formulas such as:1

Table 1

The following table shows numerically that our new formula λ n , 1 = 2 π n ( n / e ) n e M n [ 1 ] has a superiority over the the Mortici's formula μ n .

Table 2

References

  1. De Moivre A: Miscellanea Analytica de Seriebus et Quadraturis. London; 1730.

    Google Scholar 

  2. Stirling J: Methodus Differentialis. London; 1730.

    Google Scholar 

  3. Hummel PM: A note on Stirling's formoula. Am Math Monthly 1940, 42(2):97–99.

    Article  MathSciNet  Google Scholar 

  4. Cesàro E: Elements Lehrbuch der algebraischen analysis und der In-finitesimalrechnung. Leipzig, Teubner; 1922.

    Google Scholar 

  5. Uspensky JV: Introduction to Mathematical Probability. McGraw Hill, New York; 1937.

    Google Scholar 

  6. Robbins H: A remark on Stirling's formula. Am Math Monthly 1955, 62: 26–29. 10.2307/2308012

    Article  Google Scholar 

  7. Nanjundiah TS: Note on Stirling's formula. Am Math Monthly 1965, 66: 701–703.

    Article  MathSciNet  Google Scholar 

  8. Maria AJ: A remark on Striling's formula. Am Math Monthly 1965, 72: 1096–1098. 10.2307/2315957

    Article  MathSciNet  Google Scholar 

  9. Beesack PR: Improvement of Stirling's formula by elementary methods. Univ. Beograd, Publ. Elektrotenhn. Fak Ser Mat Fiz 1969, 17–21. No. 274–301

    Google Scholar 

  10. Michel R: On Stirling's formula. Am. Math. Monthly. 2002, 109(4):388–390.

    Google Scholar 

  11. Artin E: The Gamma Function (translated by M Butler). Holt, Rinehart and Winston. New York; 1964.

    Google Scholar 

  12. Impens C: Stirling's series made easy. Am Math Monthly 2003, 110(8):730–735. 10.2307/3647856

    Article  MathSciNet  Google Scholar 

  13. Mansour M, Obaid MA: Bounds of q -factorial [ n ] q ! Ars Combinatoria CII 2011, 313–319.

    Google Scholar 

  14. Andrews GE, Askey R, Roy R: Special Functions. Cambridge University Press, Cambridge; 1999.

    Chapter  Google Scholar 

  15. Mortici C: An ultimate extremely accurate formula for approximation of the factorial function. Archiv der Mathematik (Basel) 2009, 93(1):37–45. 10.1007/s00013-009-0008-5

    Article  MathSciNet  Google Scholar 

  16. Mortici C: Product approximations via asymptotic integration. Am Math Monthly 2010, 117: 434–441. 10.4169/000298910X485950

    Article  MathSciNet  Google Scholar 

  17. Mortici C: New approximations of the gamma function in terms of the digamma function. Appl Math Lett 2010, 23: 97–100. 10.1016/j.aml.2009.08.012

    Article  MathSciNet  Google Scholar 

  18. Mortici C: New improvements of the Stirling formula. Appl Math Comput 2010, 217(2):699–704. 10.1016/j.amc.2010.06.007

    Article  MathSciNet  Google Scholar 

  19. Mortici C: Best estimates of the generalized Stirling formula. Appl Math Comput 2010, 215(11):4044–4048. 10.1016/j.amc.2009.12.013

    Article  MathSciNet  Google Scholar 

  20. Mortici C: A class of integral approximations for the factorial function. Comput Math Appl 2010, 59(6):2053–2058. 10.1016/j.camwa.2009.12.010

    Article  MathSciNet  Google Scholar 

  21. Mortici C: Ramanujan formula for the generalized Stirling approximation. Appl Math Comput 2010, 217(6):2579–2585. 10.1016/j.amc.2010.07.070

    Article  MathSciNet  Google Scholar 

  22. Mortici C: A new Stirling series as continued fraction. Numer Algorithm 2011, 56(1):17–26. 10.1007/s11075-010-9370-4

    Article  MathSciNet  Google Scholar 

  23. Burnside W: A rapidly convergent series for Log N! Messenger Math 1917, 46: 157–159.

    Google Scholar 

  24. Batir N: Sharp inequalities for factorial n . Proyecciones 2008, 27(1):97–102.

    Article  MathSciNet  Google Scholar 

  25. Gosper RW: Decision procedure for indefinite hypergeometric summation. Proc Natl Acad Sci USA 1978, 75: 40–42. 10.1073/pnas.75.1.40

    Article  MathSciNet  Google Scholar 

  26. Ramanujan S: The Lost Notebook and Other Unpublished Papers (Introduced by Andrews, GE). Narosa Publishing House, New Delhi; 1988.

    Google Scholar 

Download references

Acknowledgements

This project was founded by the Deanship of Scientific Research (DSR), King Abdulaziz University, Jeddah, under grant No. 448/130/1431. The authors, therefore, acknowledge with thanks DSR support for Scientific Research.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Mansour Mahmoud.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

All authors carried out the proofs. All authors conceived of the study and participated in its design and coordination. All authors read and approved the final manuscript.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License (https://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Mahmoud, M., Alghamdi, M.A. & Agarwal, R.P. New upper bounds of n!. J Inequal Appl 2012, 27 (2012). https://doi.org/10.1186/1029-242X-2012-27

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1029-242X-2012-27

Keywords