p = poly(r), where r is a vector, returns the coefficients of the
polynomial whose roots are the elements of r.
p = poly(A), where A is an n-by-n matrix, returns the n+1 coefficients
of the characteristic polynomial of the matrix, det(λI – A).