Matrix form.
(un+1un)=(31−20)(unun−1)=An(u1u0).
Eigenvalues. det(A−λI)=(3−λ)(−λ)+2=λ2−3λ+2=(λ−1)(λ−2), so λ1=1, λ2=2.
Eigenvectors. For λ=1: v1=(1,1). For λ=2: v2=(2,1).
Diagonalization.
A=P(1002)P−1,P=(1121),P−1=(−112−1).
Power.
An=P(1002n)P−1=(1121)(1002n)(−112−1).
Compute:
An=(−1+2n+1−1+2n2−2n+12−2n).
Conclusion. un=(An)21u1+(An)22u0=(2n−1)⋅1+0=2n−1.
Verification. u0=0, u1=1, u2=3, u3=7: matches 2n−1.