Program M-File dengan Iterasi Gauss Seidle

Berikut adalah program Matlab untuk menyelesaikan SPL dengan metode iterasi Gauss Seidle :

function [hasil]=gaussd(A,b,X0,tol)

L=triu(A',1)';

c=size(A);

D=eye(c).*A;

U=triu(A,1);

hasil=[0 X0'];

k=1;

X=inv(D+L)*(b-U*X0);

while max(abs(X-X0))>tol,

    hasil=[hasil;k X'];

    X0=X;

    k=k+1;

    X=inv(D+L)*(b-U*X0);

end


Cara memanggil programnya adalah dengan mengetik 

>> A=[10 -1 2 0; -1 11 -1 3; 2 -1 10 -1; 0 3 -1 8];

     b=[6;25;-11;15];

     X0=[0 0 0 0]';

>> hasil=gaussd(A,b,X0,0.001)

Program M-File dengan Iterasi Gauss Seidle 2014-10-05T18:10:00-07:00 Rating: 4.5 Diposkan Oleh: Unknown