Doolitle Algorithm

Hello, di postingan kali ini akan dibahas mengenai algoritma doolitle. Algoritma ini digunakan untuk mencari nilai x dari Ax=b.

function [l, u, x]=doolitle(a, b)

   n=size(a,1)

   a(2:n,1)=a(2:n,1)/a(1,1)

   for k=2:n-1,

       v=a(k,1:k-1)

       for j=k:n,

           w=a(1:k-1,j)

           a(k,j)=a(k,j)-v*w

       end

       w=a(1:k-1,k)

       for j=k+1:n,

           v=a(j,1:k-1)

           a(j,k)=(a(j,k)-v*w)/a(k,k)

       end

   end

   a(n,n)=a(n,n)-a(n,1:n-1)*a(1:n-1,n)

   l=tril(a,-1)+eye(n,n)

   u=triu(a)

   y=l\b

   x=u\y

endfunction