a: 1.2
R = 0.6 (slider from 0 to 0.9999)
I_{FP}(phi) = ((a * (1 + R)^(2))) / ((1 - R^(2))^(2) + ((4 * R^(2)) * sin(phi)^(2)))