We propose an analytical five-parameter potential function for stable diatomic molecules V(x) = K/x³ - e-tx(a + bx + cx²). Here x = r/re, where r is the internuclear distance and re the equilibrium bond length. K, t, a, b, c are five parameters of the function, which are obtained from molecular spectroscopic data in terms of known molecular constants.