%A=rand(4,4);
clc
A=[1   3   5;2   4   7;1   1   0];
b=[1 4 -4];
[L U]=LU_DooLittle(A);
y=ForwardSub(L,b);
x=BackwardSub(U,y);