
function x=solve_linear_LU(A,b)

[A] = fattLU(A);
L=tril(A)-diag(diag(A))+diag(ones(size(A,1),1)); 
U=triu(A);
d=diag(U);
detA_pm=prod(d);
if prod(d) ~=0
    c=b;
    y=L\c;
    x=U\y;
else
    fprintf('\n \t SINGULAR MATRIX');
end

