[gujianlin]

各位matlab高手，
帮忙求一下下面矩阵的特征值和特征方程吧，多谢啦！


[ -d-k*I/(1+a*I^2), 0, (k*s*a*I^2-k*s)/(1+a*I^2)^2, r]
[ -d, -e-d, 0, r]
[ 0, e, -u-d, 0]
[ 0, 0, u, -d-r]
所有字母均为正常数