function area1_test ( )

%% area1_test() tests area1().
%
  fprintf ( 1, '\n' );
  fprintf ( 1, '  area1() estimates the area between the X axis\n' );
  fprintf ( 1, '  and the curve f(x)=3/4*(1-x^2) using N function values.\n' );

  area_exact = 9.0 / 16.0;

  fprintf ( 1, '\n' );
  fprintf ( 1, '     n    Area Estimate    Error\n' );
  fprintf ( 1, '\n' );
  n = 1;
  while ( n < 1024 )
    area_estimate = area1 ( n );
    area_error = abs ( area_estimate - area_exact );
    fprintf ( 1, '  %5d  %14.6g  %14.6g\n', n, area_estimate, area_error );
    n = n * 2;
  end
  fprintf ( 1, '  Exact  %14.6g\n', area_exact );

  return
end

