function cvt_uniform_test ( )

%*****************************************************************************80
%
%% cvt_uniform_test(): Lebesgue analysis for the CVT points with initialization "U".
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    22 August 2016
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'CVT_UNIFORM_TEST:\n' );
  fprintf ( 1, '  Read sets of CVT nodes initialized "U" in [-1,+1],\n' );
  fprintf ( 1, '  computed by CHEBYSHEV_1D_LUMPING.\n' );
  fprintf ( 1, '  For a set of dimension N, estimate the Lebesgue constant L,\n' );
  fprintf ( 1, '  and the ratio of L to log(N+1).\n' );

  nfun = 10001;

  k = 0;

  filename = 'cvt_005_u.txt';
  x = load ( filename );
  n = length ( x );
  xfun = linspace ( -1.0, +1.0, nfun );
  l = lebesgue_constant ( n, x, nfun, xfun );
  k = k + 1;
  n_plot(k) = n;
  l_plot(k) = l;
  l2_plot(k) = l / log ( n + 1 );

  filename = 'cvt_009_u.txt';
  x = load ( filename );
  n = length ( x );
  xfun = linspace ( -1.0, +1.0, nfun );
  l = lebesgue_constant ( n, x, nfun, xfun );
  k = k + 1;
  n_plot(k) = n;
  l_plot(k) = l;
  l2_plot(k) = l / log ( n + 1 );

  filename = 'cvt_017_u.txt';
  x = load ( filename );
  n = length ( x );
  xfun = linspace ( -1.0, +1.0, nfun );
  l = lebesgue_constant ( n, x, nfun, xfun );
  k = k + 1;
  n_plot(k) = n;
  l_plot(k) = l;
  l2_plot(k) = l / log ( n + 1 );

  filename = 'cvt_033_u.txt';
  x = load ( filename );
  n = length ( x );
  xfun = linspace ( -1.0, +1.0, nfun );
  l = lebesgue_constant ( n, x, nfun, xfun );
  k = k + 1;
  n_plot(k) = n;
  l_plot(k) = l;
  l2_plot(k) = l / log ( n + 1 );

  filename = 'cvt_065_u.txt';
  x = load ( filename );
  n = length ( x );
  xfun = linspace ( -1.0, +1.0, nfun );
  l = lebesgue_constant ( n, x, nfun, xfun );
  k = k + 1;
  n_plot(k) = n;
  l_plot(k) = l;
  l2_plot(k) = l / log ( n + 1 );

  filename = 'cvt_129_u.txt';
  x = load ( filename );
  n = length ( x );
  xfun = linspace ( -1.0, +1.0, nfun );
  l = lebesgue_constant ( n, x, nfun, xfun );
  k = k + 1;
  n_plot(k) = n;
  l_plot(k) = l;
  l2_plot(k) = l / log ( n + 1 );

  filename = 'cvt_257_u.txt';
  x = load ( filename );
  n = length ( x );
  xfun = linspace ( -1.0, +1.0, nfun );
  l = lebesgue_constant ( n, x, nfun, xfun );
  k = k + 1;
  n_plot(k) = n;
  l_plot(k) = l;
  l2_plot(k) = l / log ( n + 1 );

  r8vec2_print ( n_plot, l_plot, '  L(N) = Lebesgue constants' )
  r8vec2_print ( n_plot, l2_plot, '  L2 = L(N)/log(N+1):' )

  return
end

