“…In this paper, we study the complete monotonicity of the functions F a , then we apply these results to obtain new sharp bounds for the digamma and trigamma functions. The problem of estimating the gamma and polygamma functions has attracted the attention of many researchers, since they are close related to the theory of zeta functions [1,8,12,16,27,39], multiple gamma and related functions [7,[9][10][11]13,[24][25][26]30,33,35,38,40], gamma type distributions [15,29], or harmonic sums [23,34]. There are also many recent investigations dealing with one-sided and two-sided inequalities involving the digamma, trigamma, polygamma and other related functions, see e.g., [5,6,14,31,32].…”