function [g] = g(alpha,phi)

if alpha~=1
    g = alpha^2*(1-alpha^2)^-1*(3*phi-2);
else
    g = 0;
end