function newton_coef_1d_test ( )

%*****************************************************************************80
%
%% newton_coef_1d_test() tests newton_coef_1d().
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    04 July 2015
%
%  Author:
%
%    John Burkardt
%
  nd = 5;
  xd = [ 0.0; 1.0; 2.0; 3.0; 4.0 ];
  yd = [ 24.0; 0.0; 0.0; 0.0; 0.0 ];

  fprintf ( 1, '\n' );
  fprintf ( 1, 'newton_coef_1d_test():\n' );
  fprintf ( 1, '  newton_coef_1d() sets the coefficients for a 1D Newton interpolation.\n' );

  r8vec2_print ( xd, yd, '  Interpolation data:' );

  cd = newton_coef_1d ( nd, xd, yd );

  r8vec_print ( nd, cd, '  Newton interpolant coefficients:' );

  return
end
