function dif_antideriv_test ( )

%*****************************************************************************80
%
%% dif_antideriv_test() tests dif_antideriv();
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    28 October 2019
%
%  Author:
%
%    John Burkardt
%
  maxtab = 10;

  fprintf ( 1, '\n' );
  fprintf ( 1, 'dif_antideriv_test()\n' );
  fprintf ( 1, '  dif_antideriv() computes the difference form of the antiderivative;\n' );
  fprintf ( 1, '\n' );
%
%  Set XTAB, YTAB to X, X^3 + 2x^2 + 3x + 4.
%
  ntab = 4;
  xtab = 0 : 3;
  ytab(1:ntab) = xtab(1:ntab).^3 + 2.0 * xtab(1:ntab).^2 + 3.0 * xtab(1:ntab) + 4.0;

  ydif = data_to_dif_display ( ntab, xtab, ytab );
  dif_print ( ntab, xtab, ydif, '  The divided difference polynomial:' );
%
%  Compute a table for the antiderivative.
%
  [ ntab2, xtab2, ydif2 ] = dif_antideriv ( ntab, xtab, ydif );

  dif_print ( ntab2, xtab2, ydif2, '  The antiderivative:' );

  return
end
