% Script to test LU_Solve_gen.m
A=[  1 -1  2
    -2  1  1
    -1  2  1];
[n n]=size(A);
[L,U]=LU_factor(A,n);
L, U,
B=eye(n);
X=LU_Solve_Gen(L,U,B);
display(X);
display(X*A);