Метод блочной матричной прогонки для решения систем обыкновенных дифференциальных уравнений

Abstract

Предложен эффективный алгоритм решения краевой задачи для системы парных линейных обыкновенных дифференциальных уравнений с нелинейностью в правых частях. Рассмотрены два варианта устойчивых конечно разностных схем, приводящих к необходимости решения системы линейных алгебраических уравнений с матрицей, имеющей блочно-ленточную структуру. Для решения таких систем разработана экономичная модификация метода Гаусса с выбором главного элемента — метод f блочной матричной прогонки. Идея алгоритма заключается в реализации метода Гаусса на упакованных массивах, в которые помещаются ненулевые элементы матриц.