Next Article in Journal
Canonical Construction of Invariant Differential Operators: A Review
Next Article in Special Issue
Functional and Operatorial Equations Defined Implicitly and Moment Problems
Previous Article in Journal
Revisiting de Broglie’s Double-Solution Pilot-Wave Theory with a Lorentz-Covariant Lagrangian Framework
Previous Article in Special Issue
The Properties of Topological Manifolds of Simplicial Polynomials
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

New Accurate Approximation Formula for Gamma Function

1
Department of Mathematics, Faculty of Science, King Abdulaziz University, P.O. Box 80203, Jeddah 21589, Saudi Arabia
2
Department of Mathematics and Statistics, Faculty of Science, University of Jeddah, P.O. Box 80327, Jeddah 21589, Saudi Arabia
*
Author to whom correspondence should be addressed.
Symmetry 2024, 16(2), 150; https://doi.org/10.3390/sym16020150
Submission received: 26 December 2023 / Revised: 8 January 2024 / Accepted: 26 January 2024 / Published: 27 January 2024

Abstract

:
In this paper, a new approximation formula for the gamma function and some of its symmetric inequalities are established. We prove the superiority of our results over Yang and Tian’s approximation formula for the gamma function of order v 9 .

1. Introduction

The gamma function extends the idea of a factorial to non-integer numbers. C. F. Gauss and other mathematicians investigated it further after L. Euler introduced it for the first time in the 18th century. The gamma function is defined as follows:
Γ ( v ) = 0 e t t v 1 d t , v > 0
or by
Γ ( v ) = lim m m ! m v v ( 1 + v ) ( 2 + v ) ( m + v ) , v R { 0 , 1 , 2 , }
and satisfies the recurrence relation Γ ( v + 1 ) = v Γ ( v ) , and hence, m ! = Γ ( m + 1 ) for m N . More work has been performed by mathematicians to obtain accurate estimates of m ! and the gamma function. J. Stirling developed the following important formula:
Γ ( v + 1 ) = 2 π v v / e v 1 + O ( v 1 ) : = χ 1 ( v ) , v .
As v grows, this estimate becomes more and more accurate. For high values, Stirling’s formula, which comes from Stirling’s series expansion, provides a practical alternative to numerically integrating the defining integral for estimating the gamma function. It is crucial to remember that, while Stirling’s approximation is frequently correct for large values of v, it might not be appropriate for small values of v or situations requiring great accuracy. Other techniques can be used to obtain more precise and comprehensive approximations.
Some sharp inequalities for the ratio of gamma functions are presented by Cao and Wang [1] by using the multiple-correction method. Alzer and Jameson [2] present a new characterization of Euler’s constant γ and a concavity property of the Psi function. Yang and Tian [3] refine Windschitl’s gamma function approximation formula by providing two asymptotic expansions based on a little-known power series. The authors of [4] use several classical inequalities, such as Chebychev’s inequality for synchronous mappings, to propose some inequalities involving the extended gamma function and the Kummer confluent hypergeometric k-function. Qi and Guo [5] study the properties of the Bernstein function of the newly created ratio of finitely many gamma functions as well as the history, backgrounds, extensions, and applications of a series of ratios for finitely many gamma functions. For large values of the involved parameters, Reynolds and Stauffer [6] investigate the improved infinite sum for the incomplete gamma function. Tian and Yang [7] expand and generalize some previous results by presenting the necessary and sufficient conditions for a ratio involving q-gamma functions to be logarithmically completely monotonic using a new method. Zhang, Yin, and You [8] deduce some new inequalities and completely monotonic properties involving the generalized functions k-gamma and k-polygamma. Yildirim [9] uses the Bernstein–Widder theorem and some properties of the k-special function to present k-generalizations of some classical results and improvements to some bounds of recent results about the k-polygamma functions. Based on the incomplete gamma function, Castillo, Rojas, and Reyes [10] deduce a more flexible extension for the Fréchet distribution, along with applications. Mahmoud, Alsulami, and Almarashi [10] examine the monotonicity of some functions involving Γ ( v ) and ascertain its bounds that they were able to derive are more precise than certain inequalities that have been published previously. The research references [1,2,3,4,5,6,7,8,9,10,11] present further information regarding the gamma function and the formulas, inequalities, approximations, generalizations, and applications that go along with it.
Simple and accurate gamma function approximation formulas are essential for many applications because gamma function integrals cannot be computed directly for most non-integer values. By using approximation formulas, one may evaluate Γ ( v ) more quickly and with less processing effort. For instance, simple approximations for the gamma function can greatly increase the accuracy and efficiency of algorithms in numerical techniques and calculations. This is essential for areas like statistical analysis, optimization, and differential equation solutions [12,13]. Some intriguing estimates of gamma functions are as follows:
  • Ramanujan presents the approximation formula [14] as v
