Hamiltonian for Geodesic Flow I'm trying to prove that geodesic flow on the cotangent bundle $T^* M$ is generated by the Hamiltonian vector field $X_H$ where 
$$H = \frac{1}{2}g^{ij}p_i p_j$$
but I am stuck. Could somebody show me how to complete the calculation, or where I've made a mistake? Cheers!
I know that vector field for geodesic flow is
$$X = p^i \partial/\partial x^i - \Gamma^i_{jk}p^j p^k \partial /\partial p^i$$
so I must verify that 
$$X \ \lrcorner\  \omega = - dH$$
where $\omega = dp_i \wedge dx^i$. It's easy to check that 
$$X \ \lrcorner\  \omega = - p^i dp_i - \Gamma^l_{jk}g_{li}p^j p^k dx^i$$
$$-dH = -g^{ij}p_j dp_i -\frac{1}{2}\partial_i g^{jk}p_j p_k dx^i$$ 
Here I am stuck. Using the explicit formula for the Christoffel symbols in terms of the metric doesn't seem to work! Have I done something wrong, or am I missing something?
N.B. I'm aware of solutions that do not use coordinates, but I'd like to understand this one which does! 
 A: I've worked out my mistake - the vector field $X$ above is incorrect. The subtlety comes in raising and lowering the relevant indices. 
Firstly we can all agree that the following vector field on the tangent bundle generates geodesic flow
$$Y = p^i \partial/\partial x^i - \Gamma^i_{jk} p^j p^k \partial/\partial p^i$$
where $(x^i, p^i)$ are coordinates on $TM$. 
Now we want to use the tangent-cotangent isomorphism to determine a vector field on $T^*M$ which also gives geodesics after projection. Write out the ODEs for geodesic flow
$$\dot{x}^i = p^i $$
$$\dot{p}^i = -\Gamma^i_{jk}p^j p^k$$
Now define $p_i = g_{ij}p^j$ and calculate
\begin{align*}\dot{p}_i &= g_{ij,k}\dot{x}^k p^j + g_{ij}\dot{p}^j \\&= g_{ij,k}p^j p^k - g_{il}\Gamma^l_{jk}p^jp^k
\\ &= p^j p^k(g_{ij,k} - g_{ij,k} + \frac{1}{2}g_{jk,i}) \end{align*}
Thus the correct vector field for geodesic flow on $T^* M$ is in fact
\begin{align*}X &= g^{ij}p_j \frac{\partial}{\partial x^i}  + \frac{1}{2}\frac{\partial g_{jk}}{\partial x^i}  p^j p^k \frac{\partial }{\partial p_i}\\ &= g^{ij}p_j \frac{\partial}{\partial x^i}  - \frac{1}{2}\frac{\partial g^{jk}}{\partial x^i}  p_j p_k \frac{\partial }{\partial p_i} \end{align*}
by raising and lowering indices, and using the derivative of the inverse metric. Now this agrees with $-dH$ in the question, after contracting with $\omega$.
