function aniso_degree482_test ( )

%*****************************************************************************80
%
%% ANISO_DEGREE482_TEST compares methods of degree 4 in X, 8 in Y, 2 in Z.
%
%  Licensing:
%
%    This code is distributed under the GNU LGPL license.
%
%  Modified:
%
%    08 September 2017
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'ANISO_DEGREE482_TEST:\n' );
  fprintf ( 1, '  Plot index sets for various anisotropic methods\n' );
  fprintf ( 1, '  that reach degree 4 in X, 8 in Y and 2 in Z.\n' );
  fprintf ( 1, '\n' );

  m = 10;
  title_string = '';
%
%  Total degree <= 8.
%
  title_string = 'Total degree 2x+y+4z<=8';
  cube_display ( m, @(x,y,z)x+y+z<=-1, @(x,y,z)2*x+y+4*z<=8, title_string );
  filename = 'degree482_total.png';
  print ( '-dpng', filename );
  fprintf ( 1, '  Created "%s".\n', filename );
%
%  Maximum degree <= 8.
%
  title_string = 'Max degree (2x,y,4z)<=8';
  cube_display ( m, @(x,y,z)x+y+z<=-1, @(x,y,z)max([2*x,y,4*z])<=8, title_string );
  filename = 'degree482_max.png';
  print ( '-dpng', filename );
  fprintf ( 1, '  Created "%s".\n', filename );
%
%  Hyperbolic cross.
%
  title_string = 'Hyperbolic Cross (2x+1)(y+1)(4z+1)<=27';
  cube_display ( m, @(x,y,z)x+y+z<=-1, @(x,y,z)(2*x+1)*(y+1)*(4*z+1)<=27, title_string );
  filename = 'degree482_hyper.png';
  print ( '-dpng', filename );
  fprintf ( 1, '  Created "%s".\n', filename );
%
%  Smolyak: Log2(2*x) + Log2(y) + Log2(4*z) <= Log2(L)
%
  title_string = 'Log2(2*x) + Log2(y) + Log2(4*z) <= log2(8)';
  cube_display ( m, @(x,y,z)x+y+z<=-1, @(x,y,z)log2(max(2*x,1))+log2(max(y,1))+log2(max(4*z,1))<=log2(8.1), title_string );
  filename = 'degree482_smolyak.png';
  print ( '-dpng', filename );
  fprintf ( 1, '  Created "%s".\n', filename );

  return
end
