Smith-McMillan Form of a polynomial matrix
sage: s = var('s')
sage: G = matrix([[1,-1],[s^2+s-4,2*s^2-s-8],[s^2-4, 2*s^2-8]])
G
is a symbolic matrix, each entry of which is a polynomial of variable s
.
How can we calculate the Smith McMillan Form of the matrix G?
NOTE: G.smith_form()
doesn't work.
Thanks in advance!