This paper is showing how KL divergence can be minimized by matching the expected values of the sufficient statistics. More precisely, For any distribution p of the exponential family with pdf: $$ p_{\theta}(x) = \frac{1} {Z(\theta)}\exp(\theta^T \phi(x))$$ the distribution $$p_{θ∗} $$ which minimises the Kullback-Leibler divergence, $$ KL (p||p_{θ∗})$$ over the exponential family with natural statistic $\phi$ is implicitly given by: $$E_{p_{\theta}*(x)} [\phi(x)] = E_{p_{(x)}} [\phi(x)]$$
I need help to show that this is true for the particular case of a normal-gamma distribution. So I need to make the same proof that is on the paper but for the special case of Normal gamma.