A=[1,1,3,0;2,1,2,-1;3,6,4,1;1,-2,1,3];
b=[5;4;14;3];
% solve with backslash
% x=A\b
% solve by LU factorization
[L,U,P]=lu(A)
y=L\(P*b);
x=U\y
