Like I said, for any Hamiltonian, if you know the eigenvectors and eigenvalues, you can express an arbitrary wavefunction as a linear combination:
Given the states \left|\psi_n\right> where \hat{H}\left|\psi_n\right>=E_n \left|\psi_n\right>, you can write:
\left|\Psi\right> = \sum_n c_n \left|\psi_n\right>.
The time evolution of the eigenstates are easy:
i\hbar \frac{d\left|\psi_n(t)\right>}{dt} = \hat{H}\left|\psi_n(t)\right> = E_n \left|\psi_n(t)\right>
which has solutions: \left|\psi_n(t)\right> = e^{-iE_n t/\hbar} \left|\psi_n(0)\right>, because E_n is just a number and not an operator, so you can use the usual algebraic manipulations.
Finally, since the Schroedinger equation is linear, you have that:
\left|\Psi(t)\right> = \sum_n c_n \left|\psi_n(t)\right> = \sum_n c_n e^{-iE_n t/\hbar} \left|\psi_n(0)\right>.
Note that crucially you need to be able to find the eigenstates and eigenvalues of the Hamiltonian, which is usually quite hard.