[step:Show $I = (f)$ via Gauss's Lemma applied to divisibility]The inclusion $(f) \subset I$ is immediate from Step 2: $f \in I$, and $I$ is an ideal, so every multiple of $f$ in $k[W_1, \ldots, W_{n+1}]$ lies in $I$.
For the reverse inclusion, take any $g \in I$, i.e. $g(x_1, \ldots, x_n, y) = 0$ in $K$. View $g$ as an element of the polynomial ring $k(x_1, \ldots, x_n)[W_{n+1}]$ — that is, regard the variables $W_1, \ldots, W_n$ as substituted with $x_1, \ldots, x_n$, leaving a polynomial in the single variable $W_{n+1}$ with coefficients in the field $k(x_1, \ldots, x_n)$. Concretely, this is the image of $g$ under the homomorphism
\begin{align*}
k[W_1, \ldots, W_n][W_{n+1}] &\to k(x_1, \ldots, x_n)[W_{n+1}] \\
g(W_1, \ldots, W_n, W_{n+1}) &\mapsto g(x_1, \ldots, x_n, W_{n+1}).
\end{align*}
The image polynomial $g(x_1, \ldots, x_n, W_{n+1}) \in k(x_1, \ldots, x_n)[W_{n+1}]$ has $y$ as a root: $g(x_1, \ldots, x_n, y) = 0$ in $K$. Since $h$ is the minimal polynomial of $y$ over $k(x_1, \ldots, x_n)$ and $y$ is a root of $g(x_1, \ldots, x_n, W_{n+1})$, we have $h(W_{n+1}) \mid g(x_1, \ldots, x_n, W_{n+1})$ in $k(x_1, \ldots, x_n)[W_{n+1}]$.
Lifting back: since $f$ and $h$ are associates in $k(W_1, \ldots, W_n)[W_{n+1}]$ (Step 3) and divisibility is preserved under specialisation, $f(W_1, \ldots, W_n, W_{n+1}) \mid g(W_1, \ldots, W_n, W_{n+1})$ in $k(W_1, \ldots, W_n)[W_{n+1}]$.
Now we apply the divisibility version of *Gauss's Lemma*: if $R$ is a UFD, $F = \operatorname{Frac}(R)$, and a primitive polynomial $f \in R[T]$ divides $g \in R[T]$ in $F[T]$, then $f$ divides $g$ in $R[T]$. The hypotheses are met: $R = k[W_1, \ldots, W_n]$ is a UFD, $f$ is primitive in $R[W_{n+1}]$ (by Step 3), and $f \mid g$ in $F[W_{n+1}]$ was just established. Hence $f \mid g$ in $R[W_{n+1}] = k[W_1, \ldots, W_{n+1}]$. So $g \in (f)$, completing the inclusion $I \subset (f)$.
Combining the two inclusions, $I = (f)$.[/step]