function x=LUsolve(A, b)
% Using LU decomposition solve the system of Ax = b.
% templete for the exercise.

[m, n]= size(A);
if(m == n) 
% Step 1
  [L, U] = lualt(A);
% Initializing
   y = zeros(n, 1);
   x = zeros(n, 1);
   sum = 0;

% Fill out Step 2 to Step 5

end
% return x
x
