Памяти Юрия Еремеевича Бояринцева Текст научной статьи по специальности «Математика»



Серия «Математика» 2013. Т. 6, № 2. С. 91—103

Онлайн-доступ к журналу: http://isu.ru/izvestia

УДК 517.518

Памяти Юрия Еремеевича Бояринцева

В. Ф. Чистяков

Институт динамики систем и теории управления СО РАН

Аннотация. В работе дан краткий обзор научной и педагогической деятельности советского и российского математика Юрия Еремеевича Бояринцева.

Ключевые слова: уравнения в частных производных; уравнения переноса нейтронов; некорректные задачи; разностные схемы; устойчивость; сингулярные системы; дифференциально-алгебраические уравнения; индекс.

23 ноября 2012 года ушел из жизни известный советский и российский математик Юрий Еремеевич Бояринцев. Юрий Еремеевич родился 12 мая 1933 года в пос. Полевской Свердловской области. В 1951 году он окончил среднюю школу и поступил на физико-математический факультет Уральского государственного университета. Во время учебы в университете Юрий Еремеевич, по его воспоминаниям, специализировался в разделах математики, связанных с теорией групп.

Однако стране требовались специалисты по численным методам решения уравнений в частных производных. Поэтому в числе выпускников, показавших наилучшие успехи в учебе, после окончания университета в 1956 году он был направлен в «почтовый ящик» на Урале. В советское время это было общее название режимных организаций, в которых производились работы по оборонной тематике. Этот «почтовый ящик» сейчас известен во всем мире под названием Всесоюзный научно-исследовательский институт технической физики. Уральский период биографии Юрия Еремеевича связан с работой под руководством одного из крупнейших специалистов по численным методам решения уравнений в частных производных, будущим академиком АН СССР, Николаем Николаевичем Яненко, о котором он тепло отзывался до конца своей жизни.

1. Краткая биография

Для вычислений использовались первые советские электронно-вычислительные машины (ЭВМ). По воспоминаниям Юрия Еремеевича вычисления на первых ЭВМ требовали не только высокой математической культуры, но и таких качеств, как молниеносная реакция, знание наизусть машинных кодов, позволявших при признаках сбоя в расчетах внести коррективы в работу ЭВМ. Несмотря на жесткий контроль выполнения плановых заданий всячески поощрялись занятия научными исследованиями. Естественно, первые работы Бояринцева посвящены теории разностных схем. В числе первых в мире исследователей он изучал устойчивость разностных схем, получаемых при дискретизации уравнений в частных производных с переменными коэффициентами.

Итогом этих исследований стала защита в 1966 году кандидатской диссертации «О сходимости разностных схем с переменными коэффициентами». После защиты диссертации Юрия Еремеевича назначают заведующим лабораторией в Вычислительном центре СО АН СССР в Новосибирске.

В 1971 году Юрий Еремеевич переезжает в Иркутск и назначается заведующим кафедрой вычислительной математики в Иркутском государственном университете (ИГУ). Через некоторое время он основывает лабораторию вычислительной математики в Сибирском энергетическом институте (СЭИ) СО АН СССР и параллельно продолжает преподавательскую деятельность в ИГУ в должности доцента.

В 1975 году в рамках СЭИ был создан Отдел теории систем и кибернетики (ОТСК) с целью формирования на его базе самостоятельного научно-исследовательского института с функциями вычислительного центра коллективного пользования. Лаборатория вычислительной математики становится подразделением ОТСК с тем же названием. В 1980 году ОТСК преобразуется в отдельный институт - Иркутский вычислительный центр СО АН СССР, перемеинованный позднее в Институт динамики систем и теории управления (ИДСТУ) СО РАН. Лаборатория вычислительной математики перемеиновывается в лабораторию алгебро-дифференциальных систем.

С переездом в Иркутск Юрий Еремеевич полностью меняет тему исследований. Он начинает изучать сингулярные (вырожденные) системы дифференциальных уравнений с целью построения эффективных численных методов для их решения. Так как отсутствовали конструктивные теоремы существования, параллельно создавалась качественная теория вырожденных систем. Итогом этих исследований стала защита в 1983 году докторской диссертации «Вырожденные системы обыкновенных дифференциальных и разностных уравнений».

