Abstract
To obtain closed-form solutions for the radial Schrodinger wave equation with non-solvable potential models, we use a simple, easy, and fast perturbation technique within the framework of the asymptotic iteration method (PAIM). We will show how the PAIM can be applied directly to find the analytical coefficients in the perturbation series, without using the base eigenfunctions of the unperturbed problem. As an example, the vector Coulomb (similar to 1/r) and the harmonic oscillator (similar to r(2)) (similar to r) scalar potential parts implemented with their counterpart spin-dependent terms are chosen to investigate the meson sectors including charm and beauty quarks. This approach is applicable in the same form to both the ground state and the excited bound states and can be easily applied to other strongly non-solvable potential problems. The procedure of this method and its results will provide a valuable hint for investigating tetraquark configuration.