function blend_r_0dn_test ( )

%*****************************************************************************80
%
%% blend_r_0dn_test() tests blend_r_0dn().
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    23 October 2008
%
%  Author:
%
%    John Burkardt
%
  m = 11;
  n = 2;

  fprintf ( 1, '\n' );
  fprintf ( 1, 'blend_r_0dn_test():\n' );
  fprintf ( 1, '  BLEND_R_0DN interpolates endpoint vector data\n' );
  fprintf ( 1, '  into a list.\n' );
  fprintf ( 1, '\n' );
  fprintf ( 1, '  The list is %d rows long.\n', m );
  fprintf ( 1, '  Each vector has length %d\n', n );

  x = zeros ( m, n );

  for i = 1 : m
    r = (  i - 1 ) / ( m - 1 );
    x(i,1:n) = blend_r_0dn ( r, n, @powers_r );
  end

  r8mat_print ( m, n, x, '  Data blended by BLEND_R_0DN:' );
%
%  Load all data.
%
  for i = 1 : m
    r = (  i - 1 ) / ( m - 1 );
    for j = 1 : n
      x(i,j) = powers_r ( r, j );
    end
  end

  r8mat_print ( m, n, x, '  Exact data:' );

  return
end
