The Rayleigh–Schrodinger method or more commonly called time independent perturbation theory is a powerful method to solve a Hamiltonian system with a relatively small perturbation [W. Nolting, Theoretical Physics 7: Quantum Mechanics–Methods and Applications (Springer International, 2017); F. Schwabl, Quantum Mechanics, 4th edn. (Springer, 2007)]. In this work, we make use of the Rayleigh–Schrodinger method to formulate a general method to calculate quasibound states of static spherically symmetric black hole solutions arising from modified theories of gravity. We discover that the Schwarzschild-like term corresponds to the main Hamiltonian while the modified theory of gravity’s contribution of the spacetime metric corresponds to the perturbation terms. At the end, formula to calculate main Hamiltonian and perturbed Hamiltonian are discovered and presented.