function humps_newton_linspace ( )

  xdata = linspace ( 0.0, 2.0, 21 );
  ydata = humps ( xdata );
  [ c, z ] = data_to_newton ( xdata, ydata );
  
  xplot = linspace ( 0.0, 2.0, 101 );
  yplot = newton_value ( c, z, xplot );

  clf ( );
  hold ( 'on' );
  plot ( xplot, humps(xplot), 'b-', 'linewidth', 3 );
  plot ( xplot, yplot, 'r-', 'linewidth', 3 );
  plot ( xdata, ydata, 'k.', 'markersize', 25 );
  hold ( 'off' );
  legend ( { 'humps(x)', 'Interpolant', 'Data' } );
  grid ( 'on' );
  xlabel ( '<-- X -->' );
  ylabel ( '<-- Y(X) -->' );
  title ( 'Newton interpolant to humps() linspace data' );
  filename = 'humps_newton_linspace.png';
  print ( '-dpng', filename );
  fprintf ( 1, '  Graphics saved as "%s"\n', filename );

  return
end

