E is nonempty: 0n=0<a, so 0∈E. And a+1 is an upper bound of E: if x>a+1 then, since a>0, x>1, hence xn≥x>a+1>a, so x∈/E.
One monotonicity fact, proved directly (used in steps 1 and 4; the standalone proof is the next problem): for 0≤u<v, vn−un=(v−u)(vn−1+vn−2u+⋯+un−1)>0, because v−u>0 and the second factor is a sum of nonnegative terms including vn−1>0. Hence t↦tn is strictly increasing on [0,∞), and x≥y≥0 implies xn≥yn.
By completeness of the reals (the least-upper-bound axiom), s=supE exists. We have s≥0; we show sn=a by trichotomy: ruling out sn<a and sn>a leaves sn=a.
Suppose sn<a. For δ>0, the binomial theorem gives (s+δ)n−sn=∑k=1n(kn)sn−kδk=δ∑k=1n(kn)sn−kδk−1. For 0<δ<1, δk−1≤1, so this difference is at most δM, where M=∑k=1n(kn)(s+1)n−k depends only on s.
Picking 0<δ<min(1,(a−sn)/M) gives (s+δ)n<sn+(a−sn)=a, and s+δ≥0, so s+δ∈E — contradicting s being an upper bound.
Suppose sn>a. Since a≥0, s>0. For 0<δ<s, put t=s−δ; the factorization sn−tn=(s−t)(sn−1+sn−2t+⋯+tn−1) gives sn−(s−δ)n=δ∑j=0n−1sn−1−jtj≤nδsn−1 (each term is at most sn−1 since t≤s).
Picking 0<δ<min(s,(sn−a)/(nsn−1)) yields (s−δ)n>a. But every x≥s−δ satisfies xn≥(s−δ)n>a, so x∈/E — i.e. s−δ is also an upper bound of E, contradicting s being the least upper bound.
Hence sn=a, and since s≥0, s is the nonnegative n-th root of a. (Uniqueness is the next problem.) This proof grabs the continuity of xn by hand through the binomial expansion — an argument that works precisely because the least-upper-bound axiom guarantees a "hole-free" real line.