function f = func_eq3( theta )
%
%   Evaluate the function from equation (3') for theta in degrees

th_rad = pi*theta/180.0;  %  convert angle to radians

f = 250*cos(th_rad)*( sin(th_rad) + sqrt( sin(th_rad)^2 - 0.08 ) ) - 200;
