function ellipsoid_volume_test ( )

%*****************************************************************************80
%
%% ellipsoid_volume_test() tests ellipsoid_volume().
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    06 January 2023
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'ellipsoid_volume_test():\n' );
  fprintf ( 1, '  ellipsoid_volume_test() computes the volume of the ellipsoid\n' );
  fprintf ( 1, '    (x/a)^2+(y/b)^2+(z/c)^2=1\n' );

  fprintf ( 1, '\n' );
  fprintf ( 1, '           A           B           C        volume\n' );
  fprintf ( 1, '\n' );

  a = 1.0;
  b = 0.8;
  c = 0.625;
  volume = ellipsoid_volume ( a, b, c );
  fprintf ( 1, '  %10.4g  %10.4g  %10.4g  %10.4g\n', a, b, c, volume );

  a = 1.0;
  b = 1.0;
  c = 0.5;
  volume = ellipsoid_volume ( a, b, c );
  fprintf ( 1, '  %10.4g  %10.4g  %10.4g  %10.4g\n', a, b, c, volume );

  a = 1.0;
  b = 1.0;
  c = 1.0;
  volume = ellipsoid_volume ( a, b, c );
  fprintf ( 1, '  %10.4g  %10.4g  %10.4g  %10.4g\n', a, b, c, volume );

  a = 2.0;
  b = 1.0;
  c = 0.25;
  volume = ellipsoid_volume ( a, b, c );
  fprintf ( 1, '  %10.4g  %10.4g  %10.4g  %10.4g\n', a, b, c, volume );

  a = 2.0;
  b = 3.0;
  c = 4.0;
  volume = ellipsoid_volume ( a, b, c );
  fprintf ( 1, '  %10.4g  %10.4g  %10.4g  %10.4g\n', a, b, c, volume );

  return
end

