clear; % untuk membersihkan peubah
u = [4 2 1 0 9 1 0 5]; % Vektor baris u dari No. Stb.
v = u'; % Vektor kolom v transpose dari u
A = u*v; % Matrix A = uv [1X1]
B = v*u; % Matrix B = vu [8X8]
lambda_B = eig(B); % Nilai Eigen matrix B
k = -10:0.01:10; % Peubah k beubah dari -10 sampai 10 dengan langkah 0.01
N = length(k); % Panjangnya peubah k
I = eye(8); % Matriks Identitas [8X8]
for i = 1:N, % Untuk i dari 1 sampai N
    C = B + (k(i)*I); % Matriks C = B + kI
    eig_C = eig(C); % Nilai Eigen Matrix C
    M = length(eig_C); % Panjangnya eig_C, ada 8
    for j = 1:M,
        lambda_C(i,j)=eig_C(j);
    end
end
plot(k,lambda_C,'.'), grid on, xlabel('k'), ylabel('Nilai Eigen Matrix C')
pause
for i = 1:N, % Untuk i dari 1 sampai N
    K = [k(i) 0 0 0 0 0 0 10*k(i)]; % Matrix K diubah-ubah sesuai dengan kreasi-ta'
    E = B + (v*K); % Matriks E = B + v*K
    eig_E = eig(E); % Nilai Eigen Matrix C
    M = length(eig_E); % Panjangnya eig_C, ada 8
    for j = 1:M,
        lambda_E(i,j)=eig_E(j);
    end
end
plot(lambda_E,'d'), grid on, xlabel('k'), ylabel('Nilai Eigen Matrix E')

