function area_estimate = area1 ( n )

%% AREA1 estimates area between X axis and Y=3/4*(1-X^2).

  x_min = -1.5;
  x_max = +1.5;

  x = linspace ( x_min, x_max, n );
  y = f ( x );

  area_estimate = ( x_max - x_min ) * sum ( y ) / n;

  return
end
function value = f ( x )

  value = 0.75 * ( 1.0 - x .^ 2 );

  return
end
