Let $Y$ be a random vector of $\mathbf{R}^p.$
Assume that $Y$ is $t_p(\mu,\Sigma,\nu)$ distributed. That is the density is given by $$f_p(y,\mu,\Sigma,\nu)={\displaystyle {\frac {\Gamma \left[(\nu +p)/2\right]}{\Gamma (\nu /2)\nu ^{d/2}\pi ^{p/2}\left|{\boldsymbol {\Sigma }}\right|^{1/2}}}\left[1+{\frac {1}{\nu }}({\mathbf {y} }-{\boldsymbol {\mu }})^{T}{\boldsymbol {\Sigma }}^{-1}({\mathbf {y} }-{\boldsymbol {\mu }})\right]^{-(\nu +p)/2}}$$
Suppose that $y_1,\ldots,y_n$ denote an observed random sample from $t_p(\mu,\Sigma,\nu)$ i.e. $$y=(y_1^T,\ldots,y_n^T)^T.$$
To perform the ML estimation, we compute the log likelihood $$\log(L(\Psi))=\sum_{i=1}^n\log f_p(y_i; \mu,\Sigma,\nu)$$ $$=-\frac{1}{2}np\log(\pi\nu)+n\bigl(\log \Gamma(\frac{\nu+p}{2})-\log \Gamma(\frac{\nu}{2})\bigr)-\frac{1}{2}n\log \vert \Sigma\vert+\frac{1}{2}n(\nu+p)\log \nu-\frac{1}{2}(\nu+p)\sum_{i=1}^n \log( \nu +\delta(y_i;\mu,\Sigma))$$
where $\delta(y_i;\mu,\Sigma)$ is the Mahalanobis square distance.
Question: Is there a closed form for the solution of the likelihood equation?