  function ig = antennaint (theta)
% function ig = antennaint (theta)
% Integrand of numerical integration on page 671 in the textbook of Rao.
% Input: integration variable, angle theta.
% 
% To compute the integral equal to 0.609 on p. 671, write
% quad ('antennaint', 0, pi/2)

if abs(sin(theta)) < sqrt(eps)  ig = 0;
else ig  =  cos((pi./2).*cos(theta)).^2./sin(theta);  end
return