Γ ( v + 1 ) = π v e v 8 v 3 + 4 v 2 + v + 1 / 30 6 1 + O ( v 4 ) : = χ 2 ( v ) ,
based on some numerical calculations as conjecture.
Windschitl finds the approximation formula [15] as v
Γ ( v + 1 ) = 2 π v v e v v sinh 1 v v / 2 1 + O ( v 5 ) : = χ 3 ( v ) ,
after he noticed by coincidence the relation between some power series expansions of the extended Stirling’s formula and the hyperbolic sine function. Since it is accurate to more than eight decimal places for v > 8 , he suggested using the approximation 2 π v v / e v 1 810 v 6 + v sinh 1 v v / 2 to compute the values of the gamma function on calculators with limited program or register memory.
Smith presents the approximation formula [16] as v
Γ ( v + 1 / 2 ) = 2 π v e v 2 v tanh 1 2 v v / 2 1 + O ( v 5 ) : = χ 4 ( v ) ,
and some new representations of the gamma function and some of its related functions. In addition, he provides new continued fractions and a new formula for the beta function.
  • Mahmoud and Almuashi deduce the approximation formula [17] as v
Γ ( v + 1 ) = 2 π v v e v v 2 + 7 60 v 2 1 20 v / 2 1 + O ( v 5 ) : = χ 5 ( v ) .
Moreover, they derive some new bounds of Γ ( v ) more accurately than some of its recent ones by using Padé approximants and a new asymptotic expansion of gamma that they derived.
Nemes uses a series transformation to convert the Stirling asymptotic series approximation of the gamma function into a new one with better convergence properties and deduces the approximation formula [18] as v
Γ ( v + 1 ) = 2 π v v / e v 3 40 v 2 + 3 120 v 2 1 v 1 + O v 5 : = χ 6 ( v ) .
Yang and Chu’s deduce the approximation formula [19] as v
Γ ( v + 1 / 2 ) = 2 π v / e v e 5 v 120 v 2 + 7 1 + O v 5 : = χ 7 ( v ) ,
and some upper and lower bounds of the factorial m ! and the gamma function are presented as applications.
Nemes’s approximation formula [18] as v
Γ ( v + 1 ) = 2 π v v e v e 210 v 2 + 53 360 v ( 7 v 2 + 2 ) 1 + O ( v 7 ) : = χ 8 ( v ) .
Yang and Chu’s approximation formula [19] as v
Γ ( v + 1 / 2 ) = 2 π v / e v e 5880 v 2 1517 1440 v 98 v 2 + 31 1 + O v 7 : = χ 9 ( v ) .
Based on Windschitl’s formula, Lu, Song, and Ma derive a generated approximation of the factorial function m ! , prove several gamma function inequalities, and deduce the approximation formula [20] as v
Γ ( v + 1 ) = 2 π v v e v v sinh 1 v + 1 810 v 7 v / 2 1 + O ( v 7 ) : = χ 10 ( v ) .
Chen presents the approximation formula [21] as v
Γ ( v + 1 ) = 2 π v v / e v 168 v 3 + 48 v + 7 168 v 3 + 48 v 7 v 2 + 53 / 210 1 + O v 7 : = χ 11 ( v ) ,
and then creates an asymptotic expansion using this approximation formula.
Windschitl’s approximation formula [15] as v
Γ ( v + 1 ) = 2 π v v e v v sinh 1 v + 1 810 v 6 v / 2 1 + O ( v 7 ) : = χ 12 ( v ) .
Alzer presents the approximation formula [22] as v
Γ ( v + 1 ) = 2 π v v e v v sinh 1 v v / 2 1 1620 v 5 + 1 1 + O ( v 7 ) : = χ 13 ( v ) ,
and a double inequality of Γ ( v + 1 ) for v > 0 with the best possible constants.
Yang and Tian present the approximation formula [23] as v
Γ ( v + 1 ) = 2 π v v e v v sinh 1 v v / 2 e 7 324 v 3 35 v 2 + 33 1 + O ( v 9 ) : = χ 14 ( v )
and study the monotonicity of some functions involving Γ ( v + 1 ) .
In light of the aforementioned results, the purpose of this paper is to present the following most accurate Mahmoud and Almuashi-type approximation:
Γ ( v + 1 ) = 2 π v v e v v 2 + 7 60 v 2 1 20 v / 2 e 461 907,200 v 3 v 2 + 5197 4610 1 + O ( v 9 ) : = χ 15 ( v ) , v
which is more accurate than Yang and Tian’s approximation formula χ 14 ( v ) in [23]. Also, we proved the following approximation formula
Γ ( v + 1 ) = 2 π v v e v v 2 + 7 60 v 2 1 20 v / 2 e 461 907,200 v 5 1 + O ( v 5 ) : = χ 16 ( v ) , v .
We have used Mathematica 10 software to perform the numerical calculations throughout this work.

