1 | initial version |

Hello,

It seems that Sage is not able to solve your equation symbolically. However, you can compute the characteristic polynomial and try to analyze it

```
R = ZZ['x,sigma,gamma,omega,Q1,Q2,Q3,Q4,A1,A2,A3,mu_v']
R.inject_variables()
J = matrix(R,[[-Q1,0,0,0,omega,0,0,0,0,-A1],
[0,-Q2,0,0,0,omega,0,0,0,A1],
[0,gamma,Q3,0,0,0,omega,0,0,0],
[0,0,sigma,-Q1,0,0,0,omega,0,0],
[A2,0,0,0,-Q4,0,0,0,0,0],
[0,A2,0,0,0,-Q4,0,0,0,0],
[0,0,A2,0,0,0,-Q4,0,0,0],
[0,0,0,A2,0,0,0,-Q4,0,0],
[0,0,-A3,0,0,0,0,0,-mu_v,0],
[0,0,A3,0,0,0,0,0,0,-mu_v]])
p = J.charpoly().subs(R('x'))
p.factor()
```

As you can notice I did not use symbolic variables but polynomial variables that are more powerful (and more predictable). In particular, the factorization with **p.factor()** takes place in Z[x,sigma,gamma,Q1,...]`.

The answer of the above code shows two small factors

```
x + mu_v
(-x^2 - x*Q1 - x*Q4 - Q1*Q4 + omega*A2)^2
```

and a big factor of degree 5. The small factors give you already three eigenvalues (five if counted with multiplicity). The big factor is quite complicated and I doubt there would be any formula.

Vincent

Copyright Sage, 2010. Some rights reserved under creative commons license. Content on this site is licensed under a Creative Commons Attribution Share Alike 3.0 license.