Erst mal das char. Pol bestimmen, also
det( x*E - M ) und das ist hier x^4.
Hat also nur die Nullstelle 0, also hat M  nur einen
Eigenwert, nämlich 0.
Eigenraum dazu ist die Lösungsmenge von
M * x = 0 * x  =  Nullvektor.
Stufenform ergibt z.B.
7  10  9 15
0  1   3    -7
0  0    0     0 
0  0    0     0
also Basis des Eigenraums z.B.  (-5 ; 2 ; 0 ; 1 ) und ( 3 ; -3 ; 1 ; 0 )
und da bereits M^2 = Nullmatrix ist, ist x^2 das Min.pol.
jetzt noch 2 lin. unabhängige finden, deren Bilder diese beiden
Basisvektoren sind: 
( 4 ; -3 ; 0 ; 0 ) und  ( -1 ; 1 ; 0 ; 0 ) .
Dann ist die Transformationsmatrix also   T=
-5   4   3   -1
2   -3   -3    1
0    0     1    0
1    0     0    1
und es ist T -1 * M * T die gesuchte Jord. Normalf.