I am using Mathematica to construct a matrix for the Hamiltonian of some system. I have built this matrix already, and I have found the eigenvalues and the eigenvectors, I am uncertain if what I did next is correct: I took the normalized eigenvectors, placed them in matrix form, and did matrix multiplication with the basis set of solutions.
Let me try to be more precise since I am not sure I am using the right language when mentioning the basis solutions. In the problem we are using the set of solutions of the particle in a box model as our basis. I can increase the number of basis elements in the calculation of the matrix of the Hamiltonian (which amounts to doing $<\psi_n|H \psi_k>$ over a specified range of $n$ and $k$) in order for some of my smallest eigenvalues to begin to converge. Once I have this $H$ matrix built, and that I see that my eigenvalues are converging to some degree, I take the eigenvectors of the $H$ matrix, format them to be in matrix form, and multiply them by the set of basis solutions.
I hope that makes things clearer.