   function [x] = DepthSolv(L, b, list)
%% function [x] = DepthSolv(L, b, list)
n = size(L,1);
%%--------- Solve order according to list
x = b;
for k= 1:length(list)
      j = list(k);
      [ii, jj, vals] = find(L(:,j));
      x(ii) = x(ii) - x(j)*vals;
end
