% A2 Matrix A W/ A 2 % Da1 Eigen Values For Matrix A W/ A 1 % vB2 Eigen Vector For Matrix B W/ A 2 % Ib1 Inverse of Matrix B W/ A 1
% A2 Matrix A W/ A 2 % Da1 Eigen Values For Matrix A W/ A 1 % vB2 Eigen Vector For Matrix B W/ A 2 % Ib1 Inverse of Matrix B W/ A 1
%
%
%
A2 = matrix A w/ a = 2
dA1 = eigen values for matrix A w/ a = 1
vB2 = eigen vector for matrix B w/ a = 2
iB1 = inverse of matrix B w/ a = 1
% Matrix A for a = 1
A1 = [0 1; -1 0];
dA1 = eig(A1);
[vA1,dA1] = eig(A1);
iA1 = A1^-1;
disp('eigen values for A(1)');
disp(dA1);
disp('eigen vectors for A(1)');
disp(vA1);
disp('inverse of A(1)');
disp(iA1);
% Matrix A for a = 2
A2 = [0 1; -2 0];
dA2 = eig(A2);
[vA2, dA2] = eig(A2);
iA2 = A2^-1;
disp('eigen values for A(2)');
disp(dA2);
disp('eigen vectors for A(2)');
disp(vA2);
disp('inverse of A(2)');
disp(iA2);
% Matrix B for a = 1
B1 = [0 1 0; 0 -2 -1; 2 1 0];
dB1 = eig(B1);
[vB1, dB1] = eig(B1);
iB1 = B1^-1;
disp('eigen values for B(1)');
disp(dB1);
disp('eigen vectors for B(1)');
disp(vB1);
disp('inverse of B(1)');
disp(iB1);
% Matrix B for a = 2
B2 = [0 1 0; 0 -2 -2; 2 1 0];
dB2 = eig(B2);
[vB2, dB2] = eig(B2);
iB2 = B2^-1;
disp('eigen values for B(2)');
disp(dB2);
disp('eigen vectors for B(2)');
disp(vB2);
disp('inverse of B(2)');
disp(iB2);
0.0000 + 0.0000i
0.0000 + 0.0000i
-2.0000 + 0.0000i
0.4472 + 0.0000i
-0.8944 + 0.0000i
0.0000 + 0.0000i
inverse of B(1)
-0.5000
0
1.0000
0
-2.0000
-1.0000
0.5000
0
0
0.0000 + 0.0000i
0.0000 + 0.0000i
-2.0000 + 0.0000i
0.4472 + 0.0000i
-0.8944 + 0.0000i
-0.0000 + 0.0000i
inverse of B(2)
-0.5000
0
1.0000
0
-1.0000
-0.5000
0.5000
0
0