How can I get a "better small angle approximation"? Is there a way to improve on the small angle approximation so that I get a more accurate answer. Is there any way to use a higher order series for the $\sin$ function or to use some sort of series approximation to the elliptical solution to the equation?
For reference, the equation is $$y'' = -\sin(y).$$

Expand $-\sin(y)$ to $-y+y^{3}/6$. Assume that the solution can be presented in the form
$$
y(x) = \sum_{i = 0}^{\infty} c_{i}x^{i}
$$
Then by the Cauchy product,
$$
y(x)^2 = \sum_{i = 0}^{\infty} \left(\sum_{j = 0}^{i} c_{j}c_{i-j}\right)x^{i} =: \sum_{i=0}^{\infty} 6b_{i}x^{i}.
$$
Then 
\begin{align}
y^{3}/6 - y 
&= y\left(y^{2}/6 -1 \right) = \sum_{j=0}^{\infty}c_{j}x^{j} \left(\sum_{i=0}^{\infty} b_{i}x^{i} -1 \right)\\
&= \sum_{i=0}^{\infty} \left(\sum_{j=0}^{i} c_{j}b_{i-j}  \right)x^{i}  - \sum_{i=0}^{\infty} c_{i}x^{i} \\
&= \sum_{i=0}^{\infty} \left[\sum_{j=0}^{i} c_{j}b_{i-j}  - c_{i}  \right]x^{i}
\end{align}
Since 
$$
y''(x) = \sum_{i=0}^{\infty} (i+1)(i+2)c_{i+2}x^{i}
$$
so that equating coefficients of like powers of $x$ yields
\begin{align}
c_{i+2} 
&= \frac{1}{(i+1)(i+2)}\left[\sum_{j=0}^{i} c_{j}b_{i-j}  - c_{i}  \right]  \\
&= \frac{1}{(i+1)(i+2)}\left[\sum_{j=0}^{i} \sum_{\ell=0}^{i-j}\frac{1}{6}c_{j}c_{\ell}c_{i-j-\ell}  - c_{i}  \right]  \\
\end{align}
It's expressions like these that make people love Runge-Kutta. But in any case, 
\begin{align}
c_{2} &= -\frac{1}{2}c_{0} && c_{3} = \frac{c_{1}}{6}\left(\frac{1}{2}c_{0}^{2} -1 \right)
\end{align}
and presumably you can program the recurrence relation to get as many terms as you'd like.
But the best solution is just to type DSolve[D[y[x],{x,2}] == -Sin[y], {y,x}] into Wolfram alpha and work with the result.
