function [ g ] = g(alpha,phi)
%g(alpha,phi)=eqn 3.41 (with 3.40) in your text

g=acosd((2/5)*((5/3)*cosd(alpha)+11/6-cosd(alpha-phi)));

end

