function blend_rst_1dn_identity_test ( )

%*****************************************************************************80
%
%% BLEND_RST_1DN_IDENTITY_TEST checks for a gross error in the blend coefficients.
%
%  Licensing:
%
%    This code is distributed under the GNU LGPL license.
%
%  Modified:
%
%    25 October 2018
%
%  Author:
%
%    John Burkardt
%
  nmax = 3;

  fprintf ( 1, '\n' );
  fprintf ( 1, 'BLEND_RST_1DN_IDENTITY_TEST\n' );
  fprintf ( 1, '  Simple identity test to detect gross errors.\n' );
%
%  Set N to 3.
%
  n = 3;
%
%  Test BLEND_RST_1DN on identity.
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'Identity test for BLEND_RST_1DN:\n' );
  fprintf ( 1, '\n' );
  r = 0.0;
  s = 0.0;
  t = 0.0;
  x = blend_rst_1dn ( r, s, t, n, @identity_rst );
  fprintf ( 1, '  %8f  %8f  %8f  %8f  %8f  %8f\n', r, s, t, x(1:3) );

  r = 1.0;
  s = 0.0;
  t = 0.0;
  x = blend_rst_1dn ( r, s, t, n, @identity_rst );
  fprintf ( 1, '  %8f  %8f  %8f  %8f  %8f  %8f\n', r, s, t, x(1:3) );

  r = 0.0;
  s = 1.0;
  t = 0.0;
  x = blend_rst_1dn ( r, s, t, n, @identity_rst );
  fprintf ( 1, '  %8f  %8f  %8f  %8f  %8f  %8f\n', r, s, t, x(1:3) );

  r = 0.0;
  s = 0.0;
  t = 1.0;
  x = blend_rst_1dn ( r, s, t, n, @identity_rst );
  fprintf ( 1, '  %8f  %8f  %8f  %8f  %8f  %8f\n', r, s, t, x(1:3) );

  r = 1.0;
  s = 1.0;
  t = 1.0;
  x = blend_rst_1dn ( r, s, t, n, @identity_rst );
  fprintf ( 1, '  %8f  %8f  %8f  %8f  %8f  %8f\n', r, s, t, x(1:3) );

  r = 0.5;
  s = 0.5;
  t = 0.5;
  x = blend_rst_1dn ( r, s, t, n, @identity_rst );
  fprintf ( 1, '  %8f  %8f  %8f  %8f  %8f  %8f\n', r, s, t, x(1:3) );

  return
end
