I think this answer might essentially be the same as @orangeskid's answer, but I'm not sure, so I figured I will post it here. If it is the same or incorrect, please let me know and I will delete it.
Let $I = \langle x_1 - a_1, \dots, x_n - a_n \rangle \subset k[x_1, \dots, x_n]$ for $k$ a field. Consider any ideal $I \subsetneq J$, by the Hilbert Basis theorem $J$ is finitely generated, so we can say that $J=\langle x_1 - a_1, \dots, x_n -a_n, f_1, \dots, f_r\rangle$ for finitely many $1 \le i \le r < \infty$, with $f_i \not\in I, f_i \in k[x_1, \dots, x_n]$ and without loss of generality, all of the $f_i$ are not identically equal to zero.
Case 1: At least one of the $f_i$ is constant.
Without loss of generality, $f_r$ is constant, say $f_r \equiv c \in k, c\not=0$ so $J= \langle x_1 - a_1, \dots, c \rangle = \langle x_1 - a_1, \dots, 1 \rangle = k[x_1, \dots, x_n]$.
Case 2: None of the $f_i$ are constant.
Let's choose one of the $f_i$, which by an abuse of notation we will denote simply by $f_i(x_1, \dots, x_n)$.
Then we have that $$f_i(x_1, \dots, x_n) = f_i^m(x_1, \dots, x_n) + f_i^{m-1}(x_1 -a_1,\dots, x_n - a_n)+ \dots + f_i^1(x_1 -a_1, \dots, x_n -a_n)+ f_i(x_1-a_i,\dots,x_n-a_n) $$ where for all $1 < j \le m, \quad f_i^j(x_1, \dots, x_n) := f_i^{j-1}(x_1,\dots,x_n)-f_i^{j-1}(x_1 - a_1, \dots, x_n-a_n)$, and additionally $f_i^1(x_1, \dots, x_n):=f(x_1,\dots,x_n)-f(x_1-a_1, \dots, x_n-a_n)$, and finally that $m$ is the lowest natural number such that $f_i^m(x_1, \dots, x_n)$ is identically constant, thus for any $k >m$, $f_i^k(x_1, \dots, x_n)\equiv 0$. Such an $m$ is guaranteed to exist, because by construction, $\deg(f_i^1) < \deg(f_i)$ and $\deg(f_i^j) < \deg(f_i^{j-1})$ for all $2 \le j \le m$.
Define $c$ to be the non-zero constant such that $f_i^m(x_1, \dots, x_n) \equiv c $. Now clearly we have that $$ f_i(x_1 -a _1, \dots, x_n -a_n) \in \langle x_1 -a_1, \dots, x_n -a_n \rangle \\ f_i^j(x_1-a_1, \dots, x_n -a_n) \in \langle x_1 - a_1, \dots, x_n - a_n \rangle \quad \forall\ 1\le j \le m-1$$ Therefore we have shown that $f_i(x_1, \dots, x_n) = c + g_i (x_1, \dots, x_n)$ for some $g_i \in \langle x_1 - a_1 , \dots, x_n - a_n \rangle$. Therefore we have that $$ J = \langle x_1 - a_1, \dots, x_n - a_n, \dots, f_i, \dots, f_r \rangle \\ = \langle x_1 - a_1, \dots, x_n - a_n, \dots, c + g_i, \dots, f_r \rangle \\ = \langle x_1 - a_1, \dots, x_n - a_n, \dots, c, \dots, f_r \rangle \\ = \langle x_1 - a_1, \dots, x_n - a_n, \dots, 1, \dots, f_r \rangle = k[x_1, \dots, x_n]$$
Thus, in both cases, $I \subsetneq J \implies J = k[x_1, \dots, x_n]$, so therefore $I$ must be maximal, since clearly $I=\langle x_1 - a_1, \dots, x_n - a_n\rangle$ is proper.