A = [4 3 -1; -2 -4 5; 1 2 6];
b = [1; 6; 14;];
[L,U]=lu(A);
disp("U=");
disp(U);
disp("L=");
disp(L);
y = L\b;
disp("y=");
disp(y);
x = U\y;
disp("x=");
disp(x);