Find the eigenvalues of this matrix I was wondering if anyone wanted to give a shot at finding the eigenvalues of this $4\times 4$ system. I have tried without success, mainly because I end up having to solve a very nasty cubic polynomial that wolfram seems to be having trouble with as well. There may be some method to solve this that I didn't try, so I wanted to post this here just in case someone could help me out,
$$J(\epsilon)=\begin{bmatrix} -\mu &&0&&0 && -\beta b/\mu\\ 0&& -\sigma && 0&& \beta b/\mu\\0&& \sigma&& -(\phi+d)&& 0\\ 0 &&0&&\alpha&&-(\phi+\epsilon)\end {bmatrix}$$
Thanks!

You can do row reduction  of the matrix $$J(\epsilon)=\begin{bmatrix} -\mu &&0&&0 && -\beta b/\mu\\ 0&& -\sigma && 0&& \beta b/\mu\\0&& \sigma&& -(\phi+d)&& 0\\ 0 &&0&&\alpha&&-(\phi+\epsilon)\end {bmatrix}$$
add row 2 and row 3
$$J(\epsilon)=\begin{bmatrix} -\mu &&0&&0 && -\beta b/\mu\\
 0&& -\sigma && 0&& \beta b/\mu\\
 0&& 0&& -(\phi+d)&& \beta b/\mu\\
 0 &&0&&\alpha&&-(\phi+\epsilon)\end {bmatrix}$$
row 3 $\times \alpha /(\phi+d)+$row 4
$$J(\epsilon)=\begin{bmatrix} -\mu &&0&&0 && -\beta b/\mu\\
 0&& -\sigma && 0&& \beta b/\mu\\
 0&& 0&& -(\phi+d)&& \beta b/\mu\\
 0 &&0&&0&&-(\phi+\epsilon)+(\beta b/\mu )(\alpha /(\phi+d))\end {bmatrix}$$
now determinant will be the multiplication of the diagonal elements 
$$|J-\lambda I|=(-\mu-\lambda)(-\sigma-\lambda)(-(\phi+d)-\lambda)(-(\beta b/\mu )(\alpha /(\phi+d))-\lambda)=0$$
