Abstract
A nonperturbative but absolutely convergent algorithm is applied to the determination of the eigenfunctions of Hamiltonian with a singular potential of the form H=−d2/dx2+x2+λ/xα in the domain [0, ∞] which obey Dirichlet boundary conditons. The formal structure of the algorithm is identical to that of the Lanczos algorithm when it is formally extended to self-adjoint operators.