@Vladimir_RU

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

Задача вначале решается программированием процесса по шагам с использованием функции интегрирования rkfixed до момента, когда скорость станет отрицательно, то есть тело движется в обратную сторону. Все работает.
Теперь из этого дела нужно сделать функцию, зависящую от параметра (величины начального импульса) и вычисляющую координату в конце интегрирования. Теперь находя корень уравнения Lкон(Рнач)-Lзад через процедуру root можно найти искомую величину импульса.
Но вот тут затык. Правые части системы дифф.уравнений вычисляются и параметр в них вычисляется в самом начале и непонятно как его сделать формальным параметром функции правых частей.

С уважением к прогрессивной части чеоэчества.

0

@Symon

Нужен файл с вашими наработками.

0