Комментарии 5
Странный выбор. Вообще-то, один из самых популярных методов в современных коммерческих пакетах — метод Ньютона-Крылова. Очень быстро сходится.
Вот пример использования scipy.optimize.newton_krylov (даже без preconditioner'а) для решения аж интегро-дифференциального уравнения:
https://docs.scipy.org/doc/scipy/reference/optimize.nonlin.html
Вот пример использования scipy.optimize.newton_krylov (даже без preconditioner'а) для решения аж интегро-дифференциального уравнения:
https://docs.scipy.org/doc/scipy/reference/optimize.nonlin.html
Нехорошо как-то граничные условия Неймана аппроксимированы в одномерной задаче, в двумерной, я так полагаю, кругом Дирихле? Второе, что такое «квадратичная матрице {a}»? Ну и центральные разности для конвекции в сочетании с релаксацией… Не очень понял, кому предназначена эта публикация? Школьникам старших классов? Студентам первого курса?
У Петра Николаевича Вабищевича защищался как раз по этой теме много лет назад. Эх, ностальгия
Зарегистрируйтесь на Хабре, чтобы оставить комментарий
Численные методы решения уравнений эллиптического типа