function clenshaw_curtis_test ( )

%*****************************************************************************80
%
%% CLENSHAW_CURTIS_TEST plots Smolyak Clenshaw Curtis index sets.
%
%  Licensing:
%
%    This code is distributed under the GNU LGPL license.
%
%  Modified:
%
%    09 September 2017
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'CLENSHAW_CURTIS_TEST:\n' );
  fprintf ( 1, '  Plot Smolyak Clenshaw Curtis monomial indices.\n' );
  fprintf ( 1, '\n' );

  m = 20;
  title_string = 'Log2(X-1)+Log2(Y-1)+Log2(Z-1) <= 1';
  cube_display ( m, @cc_0, @cc_1, title_string );
  print ( '-dpng', 'cc1.png' );
  fprintf ( 1, '  Created "cc1.png".\n' );

  m = 20;
  title_string = 'Log2(X-1)+Log2(Y-1)+Log2(Z-1) <= 2';
  cube_display ( m, @cc_1, @cc_2, title_string );
  print ( '-dpng', 'cc2.png' );
  fprintf ( 1, '  Created "cc2.png".\n' );

  m = 20;
  title_string = 'Log2(X-1)+Log2(Y-1)+Log2(Z-1) <= 3';
  cube_display ( m, @cc_2, @cc_3, title_string );
  print ( '-dpng', 'cc3.png' );
  fprintf ( 1, '  Created "cc3.png".\n' );

  m = 20;
  title_string = 'Log2(X-1)+Log2(Y-1)+Log2(Z-1) <= 4';
  cube_display ( m, @cc_3, @cc_4, title_string );
  print ( '-dpng', 'cc4.png' );
  fprintf ( 1, '  Created "cc4.png".\n' );

  m = 20;
  title_string = 'Log2(X-1)+Log2(Y-1)+Log2(Z-1) <= 5';
  cube_display ( m, @cc_4, @cc_5, title_string );
  print ( '-dpng', 'cc5.png' );
  fprintf ( 1, '  Created "cc5.png".\n' );

  m = 20;
  title_string = 'Log2(X-1)+Log2(Y-1)+Log2(Z-1) <= 6';
  cube_display ( m, @cc_5, @cc_6, title_string );
  print ( '-dpng', 'cc6.png' );
  fprintf ( 1, '  Created "cc6.png".\n' );

  m = 20;
  title_string = 'Log2(X-1)+Log2(Y-1)+Log2(Z-1) <= 7';
  cube_display ( m, @cc_6, @cc_7, title_string );
  print ( '-dpng', 'cc7.png' );
  fprintf ( 1, '  Created "cc7.png".\n' );

  return
end
