clear; % untuk membersihkan peubah
u = [4 2 1 0 9 1 0 5] % Vektor baris u dari No. Stb.
pause
v = u' % Vektor kolom v transpose dari u
A = u*v % Matrix A = uv [1X1]
B = v*u % Matrix B = vu [8X8]
pause
det_B = det(B) % Determinan B = 0
k = -1:0.01:-0.25 % Peubah k beubah dari -1 sampai -0.25 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
    D = inv(C) % Matrix D adalah inverse dari Matrix C
    det_D(i) = det(D) % Determinan Matriks D
end
pause
j = 0.25:0.01:1 % Peubah j beubah dari 0.25 sampai 1 dengan langkah 0.01
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
pause
plot(k,det_D,'md'), grid on, xlabel('k'), ylabel('Determinan(D)')