Дано: прямая выходит из начала координат под углом

и приходит в точку

под углом

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

Реальные исходные данные, для которых нужно получить кривую:

Казалось бы, простая задача, да не тут-то было.
Первая попытка решения - кубическое уравнение

. Четыре неизвестных, четыре начальных уравнения. Всё бы хорошо, но для нужных исходных данных кривая имеет точку перегиба на нужном промежутке [0; 655], вначале она идет выпуклая, потом вогнутая. Это не подходит, вся кривая должна быть ниже прямой, соединяющей (0;0) и

.
Я попробовал уравнение четвертой степени, надеясь, что доп. параметр позволит варьировать место точки перегиба, но нет - ручной подбор показал, что точка перегиба остается в одном и том же месте. Параметры уравнения пятой степени подбираются плохо, не удается найти решение.
Второе, что пришло в голову - искомая кривая есть дуга эллипса с неизвестными координатами центра

и неизвестными полуосями.

Аналитически в декартовой системе координат уравнения не решаются, пробовал численно в MathCAD, точное решение не находит, только приближенное (сильно уходит от заданного угла

), а иногда только в мнимых числах.
Неужели бывают случаи, что невозможно достроить эллипс, исходя из наших начальных условий?
Интуитивно кажется, что построить можно всегда, да вот не выходит. Может искать кривую в других координатах?
Пробовал логарифмическое уравнение:

, аналитическое решение которого показывает, что я имею слишком много начальных условий и система решений не имеет. Лишний параметр

.
Может есть какие-то идеи? В каком виде поискать кривую?
Может как-то графически?
Как это делают сопроматчики и механики?