ВЫЧИСЛИТЕЛЬНАЯ МАТЕМАТИКА
УДК 519.642.5
DOI: 10.18101/2304-5728-2019-2-3-15
МНОГОШАГОВЫЕ МЕТОДЫ ДЛЯ ЧИСЛЕННОГО РЕШЕНИЯ ИНТЕГРО-АЛГЕБРАИЧЕСКИХ УРАВНЕНИЙ ИНДЕКСА ДВА1
© Будникова Ольга Сергеевна
кандидат физико-математических наук, доцент, Иркутский государственный университет Россия, 664011, г. Иркутск, ул. Нижняя Набережная, 6 E-mail: osbud@mail.ru
© Ботороева Мария Николаевна
старший преподаватель,
Иркутский государственный университет
Россия, 664011, г. Иркутск, ул. Нижняя Набережная, 6
E-mail: masha888888@mail.ru
Многие процессы, протекающие в различных экобиологических и физических системах, описываются взаимосвязанными интегральными уравнениями Вольтерра I и II рода. Их можно записать в виде системы с тождественно вырожденной матрицей перед главной частью — интегро-алгебраического уравнения. В статье исследован класс линейных интегро-алгебраических уравнений, для которого в терминах матричных полиномов сформулированы доста -точные условия существования единственного непрерывного решения. Отмечено принципиальное отличие рассматриваемых задач от интегральных уравнений Вольтерра I и II рода. Работ по качественной теории интегро-алгебраических уравнений мало, а численные методы практически не развиты. Так как многие методы, разработанные для численного решения интегральных уравнений Вольтерра, либо принципиально не применимы, либо приводят к расходящимся процессам. Для выделенного класса задач разработаны одно- и двухшаговые методы, основанные на модификации методов типа Адамса. В качестве подтверждения эффективности алгоритмов приведены результаты численных экспериментов.
Ключевые слова: многошаговые методы; интегро-алгебраические уравнения; квадратурные формулы; аппроксимация; индекс; матричный полином.
Для цитирования:
Будникова О. С., Ботороева М. Н. Многошаговые методы для численного решения интегро-алгебраических уравнений индекса два // Вестник Бурятского государственного университета. Математика, информатика. 2019. № 2. С. 3-15.
1 Исследование частично поддержано грантом РФФИ № 18-51-54001_Вьет_а
3
Введение
Интегро-алгебраические уравнения (ИАУ) — это системы интегральных уравнений типа Вольтерра с тождественно вырожденной матрицей перед главной частью. В последние годы ИАУ активно используются для моделирования развивающихся систем (см., например, [2; 9]).
В статье [11] впервые были проведены качественные исследования на предмет существования и единственности решения ИАУ и предложен численный метод их решения первого порядка точности. В публикации [17] было дано понятие индекса (характеристика сложности рассматриваемых задач), по невязке по своей сути аналогичное понятию степени некорректности для интегральных уравнений Вольтерра I рода [2]. В работах [18; 19] рассмотрено применение методов Рунге — Кутты и полиномов наилучшего приближения для численного решения полуявных ин-тегро-алгебраических уравнений. Есть работы, посвященные качественному исследованию и численному решению двумерных ИАУ (см. [15; 22]) и ИАУ со слабой особенностью в ядре (см. [16; 21]). Кроме того, ИАУ посвящено несколько разделов монографии [13]. В статье [3] рассмотрены ИАУ с переменными пределами интегрирования. Отметим, что подавляющее число работ посвящено исследованию ИАУ индекса один. Зачастую методы, разработанные для ИАУ индекса один либо принципиально не применимы для ИАУ индекса два, либо приводят к расходящимся процессам.
Данная статья посвящена построению численных методов решения ИАУ индекса два и продолжает цикл исследований по численному решению ИАУ многошаговыми методами (см. [4; 5; 8]).
Постановка задачи
Будем рассматривать систему интегральных уравнений вида
t
A(t)x(t) + J K(t, s)x(s)ds = f (t), 0 < s < t < 1. (1)
0
Предполагается, что элементы (n x n) -матриц A(t) и K(t, s) и n -мерной известной вектор-функции f (t) обладают необходимой степенью гладкости.
Непрерывную n -мерную вектор-функцию x(t) будем называть решением задачи (1), если при ее подстановке уравнение обратится в тождество.
Отметим, что при A(t) ° 0 система (1) является системой интегральных уравнений Вольтерра (ИУВ) I рода. Если det A(t) Ф 0, "t е [0,1], то (1) — это система ИУВ II рода. Если определитель матрицы A(t) вырождается на дискретном множестве точек t}- е [0,1], то (1) — это системы ИУВ III рода.
Мы рассматриваем случай
A(t) Ф O,det A(t) ° 0, (2)
где O — нулевая матрица.
Системы (1) с условием (2), продолжая классификацию, называют системами ИУВ IV рода. В первых работах такие задачи называли интегральными аналогами дифференциально-алгебраических уравнений [11] или вырожденными системами интегральных уравнений Вольтерра [6]. В настоящее время рассматриваемые задачи принято называть интегро-алгебраическими уравнениями.
При исследовании дифференциально-алгебраических уравнений В. Ф. Чистяковым было введено понятие индекса, которое можно использовать в качестве характеристики сложности выделенного класса задач.
Определение 1. [12] Минимальное число m , при котором существует дифференциальный оператор
W -W'(t ( dr>
где Wj — (n x n ) -матрицы с непрерывными элементами, такой, что
A(t )x(t )+J K (t, s) x(s)ds = A(t ) x(t ) + J K (t,
s)x(s)ds,
здесь A(t), K(t, 5) — некоторые матрицы с непрерывными элементами, причем
det A(t) Ф 0, "t е [0,1] назовем индексом системы (1).
Понятие индекса аналогично понятию степени некорректности для интегральных уравнений Вольтерра I рода, введенному А. С. Апарциным. Так, уравнение Вольтерра I рода с ядром K(t, t) Ф 0, "t е [0,l] имеет индекс один, а ИУВ I рода с K(t, t) = 0, но K\ (t, t) Ф 0, "t е [0,l], — индекс два [2].
Для получения достаточных условий существования и единственности непрерывного решения задачи (1) приведем необходимое определение. Определение 2. [14] Матричный полином
I2 B(t) + 1C (t) + D(t), где B(t), C(t), D(t) — (n x n) -матрицы, имеет простую структуру на отрезке [0,1], если для любого t е [0,1] выполнены условия:
1) rankB(t) = r = const;
2) rank k[B(t )| C (t)] = r +1 = const;
3) det(1 B(t) + 1C(t) + D(t)) = a0(t)12r+l + a1(t)12r+l-1 +...,a0(t)ф 0.
В работах [7; 14] описаны свойства матричных полиномов, имеющих простую структуру, и их приложения к исследованию различных классов уравнений. Приведем наиболее важные из них.
о
о
Лемма 1. [14] Если матричный полином ) + 1С+ ) на отрезке [0,1] имеет простую структуру и элементы матриц В^), С(^), ) принадлежат классу функций С[01], то существуют невырожденные матрицы Р^) и Q(t) с элементами из С[0 1], такие, что
Г E.
= 12
P(t )(l2 B(t) + 1C(t) + D(t ))Q(t) =
0 0 ^ (C11 (t) 0 C13 (t)| Г D11 (t) D12 (t)
0 0
0 0 0 0
+1
+
0
D21 (() D22 (() 0 0 0 E
\
n-r-l J
0 Е1 0 0 0 0
\ / \
где Сп (t), С13 (t), О11 ^), О12 (t), О21 (t), О22 (t) — некоторые блоки подходящей размерности, а Ег — единичная матрица размерности г.
Лемма 2. [7] Если матричный полином I2В(0 + 1С^) + ) на отрезке [0,1] имеет простую структуру, то матрицы Р(^В(0, Р^С(0, Р(0О(0 имеют блочный вид:
(0| ГС (0| ГО (^
P(t )в (t )=
0 0
, p (t )c (t )=
C2 (t) 0
P(t )D(t ) =
D2 (t) D3 (t)
' B1 (t Г
det C2 (t) Ф 0 "t e [0,1].
D3 (t)
V v 'J
Здесь B1 (t), C1 (t), D1 (t) — матрицы размерности (r x n), C2 (t), D2 (t) —(l x n) -матрицы, D3 (t) есть ((n - r - l)x n)-матрица.
Существование и единственность решения поставленной задачи
На основе лемм 1 и 2 можно получить условия на входные данные, при выполнении которых ИАУ (1) индекса два имеет единственное непрерывное решение.
Утверждение 1. Пусть для задачи (1) с условием (2) справедливо:
1) элементы A(t), f (t)e C^ , K(t, 5)e C[+2, A = {0 < 5 < t < 1}, p > 2;
2) матричный полином 1 A(t) + IK(t,t) + K't(t,t) имеет простую структуру на отрезке [0,1];
3) rank A(0) = rank [a(0 ) f (0)];
4) rank A ( 0) = rank [ A( 0)|( f '(0 )-( A'(0 ) + K (0,0 )) x (0 ))].
и
Тогда на отрезке [0,1] задача (1) имеет единственное непрерывное решение х() из класса Ср^.
Доказательство непосредственно вытекает из равенства
t
A(t)x'(t) + (A'(t)+K(t, t))x(t)+J K"(t, 5 )x(s )ds = f'((),
(3)
полученного из (1) дифференцированием по t, лемм 1, 2 и аналогично доказательству из [7].
Сформулированная теорема носит достаточный характер, прокомментируем ее условия. Первое условие — гладкость входных данных — является необходимым и используется при проведении доказательства. Третье и четвертое условия также являются необходимыми, так как гарантируют разрешимость системы в начальной точке. А вот второе условие носит лишь достаточный характер и гарантирует нам отсутствие сингулярных точек в решении. При нарушении этого условия исследование системы (1) затруднено, так как она может иметь множество решений. Для иллюстрации рассмотрим два простых примера.
Пример 1. Система
(t 0 У х^
0 0
t Г 3 "V x(s)>
,>(t )0
--0
2
0 s
y(s )
V" /V-' V" /У 0 имеет множество решений
где с — произвольная постоянная. Выпишем определитель матричного полинома
(.9 3
аег (л2 л^ )+лк (t, t)+, t))=^
Г 0 ^ t2
v 0
t е[0,1]
It--1
0
0
1t
=1t2 - - 12t.
Г ? "(t )\ g (()
Видно, что условие простой структуры нарушено в точке t = 0 . Пример 2. Рассмотрим еще одну систему
(t 1У хМ (-1 0 У х( )У
V0 0{у((I 0 ((- у^)у
g (0) = g'(0)= 0, t е [0,1], которая имеет множество решений
Гх ^ ) = с,
I у (t )=g"(t),
где с — произвольная постоянная.
Здесь
det (l A(()+IK (t, t)+K't (t, t)) = det
Г It — l l ^
0
1
= 1 (t — 1)— l
Условие простой структуры нарушено, так как rank A(t) = 1, а deg (det (l2 A(() + IK (t, t) + Kt(t, t ))) = ¡^ ^
В заключение данного пункта опишем сложности численного решения ИАУ. Рассматриваемые в статье задачи имеют индекс два, то есть система (1) включает в себя взаимосвязанные ИУВ II рода, ИУВ I рода, имеющих степень некорректности один и два, что хорошо демонстрирует следующий пример.
Пример 3.
Г1
о о
0 у У\ (t У
У 2 (( У Уз((у
b 1 0
c 0
(t—s у
Y У: (s У
У 2 (s У Уз(s )
Г / (t )> / (() /з (()
/2 (0) = /з (0) = /з(0) = 0, X е [0,1],
где а, Ь, и с — скалярные параметры. При гладких входных данных это уравнение имеет достаточно гладкое решение.
Для второй и третьей компоненты имеем ИУВ I рода соответственно
J У2 (s) ds = f2 (t) ,
0
t
J(t — ^Уз (s )ds = f3(()
(4)
(5)
В работе [1] показано, что простейшие методы, основанные на квадратурных формулах правых, левых, средних прямоугольников, а также метод трапеций являются сходящимися для модельного уравнения (4). Применение неявных квадратурных формул порядка выше двух точности приводит к неустойчивым процессам. Например, в [20] приведены результаты применения формул Грегори и Симпсона. Позже в [10] было показано, что ряд эффективных явных многошаговых методов для модельного уравнения (5) являются расходящимися для (4).
Обобщая все вышесказанное, можно заключить, что построение методов для численного решения ИАУ индекса два — нетривиальная задача, так как разработанные алгоритмы должны хорошо справляться с решением всех трех типов уравнений.
Численные методы
Зададим на отрезке [0,1] равномерную сетку = ¡И, г = 1,...,N, И = -1.
Опишем, прежде всего, применение явных методов типа Адамса [10; 20] к построению численных методов решения ИУВ I рода.
Для заданной функции g (), с учетом обозначения gi = g (¿г ), имеем
'¡+1 ík+l г 3+1
Ig(т)ит = |g(т^т+ 2 |g(тУТ
0 0 ] = к+1 ^
*к+1 г 0+1
14+1 (Я0,gl,...,gk,т)ит + 2 I4+1 (g^-к,g
3 - к + &]■
т ит =
3=к+1 I, к
=+ 2 И2^- -1 = И2аг+1,^1, (6)
1=0 3=к+1 1=0 1=0
где Ь[+1 (gз-к,gj-к+1,...,gj,() — интерполяционный полином степени к, зависящий от параметров , ), (g;-k+1, +1), ..., ,),
3 = к + 1,...,г .
Для к = 1,2 приведем значения коэффициентов щ+11, которые являются линейными комбинациями Ь и у1:
+1,1 = -
Г 4 Л Г 9 0 27 Л
3 3 9 5 11 23
1 3 2 3 • 1 9 5 16 7 23
2 3 2 2 3 +1,1 = - 9 5 16 12 7 23
3 2 2 23 9 5 16 12 12 7 23
V" ••0 V... ......0
Здесь при к = 1 первый нулевой столбец опущен, так как при приближенном вычислении интегралов для нечетных значений к коэффициент
юг+1,0 = 0 [10; 20].
Формула (6) носит название явной квадратурной формулы Адамса. Методы численного решения ИУВ I рода, основанные на (6), с учетом обозначений
4+1 = 4^+1 \ К+1,1 = К(¿¡+1, ^ ), /¡+1 = /(¿¡+1), х+1 » х (+1)
имеют вид:
И2т'+1.1К'+1.1х1 = ^г=k,...,N-1.
(7)
1=0
0
Вернемся к исходной задаче. Для численного решения (1) предлагается применять методы, основанные на явной квадратурной формуле Адамса, поскольку данные методы устойчивы для ИУВ I рода [10]. Но их непосредственное применение приводит к проблемам решения вырожденной СЛАУ:
г
Д+Л+1 + ^®г+1,Дг+1,Л = /г+Ъ г = ^ - 1
I=0
Приведем модификацию методов Адамса [4]. Будем вычислять первое слагаемое как значение интерполяционного полинома степени к
Пк+1(хг-к, X-к+1,..., X ) в точке ( = ^+1.
С учетом вышесказанного предлагаемые методы будут иметь вид:
к
Аг+1^а]хг-} + *1£щ+иКм,1*1 = ^ г = k,...,Я-1. ]=0 1=0 Например, одно- и двухшаговые методы имеют вид:
(8)
A+i(2X - X-i) + h^W+iA,+i,ixi = f,+i, l=0
i
A+1(3X - 3x,-1 + Xi-2)+ hYjW+1,lKi +1,lXl = fi+1.
l=0
Приведем результаты численных расчетов предложенными методами на нескольких тестовых ИАУ. Пример 4. Данная система
(
1 t t2 Y x1(t)'
t t2 t3
0 0 0
x.
(()
/
V x3(()0
Г
2tet
' et - 5 0
e"25 et+5
0 V 0
) + te- +12e
- 1)e - 3 +1 e
V-
0 e"25 ((- 5 )
Л
Ж )
X2 (5 ) V x3(5 )0
, t e [0,1]
удовлетворяет всем условиям утверждения 1 и имеет точное решение:
x(t )=(ef e - e 2 ). В таблице 2 приведены погрешности вычисления по евклидовой норме err и порядок метода p , определенный по формуле
Р = l0§2
em
err,
h/2
0
2
Таблица 2
к Одношаговый метод Двухшаговый метод
егг Р егг Р
0.2 0.7363 1.95 0.20987 2.62
0.1 0.1903 1.94 0.034186 2.92
0.05 0.0497 1.96 0.00452 2.89
0.025 0.0128 2 0.00061 2.93
Пример 5. [6] Система
(1 0У X! М 'г(0 1 У х (,)>! (/ (х)>
10 0^X2
(())
с
а(х - ^Х2 )0 ^/2 (()_
Ф 1, X е [0,1].
имеет единственное решение:
(
х(( ) =
/1 (()-г1. / (()-/'(()
1 - а
А (/()- /2(0)
1 - а
при любых / (х)е С1,/2(х)е С2, таких, что /2(0)= 0,/(0) = /2'(0) .
Результаты численных экспериментов при й = 5, / (х) = 1 - ех, /2 (х) = X2 приведены в таблице 3.
Таблица 3
к Одношаговый метод Двухшаговый метод
егг Р егг Р
0.2 0.10523 1.86 0.01667 2.76
0.1 0.02893 1.93 0.002459 2.69
0.05 0.00757 1.96 0.00038 1.97
0.025 0.00194 1.06 0.000097 2.01
Результаты численных расчетов показывают, что предложенные методы являются сходящимися. Обладают порядком точности к .
Заключение
В работе получены достаточные условия, при выполнении которых существует единственное непрерывное решение интегро-алгебраических уравнений индекса два. Разработаны численные методы решения на основе экстраполяции главной части и модификации явной квадратурной формулы типа Адамса. Эффективность данных алгоритмов подтверждена на ряде тестовых задач.
В дальнейшем планируется качественное исследование более сложного класса задач и построение методов их численного решения.
Литература
1. Апарцин А. С., Бакушинский А. Б. Приближенное решение интегральных уравнений Вольтерра I рода методом квадратур // Дифференциальные и интегральные уравнения. Иркутск: Иркутский государственный университет, 1972. Вып. 1. С. 248-258.
2. Апарцин А. С. Неклассические уравнения Вольтерра I рода: теория и численные методы. Новосибирск: Наука. Сиб. изд. фирма РАН, 1999. 193 с.
3. Ботороева М. Н., Булатов М. В. Приложения и методы численного решения одного класса интегро-алгебраических уравнений с переменными пределами интегрирования // Известия Иркутского государственного университета. Сер. Математика. 2017. № 20. С. 3-16. DOI: 10.26516/1997-7670.2017.20.3.
4. Будникова О. С., Булатов М. В. Численное решение интегро-алгебраических уравнений многошаговыми методами // Журнал вычислительной математики и математической физики. 2012. Т. 52, № 5. С. 829-839. DOI: 10/1134/S0965542512050041.
5. Будникова О. С. О модифицированных многошаговых методах для численного решения линейных интегро-алгебраических уравнений индекса два // Журнал Средневолжского математического общества. 2014. Т. 16, № 1. С. 45-54.
6. Булатов М. В. О преобразовании вырожденных систем интегральных уравнений типа Вольтерра // Вычислительные технологии. 2000. Т. 5, № 4. С. 22-30.
7. Булатов М. В., Чистякова Е. В. Об одном семействе вырожденных интег-родифференциальных уравнений // Журнал вычислительной математики и математической физики. 2011. Т. 51, № 9. С. 1665-1673. DOI: 10.1134/S0965542511090065.
8. Булатов М. В., Будникова О. С. Численное решение интегро-алгебраических уравнений со слабой особенностью в ядре k-шаговыми методами // Известия Иркутского государственного университета. Сер. Математика. 2015. № 13. С. 3-15.
9. Глушков В. М., Иванов В. В., Яненко В. М. Моделирование развивающихся систем. М.: Наука, 1983. 350 с.
10. Тен Мен Ян. Приближенное решение линейных интегральных уравнений Вольтерра I рода: дис. ... канд. физ. мат. наук. Иркутск, 1985. 215 с.
11. Чистяков В. Ф. О сингулярных системах обыкновенных дифференциальных уравнений и их интегральных аналогах // Функции Ляпунова и их применения. Новосибирск: Наука, 1987. С. 231-239.
12. Чистяков В. Ф. Алгебро-дифференциальные операторы с конечномерным ядром. Новосибирск: Наука, 1996. 280 с.
13.Brunner H. Collocation Methods for Volterra Integral and Related Functional Equations. Unversity Press. Cambridge, 2004. 612 p.
14. Bulatov M. V., Lee M. G. Application of matrix polynomials to the analysis of linear differential algebraic equations of higher order // Differential Equations. 2008. Vol. 44. P. 1353-1360.
15. Bulatov M. V., Lima P. M. Two-dimensional integral-algebraic systems: Analysis and computational methods // Journal of Computational and Applied Mathematics. 2011. Vol. 236, № 2. P. 132-140. DOI: 10.1016/j.cam.2011.06.001.
16. Bulatov M. V., Lima P. M., Wienmuller E. B. Existence and Uniqueness of Solutions to Weakly Singular Integral-Algebraic and Integro-Differential Equations // Central European Journal of Mathematics. 2014. Vol. 12, № 2. P. 308-321. DOI: 10.2478/s11533-013-0334-5.
17. Gear C. W. Differential-algebraic equations, indices, and integral algebraic equations // SIAM Journal on Numerical Analysis. 1990. Vol. 27, № 6. P. 1527-1534.
18. Hadizadeh M., Ghoreishi F., Pishbin S. Jacobi spectral solution for integral algebraic equations of index-2 // Applied Numerical Mathematics. 2011. Vol. 61, № 1. P. 131-148. DOI: 10.1016/j.apnum.2010.08.009.
19. Kauthen J. P. The numerical solution of integral-algebraic equations of index-1 by pollinomial spline collocation methods // Mathematics of Computation. 2000. Vol. 236. P. 1503-1514.
20. Linz P. Analytical and Numerical Methods for Volterra Equations. SIAM, Philadelphia, 1985. 227 p. DOI: 10.1137/1.9781611970852.
21. Pishbin S. Ghoreishi F., Hadizadeh M. The semi-explicit Volterra integral algebraic equations with weakly singular kernel: the numerical treatments // Journal of Computational and Applied Mathematics. 2013. Vol. 245, № 1. P. 121-132. DOI: 10.1016/j.cam. 2012.12.012.
22. Pishbin S. Numerical solution and structural analysis of two-dimensional integral-algebraic equations // Numerical Algorithms. 2016. Vol. 73, № 2. P. 305-322. DOI: 10.1007/s11075-016-0096-9.
MULTISTEP METHODS FOR NUMERICAL SOLUTION OF INTEGRAL ALGEBRAIC EQUATIONS OF INDEX-2
Olga S. Budnikova
Cand. Sci. (Phys. and Math.), A/Prof., Irkutsk State University
6 Nizhnyaya Naberezhnaya St., Irkutsk 664011, Russia E-mail: osbud@mail.ru
Mariya N. Botoroeva Senior Lecturer, Irkutsk State University
6 Nizhnyaya Naberezhnaya St., Irkutsk 664011, Russia E-mail: masha888888@mail.ru
Many processes in various eco-biological and physical systems are described by interconnected Volterra integral equations of the first and the second kinds. Such equations may be written as a system with an identically singular principal part, in other words, in the form of integral algebraic equation. The article studies the class of linear integral algebraic equations, for which the sufficient conditions for existence of a unique continuous solution are formulated in terms of matrix polynomials. The fundamental difference between the problems under consideration and the Volterra integral equations of the first and the second kinds is noted. There are few works on the qualitative theory of integral algebraic equations, and numerical methods for their solution are underdeveloped. Many methods for numerical solution of Volterra integral equations are not applicable or lead to divergent processes. We have proposed one- and two-step methods based on modifications of Adams-Bashforth and Adams-Moulton formulas for the selected class of problems. The results of numerical experiments confirm the efficiency of the algorithms. Keywords: multistep methods; integral algebraic equations; quadrature rule; approximation; index; matrix polynomial.
References
1. Apartsin A. S., Bakushinskii A. B. Priblizhennoe reshenie integralnykh urav-nenii Volterra I roda metodom kvadratur [Approximate Solution of Volterra Integral Equations of the First Kind by Quadrature Methods]. Differentsialnyie i integralnyie uravneniya. 1972. V. 1. Pp. 248-258.
2. Apartsin A. S. Neklassicheskie uravneniya Volterra I roda: teoriya i chis-lennye metody [Nonclassical Linear Volterra Equations of the First Kind: Theory and Numerical Methods]. Novosibirsk: Nauka Publ., 1999. 193 p.
3. Botoroeva M. N., Bulatov M. V. Prilozheniya i metody chislennogo resheniya odnogo klassa integro-algebraicheskikh uravnenii s peremennymi predelami integriro-vaniya [Applications and Methods for Numerical Solution of a Particular Class of Integral Algebraic Equations with Variable Limits of Integration]. Izvestiya Irkutskogo gosudarstvennogo universiteta. Ser. Matematika. 2017. No. 20. Pp. 3-16. DOI: 10.26516/1997-7670.2017.20.3.
4. Budnikova O. S., Bulatov M. V. Chislennoe reshenie integro-algebraicheskikh uravnenii mnogoshagovymi metodami [Numerical Solution of Integral Algebraic Equations by Multistep Methods]. Computational Mathematics and Mathematical Physics. 2012. V. 52. No. 5. Pp. 691-701. DOI: 10/1134/S0965542512050041.
5. Budnikova O. S. O modifitsirovannykh mnogoshagovykh metodakh dlya chis-lennogo resheniya lineinykh integro-algebraicheskikh uravnenii indeksa dva [On Modified Multistep Methods for Numerical Solution of Linear Integral Algebraic Equations of Index-2]. Zhurnal Srednevolzhskogo matematicheskogo obshchestva. 2014. V. 16. No. 1. Pp. 45-54.
6. Bulatov M. V. O preobrazovanii vyrozhdennykh sistem integralnykh uravnenii tipa Volterra [Transformation of Singular Systems of Volterra-Type Integral Equations]. Vychislitelnye tekhnologii. 2000. V. 5. No. 4. Pp. 22-30.
7. Bulatov M. V., Chistyakova E. V. Ob odnom semeistve vyrozhdennykh integ-rodifferentsialnykh uravnenii [On a Certain Family of Singular Integro-Differential Equations]. Computational Mathematics and Mathematical Physics. 2011. V. 51. No. 9. Pp. 1558-1566. DOI: 10.1134/S0965542511090065.
8. Bulatov M. V., Budnikova O. S. Chislennoe reshenie integro-algebraicheskikh uravnenii so slaboi osobennostyu v yadre k-shagovymi metodami [Numerical Solution of Integral Algebraic Equations with Weakly Singular Kernels by k-step Methods]. Izvestiya Irkutskogo gosudarstvennogo universiteta. Ser. Matematika. 2015. No. 13. Pp. 3-15.
9. Glushkov V. M., Ivanov V. V., Yanenko V. M. Modelirovanie razviva-yushchikhsya sistem [Modeling of Evolutionary Systems]. Moscow: Nauka Publ., 1983. 350 p.
10. Ten Men Yan. Priblizhennoe reshenie lineinykh integralnykh uravnenii Volterra I roda. Dis. ... kand.fiz.-mat. nauk [Approximate Solution of Linear Volterra Integral Equations of the First Kind. Cand. math. and phys. sci. diss.]. Irkutsk, 1985. 215 p.
11. Chistyakov V. F. Chistyakov V. F. O singulyarnykh sistemakh obyknovennykh differentsialnykh uravnenii i ikh integralnykh analogakh [On Singular Systems of Ordinary Differential Equations and Their Integral Analogues]. Funktsii Lyapunova i ikhprimeneniya. Novosibirsk: Nauka Publ., 1987. Pp. 231-239.
12. Algebro-differentsialnye operatory s konechnomernym yadrom [Algebraic Differential Operators with Finite-Dimensional Kernel]. Novosibirsk: Nauka Publ., 1996. 280 p.
13. Brunner H. Collocation Methods for Volterra Integral and Related Functional Equations. Cambridge: Unversity Press, 2004. 612 p.
14. Bulatov M. V., Lee M. G. Application of Matrix Polynomials to the Analysis of Linear Differential Algebraic Equations of Higher Order. Differential Equations. 2008. V. 44. Pp. 1353-1360.
15.Bulatov M. V., Lima P. M. Two-Dimensional Integral-Algebraic Systems: Analysis and Computational Methods. Journal of Computational and Applied Mathematics. 2011. V. 236. No. 2. Pp. 132-140. DOI: 10.1016/j.cam.2011.06.001.
16. Bulatov M. V., Lima P. M., Wienmuller E. B. Existence and Uniqueness of Solutions to Weakly Singular Integral-Algebraic and Integro-Differential Equations. Central European Journal of Mathematics. 2014. V. 12. No. 2. Pp. 308-321. DOI: 10.2478/s11533-013-0334-5.
17. Gear C. W. Differential-Algebraic Equations, Indices, and Integral Algebraic Equations. SIAM Journal on Numerical Analysis. 1990. V. 27. No. 6. Pp. 1527-1534.
18. Hadizadeh M., Ghoreishi F., Pishbin S. Jacobi Spectral Solution for Integral Algebraic Equations of Index-2. Applied Numerical Mathematics. 2011. V. 61. No. 1. Pp. 131-148. DOI: 10.1016/j.apnum.2010.08.009.
19. Kauthen J. P. The Numerical Solution of Integral-Algebraic Equations of Index-1 by Pollinomial Spline Collocation Methods. Mathematics of Computation. 2000. V. 236. Pp. 1503-1514.
20. Linz P. Analytical and Numerical Methods for Volterra Equations. Philadelphia: SIAM, 1985. 227 p. DOI: 10.1137/1.9781611970852.
21. Pishbin S. Ghoreishi F., Hadizadeh M. The Semi-Explicit Volterra Integral Algebraic Equations with Weakly Singular Kernel: The Numerical Treatments. Journal of Computational and Applied Mathematics. 2013. V. 245. No. 1. Pp. 121-132. DOI: 10.1016/j.cam. 2012.12.012.
22. Pishbin S. Numerical Solution and Structural Analysis of Two-Dimensional Integral-Algebraic Equations. Numerical Algorithms. 2016. V. 73. No. 2. Pp. 305-322. DOI: 10.1007/s11075-016-0096-9.