function r83v_jac_sl_test ( )

%*****************************************************************************80
%
%% R83V_JAC_SL_TEST tests R83V_JAC_SL.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    01 February 2016
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'R83V_JAC_SL_TEST\n' );
  fprintf ( 1, '  R83V_JAC_SL applies Jacobi iteration with an R83V matrix\n' );
  fprintf ( 1, '  to solve a linear system A*x=b.\n' );
%
%  Set A to the second difference matrix.
%
  n = 10;

  [ a, b, c ] = r83v_dif2 ( n, n );
%
%  Choose a random solution.
%
  x1 = r8vec_indicator1 ( n );
%
%  Compute the corresponding right hand side.
%
  ax = r83v_mv ( n, n, a, b, c, x1 );
%
%  Set the starting vector.
%
  x2 = zeros ( n, 1 );
%
%  Call the Jacobi routine.
%
  it_max = 25;

  for i = 1 : 3

    x2 = r83v_jac_sl ( n, a, b, c, ax, x2, it_max );

    r8vec_print ( n, x2, '  Current solution:' );

  end

  return
end
