summaryrefslogtreecommitdiff
path: root/50/CH3/EX3.27/ex_3_27.sce
blob: 242a16db7a7d9fbc1fc756ffc1a60df5f7339290 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
// example 3.27
// a) find eigenvalue and eigen vector;
// b) verify inv(S)*A*S is a diagonal matrix;

// 1) 
A=[1 2 -2 ;1 1 1;1 3 -1];

B=[1 0 0;0 1 0; 0 0 1];

 [x,lam] = geigenvectors(A,B);
 
 inv(x)*A*x
 
 // 2) 
A=[3 2 2;2 5 2;2 2 3];

B=[1 0 0;0 1 0; 0 0 1];


 [x,lam] = geigenvectors(A,B);
 
 inv(x)*A*x