У меня разница 2-го порядка. э. : y''(r)-2/r*y'(r)=(p/M**2-u**2)*y(r)+y(r)**3
К сожалению, этому уравнению удовлетворяются 2 разные области для (r<R; p=p) и (r>R; p=0), где y(R) и y'(R) должны совпадать при R для уравнения каждой области. Есть граничные условия, управляющие всем y(r), такие, что y'(0) = 0 и y'(inf) = Const.
Можно ли использовать scipy.solve_bvp для каждого региона отдельно с 3 граничными условиями для каждого региона (или 2 для одного регион и 3 для другого)? Или мне придется использовать какой-то другой численный метод для достижения решения?