function lagrange_value_1d_test ( )

%*****************************************************************************80
%
%% lagrange_value_1d_test() tests lagrange_value_1d().
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    13 January 2021
%
%  Author:
%
%    John Burkardt
%
  nd = 4;
  ni = 21;

  xd = [   0.0;   2.0;  5.0;   10.0 ];
  yd = [ -28.0; +10.0; -8.0; +162.0 ];

  fprintf ( 1, '\n' );
  fprintf ( 1, 'lagrange_value_1d_test():\n' );
  fprintf ( 1, '  lagrange_value_1d() evaluates a Lagrange interpolant.\n' );

  x_min = 0.0;
  x_max = 10.0;
  xi = linspace ( x_min, x_max, ni );

  yi = lagrange_value_1d ( nd, xd, yd, ni, xi );

  r8vec2_print ( xi, yi, '  Table of interpolant values:' );

  return
end
