Solution (source code)

= Solution

Substituting $\widetilde\theta=\tau^2y/(\sigma^2+\tau^2)$ and $\sigma_\theta^2=\sigma^2\tau^2/(\sigma^2+\tau^2)$ into parts c and d and simplifying gives
$$
\mathbb E_{\theta\mid y}\log L(\theta)
-D_{\mathrm{KL}}\{p(\theta\mid y)\Vert\pi\}
=-\frac12\log\{2\pi(\sigma^2+\tau^2)\}
-\frac{y^2}{2(\sigma^2+\tau^2)}
=\log Z.
$$
Thus the equality holds.