function [ gbad ] = gbad(alpha,phi)
%gbad(alpha,phi)=eqn 3.32 (with 3.31) in your text

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

end