После защиты докторской диссертации Юрий Еремевич продолжает много и плодотворно трудиться по развитию теории ДАУ и построению численных методов для их решения. В 2002 году он передает руководство лабораторией автору этих строк, переходит на должность

главного научного сотрудника и работает в этой должности до самой своей смерти.

2. Работа на Урале и в Новосибирске

Как выше уже говорилось, в этот период своей деятельности Юрий Еремеевич исследует численные методы решения систем дифференциальных уравнений в частных производных. Результаты исследований разностных схем были опубликованы в Докладах АН СССР [1], [2] и подитожены в развернутой статье [3].

На примере начально-краевой задачи для волнового уравнения

д2п д , ,дп , . , .

= ^ 0<Ко< х) <Х< со, (2.1)

д , \

< L, 0 < t < T < ж,

где к0, х, L, T—фиксированные константы,

д

и(х, 0) = <р(х), 0) = 1р(х), и{0, t) = и{ 1, t) = 0, (2.2)

демонстрируется общая методика. При определенных условиях на гладкость в области [0,1] х [0, T] вводится сетка c шагами т и h. Предполагается, что при измельчении шагов т/h = £ = const. Строится разностная схема, допускающая матричное представление

(E + alBn+l)Un+1 + (—2E + aoBn)Un + (E + a-iBn-1)Un-1 = 0, (2.3) где E—единичная матрица размерности N, (N + 1)h = 1,

Bn = У — б^тП-г + 6j(rn-i + гП) + Sj+1 rn\\N, 0 < пт < T, rn = K((i + 1/2)h, пт )£2,

5j —символ Кронекера, Un — N-мерный вектор-столбец: Un = (U'n)i . Матрицей Bn и вектором Un учтены краевые условия

Uon = UN+1 = 0, U° = p(ih), Ui — Ui = ф(ih).

Центральной идеей доказательства устойчивости и сходимости разностных схем является предположение: если оператор перехода устойчив (например, удовлетворяет критерию устойчивости Неймана) на каждом шагу процесса, то и разностная схема устойчива. На основе изучения свойств корней матричного уравнения

(E + aiBn)X2 + (—2E + aoBn)X + (E + a-iBn) = 0,

которое является некоторым аналогом характеристического уравнения, Бояринцевым доказано следующее утверждение.

94 в. ф. чистяков

Теорема 1. Пусть величина z = 4(а0 — 4aiа_1)_1 не принадлежит промежутку (0, 4R + е), где R = sup{K(x,t)(2, x £ [0,1],t £ [0,Т]}7 е > 0, ив (3) \\Ul2 II < 1 + Ct, C = const.

Тогда решения Ui из схемы (2.3) сходятся к решению задачи (2.1), (2.2) при стремлении t к нулю.

Далее показано, что результаты можно распространить на краевые задачи для систем общего вида

д2п 1 dj . .dj u

0 < к0 < Ki,j (t,x) < x < ж.

Второе значительное достижение сввязано с изучением метода сферических гармоник для уравнения переноса нейтронов (статья [4], написанная совместно с О. П. Узнадзе). Уравнение переноса нейтронов записывается в цилиндрической системе координат (r,z,ф):

П 1

/#// Wi^T'M.^.yi'j^fv.iji), (2.4) о_1

где обозначения взяты из [5]. Функция ф = r,z,t,Y,ф) описывает пространственно-угловое распределение в цилиндре D. Предполагается, что в момент t = 0 известна функция s(r,z,Y,ф) = r,z, 0^,ф), а на поверхность цилиндра Г со стороны вакуума падают нейтроны с распределением f (Г, Q,t), Q—единичный вектор в направлении переноса нейтронов. Для простоты предполагается, что f = 0.

Решение и свободный член представляются в виде ряда

ф, z, t, Y, ф) = J2(2n + 1)Yn(r, z, t, Y, ф), (2.5)

n=0

n

Yn(r,z,t,Y,Ф ) = £ ^m(r,z,t)pnm (Y )

m=0

<x

S(r, z, t, y, ф) =J2(2n + 1)Sn(r, z, t, y, ф), (2.6)

