function x = r8mat_cholesky_solve ( n, l, b )

%*****************************************************************************80
%
%% R8MAT_CHOLESKY_SOLVE solves a Cholesky factored linear system A * x = b.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    30 November 2017
%
%  Author:
%
%    John Burkardt
%
%  Input:
%
%    integer N, the number of rows and columns of the matrix A.
%
%    real L(N,N), the N by N Cholesky factor of the
%    system matrix A.
%
%    real B(N,1), the right hand side of the linear system.
%
%  Output:
%
%    real X(N,1), the solution of the linear system.
%

%
%  Destroy all row vectors!
%
  b = b(:);
%
%  Solve L * y = b.
%
  x = r8mat_l_solve ( n, l, b );
%
%  Solve L' * x = y.
%
  x = r8mat_lt_solve ( n, l, x );

  return
end
