[step:Clear denominators to obtain $f^N \in I$]Each $h_i(X, T)$ is a polynomial in $T$ with coefficients in $k[X_1, \ldots, X_n]$; write
\begin{align*}
h_i(X, T) = \sum_{j=0}^{d_i} h_{i,j}(X)\, T^j
\end{align*}
with $h_{i,j} \in k[X_1, \ldots, X_n]$ and $d_i \ge 0$. Substituting $T = 1/f$,
\begin{align*}
h_i(X, 1/f) = \sum_{j=0}^{d_i} h_{i,j}(X)\, f^{-j} = f^{-d_i} \sum_{j=0}^{d_i} h_{i,j}(X)\, f^{d_i - j}.
\end{align*}
Let $D = \max_{1 \le i \le r} d_i$. Multiplying $(\ast\ast)$ through by $f^D \in k[X_1, \ldots, X_n]$,
\begin{align*}
f^D = \sum_{i=1}^r f^{D - d_i} \cdot f^{d_i} \cdot h_i(X, 1/f) \cdot g_i(X) = \sum_{i=1}^r \tilde{h}_i(X) \cdot g_i(X),
\end{align*}
where
\begin{align*}
\tilde{h}_i(X) := f^{D - d_i}(X) \cdot \sum_{j=0}^{d_i} h_{i,j}(X)\, f(X)^{d_i - j} = \sum_{j=0}^{d_i} h_{i,j}(X)\, f(X)^{D - j} \in k[X_1, \ldots, X_n].
\end{align*}
The polynomial $\tilde{h}_i$ has no denominators and lies in $k[X_1, \ldots, X_n]$. We have written
\begin{align*}
f^D = \sum_{i=1}^r \tilde{h}_i \cdot g_i \in (g_1, \ldots, g_r) = I.
\end{align*}
We now choose $N \ge 1$. If $D \ge 1$, set $N = D$: then $f^N \in I$ directly, so $f \in \sqrt{I}$.
If $D = 0$ — meaning every $h_i$ is constant in $T$, so $h_i(X, 1/f) = h_{i,0}(X)$ — the identity $(\ast\ast)$ reads
\begin{align*}
1 = \sum_{i=1}^r h_{i,0}(X) \cdot g_i(X),
\end{align*}
which is already an identity in $k[X_1, \ldots, X_n]$ and exhibits $1 \in I$. Hence $I = k[X_1, \ldots, X_n]$, so $f \cdot 1 = f \in I$ and therefore $f^1 \in I$, giving $f \in \sqrt{I}$ with $N = 1$. (Geometrically, $I = k[X_1, \ldots, X_n]$ means $V(I) = \varnothing$, and the equality $I(V(I)) = \sqrt{I}$ reduces to $I(\varnothing) = k[X_1, \ldots, X_n] = \sqrt{I}$, which holds directly from the definitions.)
As $f \in I(V(I))$ was arbitrary, $I(V(I)) \subset \sqrt{I}$.
Combining with Step 1, $I(V(I)) = \sqrt{I}$, completing the proof.[/step]