1) Вызывают сомнения символы Кристоффеля. Вы напрасно их записали в виде матрицы. Так как у них три индекса, при двух координатах всего 8 компонент. С учетом симметрии

— 6 независимых компонент. Поэтому трудно понять, где у Вас какие компоненты, относятся ли обе матрицы к кристоффелям второго рода (которые, собственно, и входят в уравнение геодезической) или первая к первому роду, вторая ко второму. Потом, в нужных кристоффелях

в знаменателе, но по Вашим формулам трудно понять, так у Вас получилось или нет.
У меня получилось следующее:








2) Уравнение геодезической записывается вовсе не для произвольного параметра, поэтому просто так вместо

подставить

или

нельзя. Такой вид оно имеет лишь для канонического параметра, в нашем случае это натуральный (равный длине дуги кривой от выбранного начала) или полученный из такового линейной заменой. Нам же надо получить связь между

и

, поэтому

не нужен абсолютно, и от него надо избавиться. Но при переходе к другому параметру уравнение изменяет свой вид:

,
где

. Это уравнение ещё вывести надо, хоть это и несложно.
3) Надо догадаться из геометрических соображений, что в качестве независимого параметра надо взять именно

. Далее всё подставляем, упрощаем, решаем. В конечном счете всё сводится к одному уравнению

, решением которого будет полином второй степени — заветная парабола. Это писанины ещё где-то на страницу.
В общем, на каждом этапе почти наверняка понадобятся ещё подсказки.