function [ x ] = LUsolve( L,U,b )
%LUsolve: solve LUx=b where (L,U) are (lower,upper) triag
  y = Lsolve(L,b);
  x = Usolve(U,y);
end
