function [phiu] = phi(alpha3)

if alpha3~=1
    phiu = alpha3*(1-alpha3^2)^(-3/2)*(acos(alpha3)-alpha3*sqrt(1-alpha3^2));
else
    phiu = 0;
end