2. Main Results

Recall that a real valued function g, which is defined and infinitely differentiable on v ( 0 , ) , is said to be completely monontonic (CM) if for all p 0 that ( 1 ) p g ( p ) ( v ) 0 on v ( 0 , ) . We refer to [24,25,26,27] for further information on CM functions and their applications. The function g is CM, according to Bernstein’s theorem [28], if and only if g ( v ) = 0 e v t d ς ( t ) , where ς ( t ) is a non-negative measure on t ( 0 , ) such that the integral converges for v ( 0 , ) . If the function g ( v ) is CM for v ( 0 , ) , lim v g ( v ) g ( ) and the function v ϖ g ( x ) g ( ) is CM for v ( 0 , ) if and only if ϖ [ 0 , δ ] , then the real number δ is called the completely monotonic degree of g ( v ) with respect to v ( 0 , ) (see [29,30]) and is denoted by deg cm v [ g ( v ) ] = δ . This concept can aid in more accurate measurements of CM functions.
To obtain our first bounds of gamma functions, we first give the following theorem:
Theorem 1.
Q 1 ( v ) = 1 v + 461 181,440 ( v + 1 ) 6 + 1 2 v + 2 + 1 20 v ( v + 2 ) + 19 + 7 60 v ( v + 2 ) + 67
1 2 ln ( v + 1 ) 2 ( v + 1 ) 2 + 7 60 ( v + 1 ) 2 1 20 + ψ ( v ) , v 0
is CM function with
deg cm v [ Q 1 ( v ) ] = 0 ,
where ψ ( v ) is the digamma function.
Proof. 
Using Gauss integral form of the digamma function [31]
ψ ( v ) = 0 1 t e t + e v t e t 1 d t , v > 0
we obtain
Q 1 ( v ) = 0 e t W 1 ( t ) t e v t d t ,
where
W 1 ( t ) = 461 t 6 21,772,800 + e t 2 5 t 4 5 t 2 e t 2 5 t 4 5 + e t t e t 1 + 1 2 e t 2 5 + 1 2 e t 2 5 1 2 7 15 t sin 1 2 7 15 t + cos 1 2 7 15 t 1 .
Consider the following functions for t 0
W 2 ( t ) = 461 t 6 21,772,800 + e t 2 5 t 4 5 t 2 e t 2 5 t 4 5 + 1 2 e t 2 5 + 1 2 e t 2 5 1 ,
W 3 ( t ) = 4799 t 6 217,728,000 t 4 3200 t 2 40 t 2 ,
W 4 ( t ) = e t t e t 1 + 1 2 7 15 t sin 1 2 7 15 t + cos 1 2 7 15 t ,
W 5 ( t ) = 1 2 7 15 t e t t e t 1 1 ,
and
W 6 ( t ) = 7 7 t 2 7 t 2 300 t 2 + 324,000 77,760,000 t 2 9,331,200,000 e t t e t 1 + 1 .
For t 0 , we have
W 2 ( t ) W 3 ( t ) = t 2 + 120 t 2 + 240 t 2 1,152,000 t sinh t 2 5 2 5 + cosh t 2 5 1 = p = 4 20 p t 2 p 2 p ( 2 p 2 ) ! < 0 .
Then, we have
W 2 ( t ) < W 3 ( t ) , t 0 .
Using the inequality
1 2 7 15 t sin 1 2 7 15 t + cos 1 2 7 15 t > 1 2 7 15 t 1 , t ,
we obtain
W 4 ( t ) > W 5 ( t ) , t .
Now,
W 4 ( t ) W 3 ( t ) > W 5 ( t ) W 3 ( t ) = 1 217,728,000 e t 1 W 7 ( t )
with
W 7 ( t ) = 4799 e t t 6 4799 t 6 + 68,040 e t t 4 68,040 t 4 + 5,443,200 e t t 2 5,443,200 t 2 7,257,600 105 e t t 108,864,000 e t t + 7,257,600 105 t 108,864,000 t 217,728,000 e t + 217,728,000 = W 8 ( t ) + n = 7 a n e 8 ( t 8 ) n n ! > 0 , t 8
where
a n = 4799 n 6 + 158,367 n 5 + 2,779,475 n 4 + 30,251,105 n 3 + 213,687,206 n 2 + 814,213,960 7,257,600 105 n + 1024 777,779 56,700 105 > 0 , n = 7,8,9 ,
and the polynomial
W 8 ( t ) = 9.90552 × 10 10 t 6 4.3738 × 10 12 t 5 + 8.10277 × 10 13 t 4 8.05282 × 10 14 t 3 + 4.52427 × 10 15 t 2 1.36148 × 10 16 t + 1.71352 × 10 16
is positive and has no real zeros for t 8 . Then,
W 4 > W 3 , t 8 .
Using the two inequalities [32]
cos ( t ) > t 6 720 + t 4 24 t 2 2 + 1 , t > 0
and
sin ( t ) < t 7 5040 + t 5 120 t 3 6 + t , t > 0
we obtain
W 4 ( t ) > W 6 ( t ) , t > 0 .
Then,
W 4 ( t ) W 3 ( t ) > W 4 ( t ) W 6 ( t ) = 2401 t 8 + 2,160,000 t 6 90,720,000 t 4 + 5,443,200,000 t 2 + 65,318,400,000 65,318,400,000 1 2 t coth t 2 .
Using
t 3 42 + t e t + 1 t 4 840 + 3 t 2 14 + 2 e t 1 = p = 9 ( p 8 ) ( p 7 ) ( p 6 ) ( p 5 ) 840 p ! t p ,
we obtain
coth t 2 < t 4 840 + 3 t 2 14 + 2 t 3 42 + t , t > 0
and hence,
W 4 ( t ) W 3 ( t ) > W 4 ( t ) W 6 ( t ) > t 8 ( 2,059,158 2401 t 2 ) 65,318,400,000 t 2 + 42 > 0 , 0 < t < 2,059,158 49 .
Then,
W 4 ( t ) > W 3 ( t ) , 0 < t 8 .
From inequalities (8)–(10), we obtain
W 1 ( t ) < 0 , t 0
and therefore, Q 1 ( v ) is CM function.
Now, if we suppose that v α Q 1 ( v ) is a CM function for v > 0 , then it is a decreasing function, that is
α v Q 1 ( v ) Q 1 ( v ) .
Using the following asymptotic expansion and its derivative
ψ ( v ) 1 v γ + π 2 v 6 ζ ( 3 ) v 2 + π 4 v 3 90 ζ ( 5 ) v 4 + , v 0
where ζ ( v ) is the Riemann zeta function for v > 1 , we obtain
α 81,155,027,309 49,004,796,960 π 2 6 v + v 2 2 ζ ( 3 ) 22,195,771,076,117 8,911,872,361,440 + O ( v 3 ) 152,361,413 230,973,120 γ 1 2 log 67 57 + π 2 6 81,155,027,309 49,004,796,960 v + O ( v 2 ) 0 as v 0 .
Hence, α = 0 or
deg cm v [ Q 1 ( v ) ] = 0 .
Theorem 2.
The function
Q 2 ( v ) = Γ ( v + 1 ) 2 π v v e v v 2 + 7 60 v 2 1 20 v / 2 e 461 907200 v 5 , v 1
is increasing. Furthermore, we have the following symmetric inequality:
a 1 2 π v v e v v 2 + 7 60 v 2 1 20 v / 2 e 461 907,200 v 5 < Γ ( v + 1 ) < a 2 2 π v v e v v 2 + 7 60 v 2 1 20 v / 2 e 461 907,200 v 5 , v 1
with best possible constants a 1 = Q 2 ( 1 ) 0.999733 and a 2 = 1 .
Proof. 
Using the relation
Q 1 ( v ) = d d v ln Q 2 ( v + 1 ) , v > 0
and with Q 1 ( v ) a CM function for v > 0 , then Q 2 ( v ) is an increasing function for v 1 . Also, using the asymptotic expansion [17]
Γ ( v + 1 ) 2 π v ( v / e ) v v 2 + 7 60 v 2 1 20 v / 2 1 + r = 1 λ r v r , v
where
λ r = 1 r j = 1 r j λ r j ρ j ρ r = B r + 1 r ( r + 1 ) 1 2 ω r + 1 ω r = ξ r 1 r j = 1 r 1 j ω j ξ r j , r = 1 , 2 , 3 ,
with
ξ 0 = 1 , ξ 2 s = 10 3 ( 20 ) s , ξ 2 s 1 = 0 , s = 1 , 2 , 3 ,
we obtain
Q 2 ( v ) e 461 907,200 v 5 1 + 461 907,200 v 5 5197 9,072,000 v 7 + 1,436,249 1,710,720,000 v 9 + , v
and
a 2 = lim v Q 2 ( v ) = 1 .
Remark 1.
From the expansion (13), we conclude that
Γ ( v + 1 ) = 2 π v v e v v 2 + 7 60 v 2 1 20 v / 2 e 461 907,200 v 5 1 + O ( v 5 ) : = χ 16 ( v ) , v .
Now, we present our second bounds of gamma function.
Theorem 3.
The function
G 1 ( v ) = Γ ( v + 1 ) 2 π v v e v v 2 + 7 60 v 2 1 20 v / 2 e 461 907,200 v 3 v 2 + 5197 4610 , v 13 10
is decreasing. Furthermore, we have the following symmetric inequality:
a 3 < Γ ( v + 1 ) 2 π v v e v v 2 + 7 60 v 2 1 20 v / 2 e 461 907,200 v 3 v 2 + 5197 4610 < a 4 , v 13 10
with best possible constants a 3 = 1 and a 4 = G 1 13 10 1.00000004 .
Proof. 
Consider the function
K 1 ( v ) = d d v ln G 1 v + 13 10 d d v ln G 1 v + 23 10 , v > 0
then,
K 1 ( v ) = 5 10 v + 13 2,449,304,525 61,255,978,812 ( 10 v + 13 ) 2 + 26,565,125 1,964,466 ( 10 v + 13 ) 4 5 10 v + 23 + 2,449,304,525 61,255,978,812 ( 10 v + 23 ) 2 26,565,125 1,964,466 ( 10 v + 23 ) 4 + 5 4 ( 5 v ( 5 v + 13 ) + 41 ) + 35 60 v ( 5 v + 13 ) + 542 + 1,129,129,386,025 61,255,978,812 ( 9220 v ( 5 v + 13 ) + 129,879 ) 5,645,646,930,125 2,946,699 ( 9220 v ( 5 v + 13 ) + 129 , 879 ) 2 5 20 v ( 5 v + 23 ) + 524 35 2 ( 30 v ( 5 v + 23 ) + 811 ) 1,129,129,386,025 61,255,978,812 ( 9220 v ( 5 v + 23 ) + 295,839 ) + 5,645,646,930,125 2,946,699 ( 9220 v ( 5 v + 23 ) + 295 , 839 ) 2 log v + 13 10 + log v + 23 10 1 2 log 25 6 ( 5 v ( 5 v + 13 ) + 41 ) + 1 + 1 2 log 25 30 v ( 5 v + 23 ) + 786 + 1
and
K 1 ( v ) = K 2 ( v ) K 3 ( v ) > 0 , v > 0
where
K 2 ( v ) = 9.76563 × 10 6 ( 1.40931 × 10 47 v 26 + 6.59557 × 10 48 v 25 + 1.48138 × 10 50 v 24 + 2.12565 × 10 51 v 23 + 2.1883 × 10 52 v 22 + 1.72065 × 10 53 v 21 + 1.0741 × 10 54 v 20 + 5.46199 × 10 54 v 19 + 2.30366 × 10 55 v 18 + 8.16187 × 10 55 v 17 + 2.45138 × 10 56 v 16 + 6.28063 × 10 56 v 15 + 1.37817 × 10 57 v 14 + 2.59529 × 10 57 v 13 + 4.19512 × 10 57 v 12 + 5.81144 × 10 57 v 11 + 6.8751 × 10 57 v 10 + 6.90658 × 10 57 v 9 + 5.84297 × 10 57 v 8 + 4.1146 × 10 57 v 7 + 2.37302 × 10 57 v 6 + 1.09554 × 10 57 v 5 + 3.91652 × 10 56 v 4 + 1.03052 × 10 56 v 3 + 1.83331 × 10 55 v 2 + 1.87037 × 10 54 v + 7.09549 × 10 52 )
and
K 3 ( v ) = 378 ( 10 v + 13 ) 5 ( 10 v + 23 ) 5 25 v 2 + 65 v + 41 2 25 v 2 + 115 v + 131 2 150 v 2 + 390 v + 271 2 150 v 2 + 690 v + 811 2 46,100 v 2 + 119,860 v + 129,879 3 46,100 v 2 + 212,060 v + 295,839 3 .
Then, K 1 ( v ) is increasing function with lim v K 1 ( v ) = 0 since
K 1 ( v ) 33,094,699,909,052,147 9,022,995,596,160,000 v 13 + 213,922,547 619,584,000 v 12 213,922,547 12,267,763,200 v 11 + , v .
Then, K 1 ( v ) < 0 or
d d v ln G 1 v + 13 10 < d d v ln G 1 v + 23 10 , v > 0 .
Using expansion (12), we obtain
d d v ln G 1 v 212,521 23,050 v 2 + 15,591 90,720 v 4 4610 v 2 + 5197 2 461 181,440 v 6 + 5197 1,296,000 v 8 1,436,249 190,080,000 v 10 + 26,863,154,077 1,273,708,800,000 v 12 + , v
and then lim v d d v ln G 1 v + 13 10 = 0 . But, if a function g : ( v 0 , ) R satisfies lim v g ( v ) = 0 and g ( v ) < g ( v + k ) , k N , then g ( v ) < 0 for v > v 0 (see [33]). Then, d d v ln G 1 v + 13 10 < 0 or ln G 1 v + 13 10 is decreasing for v > 0 . Hence, G 1 ( v ) is decreasing for v > 13 10 . Using the expansion
G 1 ( v ) e 461 907,200 5197 4610 v 2 + 1 v 5 1 + 461 907,200 v 5 5197 9,072,000 v 7 + 1,436,249 1,710,720,000 v 9 + , v
we obtain lim v G 1 ( v ) = 1 . Then, the function G 1 ( v ) is bounded by the two best possible constants a 3 = 1 and a 4 = G 1 ( 13 10 ) . □
Remark 2.
From the expansion (15), we conclude that
Γ ( v + 1 ) = 2 π v v e v v 2 + 7 60 v 2 1 20 v / 2 e 461 907,200 v 3 v 2 + 5197 4610 1 + O ( v 9 ) = χ 15 ( v ) , v .

