Deriving effective Hamiltonian models plays an essential role in quantum theory, with particular emphasis in recent years on control and engineering problems. In this work, we present two symbolic methods for computing effective Hamiltonian models: the Non-perturbative Analytical Diagonalization (NPAD) and the Recursive Schrieffer-Wolff Transformation (RSWT). NPAD makes use of the Jacobi iteration and works without the assumptions of perturbation theory while retaining convergence, allowing to treat a very wide range of models. In the perturbation regime, it reduces to RSWT, which takes advantage of an in-built recursive structure where remarkably the number of terms increases only linearly with perturbation order, exponentially decreasing the number of terms compared to the ubiquitous Schrieffer-Wolff method. Both methods consist of elementary algebra and can be easily automated to obtain closed-form expressions. To demonstrate the application of the methods, we study the ZZ interaction of superconducting qubits systems in two different parameter regimes. In the near-resonant regime, the method produces an analytical expression for the effective ZZ interaction strength, with an improvement of at least one order of magnitude compared to neglecting the comparably further detuned state. In the dispersive regime, we calculate perturbative corrections up to the 8th order and accurately identify a regime where ZZ interaction is suppressed in a simple model with only resonator-mediated coupling.