(I((Subscript[k, x]),(Subscript[k, y]))≈((\[LeftBracketingBar](Subscript[E, R])\[RightBracketingBar])^2)[1+(a^2)+2a cos  (4 π(((Subscript[ε, 1])/( 1/(Subscript[P, 1])))-((Subscript[ε, 2])/( (Subscript[P, 2])))))]Equation   6)