3. Comparison among Some Approximation Formulas of the Gamma Function

The approximation formula χ 14 ( v ) is of course better than the other mentioned formulas from χ 1 ( v ) to χ 13 ( v ) based on the rate of convergence. Also, the two approximation formulas χ 14 ( v ) and χ 15 ( v ) are converging to Γ ( v + 1 ) with a rate like v 9 as v . To compare between the formulas χ 14 ( v ) and χ 15 ( v ) , consider the function
K 4 ( v ) = ln 3 v 20 v 2 1 sinh 1 v 60 v 2 + 7 + 1,597,365 v 2 + 3,172,927 45,360 v 4 35 v 2 + 33 4610 v 2 + 5197 , v 1
and then
K 4 1 v = 14 81 33 v 2 + 35 212,521 11,340 5197 v 2 + 4610 v 5 1 v + coth ( v ) + 245 81 35 v 2 + 33 2 97,972,181 2268 4610 v 2 + 5197 2 400 7 v 4 + 80 v 2 + 1200 v , 0 < v 1 .
Using the expansion
21 e 2 v + 1 v v 4 + 60 v 2 + 495 2 e v v 6 + 210 v 4 + 4725 v 2 + 10,395 sinh ( v ) 10,395 = p = 13 2 p 6 ( p 12 ) ( p 11 ) ( p 10 ) ( p 9 ) ( p 8 ) ( p 7 ) 10,395 p ! v p > 0 , v > 0
we obtain
coth ( v ) < v 6 10,395 + 2 v 4 99 + 5 v 2 11 + 1 v 5 495 + 4 v 3 33 + v , v > 0 .
Then,
K 4 1 v < v 9 K 5 ( v ) K 6 ( v ) < 0 , 0 < v 1
where
K 5 ( v ) = 3,809,121,073,989 v 10 + 308,281,165,523,325 v 8 + 19,393,944,515,758,250 v 6 + 272,405,703,557,635,250 v 4 + 603,315,377,277,995,625 v 2 + 344,328,389,211,515,625
and
K 6 ( v ) = 11 , 340 v 2 20 7 v 2 + 60 33 v 2 + 35 2 5197 v 2 + 4610 2 v 4 + 60 v 2 + 495 .
Hence, K 4 ( v ) is an increasing function for v 1 . Also,
K 4 ( v ) 758 , 772 , 059 3,864,345,408,000 v 10 + 2 , 050 , 016 , 121 , 570 , 689 4,863,394,626,330,240,000 v 12 + , v
and hence, lim v K 4 ( v ) = 0 . Then, K 4 ( v ) < 0 for v 1 or
e 7 324 v 3 35 v 2 + 33 v sinh 1 v v / 2 < v 2 + 7 60 v 2 1 20 v / 2 e 461 907,200 v 3 v 2 + 5197 4610 , v 1 .
Then, using the inequalities (14) and (16), we obtain
χ 14 ( v ) < χ 15 ( v ) < Γ ( v ) , v > 13 / 10 .
Therefore, the approximation formula χ 15 ( v ) is better than χ 14 ( v ) for v 13 10 .

