1. The problem statement, all variables and given/known data This problem is from Zetelli 3.21 http://imgur.com/wYTNVwz http://imgur.com/wYTNVwz 2. Relevant equations Just the standard probability via product between the eigenfunction and the wavefunction 3. The attempt at a solution I've found the eigenvectors for the Hamiltonian (-i,5,3) (-i,-2,3) (3i,0,1) I can orthonormalize them, as well as normalize the initial state |y(0)> however I run into a problem when actually calculating the probabilities. I only get complex values and can't seem to get rid of them, especially since the normaliztion factor for |y(0)> is 1/sqrt(-1-16i). Can anyone help?