For a large family of real-valued Radon measures m on R^d, including the Kato
class, the operators -\Delta + C^2 \Delta^2 + m tend to the Schrodinger
operator -\Delta +m in the norm resolvent sense as C tends to zero. If the
measure is moreover finite and the dimension smaller than four, the former
operator can be approximated by a sequence of operators with point measures in
the norm resolvent sense. The combination of both convergence results leads to
an efficient method for the numerical computation of the eigenvalues in the
discrete spectrum and corresponding eigenfunctions of Schrodinger operators. We
illustrate the approximation by numerical calculations of eigenvalues for one
simple example of measure m.