We present arguments demonstrating that the application of the Nikiforov-Uvarov polynomial method to solve the Schrödinger equation with the Tietz-Hua potential is valid only when e −b h re ≤ c h < 1 and r0 < r < +∞. In particular, it is point out that the numerical results with c h = 0 for the diatomic molecules HF, N2, I2, H2, O2 and O + 2 given in Tables 3-5 by Hamzavi and co-workers are wrong. When −1 < c h < 0 or 0 < c h < e −b h re , this approach is not suitable. In both cases, it is shown that the solutions of the Schrödinger equation are expressed in terms of the generalized hypergeometric functions 2F1(a, b, c; z). The determination of the energy levels requires the solution of transcendental equations involving the hypergeometric function by means of the numerical procedure.