function r8cb_det_test ( )

%*****************************************************************************80
%
%% r8cb_det_test() tests r8cb_det().
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    25 July 2016
%
%  Author:
%
%    John Burkardt
%
  m = 5;
  n = 5;
  ml = 1;
  mu = 1;
  seed = 123456789;

  fprintf ( 1, '\n' );
  fprintf ( 1, 'r8cb_det_test()\n' );
  fprintf ( 1, '  r8cb_det() computes the determinant of a matrix\n' );
  fprintf ( 1, '  that has been factored by R8CB_NP_FA.\n' );
  fprintf ( 1, '\n' );
  fprintf ( 1, '  Matrix order N     = %d\n', n );
  fprintf ( 1, '  Lower bandwidth ML = %d\n', ml );
  fprintf ( 1, '  Upper bandwidth MU = %d\n', mu );
%
%  Set the matrix.
%
  a = r8cb_dif2 ( m, n );

  r8cb_print ( m, n, ml, mu, a, '  The compact band matrix:' );
%
%  Copy the matrix into a general array.
%
  a2 = r8cb_to_r8ge ( m, n, ml, mu, a );
%
%  Factor the matrix.
%
  [ a_lu, info ] = r8cb_np_fa ( n, ml, mu, a );
%
%  Compute the determinant.
%
  determ = r8cb_det ( n, ml, mu, a_lu );

  fprintf ( 1, '\n' );
  fprintf ( 1, '  R8CB_DET computes the determinant = %f\n', determ );
%
%  Factor the general matrix.
%
  [ a2_lu, pivot, info ] = r8ge_fa ( n, a2 );
%
%  Compute the determinant.
%
  determ = r8ge_det ( n, a2_lu, pivot );

  fprintf ( 1, '  R8GE_DET computes the determinant = %f\n', determ );
%
%  Compute the determinant.
%
  determ = det ( a2 );

  fprintf ( 1, '  DET computes the determinant      = %f\n', determ );

  return
end
