
function x=solve_linear_LUP(A,b)

[L,U,P]=lu(A);
d=diag(U);
detA_pm=prod(d);
if prod(d) ~=0
    c=P*b;
    y=L\c;
    x=U\y;
else
    fprintf('\n \t SINGULAR MATRIX');
end

