Scribd
Upload
38 views

Lu Factorization Matlab Code:: If Else For

This MATLAB code demonstrates LU factorization using Doolittle and Crout's methods. It takes a 6x6 matrix A and factors it into lower and upper triangular matrices (L and U) using both methods. For Doolittle, L has ones along the diagonal and U is the upper triangular portion of A. For Crout's method, L is the lower triangular portion of A and U has ones along the diagonal.

Uploaded by

Engr. Naveed Mazhar
Copyright
© © All Rights Reserved
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
38 views

Lu Factorization Matlab Code:: If Else For

This MATLAB code demonstrates LU factorization using Doolittle and Crout's methods. It takes a 6x6 matrix A and factors it into lower and upper triangular matrices (L and U) using both methods. For Doolittle, L has ones along the diagonal and U is the upper triangular portion of A. For Crout's method, L is the lower triangular portion of A and U has ones along the diagonal.

Uploaded by

Engr. Naveed Mazhar
Copyright
© © All Rights Reserved
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 2

LU FACTORIZATION

MATLAB Code:
clc;clear all;
A=[ 1 2 3 4 0 1; 0 1 3 2 1 4;
1 2 -1 -2 4 7;7 1 3 8 1 5;
-1 -8 3 2 1 0;1 1 2 2 3 1];
[n,m]=size(A);
if (n~=m)
disp('A should be a square matrix.');
else
for i=1:n-1
A(i+1:n,i)=A(i+1:n,i)/A(i,i);
A(i+1:n,i+1:n)=A(i+1:n,i+1:n)-A(i+1:n,i)*A(i,i+1:n);
end
%%%% DOOLITTLE METHOD %%%%
L_doo=eye(n,n)+tril(A,-1);
U_doo=triu(A);
A,L_doo,U_doo,
%%%% CROUTS METHOD %%%%
L_cr=tril(A);
U_cr=eye(n,n)+triu(A,1);
L_cr,U_cr,
end

OUTPUT:
A=
1.0000 2.0000 3.0000 4.0000

0 1.0000

0 1.0000 3.0000 2.0000 1.0000 4.0000


1.0000

0 -4.0000 -6.0000 4.0000 6.0000

7.0000 -13.0000 -5.2500 -25.5000 35.0000 81.5000


-1.0000 -6.0000 -6.0000 0.7059 6.2941 3.4706
1.0000 -1.0000 -0.5000 0.1176 0.2991 -3.6262
L_doo =
1.0000

0 1.0000

1.0000

0 1.0000

7.0000 -13.0000 -5.2500 1.0000

-1.0000 -6.0000 -6.0000 0.7059 1.0000

0
0

1.0000 -1.0000 -0.5000 0.1176 0.2991 1.0000

U_doo =
1.0000 2.0000 3.0000 4.0000

0 1.0000

0 1.0000 3.0000 2.0000 1.0000 4.0000


0

0 -4.0000 -6.0000 4.0000 6.0000

0 -25.5000 35.0000 81.5000

0 6.2941 3.4706

0 -3.6262

L_cr =
1.0000

0 1.0000

1.0000

0 -4.0000

7.0000 -13.0000 -5.2500 -25.5000

-1.0000 -6.0000 -6.0000 0.7059 6.2941

0
0

1.0000 -1.0000 -0.5000 0.1176 0.2991 -3.6262


U_cr =
1.0000 2.0000 3.0000 4.0000

0 1.0000

0 1.0000 3.0000 2.0000 1.0000 4.0000


0

0 1.0000 -6.0000 4.0000 6.0000

0 1.0000 35.0000 81.5000

0 1.0000 3.4706

0 1.0000

You might also like