n=0

n

Sn(r,z,t,Y,ф) = £ Sm(r,z,t)Pm(Y),

m=0

где Уп, Бп суть сферические функции, зависящие от г,г,Ь, РП^О)—полиномы Лежандра, а индикатриссу рассеяния запишем в виде

1 те

д(П ^ П1) = «7[(П'П)] = дЫ = — ]Г(2п + 1)дпРпЫ, (2-7)

п=0

где

1

9п — ! д(^0 )Рп -1

С учетом соотношений (2.5)-(2.7) и ортогональности полиномов Лежандра, интегральный член в (2.4) приобретает вид

п 1 ^

V Г Г V те

Js

2тг

/С V ^^

J -»• п')<р(г, г, I, О!) = ^ + 1 )¥п(г, г, 1,7, <ф).

о -1 п=0

Таким образом, относительно коэффициентов разложения ^(г, х-> ¿) вместо интегро-дифференциального уравнения получается бесконечная система дифференциальных уравнений

дУп 1 9 тш+ил

Ш 2г(т+1) дг

(п + т + 1)(п + т + 2) т+1 (п - т - 1)(п - т - 2) т+1 2^+1 ^1 + 2^+1

(1 + $1,т)гГП 1 д \ т-1(, „т-1 ,т-1\ 1 I

2(2п + 1) д~г[г +

д_

дг

п + т + 1 т п — т

1 + ^ггт^-1

I у1 т _ от

+ Vn^n —

2п + 1 Т п^ ' 2п + Г п-п — 0,1, 2, ••• ; т — 0,1.2, ■■■ п, (2.8)

где ¿1;т—символ Кронекера, V п - V Vsgn.

Если ограничиться конечным числом членов разложения, то получается Р^—приближение сферических гармоник (п — 0,1, 2, ■■■ ,И). Чтобы замкнуть задачу к системе (2.8) добавляются начальные граничные условия. Они формируются расщеплением начальных и краевых условий исходной задачи. Доказана соответствующая теорема. Краевая задача для (2.8) решается методом расщепления. Вопрос о сходимости рассмотрен при условиях постоянства шагов дискретизации Дг, Дг и отношения т/Дг, т/Дг фиксированы, где т—шаг дискретизации по времени, все гармоники ^т = 0 на внешней границе цилиндра, а также при предположении, что решение краевой задачи (2.8) в конечной окрестности цилиндрической оси известно.

Юрий Еремеевич всегда интересовался некорректными задачами. Им выполнена (совместно с В. Г. Васильевым) работа [6] по применению метода квазиобращения к задаче теплопроводности с обратным временем. Рассматривалось уравнение

f = -а2§, (ж, t) £ [0,1] х [О, Т], (2.9)

где a = const > 0, T = const > 0, с начально-краевыми условиями,

u(x, 0) = p(x), u(0, t) = u(1,t) = 0. (2.10)

Задача (2.9), (2.10) некорректна. Сколь угодно малым возмущениям входных данных могут соответствовать сколь угодно большие возмущения решения. Метод квазиобращения является некоторым методом регуляризации. Исходной задаче в соответствие ставится уравнение

dv д2 v d4v

а>0' 0М)е[оД]х[о,т], (2.11)

с начально-краевыми условиями,

, ч , ч , ч d2v(0,t) d2v(1,t) v(x, 0) = ф), v(0, t) = v(l,t) = 0, ' = gyJ ' = 0. (2.12)

Известно, что при малом а решения задач (2.9), (2.10) и (2.11), (2.12) близки. Для решения регуляризованной задачи применяется разностная схема

vn+1 - vn _ 2 A2vn+1 A4vn+1

= —а —гт;--а-

Н2 Н4 '

где Дvгra+1 = - = 1,2, ••• (Ж + 1)Н = 1, 0 < пт < Т, т <

4а,

V0 = фh), = V-1 + у1п = + УПы+2-

Оценка оператора шага в Ь2 имеет вид

\\С|| = шах{М(т, к), к = 1,2, ••• ,Щ<1 + 0(т),

где

4Т „ КПП / 4а „ КПП „ \ 1

м(г, к) =

4т 2 КпП (4а 2 Knh 2

1+VSm — USm I" "а

Эта оценка говорит о том, что будут накапливаться ошибки округления. Но от этого дефекта можно избавиться.

Теорема 2. Если На/2а1/2 < 1, то

max{M(т, К), а/2а1/2 < К < N} < 1.

3. Иркутский период

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

где ^ : И2п х Т ^ И, ' — й/дЛ,, х = ж(£) — п-мерная искомая вектор-функция, и удовлетворяющие условию

во всей области определения F. Если v = n, то условие (3.2) эквивалентно условию det (dF/dx) = 0.

Для обозначения объекта исследования в нашей и зарубежной литературе использовались и другие термины: «вырожденные системы ОДУ» «дескрипторные системы», «алгебро-дифференциальные системы (АДС)», «сингулярные системы», «системы обыкновенных дифференциальных уравнений (ОДУ), неразрешенные относительно производных», «вырожденные системы ОДУ» и даже «неявные (implicit) системы ОДУ». В настоящее время термин ДАУ потеснил все другие названия. Это вызвано тем, что в приложениях математическая модель объекта исследования часто имеет вид взаимосвязанных дифференциальных и алгебраических уравнений. В своих работах Бояринцев остановился на термине «АДС».

Интерес к ДАУ и численным методам их решения возник не только как к интересному математическому объекту. Он стимулировался проблемами математического моделирования в прикладных областях, в частности, в теориях электронных схем и электрических цепей, механике и химической кинетике, гидродинамике и теплотехнике. В начале 70-х годов прошлого века в зарубежной литературе были уже опубликованы работы по применению для численного решения ДАУ разностных методов, основанных на формуле дифференцирования назад (ФДН). Были также опубликованы ряд работ в СССР и за рубежом исследователей в прикладных областях с применением разностных схем, используемых для решения ОДУ в нормальной форме. Но эти работы носили экспериментальный характер. Решали системы вида (3.1) с применением известных разностных схем и получали (не получали) приближение к решению. Положение осложнялось тем, что не существовало конструктивных теорем существования.

Юрий Еремеевич начал исследования с простейшего случая - линейного ДАУ с постоянными коэффициентами

F(x, x, t) = 0, t e T =[a, в] С R1,

(3.1)

rank (dF/dx) < min (n, v)

(3.2)

F(x, x, t) = Ax + Bx - =0, t e T,

(3.3)

и простейшей разностной схемы

АУг+1~Уг + Вуг+1 - ф(1г+1) = 0, (3.4)

т

г = 0,1, ••• , N, т = (в - а)/Я, и = а + гт,

и ее первого «дифференциального приближения», полученного при применении к схеме (3.4) методики из [7]

[А + еВ]у£ + В(г)уе - ф(г) = 0, г е Т, е> 0, (3.5)

где А, В -(п х п)-постоянные матрицы, ф(г) —заданная вектор-функция из Сп(Т). Предполагалось, что пучок матриц коэффициентов ХА + В регулярен, т.е. существует число Хо (в общем случае комплексное) и такое, что (ХоА + В) = 0. Краевые задачи ставились в предельно общем виде. Предполагалось, что на решение системы (3.3) наложено условие в виде интеграла Стильтьеса

в

У С(s)x(s)da(s) = а,

а

где С^) - (д х п)-матрица, а-заданный вектор. В частности, в таком виде можно записать любые многоточечные краевые условия.

Выбор был не случайным. Для ДАУ с постоянными матрицами коэффициентов было известно [8], что для любой системы АХ + Вх -ф(г) = 0 с матрицами коэффициентов размерности (и х п) существуют постоянные неособенные матрицы Р, Q подходящей размерности со свойством

Р[А^у + BQy - ф(г)] =

= diag{Ed, Ж, М, Ь, 0}у + 3, Е3, М1 ,ЬЪ 0}у - Рф(г) = 0, (3.6)

где х = Qy, Ed, Еа-единичные матрицы размерности d и s соответственно, М, Ь-прямоугольные матрицы полного ранга и подходящей размерности, начиная с некоторого натурального к: Nк = 0. Выражение

Р(ХА + В^ = Хdiag{Ed, Ж, М, Ь, 0} + 3, Еа, Мх,ЬЪ 0}, (3.7)

называется канонической структурой пучка ХА + В в смысле Веейр-штрасса-Кронекера. Если пучок матриц ХА + В регулярен, то правая часть равенства (3.7) имеет вид

Х^Е Ж} + 3, Еп^}, (3.8)

