Faktorisasi LU Metode Doolittle dan Crout

Doolittle

mprintf("               Tugas Akhir Praktikum Metode Numeris                   \n")
mprintf("\n")
mprintf("\n")
mprintf("Program untuk mencari nilai x persamaan Ax=B dengan metode doolitle \n")
mprintf("Damar Rizki Ananda \n")
mprintf("17/413505/TK/45945 \n")
mprintf("Inputkan matriks A(3x3) dan matriks B(3x1) yang akan dicari nilai x nya \n")
  
function [x, y, u, l]=doolittle(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)
                if A(k,k)==0 then
        error("nilai x tidak terdefinisi karena pembaginya 0(nol)")
                end
        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
        mprintf("nilai x dari persamaan Ax=B adalah\n")
    endfunction

Crout

mprintf("               Tugas Akhir Praktikum Metode Numeris                   \n")
mprintf("\n")
mprintf("\n")
mprintf("Program untuk mencari nilai x persamaan Ax=B dengan metode crout \n")
mprintf("Damar Rizki ananda \n")
mprintf("17/413505/TK/45945 \n")
mprintf("Inputkan matriks A(3x3) dan matriks B(3x1) yang akan dicari nilai x nya \n")
  function [x, y, u, l]=crout(A, B)
    n=size(A,1)
    A(1,2:n)=A(1,2:n)/A(1,1)
    for j=2:n,
        for k=2:j,
            v=A(j,1:k-1)
            w=A(1:k-1,k)
            A(j,k)=A(j,k)-v*w
        end,
        v=A(j,1:j-1)
        for k=j+1:n,
            w=A(1:j-1,k)
                    if A(j,j)==0 then
        error("nilai x tidak terdefinisi karena pembaginya 0(nol)")
    end
            A(j,k)=(A(j,k)-v*w)/A(j,j)
         end
    end,
    l=tril(A)
    u=triu(A,1)+eye(n,n)
    y=l\B
    x=u\y
    mprintf("nilai x dari persamaan Ax=B adalah\n")
    endfunction