Here's a proof without using complex analysis, differential geometry or the projective plane, and only using results from general topology.
Consider the equivalence relation on $\mathbb{C}^n$ given by $x\sim y$ iff $y = sx := (x_{s(1)}, ..., x_{s(n)})$ for some $s\in S_n$. Give $\mathbb{C}^n/\sim$ the quotient topology and let $q:\mathbb{C}^n\to\mathbb{C}^n/\sim$ be the quotient map. Note that $q^{-1}q[x+B_r] = \bigcup_{s\in S_n} (sx+B_r)$ is open, where $B_r = \{y : |y| < r\}$, so that $q[x+B_r]$ is open. If $U\subseteq \mathbb{C}^n/\sim$ is open and $q(x)\in U$ then by taking $r > 0$ such that $x+B_r\subseteq q^{-1}[U]$, it follows that $\{q(x+B_r) : x\in \mathbb{C}^n, r > 0\}$ forms a basis for $\mathbb{C}^n/\sim$. Let $$d(q(x), q(y)) = \min\{|x'-y'| : x'\sim x, y'\sim y\}.$$ Clearly $d(q(x), q(y)) = 0$ iff $q(x) = q(y)$ and $d(q(x), q(y)) = d(q(y), q(x))$. Given $x, y, z\in\mathbb{C}^n$, let $x'\sim x$ be such that $d(q(x), q(y)) = |x'-y|$. Then $$d(q(x), q(z))\leq |x'-z'| \leq |x'-y|+|y-z'| = d(q(x), q(y))+|y-z'|$$ for $z'\sim z$. Taking minimum over all $z'\sim z$, we obtain triangle inequality for $d$. So $d$ is a metric. Notice that $q[x+B_r] = \{q(y) : d(q(x), q(y)) < r\}$ so that $\mathbb{C}^n/\sim$ is metrizable with metric $d$.
Let $\sigma:\mathbb{C}^n\to\mathbb{C}^n$ map $(\lambda_1, ..., \lambda_n)$ to coefficients of $(x-\lambda_1)...(x-\lambda_n) = x^n+a_1x^{n-1}+...+a_n$. From Viete's formulas, $\sigma$ is continuous. Note the following theorem:
Theorem 1. Let $f:X\to Y$ be a continuous map between metrizable spaces. Then $f$ is perfect iff for all $x_n\in X$ such that $f(x_n)$ is convergent, it follows that $x_n$ has a convergent subsequence.
For reference of this see chapter 2 of Geometric Aspects of General Topology by K. Sakai.
Suppose that $\lambda^{(m)} = (\lambda_1^{(m)}, ..., \lambda_n^{(m)})$ is such that $\sigma(\lambda^{(m)}) = (a_1^{(m)}, ..., a_n^{(m)})$ is convergent. Note that $$|\lambda_i^{(m)}|^n = \left|\sum_{j=1}^n a_j^{(m)}(\lambda_i^{(m)})^{n-j}\right|\leq \sum_{j=1}^n |a_j^{(m)}||\lambda_i^{(m)}|^{n-j}.$$ If $|\lambda_i^{(m)}|\geq 1$ then by dividing above inequality by $|\lambda_i^{(m)}|^{n-1}$, it follows that $|\lambda_i^{(m)}|\leq \sum_{j=1}^n |a_j^{(m)}|$, thus for arbitrary $\lambda_i^{(m)}$ we have $|\lambda_i^{(m)}|\leq \max(1, \sum_{j=1}^n |a_j^{(m)}|)$ so that $\lambda_i^{(m)}$ are bounded. Thus we can take convergent subsequence $\lambda^{(m_k)}$ of $\lambda^{(m)}$. It follows that $\sigma$ is a perfect map.
Theorem 2. Let $f:X\to Y$, $q:X\to Z$ a quotient map and $g = f\circ q^{-1}:Z\to Y$ be well-defined. If $f$ is continuous/open/closed, then so is $g$.
See e.g. chapter 6 of Dugundji's Topology.
Thus since $\sigma$ is a perfect map, it follows that $\sigma\circ q^{-1}:\mathbb{C}^n/\sim\to \mathbb{C}^n$ is a closed continuous bijection, thus a homeomorphism. As a corollary we can state:
Corollary 1. Let $a_1, ..., a_n\in\mathbb{C}$ and $\lambda_1, ..., \lambda_n$ be roots of $x^n+a_1x^{n-1}+...+a_n$, $\varepsilon > 0$, then there exists $\delta > 0$ such that if $|a_i-b_i| < \delta$ then there exists an arrangement $\mu_1, ..., \mu_n$ of roots of $x^n+b_1x^{n-1}+...+b_n$ such that $|\lambda_i-\mu_i| < \varepsilon$ for $i = 1, ..., n$.
Corollary 2. If $a_1, ..., a_n\in\mathbb{R}$, let $\lambda_1, ..., \lambda_n$ be roots of $x^n+a_1x^{n-1}+...+a_n$ ordered so that $\text{Re}(\lambda_1) \leq ... \leq \text{Re}(\lambda_n)$. Define $f:\mathbb{R}^n\to\mathbb{R}^n$ as $f(a_1, ..., a_n) = (\text{Re}(\lambda_1), ..., \text{Re}(\lambda_n))$. Then $f$ is continuous.
Proof: The map $g:\mathbb{R}^n\to\mathbb{R}^n$ ordering elements of $\mathbb{R}^n$ increasingly is continuous, since if $a\in \mathbb{R}^n$ and $r_1 < ... < r_m$ are distinct among elements in $g(a)$, then for any $\varepsilon > 0$, if $0 < \delta < \min(\varepsilon, \frac{r_2-r_1}{2}, ..., \frac{r_m-r_{m-1}}{2})$ and $|a-b| < \delta$, then $|g(a)-g(b)| = |a-b| < \varepsilon$. The map $f$ is then composition of $(\sigma\circ q^{-1})^{-1}\restriction_{\mathbb{R}^n}$, the map $\mathbb{C}^n/\sim\to \mathbb{R}^n/\sim$ induced by the map $(z_1, ..., z_n)\to (\text{Re}(z_1), ..., \text{Re}(z_n))$, and the map $\mathbb{R}^n/\sim\to\mathbb{R}^n$ induced by $g$, and as such is continuous. $\square$