An exactly separable version of the Bohr Hamiltonian is developed using a potential of the form u(β) + u(γ)/β 2 , with the Davidson potential u(β) = β 2 + β 4 0 /β 2 (where β 0 is the position of the minimum) and a stiff harmonic oscillator for u(γ) centered at γ = 0 • . In the resulting solution, called exactly separable Davidson (ES-D), the ground state band, γ band and 0 + 2 band are all treated on an equal footing. The bandheads, energy spacings within bands, and a number of interband and intraband B(E2) transition rates are well reproduced for almost all well-deformed rare earth and actinide nuclei using two parameters (β 0 , γ stiffness). Insights regarding the recently found correlation between γ stiffness and the γ-bandhead energy, as well as the long standing problem of producing a level scheme with Interacting Boson Approximation SU(3) degeneracies from the Bohr Hamiltonian, are also obtained.