function output = struveh1approx(x)
% Struve H1 approximation obtained from paper by Aarts and Janssen (JASA,
% 2003)
%
% Robert J. McGough
% Michigan State University
% 18 Dec 2005

output = 2/pi - besselj(0, x) + (16/pi - 5) * sin(x) ./ x + (12 - 36/pi) .* (1 - cos(x)) ./ x .^2;