  We discuss the application of the Discrete Variable Representation to
Schr\"odinger problems which involve singular Hamiltonians. Unlike recent
authors who invoke transformations to rid the eigenvalue equation of
singularities at the cost of added complexity, we show that an approach based
solely on an orthogonal polynomial basis is adequate, provided the
Gauss-Lobatto or Gauss-Radau quadrature rule is used. This ensures that the
mesh contains the singular points and by simply discarding the DVR functions
corresponding to those points, all matrix elements become well-behaved, the
boundary conditions are satisfied and the calculation is rapidly convergent.
The accuracy of the method is demonstrated by applying it to the hydrogen atom.
We emphasize that the method is equally capable of describing bound states and
continuum solutions.
