# 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!