@Patrick87 gives a great answer for your specific case, I thought I would give a tip of how to find $s_L(n)$ in the more general case of any language $L$ that can be represented by an irreducible DFA (i.e. if it is possible to get to any state from any state). Note that your language is of this type.
Proof of theorem for irreducible DFAs
Let $D$ be the transition matrix of your $m$-state DFA, since it is irreducible, the matrix is normal and has a full eigenbasis $|\lambda_1\rangle ... |\lambda_m\rangle$. Let $|A\rangle$ be the accept vector: i.e. $\langle i | A \rangle$ is 1 if $i$ is an accept state, and 0 otherwise. WLOG assume that $|1\rangle$ is the initial state, and since we have a complete eigenbasis, we know that $|1\rangle = c_1|\lambda_1\rangle + ... + c_m|\lambda_m\rangle$ for some coefficients $c_1 ... c_m$ (note that $c_i = \langle \lambda_i | i \rangle$).
Now we can prove a restricted case of the theorem in the question (restricted to irreducible DFAs; as an exercise generalize this proof to the whole theorem). Since $D$ is the transition matrix $D|1\rangle$ is the vector of states reachable after reading any one character, $D^2|1\rangle$ is the same for two characters, etc. Given a vector $|x\rangle$, $\langle A|x\rangle$ is simply the sum of the components of $|x\rangle$ that are accept states. Thus:
$$
\begin{align}
s_L(n) & = \langle A |D^n| 1 \rangle \\
& =\langle A | D^n (c_1 |\lambda_1\rangle ... c_m |\lambda_m\rangle) \\
& = c_1 \lambda_1^n \langle A | \lambda_1 \rangle + ... + c_m \lambda_m^n \langle A | \lambda_m \rangle \\
& = \langle A | \lambda_1 \rangle\langle\lambda_1|1\rangle \lambda_1^n + ... + \langle A | \lambda_m \rangle\langle \lambda_m | 1 \rangle \lambda_m^n \\
& = p_1\lambda_1^n + ... + p_m \lambda_m^m
\end{align}
$$
Now we know that for an irreducible m-state DFA, $p_1 ... p_m$ will be zero order polynomials (i.e. constants) that depends on the DFA and $\lambda_1 ... \lambda_m$ will be eigenvalues of the transition matrix.
Generality note
If you want to prove this theorem for arbitrary DFA, then you will need to look at the Schur decomposition of $D$ and then polynomials of non-zero degree will pop up because of the nilpotent terms. It is still enlightening to do this, since it will let you bound the max degree of the polynomials. You will also find a relationship between how complicated the polynomials are and how many $\lambda$s you will have.
Application to specific question
For your language $L$ we can select the DFA with transition matrix:
$$ D = \begin{pmatrix} 0 & 1 \\ 1 & 0 \end{pmatrix} $$
and accept vector:
$$ A = \begin{pmatrix} 1 \\ 0 \end{pmatrix}$$
Find the eigenvectors and their eigenvalues $\lambda_1 = 1$ with $| \lambda_1 \rangle = \frac{1}{\sqrt{2}}\begin{pmatrix} 1 \\ 1 \end{pmatrix}$ and $\lambda_2 = -1$ with $| \lambda_2 \rangle = \frac{1}{\sqrt{2}}\begin{pmatrix} 1 \\ -1 \end{pmatrix}$. We can use this to find $p_1 = 1/2$ and $p_2 = 1/2$. To give us:
$$s_L(n) = \frac{1}{2} + \frac{1}{2}(-1)^n$$