Metode Numerik Eliminasi Gauss MATLAB Script

Dalam postingan kali ini, saya akan memberikan Script M-File untuk melakukan Eliminasi Gauss pada Sistem Persamaan Linear dalam perhitungan analisis numerik atau metode numerik. Berikut ini coding script yang dapat dijalankan pada program MATLAB dan script ini dapat dituliskan di Blank M-File MATLAB yang selanjutnya dapat dirun dengan nama GAUSS_ELIM:





function x = GAUSS_ELIM(A, b)
A=[2 0 0 0; 1 3 0 0; 3 4 2 0; -1 3 -6 -1];
b=[-10; 4; 2; 8];
n = size(A, 1);
r = zeros(n, 1);
for i = 1 : 1 : n
    r(i) = i;
end

x=zeros(n,1);
for k = 1 : 1 : n
    max = abs(A(r(k), r(k)));
    max_pos = k;   
    for l = k : 1 : n
        if abs(A(r(l), r(k))) > max;
            max=abs(A(r(l), r(k)));
            max_pos = l;               
        end
    end
    temp_r=r;
    r(k)=temp_r(max_pos);
    r(max_pos)=temp_r(k);
    for i = 1:1:n
        if i ~= k
            zeta = A(r(i),k)/A(r(k),k);
            for j = k : 1 : n
                A(r(i),j)=A(r(i),j)-A(r(k),j)*zeta;
            end
            b(r(i)) = b(r(i)) - b(r(k))*zeta;
        end
    end
end

for i = 1 : 1 : n
    x(i) = b(r(i))/A(r(i),i);
end

end

Komentar

Postingan Populer