We consider the question of whether the high-energy eigenfunctions of certain Schrödinger operators on the d-dimensional hyperbolic space of constant curvature −κ2 are flexible enough to approximate an arbitrary solution of the Helmholtz equation Δh + h = 0 on Rd, over the natural length scale O(λ−1/2) determined by the eigenvalue λ ≫ 1. This problem is motivated by the fact that, by the asymptotics of the local Weyl law, approximate Laplace eigenfunctions do have this approximation property on any compact Riemannian manifold. In this paper we are specifically interested in the Coulomb and harmonic oscillator operators on the hyperbolic spaces Hd(κ). As the dimension of the space of bound states of these operators tends to infinity as κ ↘ 0, one can hope to approximate solutions to the Helmholtz equation by eigenfunctions for some κ > 0 that is not fixed a priori. Our main result shows that this is indeed the case, under suitable hypotheses. We also prove a global approximation theorem with decay for the Helmholtz equation on manifolds that are isometric to the hyperbolic space outside a compact set, and consider an application to the study of the heat equation on Hd(κ). Although global approximation and inverse approximation results are heuristically related in that both theorems explore flexibility properties of solutions to elliptic equations on hyperbolic spaces, we will see that the underlying ideas behind these theorems are very different.