clear,clc
A=[1 1 0 3;
   2 1 -1 1;
   3 -1 -1 2;
  -1 2 3 -1]

[L1,U1]=lu_decompose(A)

B=[8;7;14;-7]
    
V=inv(L1)*B

X=inv(U1)*V