Я пытаюсь решить систему двух дифференциальных уравнений второго порядка в двух разных областях пространства:
y1 '' = f (x, y1) for -x1 <= x <= 0 </p>
y2 '' = g (x, y2) для 0 <= x <= x2 </p>
с известными граничными условиями y1 (-x1) и y2 (x2), и я также знаю, что y1 ' (x = 0) = y2 '(x = 0) и y1 (x = 0) = y2 (x = 0), т.е. функции y и их производные должны быть непрерывными по x = 0. (На самом деле должна быть только одна функция y, но для ясности я разделил ее на две части).
В Matlab есть очень простой способ решить эту проблему; решатель matlab bv4 c включает эту возможность.
Я не могу понять, как (если это возможно) решить эту проблему с помощью решателя scipy bv4pv. Кто-нибудь может мне помочь с этим вопросом? Спасибо.