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

Исследуется линейная дифференциально-алгебраическая система уравнений в частных производных первого порядка с произвольным индексом. Понятие индекса системы определяется максимальной в области определения степенью элементарных делителей, соответствующих нулевым и бесконечным корням характеристического многочлена матричного пучка, построенного по коэффициентам дифференциально-алгебраической системы. Численное решение дифференциально-алгебраических систем высокого индекса до сих пор остается нерешенной проблемой. Настоящая работа направлена на ее решение. Предлагается итерационный алгоритм численного решения, основанный на специальном расщеплении матричного пучка и применении к расщепленной дифференциально-алгебраической системе метода последовательных приближений. Доказывается устойчивость алгоритма. На тестовых примерах демонстрируется его эффективность. Библ. 15. Фиг. 1.