OK, so far so good. I asked for the eigenvalues of the matrix A and it shows the 4 eigenvalues that I expected. So then for the sake of understanding the solver, I though I would try to compute the eigenvalues by setting the determinant of (sI-A) equal to zero:
>> solve('det(s*eye(4)-A) = 0',s) Warning: Explicit solution could not be found.
Why doesn't this work? After some experimentation I found this works:
So as you can see the solver can compute the eigenvalues correctly, but this form of having to convert the determinant to a character array seems kind of awkward. Is that necessary, or is there a simpler way of entering this equation into the solver?