function newton_to_power_test ( )

%*****************************************************************************80
%
%% newton_to_power_test tests newton_to_power.
%
%  Licensing:
%
%    This code is distributed under the GNU LGPL license.
%
%  Modified:
%
%    16 November 2019
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'newton_to_power_test\n' );
  fprintf ( 1, '  newton_to_power converts a Newton polynomial\n' );
  fprintf ( 1, '  to the standard power form\n' );

  r = 3.0;
  [ c, z ] = roots_to_newton ( r );
  newton_print ( c, z, '  The Newton polynomial:' );
  c = newton_to_power ( c, z );
  poly_print ( c, '  The power polynomial:' );

  r = [ 3.0, 1.0 ];
  [ c, z ] = roots_to_newton ( r );
  newton_print ( c, z, '  The Newton polynomial:' );
  c = newton_to_power ( c, z );
  poly_print ( c, '  The power polynomial:' );

  r = [ 3.0, 1.0, 2.0 ];
  [ c, z ] = roots_to_newton ( r );
  newton_print ( c, z, '  The Newton polynomial:' );
  c = newton_to_power ( c, z );
  poly_print ( c, '  The power polynomial:' );

  r = [ 3.0, 1.0, 2.0, 4.0 ];
  [ c, z ] = roots_to_newton ( r );
  newton_print ( c, z, '  The Newton polynomial:' );
  c = newton_to_power ( c, z );
  poly_print ( c, '  The power polynomial:' );

  return
end
