Wenn a∈K ein Eigenwert von M ist, dann gibt es ein v∈ℝ^n \ {0} mit
M*v = a*v
Dann ist auch M^2*v = M( M(v)) = M(a*v)
wegen der Linearität = a*M(v) = a*a*v = a^2 * v.
Entsprechend M^n*v=a^n*v.
Ist nun q = bkX^k + bk-1X^(k-1) + ... + b0 aus K[x]
dann q(M)*v = bkM^k*v + bk-1M^(k-1)*v + ... + b0*v
= bka^k*v + bk-1a^(k-1)*v + ... + b0*v
= q(a) * v q.e.d.