function r = phi(alpha,gamma)

global nu sigma k kt L q h

r = (alpha*(1+alpha)*(1+gamma)^2)/(L + ((1+alpha)^2)*((1+gamma)^2));
