Abstract
The adaptive perturbation chooses a non-standard decomposition. The Hamiltonian becomes a sum of solvable and perturbation parts. We calculate the spectrum using the adaptive perturbation method at the leading-order to compare to numerical solutions. The maximum deviation is around 5% for different coupling regions. A perturbation study relies on whether a choice of leading-order is suitable. Our result with different parameters should show that the adaptive perturbation method provides appropriate saddle points to all coupling regions. In the end, we show that the perturbation parameters should not be a coupling constant.