cos⁡ψ=(cos⁡(2⁢θ)+k2⁢sin4⁡θ)/(1-k2⁢sin4⁡θ),