clear;
% PRAKTIKUM ALJABAR LINIER 2011 (TAMBAHAN)
'PRAKTIKUM ALJABAR LINIER 2011 (TAMBAHAN)' % Text yang akan ditampilkan di "Command Window"
'Hubungan nilai "k" dengan determinan matrix'
% Meng-INPUT DATA
A = input('Ketik No. Stb. anda tanpa huruf D dalam bentuk vektor baris ')
pause
% Men-TRANSPOSE MATRIX
'B adalah A-transpose'
B = A'
pause
% Mengalikan MATRIX
'C adalah A*B'
C = A*B
pause
'D adalah B*A'
D = B*A
pause
% Determinan MATRIX
'Menghitung Determinan Matrix D untuk k = 0'
detD0 = det(D)
pause
'Membentuk Matrix Identitas I  [8X8]'
I = eye(8)
pause
'Menghitung Determinan Matrix E = D + k*I untuk k = -1'
k = -1
E = D + k*I
det_E = det(E)
'Determinan E tidak lagi nol seperti D'
pause
'Menghitung Determinan Matrix F = D + k*I untuk k = 1'
k = 1
F = D + k*I
det_F = det(F)
'Determinan F tidak lagi nol seperti D'
pause
'Sekarang bagaimana kalau k berubah dari -1 s/d 1 ???'
k = [-1:0.01:1];
'k berubah dari -1 sampai 1 dengan step 0.01'
'k = [-1:0.01:1]'
'N adalah panjangnya peubah k'
N = length(k)
'Matrix G = D + k*I dihitung determinannya untuk -1 <= k <= 1'
for i = 1:N
    G = D + k(i)*I; 
    clf;
    plot(G), grid on, xlabel('matrix G = D + k*I')
    det_G(i) = det(G);
end
'LIHAT WINDOWS Figure sebelom menekan ENTER!!!'
pause
plot(k,det_G, '*'), grid on, xlabel('k'), ylabel('det matrix G = D + k*I')
pause
'Sekarang bagaimana dengan inverse dari G ?'
k1 = [-1:0.01:-0.2]
'k1 berubah dari -1 sampai -0.2 dengan step 0.01'
'k1 = [-1:0.01:-0.2]'
'N1 adalah panjangnya peubah k1'
N1 = length(k1)
'Matrix H = inv(G) dihitung determinannya untuk -1 <= k1 <= -0.2'
for i = 1:N1
    G = D + k1(i)*I;
    H = inv(G);
    clf;
    plot(H), grid on, xlabel('matrix H = inv(G)')
    det_H(i) = det(H);
    kH(i) = k1(i);
end
'LIHAT WINDOWS Figure sebelom menekan ENTER!!!'
pause
k2 = [0.2:0.01:1]
'k2 berubah dari 0.2 sampai 1 dengan step 0.01'
'k2 = [0.2:0.01:1]'
'N2 adalah panjangnya peubah k2'
N2 = length(k2)
'Matrix H = inv(G) dihitung determinannya untuk 0.2 <= k2 <= 1'
for i = 1:N2
    G = D + k2(i)*I;
    H = inv(G);
    clf;
    plot(H), grid on, xlabel('matrix H = inv(G)')
    det_H(N1+i) = det(H);
    kH(N1+i) = k2(i);
end
'LIHAT WINDOWS Figure sebelom menekan ENTER!!!'
pause
plot(kH,det_H,'*'), grid on, xlabel('k'), ylabel('det matrix H = inv(G)')
pause

