Ask Your Question

Revision history [back]

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