function newton_print ( c, z, title )

%*****************************************************************************80
%
%% newton_print prints a Newton polynomial.
%
%  Licensing:
%
%    This code is distributed under the GNU LGPL license.
%
%  Modified:
%
%    16 November 2019
%
%  Author:
%
%    John Burkardt
%
%  Input:
%
%    real c(n), the coefficients.
%
%    real z(n), the centers.
%
%    string TITLE, a title.
%
  n = length ( z );

  fprintf ( 1, '\n' );
  fprintf ( 1, '%s\n', title );
  fprintf ( 1, '\n' );
  fprintf ( 1, '  p(x) =                           %14f\n', c(1) );

  for i = 2 : n
    fprintf ( 1, '       + ( x - %14f) * ( %14f\n', z(i-1), c(i) );
  end

  fprintf ( 1, '  ' );
  for i = 1 : n - 1
    fprintf ( 1, ')' );
  end
  fprintf ( 1, '\n' );

  return
end
