For a free particle in curved spacetime with signature $(-,+,+,+)$, consider the Hamiltonian $$H = \frac{1}{2}\left(g^{\mu\nu}p_\mu p_\nu + m^2\right).\tag{0}$$ Hamilton's equations are given by:
$$\dot{x}^\mu = \frac{\partial H}{\partial p_\mu} \qquad \dot{p}^\mu = - \frac{\partial H}{\partial x^\mu}~. \tag{1}$$
I want to show that this is equivalent to the geodesic equation
$$\frac{d^2 x^\mu}{d\tau^2} + \Gamma^\mu_{\alpha \beta}\dot{x}^\alpha \dot{x}^\beta = 0 =\dot{p}^\mu + \Gamma^\mu_{\alpha \beta}p^\alpha p^\beta~. \tag{2}$$
Direct computation from the Hamiltonian gives $\dot{x}^\mu = p^\mu$. The other equation gives me
$$\dot{p}^\mu = g^{\mu\nu}\dot{p}_\nu = g^{\mu\nu}\Gamma^{\rho}_{\nu\beta}p^\beta p_\rho~. \tag{3}$$
This must be wrong, unless $g^{\mu\nu}\Gamma^{\rho}_{\nu\beta}p^\beta p_\rho = -\Gamma^\mu_{\alpha \beta} p^\beta p^\alpha$, but I am unable to prove this.