function x = FwdSubRow(L,b)
%  x = FwdSubRow(L,b)
%  
% Inputs:
%  L   Lower triangular matrix
%  b   Right side vector
%  
%  Solves Lx = b by forward substitution (row-oriented version)
% Output:
%  x   solution vector
%  
[n,n] = size(L);
x = zeros(n,1);
x(1) = b(1)/L(1,1);
for k = 2:n
  x(k) = (b(k)-L(k,1:k-1)*x(1:k-1))/L(k,k);
end