# Smith normal form of a polynomial matrix.

I have the following matrix: $P(s) :=$

\begin{bmatrix} s^2 & s-1 \\ s & s^2 \end{bmatrix}

How does one compute the Smith normal form of this matrix? I can't quite grasp the algorithm.

-