function r8ncf_dif2_test ( )

%*****************************************************************************80
%
%% R8NCF_DIF2_TEST tests R8NCF_DIF2.
%
%  Licensing:
%
%    This code is distributed under the MIT license. 
%
%  Modified:
%
%    20 July 2016
%
%  Author:
%
%   John Burkardt
%
  m = 7;
  n = 5;

  fprintf ( 1, '\n' );
  fprintf ( 1, 'R8NCF_DIF2_TEST\n' );
  fprintf ( 1, '  R8NCF_DIF2 sets up a R8NCF second difference matrix;\n' );
  fprintf ( 1, '\n' );
  fprintf ( 1, '  Matrix rows M =    %d\n', m );
  fprintf ( 1, '  Matrix columns N = %d\n', n );

  nz_num = r8ncf_dif2_nz_num ( m, n );

  fprintf ( 1, '  Matrix nonzeros =  %d\n', nz_num );

  rowcol = r8ncf_dif2_rowcol ( m, n, nz_num );

  a = r8ncf_dif2 ( m, n, nz_num, rowcol );

  r8ncf_print ( m, n, nz_num, rowcol, a, '  The R8NCF matrix:' );

  return
end
