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]
det_B = det(B); % Determinan B = 0
k = -0.25:0.005:-0.05; % Peubah k beubah dari -0.25 sampai -0.05 dengan langkah 0.005
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
    D = inv(C); % Matrix D adalah inverse dari Matrix C
    det_D(i) = det(D); % Determinan Matriks D
end
j = 0.05:0.005:0.25; % Peubah j beubah dari 0.05 sampai 0.25 dengan langkah 0.005
M= length(j); % Panjangnya peubah j
I = eye(8); % Matriks Identitas [8X8]
for i = 1:M, % Untuk i dari 1 sampai N
    C = B + (j(i)*I); % Matriks C = B + kI
    D = inv(C); % Matrix D adalah inverse dari Matrix C
    det_D(N+i) = det(D); % Determinan Matriks D
    k(N+i) = j(i); % untuk nilai k yang positif
end
plot(k,det_D,'md'), grid on, xlabel('k'), ylabel('Determinan(D)')