Motode Substitusi Mundur

Algoritma

Input   : Masukkan matriks Aij,i=1:n, j=1:m, nilai toleransi, 
Output : X(i), i=1:n
Langkah:
1.  jika |Ann|<eps, maka gagal
2.  X(n)=A(n,m)/A(n,n)
2.  untuk i=n-1:1
             S=0
            untuk j=i-1:n 
               S=S+Aij*Xj

                X(i)=(Ai,j+1-S)/Aii
4. jika |Aii| <eps berhenti
5. jika tidak berhenti ulangi lagi pada langkah 3

Algoritma pada MATLAP

A=input('Matriks A=');
eps=input('nilai toleransi=')
[n,m]=size(A);
X=zeros(n,1);

X(b)=A(b,m)/A(b,b)
if A(b,b)<eps
    disp('metode gagal');
    break
end
for i=n-1:-1:1
    S=0;
    for j=i+1:n   
      S=S+A(i,j)*X(j)
      X(i)=(A(i,j+1)-S)/A(i,i)
    end
    if A(i,i)<eps
        disp('metode gagal');
            break
    end
end


Tanks
 salim_lombe