In this work we have applied ansatz method to solve for the approximate ℓ-state solution of time independent Schrödinger wave equation with modified Möbius squared potential plus Hulthén potential to obtain closed form expressions for the energy eigenvalues and normalized radial wave-functions. In dealing with the spin-orbit coupling potential of the effective potential energy function, we have employed the Pekeris type approximation scheme, using our expressions for the bound state energy eigenvalues, we have deduced closed form expressions for the bound states energy eigenvalues and normalized radial wave-functions for Hulthén potential, modified Möbius square potential and Deng-Fan potential. Using the value 0.976865485225 for the parameter ω, we have computed bound state energy eigenvalues for various quantum states (in atomic units). We have also computed bound state energy eigenvalues for six diatomic molecules: HCl, LiH, TiH, NiC, TiC and ScF. The results we obtained are in near perfect agreement with numerical results in the literature and a clear demonstration of the superiority of the Pekeris-type approximation scheme over the Greene and Aldrich approximation scheme for the modified Möbius squares potential plus Hulthén potential.