[step:Show $(X_1 - a_1, \ldots, X_n - a_n)$ is itself maximal, forcing equality]
Define the evaluation $k$-algebra homomorphism
\begin{align*}
\mathrm{ev}_a : k[X_1, \ldots, X_n] &\to k \\
f &\mapsto f(a_1, \ldots, a_n).
\end{align*}
This is a ring homomorphism (polynomial addition and multiplication commute with evaluation), and surjective onto $k$ (every $c \in k$ is the image of the constant polynomial $c$). By the First Isomorphism Theorem,
\begin{align*}
k[X_1, \ldots, X_n] / \ker(\mathrm{ev}_a) \cong k.
\end{align*}
Since $k$ is a field, the right-hand side is a field, hence $\ker(\mathrm{ev}_a) \subset k[X_1, \ldots, X_n]$ is a maximal ideal.
We claim $\ker(\mathrm{ev}_a) = (X_1 - a_1, \ldots, X_n - a_n)$.
\textbf{Inclusion $(X_1 - a_1, \ldots, X_n - a_n) \subseteq \ker(\mathrm{ev}_a)$.} Each generator $X_i - a_i$ evaluates to $a_i - a_i = 0$ at the point $a = (a_1, \ldots, a_n)$, so $X_i - a_i \in \ker(\mathrm{ev}_a)$ for each $i$. The kernel is an ideal, hence contains the ideal generated by the $X_i - a_i$.
\textbf{Inclusion $\ker(\mathrm{ev}_a) \subseteq (X_1 - a_1, \ldots, X_n - a_n)$.} Let $f \in \ker(\mathrm{ev}_a)$, so $f(a_1, \ldots, a_n) = 0$. Substitute $X_i = Y_i + a_i$ (equivalently $Y_i := X_i - a_i$) — this defines a $k$-algebra automorphism
\begin{align*}
\Phi : k[X_1, \ldots, X_n] &\to k[Y_1, \ldots, Y_n] \\
f(X_1, \ldots, X_n) &\mapsto f(Y_1 + a_1, \ldots, Y_n + a_n)
\end{align*}
(the inverse is $Y_i \mapsto X_i - a_i$). Expand $\Phi f$ as a $k$-linear combination of monomials in $Y$:
\begin{align*}
(\Phi f)(Y) = \sum_{\alpha \in \mathbb{N}^n} c_\alpha Y^\alpha = c_0 + \sum_{|\alpha| \geq 1} c_\alpha Y^\alpha.
\end{align*}
Evaluating both sides at $Y = 0$ gives $(\Phi f)(0) = c_0$ on one hand, and $f(a_1, \ldots, a_n) = 0$ on the other, so $c_0 = 0$. Hence
\begin{align*}
(\Phi f)(Y) = \sum_{|\alpha| \geq 1} c_\alpha Y^\alpha,
\end{align*}
and every term has $|\alpha| \geq 1$, so is divisible by some $Y_i$ — i.e.\ lies in the ideal $(Y_1, \ldots, Y_n) \subset k[Y_1, \ldots, Y_n]$. Reversing the substitution: $\Phi f \in (Y_1, \ldots, Y_n)$ corresponds to $f \in \Phi^{-1}((Y_1, \ldots, Y_n)) = (X_1 - a_1, \ldots, X_n - a_n)$. This is the desired inclusion.
Combining both inclusions,
\begin{align*}
\ker(\mathrm{ev}_a) = (X_1 - a_1, \ldots, X_n - a_n).
\end{align*}
Hence $(X_1 - a_1, \ldots, X_n - a_n)$ is a maximal ideal of $k[X_1, \ldots, X_n]$.
[/step]