Hi blue_leaf! Thanks for that, I just realized that I interpreted the question wrong. I thought the matrix given was only the unperturbed Hamiltonian.
So I have separated it into its un/pertrubed matrices. I just don't think I understand what the question is asking, when it says the change in eigenvalues.
##E_1\approx E_1^{(0)}+E_1^{(1)}+E_1^{(2)}##
##~~~~~=E_0+0-\frac{4}{7}E_0\lambda^2##
##E_2\approx E_2^{(0)}+E_2^{(1)}+E_2^{(2)}##
##~~~~~=8E_0+0+\frac{4}{7}E_0\lambda^2##
##E_3\approx E_3^{(0)}+E_3^{(1)}+E_3^{(2)}##
##~~~~~=3E_0+E_0\lambda+0##
##E_4\approx E_4^{(0)}+E_4^{(1)}+E_4^{(2)}##
##~~~~~=7E_0+0+0##
So are the changes in the eigenvalues ##-\frac{4}{7}E_0\lambda^2,~\frac{4}{7}E_0\lambda^2,~E_0\lambda \text{ and } 0## respectively ?