A Schauder basis of a Hilbert space (or more generally a Banach space) is a linearly independent set of vectors such that every vector can be uniquely written as the sum of a norm convergent series where the individual terms of the series are multiples of the basis vectors.
Any polynomial is a linear combination of a finite number of our basis vectors.
Uniform convergence on a finite measure space (such as a bounded interval with Lebesgue measure) implies convergence in quadratic mean; in fact, for any two essentially bounded functions $f$ and $g$
$$\|f-g\|_2\leq\|f-g\|_\infty\cdot\sqrt{b-a}.$$
Finally note that the indicator function of any interval, and therefore any block function, can be approximated in $L^2$ norm by continuous functions.
From all of this it follows that the span of our normalised set of vectors is dense in $L^2.$ So an arbitrary $f\in L^2$ can be written as the limit of a sequence of polynomials $f_n.$ We shall now argue that it can even be written as the sum of the series $\sum_{i=1}^\infty\langle f,e_i\rangle e_i.$
If $f_n$ is a polynomial of degree $m(n),$ we obtain a better approximation replacing $f_n$ with the orthogonal projection of $f$ on the $m(n)$-dimensional subspace generated by the first $m(n)$ vectors of our Gram-Schmidt basis:
$$\overline f_n=\sum_{i=1}^{m(n)}\langle f,e_i\rangle e_i.$$
Thus at least some partial sums of the series $\sum_{i=1}^\infty\langle f,e_i\rangle e_i$ converge to $f.$ But the other partial sums are never worse approximations than the previous ones, i.e., $\sum_{i=1}^{k+1}\langle f,e_i\rangle e_i$ is always at least as close to $f$ as $\sum_{i=1}^{k}\langle f,e_i\rangle e_i$ is, by orthonormality.
The arbitrary nature of $f$ proves that $\{e_1,e_2,\ldots\}$ is a Schauder basis.