function lebesgue_equidistant3_test ( )

%*****************************************************************************80
%
%% LEBESGUE_EQUIDISTANT3_TEST looks at Equidistant3 points.
%
%  Licensing:
%
%    This code is distributed under the GNU LGPL license.
%
%  Modified:
%
%    17 January 2015
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'LEBESGUE_EQUIDISTANT3_TEST:\n' );
  fprintf ( 1, '  Analyze Equidistant3 points.\n' );

  n_max = 11;
  l = zeros ( n_max, 1 );
  xfun = linspace ( -1.0, +1.0, 501 );
  xfun = xfun';

  for n = 1 : n_max
    x = equidistant3 ( n );
    l(n) = lebesgue_constant ( n, x, xfun );
  end

  r8vec_print ( n_max, l, '  Equidistant3 Lebesgue constants for N = 1 to 11:' )
%
%  Examine one case more closely.
%
  n = 11;
  x = equidistant3 ( n );
  r8vec_print ( n, x, '  Equidistant3 points for N = 11' );

  label = 'Equidistant3 points for N = 11';
  filename = 'equidistant3.png';
  lebesgue_plot ( n, x, xfun, label, filename );
  fprintf ( 1, '\n' );
  fprintf ( 1, '  Plot file saved as "%s"\n', filename );

  return
end

