function humps_deriv2_test ( )

%*****************************************************************************80
%
%% humps_deriv2_test() tests humps_deriv2().
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    04 August 2019
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'humps_deriv2_test():\n' );
  fprintf ( 1, '  MATLAB version\n');
  fprintf ( 1, '  humps_deriv2() evaluates the humps second derivative.\n' );

  a = 0.0;
  b = 2.0;
  x = linspace ( a, b, 101 );
  y = humps_deriv2 ( x );

  clf ( );

  hold ( 'on' )

    plot ( x, y, 'r-', 'Linewidth', 2 );

    if ( a <= 0.0 && 0.0 <= b )
      plot ( [a,b], [0,0], 'k-', 'linewidth', 2 );
    end

    ymin = min ( y );
    ymax = max ( y );
    if ( ymin <= 0.0 && 0.0 <= ymax )
      plot ( [0, 0], [ymin,ymax], 'k-', 'linewidth', 2 );
    end

  hold ( 'off' );

  grid ( 'on' );
  xlabel ( '<--- X --->' );
  ylabel ( '<--- Y --->' );
  title ( 'y = humps"(x)' );
  pause ( 5 );
  filename = 'humps_deriv2.png';
  print ( '-dpng', filename );
  fprintf ( 1, '  Graphics saved as "%s"\n', filename );
%
%  Terminate.
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'humps_deriv_test:\n' );
  fprintf ( 1, '  Normal end of execution.\n' );

  return
end

