function r83v_cr_sls_test ( )

%*****************************************************************************80
%
%% R83V_CR_SLS_TEST tests R83V_CR_SLS.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    01 February 2016
%
%  Author:
%
%    John Burkardt
%
  n = 5;
  debug = 0;

  fprintf ( 1, '\n' );
  fprintf ( 1, 'R83V_CR_SLS_TEST\n' );
  fprintf ( 1, '  R83V_CR_SLS solves multiple linear systems A*x1:xn=b1:bn\n' );
  fprintf ( 1, '  after R83V_CR_FA has factored it.\n' );
  fprintf ( 1, '\n' );
  fprintf ( 1, '  Matrix order N = %d\n', n );
  fprintf ( 1, '  Demonstrate multiple system solution method.\n' );
  fprintf ( 1, '\n' );
%
%  Set the matrix values.
%
  [ a, b, c ] = r83v_dif2 ( n, n );
%
%  Print the matrix.
%
  r83v_print ( n, n, a, b, c, '  Input matrix A:' );
%
%  Factor the matrix once.
%
  a_cr = r83v_cr_fa ( n, a, b, c );
%
%  Set the number of right hand sides.
%
  nb = 2;
  ax = zeros ( n, nb );

  ax(1:n-1,1) = 0.0E+00;
  ax(n,1) = n + 1;

  ax(1,2) = 1.0E+00;
  ax(2:n-1,2) = 0.0E+00;
  ax(n,2) = 1.0E+00;

  r8ge_print ( n, 2, ax, '  Right hand sides b1:b2' );
%
%  Solve the linear systems.
%
  x = r83v_cr_sls ( n, a_cr, nb, ax );

  r8ge_print ( n, 2, x, '  Solutions x1:x2' );

  return
end