и число к называется индексом пучка. Иначе говоря формула (3.6) позволяет выписать общее решение ДАУ в конечной форме, сформулировать условия совместности ДАУ и разрешимости начальных и краевых задач.

В работе [9] впервые в известной автору литературе было исследовано влияние структуры пучка матриц коэффициентов \Л + В на поведение численных методов и обнаружено интересное явление, возникающее при применении разностных методов для решения ДАУ, названное позже «пограничным слоем ошибок». Суть этого явления состоит в том, что в первых точках сетки отклонение разностного приближения (3.4) от решения ДАУ стремится к бесконечности с порядком О(1/тк-2), при стремлении шага сетки т к нулю. Здесь к индекс пучка матриц коэффициентов. Статья [9] прояснила ряд принципиальных вопросов и сыграла большую роль в выборе направления исследований всей лаборатории.

В монографии [10] для случая регулярного пучка матриц коэффициентов впервые обоснован «метод возмущения» (3.5). Показано, что «пограничному слою ошибок» соответствует пограничный слой в окрестности г = а вида

Ш) - х(тс{т) = 0(е^/е2к~2), г<ЕТе = [а,а + 0(е)}.

Ь(г) -х(1)\\с(т) = О(е), г / Т.

Практически одновременно были построены общие решения системы (3.1) в конечной форме. Было показано, что

х(г) =

г

еА°гЛ°Ах(0) + [ еА°(—)А°$(з)й8 + Тк-1(Еп - Л°Л)ф(г)

= eXot

где А = (А0А + B) 1 A, ip(t) = e Xot^(t), T—оператор действующий на вектор-функцию ф = p(t) по правилу

d

TV = <р, т-1 = О,Т<р = ~{Аф) + {Еп - ADA)<p,

MD—матрица Дразина, единственным образом определяемая для произвольной квадратной матрицы M по правилу

MMD = MD M, MD MMD = MD, MD Mk+1 = Mk

где k—индекс матрицы M: наименьшее из натуральных чисел, для которого rank Mk = rank Mk+1. Отметим, что индексы матрицы А и пучка АА + B совпадают.

В ходе этих исследований сформировались важнейшие понятия теории ДАУ, такие как индекс ДАУ, допустимые начальные данные, пр-стая структура, особые точки и т.д.

Юрий Еремеевич получил обобщение уравнения Ляпунова.

и

Лемма 1. Для того чтобы, все отличные от нуля собственные собственные числа матрицы А имели отрицательную вещественную часть необходимо и достаточно, чтобы уравнение

А*у + га = -(А° А)*А° А

имело решение У являющейся эрмитовой матрицей, положительно определенной на образе матрицы А° А.

На следующем этапе исследований рассматривались линейные системы уравнений с переменными коэффициентами

Х[А(г)х] + В(г)х = ф(г), г е Т,

А(г)\х + В(г)х = ф(г), г е т,

где А(г), В (г) - (V х п)-произвольные матрицы, А-скалярный оператор, например, дифференциальный, разностный, интегральный. Были построены общие решения. Техника доказательств очень сложна, выражения очень громоздкие. Поэтому мы не можем привести их здесь. Для желающих ознакомиться наилучшим руководством будет докторская диссертация Бояринцева «Вырожденные системы обыкновенных дифференциальных и разностных уравнений», доступ к которой можно получить через Интернет.

Скажем следующее. Кронекерова структура пучка матриц ХА(г) + В (г) неинвариантна относительно преобразований, использующих замену переменных в ДАУ. А именно, пусть мы произвели замену х = Q(г)y. Структуры пучков

ХА(г) + В (г), ХА(г^(г) + [В (г) + А(г^(г)]

могут совершенно различаться. Структура общего решения ДАУ для систем с постоянными коэффициентами полностью определяется кро-некеровой канонической формой пучка. Бояринцевым были выделены классы ДАУ со свойством совершенства и свойством О, для которых неособенные преобразования не меняют типа кронекеровой структуры пучка матриц исходной системы. Для ДАУ, обладающих этими свойствами, и их разностных аналогов выписаны формулы общих решений.

В последние годы жизни Юрий Еремеевич много работал над способами вычисления так назваемых базовых матриц Со, С1, • • • ,Ск. Любое решение системы (1.3) с регулярным пучком матриц при достаточно гладкой ф (г) при достаточно гладкой ф (г) представимо в виде

к / d V-1

х(1) = С0г(1) + Е ( л ) = + Ш (3-9)

г=1 ^ '

В формуле (3.9) дифференциальный оператор можно заменять на разностный и это свойство можно применять для построения численных методов. Базовые матрицы использовались Бояринцевым для исследования нелинейных ДАУ с выделенной линейной частью

АХ + Вх = Ф(х,г).

Сейчас исследователи в области автоматического регулирования начинают использовать этот аппарат для своих целей.

В заключение параграфа отметим еще одно обстоятельство. В монографиях содержатся результаты исследований самого Бояринцева и других авторов по теории обобщенных обратных матриц (включая матрицу Дразина и ее обобщения). Эти монографии являются хорошим и пока единственным в России введением в этот предмет.

4. Педагогическая деятельность и научное руководство

Юрий Еремеевич был широко эрудированным ученым и прекрасным педагогом. В рамках тематики, начало которой он положил, защищено три докторских и несколько кандидатских диссертаций. Бояринцев задавал общее направление исследований и вмешивался только для того, чтобы поправить или дать консультацию. При чтении лекций Бояринцев очень тщательно подходил к отбору учебного материала. Доказывалось немного теорем, но каждая теорема разбиралась «по косточкам». Автору особенно запомнились лекции посвященные теореме Лакса и критерию Неймана.

В иркутский период научной деятельности Юрий Еремеевич предпочитал публиковать результаты своих исследований в виде монографий, где можно подробно и с большим количеством поучительных примеров изложить материал. Им написано (3 из них в соавторстве) и опубликовано в издательстве «Наука» 7 монографий. Книга [11] переиздана на английском языке [12].

По инициативе и участии Бояринцева на математическом факультете ИГУ читались спецкурсы по новым тогда разделам вычислительной математики. Таисия Ивановна Назаренко прочитала курсы лекций «Интервальный анализ», «Численные методы решения ОДУ типа Адамса». Валентина Никитична Васильева прочитала курсы лекций «Теория сплайнов», «Теория некоррекктных задач», включая метод квазиобращения Лионса. Анатолий Соломонович Апарцин прочитал курсы лекций «Численные методы решения интегральных уравнений Вольтерра первого рода». Сам Юрий Еремеевич Бояринцев читал спецкурс «Теория разностных схем (теорема Лакса и критерий устойчивости Неймана, численное решение ДАУ)». Лекции оказались настолько увлекательными, что автор этих строк решил специализироваться в

этой области математики и стал доктором физико-математических наук. Некотрые курсы были настолько полными, что вполне заслуживали издания в виде отдельной книги.

Заканчивая этот обзор, скажем о прекрасных человеческих качествах Юрия Еремеевича Бояринцева. За тридцать пять лет работы, теснейшим образом общаясь с ним, автор ни разу (ни разу!) не слышал грубого слова. Неизменно выдержанный, вежливый, строго одетый, внимательный Юрий Еремеевич был образцом русского интеллигента. Математика была для него жизнью. Автору известно только одно нематематическое увлечение Юрия Еремеевича: классическая музыка. Он был тонким ее ценителем.

Список литературы

1. Janenko N. N. Convergence of difference schemes for the heat equation with variable coefficients / N. N. Janenko, Ju. E. Bojarincev (Russian) // Dokl. Akad. Nauk SSSR. - 1961. - Vol. 139. - P. 1322-1324.

2. Bojarincev Ju. E. On convergence of difference schemes for the waveequation with variable coefficients / Ju. E. Bojarincev (Russian) // Dokl. Akad. Nauk SSSR. -1965. - Vol. 165. - P. 474-475.

3. Бояринцев Ю. Е. О сходимости разностных схем для уравнений с переменными коэффициентами / Ю. Е. Бояринцев // Тр. МИАН СССР. - 1966. - Т. 74. -C. 16-37

4. Бояринцев Ю. Е. Об одном способе решения нестационарного уравнения переноса / Ю. Е. Бояринцев, О. П. Узнадзе // Журн. вычисл. математики и мат. физики. - 1967. - Т. 7, № 6. C. 1406-1413.

5. Марчук Г. И. Методы расчета ядерных реакторов / Г. И. Марчук. - М. : Гос. изд-во лит. в обл. атом. науки и техники. - 666 с.

6. Бояринцев Ю. Е. Об устойчивости метода квазиобращения при решении некорректных эволюционных уравнений / Ю. Е. Бояринцев, В. Г. Васильев // Журн. вычисл. математики и мат. физики. - 9:4 (1969). C. 951-952.

7. Шокин Ю. И. О связи корректности первых дифференциальных приближений и устойчивости разностных схем для гиперболических систем уравнений / Ю. И. Шокин, Н. Н. Яненко // Мат. заметки. - 1968. - Т. 4, № 5. - C. 493-502.

8. Гантмахер Ф. Р. Теория матриц / Ф. Р. Гантмахер. - М. : Наука, 1967. - 575 c.

9. Бояринцев Ю. Е. Применение разностных методов к решению регулярных систем обыкновенных дифференциальных уравнений / Ю. Е. Бояринцев, B. M. Корсуков // Вопросы прикладной математики. - Иркутск : Изд. СЭИ СО АН СССР, 1975.

10. Бояринцев Ю. Е. Регулярные и сингулярные системы линейных обыкновенных дифференциальных уравнений / Ю. Е. Бояринцев. - Новосибирск : Наука, Сиб. отд-ние, 1980. - 222 с.

11. Бояринцев Ю. Е. Методы решения вырожденных систем обыкновенных дифференциальных уравнений / Ю. Е. Бояринцев. - Новосибирск : Наука, 1988. -158 с.

12. Boyarintsev Y. Methods of solving singular systems of ordinary differential equations / Yuri Boyarintsev ; Translated from the 1988 Russian original by Vassily Michalkowski. Pure and Applied Mathematics (New York). A Wiley-Interscience Publication. John Wiley and Sons, Ltd., Chichester, 1992. - 163 p.

13. Численные методы решения сингулярных систем / Ю. Е. Бояринцев, В. А. Данилов, А. А. Логинов, В. Ф. Чистяков. - Новосибирск : Наука. Сиб. отделение, 1989. - 261 с.

14. Бояринцев Ю. Е. Методы решения непрерывных и дискретных задач для сингулярных систем уравнений / Ю. Е. Бояринцев. - Новосибирск : Наука, Сиб. отд-ние, 1996. - 261 с.

15. Бояринцев Ю. Е. Алгебро-дифференциальные системы: Методы решения и исследования / Ю. Е. Бояринцев, В. Ф. Чистяков. - Новосибирск : Наука, 1998. -224 с.

16. Бояринцев Ю. Е. Линейные и нелинейные алгебро-дифференциальные системы / Ю. Е. Бояринцев. - Новосибирск : Наука, 2000. - 222 с.

17. Бояринцев Ю. Е. Пучки матриц и алгебро-дифференциальные системы / Ю. Е. Бояринцев, И. В. Орлова. - Новосибирск : Наука, 2006. - 124 с.

V. F. Chistyakov

In the memory of Bojarintsev Yuriy Yeremeyevich

Keywords: partial differential equations, neutron transport equation, ill-posed problems, difference schemes, stability, singular systems, differential algebraic equations, index.

Abstract. In this paper we give a brief review of the scientific work by Soviet and Russian mathematician Bojarintsev Yuriy Yeremeyevich.

Чистяков Виктор Филимонович, доктор физико-математических наук, главный научный сотрудник, Институт динамики систем и теории управления СО РАН, chist@icc.ru

Chistyakov Viktor Filimonovich, Doctor of Sciences, chief research fellow, Institute for System Dynamics and Control Theory SB RAS