``` saturate((cos(alpha)-cos(phi))/(cos(theta)-cos(phi))^P) ```