% A1 = [2 -1 1 4; 4 -3 7 14; 0 -3 18 19; 6 -2 7 14];
% b1 = [0; 0; 2; 1];
% [x1,L1,U1] = lusolve(A1,b1,0);
% A1\b1;


% 
A2 = [2 -2 -4 4; 4 -3 -6 14; -8 3 6 19; 6 -2 7 14];
b2 = [2; 36; 92; 79];
[x2,L2,U2] = lusolve(A2,b2,1)
% A2\b2
% [L,U,P] = lu(A2)