Решение граничного уравнения второго порядка на неоднородной сетке

#python #numpy #scipy #numerical-methods #differential-equations

#python #numpy #scipy #численные методы #дифференциальные уравнения

Вопрос:

У меня есть уравнение вида

y'' a(x) y' b(x) y = f(x)
y(0) = y(1) = 1

где x неравномерно разнесено.

Как я могу решить этот тип краевой задачи второго порядка в python?

Комментарии:

1. Знаете ли вы, с чего начать, если сетка была однородной? Если у вас нет опыта в численных методах, возможно, было бы лучше указать вам на онлайн-ресурс для обучения.

2. В качестве альтернативы, знаете ли вы, как это решить, если бы это была проблема с начальным значением y(0) = 1 y'(0) = s для некоторой константы s ? В этом случае вы хотите посмотреть метод съемки .

Ответ №1:

Вы ищете семейство функций [ y x (t) ], где t находится переменная, на которую вы выводите y ?

Верны ли ваши граничные условия для всех x?

Возможно, вы захотите взглянуть на http://docs.scipy.org/doc/scipy/reference/generated/scipy.integrate.ode.html