LU decomposition using Doolittle's Method with MATLAB code

The contents of this video lecture are:
📜Contents 📜
📌 (0:03​​​​​​) LU Decomposition
📌 (2:55​​​) Doolittle's Method
📌 (7:47​​​) MATLAB code of Doolittle's Method
📌 (29:38) MATLAB code of Forward substitution
📌 (33:50​​​) MATLAB code of Back substitution

#LUdecomposition
#doolittlesmethod
#forwardsubstitution
#backsubstitution
#linearsystems​​​
#numericalanalysis​​​​​​
#numericalcomputation​

[usernameisamyth]

code:



clc;clear all;

A=input('Enter the coefficient matrix: ');
b=input('Enter the vector of constants: ');

N=length(A);
L=zeros(N, N);
U=zeros(N, N);

for i=1:N
L(i, i)=1;
end

U(1, :)=A(1, :);
L(:, 1)=A(:, 1)/U(1, 1);

for i=2:N
for j=i:N
U(i, j)=A(i, j)-L(i, 1:i-1)*U(1:i-1, j);
end
for k=i+1:N
L(k, i)=(A(k, i)-L(k, 1:i-1)*U(1:i-1, i))/U(i, i);
end
end
L, U
y=L\b
x=U\y

[eedeatlikkle]

You are a king! Best explanations I have ever come across!

I was wondering if you will be covering initial value problems and boundary value problems using matlab?
Since you explain concepts very well, it would be very helpful.

[m.eminerdogan]

firstly thx for video but ı have problem with NAN and İnf, how can ı fix that ?

[mrigankchhabra]

Greetings Mr.Iqbal, I am watching your video lectures of coding matrix problems, and I am finding it incredibly useful. But have a small query too, may you please explain it again the reason behind of taking zeros, Is this a function and what does it actually signify? Regards

[yonathanbomisso]

Puis-je avoir le pdf de cette méthode ?

[siddhanthmadhavan]

sir, i am getting error like "Index in position 1 exceeds array bounds (must not exceed 1)" in line number 13 in my program when i typed the same like u have typed sir. pls help in resolving my error asap.

[omarsoliman]

Unrecognized function or variable 'B'.

Error in LU_Decomposition (line 20)
Y(1) = B(1)/L(1, 1);

[tirthakgayen]

can you make a video for cholesky decomposition, please

[daniyal]

this is a good method to solve lud composition, but there is a built in function in matlab to solve lud coposition.
clear;
clc;

A = [2 4 -6; 1 5 3; 1 3 2];
b = [-4; 10; 5];

[L, U] = lu(A);
y = L\b;
x = U\y;

disp("Solution:");
disp(x);