I've been experimenting with the Beltrami-Klein model of hyperbolic geometry. According to Wikipedia, in this model the geodesics are given by Euclidean line segments, and I want to prove this. Let $x^{\alpha}$ denote a geodesic with parameter $t$, then the geodesic is given by parallel transporting the tangent along itself, i.e.;
$$\frac{d^2x^{\alpha}}{dt^2}+\Gamma^{\alpha}_{\ \beta \gamma}\frac{dx^{\beta}}{dt}\frac{dx^{\gamma}}{dt}=0$$
Where $\Gamma^{\alpha}_{\ \beta \gamma}$ are the Christoffel symbols. For this model we have that: \begin{align*} &\Gamma^{1}_{\ 11}=\frac{2x}{1-x^2-y^2}\\ &\Gamma^{1}_{\ 12}=\frac{y}{1-x^2-y^2}\\ &\Gamma^{2}_{\ 12}=\frac{x}{1-x^2-y^2}\\ &\Gamma^{2}_{\ 22}=\frac{2y}{1-x^2-y^2}\\ \end{align*}
With $\Gamma^{\alpha}_{\ \beta \gamma}$ symmetric on the lower indices and all other symbols zero. This yields 2 second order ODES for the components of the geodesic:
$$\frac{d^2x}{dt^2}+\frac{2}{1-x^2-y^2}\left(x\left(\frac{dx}{dt}\right)^2+y\frac{dx}{dt}\frac{dy}{dt}\right)=0$$
$$\frac{d^2y}{dt^2}+\frac{2}{1-x^2-y^2}\left(y\left(\frac{dy}{dt}\right)^2+x\frac{dy}{dt}\frac{dx}{dt}\right)=0$$
I have no idea how to solve these equations and any help yould be much appreciated. I tried a simple case however: Consider the geodesic joining $(x,y)=(0,0)$ to $(x,y)=(\frac{1}{2},0)$. If this is a Euclidean line as stated on wikipedia, it should be given by $x=t$, $y=0$, for $0\leq t\leq \frac{1}{2}$. Plugging this into the system of ODES yields:
$$\frac{2t}{1-t^2}=0$$
$$0=0$$
So clearly this is not a solution! I thought that if Euclidean line segments were geodesics then they would satisfy the geodesic equation.
I would really appreciate any help in solving the general ODES for the geodesic joining two points, as well as in understanding why the Euclidean line given above doesn't appear to be a geodesic.