Hi! This is my first post here. If I make any mistake, please tell me. I have a question: if I want to find a pole of a polynomial, ex: " a*s^2+b*s+c ", a,b,and c are constant given. I just use simple command: roots([a, b,c]), I can get the answer. However, if now I meet a situation is " A*s^2+C=0 ", A and C are all 3*3 constant matrices, I used the " sym and solve " to find the equation. It's OK. However, except this method, is there any simple command, just like " roots " to solve the polynomial with matrix coefficient? Thanks!