what you are trying to do is essentially the same as proving a real symmetric matrix is diagonalizable. i know "a" proof of that, how "easy" it is depends on your perspective.
the idea is to enlarge our scope, and consider complex matrices instead (so that our eigenvalues are always in our underlying field). furthermore, we need to consider the complex "analogues" of orthogonal matrices, and symmetric matrices:
unitary matrices, and hermetian matrices.
to definite that, we need to define the analogue of the transpose for complex matrices: the conjugate-transpose:
$A^{\ast} = (\overline{a_{ji}})$ for $A = (a_{ij})$
so a unitary matrix is one for which:
$U^{\ast} = U^{-1}$
and a hermetian matrix is one for which:
$A^{\ast} = A$.
now, on to the business at hand. the nail we're going to use to crack this nut is this lemma:
for every square complex matrix A, there exists a *unitary* matrix U, with:
UAU* = UAU-1 = T, where T is upper-triangular.
we'll use induction on n. for a 1x1 matrix, this is trivial.
now we assume this is true for any (n-1)x(n-1) complex matrix. let A be an nxn complex matrix. now det(xI - A) is a complex polynomial of degree n, which has a root, $\lambda$, and thus A has an eigenvector, u belonging to $\lambda$ (there is always at least non-zero one solution to (A - λI) because it is singular). we can normalize this vector, giving a unit eigenvector, so we lose nothing by assuming u is a unit eigenvector to begin with.
using gram-schmidt, we can extend {u} first to a basis for $\Bbb C^n$ and then to an orthogonal basis for $\Bbb C^n$. the orthogonality of this basis (under the standard inner product for $\Bbb C^n$) means that if we form a matrix U1 whose first column is u, and the other columns are the other orthogonal basis vectors, then this matrix is unitary. let e1 be the complex vector (1,0,...,0).
now note that U1-1AUe1 = U1Au = U1-1(λu) = λU1-1u = λe1.
thus in block form, we have:
$U_1^{-1}AU_1 = \begin{bmatrix}\lambda&*\\0&B \end{bmatrix}$
where B is an (n-1)x(n-1) matrix.
applying our induction hypothesis, there is an (n-1)x(n-1) unitary matrix U' such that:
U'-1BU = U'*BU = T', is upper-triangular.
define U2 to be the matrix (in block form):
$U_2 = \begin{bmatrix}1&0\\0&U' \end{bmatrix}$
then:
$U_2^{\ast} = \begin{bmatrix}1&0\\0&(U')^{\ast} \end{bmatrix}$
so $U_2U_2^{\ast} = \begin{bmatrix}1&0\\0&U'(U')^{\ast} \end{bmatrix} = I$
since U' is unitary, thus U2 is likewise unitary. ok, so now we have:
(U1U2)-1A(U1U2) = U2-1(U1-1AU1)U2 =
$\begin{bmatrix}1&0\\0&(U')^{\ast} \end{bmatrix} \begin{bmatrix}\lambda&\ast \\0&B \end{bmatrix} \begin{bmatrix}1&0\\0&U' \end{bmatrix}$
$= \begin{bmatrix}\lambda&\ast \\0&(U')^{\ast}BU' \end{bmatrix} = \begin{bmatrix}1&\ast \\0&T' \end{bmatrix}$
which is upper-triangular, since T' is upper-triangular.
so let U = U1U2. then this is the desired unitary matrix, since:
U-1 = (U1U2)-1 = U2-1U1-1 = U2*U1* = (U1U2)* = U*
(the proof of the 4th equality is easy and analogous to the similar proof for the transpose, left to the interested reader).
this concludes the proof of the lemma (which actually does most of the heavy lifting for us). it's pretty much downhill from here.
observation #1: if the matrix A is hermetian, so is the upper-triangular matrix T which is unitarily similar to A.
why? because (U*AU)* = U*A*(U*)* = U*AU (since A* = A, if A is hermetian).
but...if T is upper-triangular, AND hermetian, it must be LOWER triangular, as well. which means when A is hermetian, T is DIAGONAL. furthermore, since the diagonal elements of T equal their own complex-conjugates, T is a REAL matrix.
observation #2: any real symmetric matrix is hermetian. thus any real symmetric matrix is unitarily similar to a real diagonal matrix, D.
observation #3: the eigenvalues of A are the eigenvalues of D, and vice-versa.
proof: det(xI - A) = det(U-1(xI - A)U) = det(U-1(xI)U - U-1AU) = det(xI - D)
observation #4: since the eigenvalues of A (a real symmetric matrix) are real, the eigenvectors are likewise real. thus we may take U to be a real unitary matrix, that is, an orthogonal one.
now suppose that a real, symmetric matrix A has an eigenvalue of (algebraic) multiplicity m. since A is orthogonally similar to a diagonal matrix D with the same eigenvalues (via an orthogonal matrix Q), D has m diagonal entries that are the same.
note that Q-1AQ = D, means AQ = QD. since Q is invertible, its columns are linearly independent (it's square of full rank).
suppose $\{u_{j_1},\dots,u_{j_m}\}$ are the columns of Q corresponding to the diagonal entry λ of D repeated m times (the eigenvalue of algebraic multiplicity m). then:
$A(u_{j_k}) = (u_{j_k})^TD = \lambda u_{j_k}\ k = 1,\dots,m$
giving m linearly independent eigenvectors of the eigenspace Eλ.