Accurate analytic presentation of solution of the Schrödinger equation with arbitrary physical potential

High precision approximate analytic expressions of the ground state energies and wave functions for the arbitrary physical potentials are found by first casting the Schrödinger equation into the nonlinear Riccati form and then solving that nonlinear equation analytically in the first iteration of the quasilinearization method (QLM). In the QLM the nonlinear differential equation is treated by approximating the nonlinear terms by a sequence of linear expressions. The QLM is iterative but not perturbative and gives stable solutions to nonlinear problems without depending on the existence of a smallness parameter. The choice of zero iteration is based on general features of exact solutions near the boundaries. The approach is illustrated on the examples of the Yukawa, Woods-Saxon and funnel potentials. For the latter potential, solutions describing charmonium, bottonium and topponium are analyzed. Comparison of our approximate analytic expressions for binding energies and wave functions with the exact numerical solutions demonstrates their high accuracy in the wide range of physical parameters. The accuracy ranging between 10^-4 and 10^-8 for the energies and, correspondingly, 10^-2 and 10^-4 for the wave functions is reached. The derived formulas enable one to make accurate analytical estimates of how variation of different interactions parameters affects correspondent physical systems.