Faktorisasi LU dengan metode doolittle

Algoritma faktorisasi LU dengan metode doolittle
  1. Tentukan jumlah baris dari matriks A
  2. Untuk k=1,...,n


  3. Untuk j=2,...,n

  4. Untuk k=2,...,n
    1. Untuk j=2,...,k
    2. Untuk j=k+1,...,n


Source program MATLABnya adalah seperti berikut
function [l, u]=doolittle(a)
    n=size(a,1)
    for k=1:n,
        u(1,k)=a(1,k);
        m(k,k)=1;
    end
    for j=2:n,
        m(j,1)=a(j,1)/u(1,1);
    end
    for k=2:n,
        for j=2:k,
            u(j,k)=a(j,k)-m(j,1:j-1)*u(1:j-1,k);
        end
        for j=k+1:n,
            m(j,k)=(a(j,k)-m(j,1:k-1)*u(1:k-1,k))/u(k,k);
        end
    end
endfunction
Share on Google Plus

About Izzuddin Fathin Azhar

A young man with awesome spirit.
    Blogger Comment

1 komentar: