blob: 020ed5ad98a4e819abfac23f4a4fa71129cf9a37 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
|
//Example No. 14_03
//Eigen Vectors
//Pg No. 473
clear ; close ; clc ;
A = [8 -4 ; 2 2 ] ;
lamd = poly(0,'lamd')
p = det(A - lamd*eye())
root = roots(p)
mprintf('\n The roots are \n lamda1 = %f \n lamda2 = %f \n ',root(1),root(2))
A1 = A - root(1)*eye()
X1 = [-1*A1(1,2)/A1(1,1) ; 1]
disp(X1,'X1 = ')
A2 = A - root(2)*eye()
X2 = [-1*A2(1,2)/A2(1,1) ; 1]
disp(X2,'X2 = ')
|
