a=[ 3 2 1
      2 0 1
      -1 1 0];
b=[6 3 0]';

[l,u,x,my_det] = my_gauss(3,a,b);
% [l,u,x,my_det] = my_gaussb(3,a,b);