eq:mkappa is similar to what we see in vMF maximum-likelihood estimation, therefore, we can use the maximum-likelihood derivations from now on (See Appendix of [1]} equations A.7 to A.8) which leads to maximum-likelihood estimation
\(\kappa = A^{-1}_d(\bar{r})~,\) 
