function [ fprime ] = fprime(alpha,phi)
%fprime(alpha,phi)=eqn 3.63 in your text
%Assuming alpha is constant, i.e., 40 degrees

fprime=(5/2)*sind(phi)-sind(alpha-phi);

end