4. Conclusions

Applications of the gamma function may be found in a wide range of real-world fields, from finance and economics (e.g., option pricing models) to medical research (e.g., modeling disease spread). Approximation formulae can improve these studies’ precision and speed. The main conclusions of this paper are stated in Theorems (2) and (3). Concretely speaking, based on Mahmoud and Almuashi’s formula, the authors studied the monotonicity and complete monotonicity of some functions related to Γ ( v ) to present the Formulas (4) and (5) and some symmetric inequalities for Γ ( v ) . Our new approximation formula χ 15 ( v ) and Yang and Tian’s approximation χ 14 ( v ) are of the same order v 9 but the superiority of our results is proven.

Author Contributions

Writing of original draft, M.M. and H.A. All authors contributed equally to the writing of this paper. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Cao, X.; Wang, R. Some inequalities for the ratio of Gamma functions. J. Inequal. Appl. 2015, 2015, 178. [Google Scholar] [CrossRef]
  2. Alzer, H.; Jameson, G. A harmonic mean inequality for the Digamma function and related results. Rend. Sem. Mat. Univ. Padova 2017, 137, 203–209. [Google Scholar] [CrossRef] [PubMed]
  3. Yang, Z.-H.; Tian, J. Two asymptotic expansions for Gamma function developed by Windschitl’s formula. Open Math. 2017, 16, 1048–1060. [Google Scholar] [CrossRef]
  4. Nisar, K.S.; Qi, F.; Rahman, G.; Mubeen, S.; Arshad, M. Some inequalities involving the extended Gamma function and the Kummer confluent hypergeometric k-function. J. Inequal. Appl. 2018, 2018, 135. [Google Scholar] [CrossRef] [PubMed]
  5. Qi, F.; Guo, B.-N. From inequalities involving exponential functions and sums to logarithmically complete monotonicity of ratios of Gamma functions. J. Math. Anal. Appl. 2021, 493, 124478. [Google Scholar] [CrossRef] [PubMed]
  6. Reynolds, R.; Stauffer, A. A Note on the summation of the incomplete Gamma function. Symmetry 2021, 13, 2369. [Google Scholar] [CrossRef]
  7. Tian, J.-F.; Yang, Z. Logarithmically complete monotonicity of ratios of q-gamma functions. J. Math. Anal. Appl. 2022, 508, 125868. [Google Scholar] [CrossRef]
  8. Zhang, J.-M.; Yin, L.; You, H.-L. Complete monotonicity and inequalities involving the k-Gamma and k-Polygamma functions. Math. Slovaca 2023, 73, 1217–1230. [Google Scholar] [CrossRef]
  9. Yildirim, E. Complete monotonicity of functions involving k-Trigamma and k-Tetragamma functions with related inequalities. Turkish J. Ineq. 2023, 7, 12–21. [Google Scholar]
  10. Castillo, J.S.; Rojas, M.A.; Reyes, J. A more flexible extension of the Fréchet distribution based on the incomplete Gamma function and applications. Symmetry 2023, 15, 1608. [Google Scholar] [CrossRef]
  11. Mahmoud, M.; Alsulami, S.M.; Almarashi, S. On Some Bounds for the Gamma Function. Symmetry 2023, 15, 937. [Google Scholar] [CrossRef]
  12. Ismail, G.M.; El-Moshneb, M.M.; Zayed, M. Analytical technique for solving strongly nonlinear oscillator differential equations. Alex. Eng. J. 2023, 74, 547–557. [Google Scholar] [CrossRef]
  13. Wang, K.J.; Wang, G.D. Gamma function method for the nonlinear cubic-quintic Duffing oscillators. J. Low Freq. Noise Vib. Active Control 2022, 41, 216–222. [Google Scholar] [CrossRef]
  14. Andrews, G.E.; Berndt, B.C. Ramanujan’s Lost Notebook: Part IV; Springer Science+ Business Media: New York, NY, USA, 2013. [Google Scholar]
  15. Programmable Calculators. 2002. Available online: http://www.rskey.org/CMS/the-library/11 (accessed on 20 April 2020).
  16. Smith, W.D. The Gamma Function Revisited. 2006. Available online: http://schule.bayernport.com/gamma/gamma05.pdf (accessed on 20 April 2020).
  17. Mahmoud, M.; and Almuashi, H. On Some Asymptotic Expansions for the Gamma Function. Symmetry 2022, 14, 2459. [Google Scholar] [CrossRef]
  18. Nemes, G. New asymptotic expansion for the Gamma function. Arch. Math. 2010, 95, 161–169. [Google Scholar] [CrossRef]
  19. Yang, Z.-H.; Chu, Y.-M. Asymptotic formulas for gamma function with applications. Appl. Math. Comput. 2015, 270, 665–680. [Google Scholar] [CrossRef]
  20. Lu, D.; Song, L.; Ma, C. A generated approximation of the gamma function related to Windschitl’s formula. J. Number Theory 2014, 140, 215–225. [Google Scholar] [CrossRef]
  21. Chen, C.-P. A more accurate approximation for the gamma function. J. Number Theory 2016, 164, 417–428. [Google Scholar] [CrossRef]
  22. Alzer, H. Sharp upper and lower bounds for the Gamma function. Proc. R. Soc. Edinb. Sect. Math. 2009, 139, 709–718. [Google Scholar] [CrossRef]
  23. Yang, Z.-H.; Tian, J.-F. An accurate approximation formula for Gamma function. J. Inequal. Appl. 2018, 2018, 56. [Google Scholar] [CrossRef]
  24. Alzer, H. On some inequalities for the gamma and psi function. Math. Comput. 1997, 66, 373–389. [Google Scholar] [CrossRef]
  25. Burić, T.; Elezović, N. Some completely monotonic functions related to the psi function. Math. Inequal. Appl. 2011, 14, 679–691. [Google Scholar]
  26. Qi, F.; Wang, S.-H. Complete monotonicity, completely monotonic degree, integral representations, and an inequality related to the exponential, trigamma, and modified Bessel functions. Glob. J. Math. Anal. 2014, 2, 91–97. [Google Scholar] [CrossRef]
  27. Qi, F.; Guo, S.; Guo, B.-N. Complete monotonicity of some functions involving polygamma functions. J. Comput. Appl. Math. 2010, 233, 2149–2160. [Google Scholar] [CrossRef]
  28. Widder, D.V. The Laplace Transform; Princeton University Press: Princeton, NJ, USA, 1946. [Google Scholar]
  29. Guo, B.-N.; Qi, F. A completely monotonic function involving the tri-gamma function and with degree one. Appl. Math. Comput. 2012, 218, 9890–9897. [Google Scholar] [CrossRef]
  30. Koumandos, S. Monotonicity of some functions involving the gamma and psi functions. Math. Comp. 2008, 77, 2261–2275. [Google Scholar] [CrossRef]
  31. Abramowitz, M.; Stegun, I.A. (Eds.) Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables, National Bureau of Standards, Applied Mathematics Series 55; 9th printing; Dover Publications: New York, NY, USA; Washington, DC, USA, 1972. [Google Scholar]
  32. Qi, F.; Niu, D.W.; Guo, B.N. Refinements, generalizations, and applications of Jordan’s inequality and related problems. J. Inequal. Appl. 2009, 2009, 271923. [Google Scholar] [CrossRef]
  33. Elbert, A.; Laforgia, A. On some properties of the Gamma function. Proc. Am. Math. Soc. 2000, 128, 2667–2673. [Google Scholar] [CrossRef]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Mahmoud, M.; Almuashi, H. New Accurate Approximation Formula for Gamma Function. Symmetry 2024, 16, 150. https://doi.org/10.3390/sym16020150

AMA Style

Mahmoud M, Almuashi H. New Accurate Approximation Formula for Gamma Function. Symmetry. 2024; 16(2):150. https://doi.org/10.3390/sym16020150

Chicago/Turabian Style

Mahmoud, Mansour, and Hanan Almuashi. 2024. "New Accurate Approximation Formula for Gamma Function" Symmetry 16, no. 2: 150. https://doi.org/10.3390/sym16020150

APA Style

Mahmoud, M., & Almuashi, H. (2024). New Accurate Approximation Formula for Gamma Function. Symmetry, 16(2), 150. https://doi.org/10.3390/sym16020150

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop