function blend_ij_0d1_test ( )

%*****************************************************************************80
%
%% blend_ij_0d1_test() tests blend_ij_0d1().
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    23 October 2008
%
%  Author:
%
%    John Burkardt
%
  m1 = 5;
  m2 = 4;

  fprintf ( 1, '\n' );
  fprintf ( 1, 'blend_ij_0d1_test():\n' );
  fprintf ( 1, '  BLEND_IJ_0D1 interpolates data in a table,\n' );
  fprintf ( 1, '  from corner data.\n' );
  fprintf ( 1, '\n' );
  fprintf ( 1, '  The table is %d rows by %d columns.\n', m1, m2 );

  x = zeros ( m1, m2 );
%
%  Load data in the corners only.
%
  i = 1;
  j = 1;
  r = (  i - 1 ) / ( m1 - 1 );
  s = (  j - 1 ) / ( m2 - 1 );
  x(i,j) = cubic_rs ( r, s, 1 );

  i = m1;
  j = 1;
  r = (  i - 1 ) / ( m1 - 1 );
  s = (  j - 1 ) / ( m2 - 1 );
  x(i,j) = cubic_rs ( r, s, 1 );

  i = 1;
  j = m2;
  r = (  i - 1 ) / ( m1 - 1 );
  s = (  j - 1 ) / ( m2 - 1 );
  x(i,j) = cubic_rs ( r, s, 1 );

  i = m1;
  j = m2;
  r = (  i - 1 ) / ( m1 - 1 );
  s = (  j - 1 ) / ( m2 - 1 );
  x(i,j) = cubic_rs ( r, s, 1 );

  x = blend_ij_0d1 ( x, m1, m2 );

  fprintf ( 1, '\n' );
  fprintf ( 1, '  Values interpolated by BLEND_IJ_0D1:\n' );
  fprintf ( 1, '\n' );

  for i = 1 : m1
    fprintf ( 1, '  %10f  %10f  %10f  %10f\n', x(i,1:m2) );
  end

  return
end
