У нас есть поверхность, полученная склейкой двух поверхностей. Необходимо найти геодезическую для двух заданных краевых точек.
Отдельно находим геодезическую для обеих частей поверхности, подставляя соответствующие начальные данные и получаем два уравнения вида x=f(y,x*), где x* - неизвестная координата по иксу на стыке двух частей поверхности и, соответственно, на стыке двух геодезических.
Получается, нужно минимизировать длину относительно x* для обеих геодезических сразу. а?
Дана поверхность:
![$ x=x, y=y, x \in [0,1], y \in [0,2] $ $ x=x, y=y, x \in [0,1], y \in [0,2] $](https://dxdy-01.korotkov.co.uk/f/8/d/0/8d015b85b9f6ef7d917523fdebc08cce82.png)
, при
![$ y \in [0,1] , z = \sqrt {y}$ $ y \in [0,1] , z = \sqrt {y}$](https://dxdy-02.korotkov.co.uk/f/5/c/a/5ca5877ca7907d718d548cb476442eaf82.png)
, при
![$ y \in [1,2] , z = \sqrt {2-y}$ $ y \in [1,2] , z = \sqrt {2-y}$](https://dxdy-03.korotkov.co.uk/f/a/2/d/a2d281566e33a5dfd8c5e8e2f79c489a82.png)
.
Точки на поверхности, между которыми нужно найти геодезическую:

, точка

находится на стыке поверхностей и геодезических, координата по иксу неизвестна.
Для каждой поверхности в отдельности выписали уравнение Эйлера и нашли геодезические:
![$x= \left ( \dfrac {a-x^*}{1.3} \right ) \left [- \ln{ \left ( \sqrt { \dfrac{4y+1}{y} } +2 \right ) } + \dfrac{1}{8} \ln{ \left ( \sqrt { \dfrac{4y+1}{y} } -2 \right ) } - \dfrac{ \sqrt { \dfrac{4y+1}{y} } }{\dfrac{2(4y+1)}{y} -8} \right ] + 2.4a-1.4x^*$ $x= \left ( \dfrac {a-x^*}{1.3} \right ) \left [- \ln{ \left ( \sqrt { \dfrac{4y+1}{y} } +2 \right ) } + \dfrac{1}{8} \ln{ \left ( \sqrt { \dfrac{4y+1}{y} } -2 \right ) } - \dfrac{ \sqrt { \dfrac{4y+1}{y} } }{\dfrac{2(4y+1)}{y} -8} \right ] + 2.4a-1.4x^*$](https://dxdy-02.korotkov.co.uk/f/9/e/0/9e0bc10667d205dfb7d4b50b092fd22d82.png)

Соответственно, для первого куска поверхности и для второго. Здесь уже константы при игтенрировании выражены и подставленны из начальных условий. Если брать общую длинну этой кривой - получится сумма интегралов, в каждом из которых можно будет отделить терминант, не зависящий от

и проинтегрировать. Но как минимизировать сумму интегралов с разными пределами, где

, в сущности, константа - я не знаю.
Единственная мысль - приравнять функции и подставить

- тогда можно будет найти

- с той точки зрения, что длинна всей кривой - суммы интегралов - будет минимальна, когда длинна каждого куска геодезической - каждого интеграла - будет минимальной, а экстремум здесь, по идее, достигается если приравнять производную к 0 - т.е. саму функцию ( и первую и вторую). Координата по

, по идее, должна сохраниться. Но если по-отдельности приравнивать к 0 каждую из геодезических, могут получится разные

- да они и получатся разными, поэтому, возможно, стоит просто приравнять их друг к другу и выразить

при

?
И еще одна мысль: куски поверхности - симметричные, геодезические, соответственно, тоже будут симметричными, а в координатах точек, между которыми нужно найти геодезическую, идет сдвиг на

- так может быть искомый минимум

будет просто-напросто серединой отрезка
![$[a,b]$ $[a,b]$](https://dxdy-04.korotkov.co.uk/f/f/e/4/fe477a2781d275b4481790690fccd15f82.png)
? Т.е.
А если найти

при развертке как пересечение прямой, соединяющей a и b, и прямой

? Найденная

будет удовлетворять уравнениям геодезических? Т.е. это точка, конечно, обеспечивает по логике вещей кратчайшее расстояние между a и b, но будет ли она минимизировать геодезические? И как проверить, есть ли там минимум?