I(k)=A(k)(∑_(i=1)^p)[((Subscript[I, ref])/(2p))+((Subscript[I, i])/2)+(Sqrt[(Subscript[I, ref])(Subscript[I, i])])cos(2(Subscript[kz, i]))]