summaryrefslogtreecommitdiff
path: root/26/CH5/EX5.3.1/5_3_1.sce
blob: 101f3dbc2e96798450b71ab5c5b162c1e896477c (plain)
1
2
3
4
5
6
7
8
9
10
11
disp('The given eigenvector matrix is:')
p=[5 7;2 3]
disp(p,'P=')
disp('The diagonal matrix is:')
d=[2 0;0 1]
disp(d,'D=')
disp('Therefore, matrix A=PD(p^-1)')
s=inv(p)
disp(p*d*s)
disp('Hence, A^4=P(D^4)(P^-1)')
disp(p*(d^4)*s)