Диффуры 2го порядка и Метод Рунге-Кутта

[hou]

Отделено от чужой темы. //AKM

k1 = h*p.Fyl(l_temp, yp_temp);
k2 = h*p.Fyl(l_temp + h / 2d, yp_temp + h * k1 / 2d);
k3 = h*p.Fyl(l_temp + h / 2d, yp_temp + h * k2 / 2d);
k4 = h*p.Fyl(l_temp + h, yp_temp + h * k3);
d = (k1 + 2 * k2 + 2 * k3 + k4) / 6d;

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

как лучше выйти из ситуации и какой метод применить вторым.

[AKM]

Тема перемещена из форума «Помогите решить / разобраться (М)» в форум «Карантин»
Тема перемещена в Карантин по следующим причинам: невнятное описание задачи, не оформлены формулы.

Исправьте все Ваши ошибки и сообщите об этом в теме Сообщение в карантине исправлено.
Настоятельно рекомендуется ознакомиться с темами Что такое карантин и что нужно делать, чтобы там оказаться и Правила научного форума.