function polyinterp_test ( )

%*****************************************************************************80
%
%% polyinterp_test tests polyinterp.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    07 January 2020
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'polyinterp_test\n' );
  fprintf ( 1, '  Test polyinterp.\n' );

  x = 1:6;
  y = [ 16, 18, 21, 17, 15, 12 ];
  disp ( [ x; y ] )
  u = 0.75 : 0.05 : 6.25;
  v = polyinterp ( x, y, u );
  plot ( x, y, 'o', u, v, 'r-' );
  grid ( 'on' );
  title ( 'Data interpolation by polyinterp' );
  filename = 'polyinterp_test.png';
  print ( '-dpng', filename );
  fprintf ( 1, '  Graphics saved as "%s"\n', filename );

  return
end

