I came across the following problem, let $N![z^N]A(z)$ denote the coefficient of an exponential generating function (EGF) $A(z)$. The EGF is similar to an ordinary generating function (OGF) $A'(z)$ except that instead of the series $A'(z)=\sum_0^Na_Nz^N$ for an OGF, we have $A(z)=a_Nz^N/N!$ for an EGF $A(z)$.
For instance, for EGF $A(z)=e^z$, we have $N![z^N]A(z)=N![z^N]e^z=1$, i.e. the EGF coefficients of $e^z$ are $1$ for all $N \in \mathbb{N}$, i.e. the coefficients of the set $ \{z^0/0!,z^1/1!,z^2/2!,...,z^N/N! \}$ are all $1$ for any $N$ given $e^z$. Similarly, for EGF $A(z)=1/(1-z)$, we have $N!z[^N](1/(1-z))=N!$.
Now, given the following EGF $A(z)$:
$$ A(z)=e^z\int^z_0\frac{1-e^{-t}}{t}dt $$
We are supposed to get $N![z^N]A(z)=H_N$, where $H_N$ is the $N$th harmonic number, i.e.
$$ N![z^N]e^z\int^z_0\frac{1-e^{-t}}{t}dt = H_N $$
I could not think of a way to prove the above statement. The problem gave a hint that proving this statement involves forming a differential equation for the EGF $H(z)=\sum_{N \geq 0}H_Nz^N/N!$...
Any help?