Abstract
In this work the Grimm and Storer method is employed to calculate all the higher excited states wave functions and their corresponding energy eigenvalues , for the first time by modifying the iterative numerical procedure. Using symmetry arguments of an even Hamiltonian, the even and odd excited states can be calculated separately