Motion of a pendulum The equations of motions for a simple pendulum is given by 
$$\ddot{\theta} ~=~ -\frac{g}{\ell}\sin(\theta),$$
where $g$ is acceleration due to gravity and $\ell$ is the length of the pendulum's string. Notice that the differential equation is of second order, does this mean that if I solve this equation numerically, the numbers that I get refers to the change in the velocity of the pendulum?
 A: This is really the same as Mark's answer but phrased in a different way. If you write your equation in Leibniz's notation it is:
$$ \frac{d^2\theta}{dt^2} = -\frac{g}{\ell}\sin(\theta) $$
This makes it clearer that the solution is going to be the function $\theta(t)$ i.e. the angle as a function of time.
A: If you solve the equation for $\theta$, whether numerically or otherwise, you will have $\theta$, not $\dot\theta$ or $\ddot{\theta}$. You can get those by differentiating, of course, or you may find them as part of your numerical algorithm along the way. The fact that the equation is second order simply means you will need two initial conditions to find a solution.
A: This might be helpful if you want to check the accuracy of your
numerical integration scheme.
Normalize time using $\tau=t\sqrt{g/l}$ giving
$$\frac{d^2\theta}{d\tau^2}=-\sin\theta$$
which can be written as
$$
\frac12\frac{d}{d\theta}\left(\frac{d\theta}{d\tau}\right)^2=-\sin\theta
$$
Integrating this gives,
$$
\frac{d\theta}{d\tau}=\sqrt{c+2\cos\theta}
$$
Put $\theta=2\lambda$ so
$$
d\tau=\frac{2d\lambda}{\sqrt{c+2-4\sin^2\lambda}}
$$
and thus
$$
\tau(\theta)=\tau\left(\theta_0\right)+\frac{2}{\sqrt{2+c}}\int_{\theta_0/2}^{\theta/2}\frac{d\lambda}{\sqrt{1-\frac{4}{2+c}\sin^2\lambda}}
$$
which is the tabulated form of the standard elliptic integral with modulus $4/(2+c)$. Numerical solution requires choosing initial values $\theta_0$ and $\left(\frac{d\theta}{d\tau}\right)_0$ and thus
$$
c=\left(\frac{d\theta}{d\tau}\right)_0^2-2\cos\theta_0
$$
