function [ f ] = f(alpha,phi)
%f(alpha,phi)=eqn 3.5 in your text

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

end

