function q = rectangle_integral2 ( )

%*****************************************************************************80
%
%% rectangle_integral2 estimates an integral over a rectangle using integral2().
%
%  Output:
%
%    real q: the integral estimate.
%
  a = 1.0;
  b = 5.0;
  c = 2.0;
  d = 3.0;

  q = integral2 ( @(x,y)sinexp(x,y), a, b, c, d );
  
  fprintf ( 1, '  Using integral2(), integral estimate is %g\n', q );

  return
end

