Hello,

My question is regarding the Hamiltonian of the cooper pair box. The Hamiltonian looks like $$H=E_C(n-n_g)^2\left|n\right\rangle\left\langle n \right|-\frac{E_J}{2}\left[\left|n\right\rangle\left\langle n+1\right|+\left|n+1\right\rangle\left\langle n \right|\right]$$ Which in matrix form looks like $$\begin{pmatrix} E_C(0-n_g)^2 & -\frac{E_J}{2} & 0 & 0 & 0 & 0\\ -\frac{E_J}{2} & E_C(1-n_g)^2 & -\frac{E_J}{2} & 0 & 0 & 0\\ 0 & -\frac{E_J}{2} & E_C(2-n_g)^2 & -\frac{E_J}{2} & 0 & 0\\ 0 & 0 & -\frac{E_J}{2} & E_C(3-n_g)^2 & -\frac{E_J}{2} & 0\\ 0 & 0 & 0 & -\frac{E_J}{2} & E_C(4-n_g)^2 & -\frac{E_J}{2} \\ 0 & 0 & 0 & 0 & -\frac{E_J}{2} & E_C(5-n_g)^2 \end{pmatrix}$$ For states 0 to 5. I entered this matrix into mathematica like

matr = {{Ec*(0 - ng)^2, -Ej/2, 0, 0, 0, 0}, 
       {-Ej/2, Ec*(1 - ng)^2, -Ej/2, 0, 0, 0}, 
       {0, -Ej/2, Ec*(2 - ng)^2, -Ej/2, 0, 0}, 
       {0, 0, -Ej/2, Ec*(3 - ng)^2, -Ej/2, 0}, 
       {0, 0, 0, -Ej/2, Ec*(4 - ng)^2, -Ej/2}, 
       {0, 0, 0, 0, -Ej/2, Ec*(5 - ng)^2}}

Where Ec is $E_C$, Ej is $E_J$ and ng is $n_g$. Then I used the Eigenvalues function with the following code

Eigenvalues[matr]

This produced a lot of lines of output which was fairly indecipherable. Is this the correct way to use this function to return the eigenvalues of this matrix? I would be grateful for any help.