A method for calculating vibrational bound states: Iterative solution of the collocation equations constructed from localized basis sets