Solution (source code)

= Solution

The <log-likelihood> has derivative
$$
\ell'(\theta)=\sum_iY_i-nA'(\theta).
$$
Differentiating the normalizing identity for the family gives
$$
A'(\theta)=\mathbb E_\theta Y_1=\mu(\theta),
\qquad
A''(\theta)=\operatorname{Var}_\theta(Y_1)\geq0.
$$
Hence an interior maximum satisfies, and under nondegeneracy uniquely satisfies,
$$
\boxed{\mu(\widehat\theta_{\rm MLE})=\overline Y}.
$$