function clenshaw_curtis_test01 ( )

%% TEST01 simply calls CC_GRID once.
%
%  Licensing:
%
%    This code is distributed under the GNU LGPL license.
%
%  Modified:
%
%    12 October 2006
%
%  Author:
%
%    John Burkardt
%
  dim_num = 3;
  order_1d  = [ 3, 4, 2 ];

  fprintf ( 1, '\n' );
  fprintf ( 1, 'TEST01:\n' );
  fprintf ( 1, '  CC_GRID returns a grid of Clenshaw-Curtis points.\n' );
  fprintf ( 1, '  Here, we simply call for a specific grid.\n' );
  fprintf ( 1, '\n' );
  fprintf ( 1, '  Spatial dimension of grid = %d\n', dim_num );

  order_nd = prod ( order_1d(1:dim_num) );

  fprintf ( 1, '\n' );
  fprintf ( 1, '  Total number of points in the grid = %d\n', order_nd );
  fprintf ( 1, '\n' );

  grid_point = cc_grid ( dim_num, order_1d, order_nd );

  j = 1;
  q = sum ( order_1d(1:dim_num) );

  fprintf ( 1, '\n' );
  fprintf ( 1, '         I         Q          Grid orders:\n' );
  fprintf ( 1, '\n' );

  fprintf ( 1, '  %8d  %8d', j, q );
  for j = 1 : dim_num
    fprintf ( 1, '  %8d', order_1d(j) );
  end
  fprintf ( 1, '\n' );

  fprintf ( 1, '\n' );
  fprintf ( 1, '  Grid points:\n' );
  fprintf ( 1, '\n' );

  for j = 1 : order_nd
    fprintf ( 1, '  %8d', j );
    for dim = 1 : dim_num
      fprintf ( 1, '  %12f', grid_point(dim,j) );
    end
    fprintf ( 1, '\n' );
  end

  return
end