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) k(i) k(i) k(i)k(i 
    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')