Быстрое решение модельной задачи для бигармонического уравнения Текст научной статьи по специальности «Математика»

УДК 519.63

DOI: 10.14529/mmph190105

БЫСТРОЕ РЕШЕНИЕ МОДЕЛЬНОЙ ЗАДАЧИ ДЛЯ БИГАРМОНИЧЕСКОГО УРАВНЕНИЯ

А.Л. Ушаков

Южно-Уральский государственный университет, г. Челябинск, Российская Федерация E-mail: ushakoval@susu.ru

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

Ключевые слова: фиктивное продолжение; итерационные факторизации.

Введение

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

Непрерывная задача

Исходная задача:

U е Zi: Л(«1, Vi) = ^ф) "V е Zb g е (1)

где

Zi = Zi (W) = jvi е Ж,2(W): v|Г = 0, г2 = о};

пространство Соболева из функций на области прямоугольного вида:

W = (0;bi)x(0;b2), с Г = {bi}x(0;b2) U(0;bi)x{b,}, Г, ={0}x(0;b2)U (0;bi)x{0},

а

L(Ui, Vi) = J AuiAVid W

W

Ушаков А.Л. Быстрое решение модельной задачи

для бигармонического уравнения

- билинейная форма, где заданные константы Ъ1, Ъ2 > 0. Решение задачи (1) существует и единственно [2, 3]. При

&1 (V ) = | П,

п

где /1 - заданная функция, суммируемая с квадратом, задача из (1) записывается в виде:

Л^ = /1, щ\Г1 =Ащ |Г1 = 0, Ц|Г2 ^Г2 = 0, (2)

где п - внешняя нормаль на границе области ЭП .

Дискретная задача и ее продолжение

Рассматривается линейная система алгебраических уравнений, получаемая при дискретизации (1), (2) на основе метода сумматорных тождеств:

и1 е : Щ = /1, / е , (3)

здесь векторы:

V е Мы: V = (v11,...,у1 ы)', N = т • п, т , п е N,

полагается, что

^,т(у-1)+, = v1,i,у , * = 1,...,т , ] = 1,...,П ,

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

(хг, у у) = ( (/ - 0,5)Л1,( у - 0,5)^2), когда шаги сетки выбираются следующие:

Н1 = Ъ1/(п + 0,5), Л2 = Ъ2/(п + 0,5), сетка состоит из узлов указанных выше, а матрица Л размерности NX N определяется так:

т п

(М^V!) = ХЕ ((2иМ,у - ы1,»-1,у - ы1,»+1,у )А1-2 + (2иМ,у - ы1,»,у-1 -

;=1 у=1

-и1,»,у +1 )((2^,/,у - v1,i-1,у - v1,i+1,у )А1-2 + (2\,,у - v1,i,у-1 - v1,i,у+1 ,

и1,г,и+1 = v1,i,и+1 = 0 и1,1,0 = u1,i,1, ^'1,1,0 = VU,1, ' = 1,...,^ и1,т+1,у = ^1^+1,у = 0 и1,0,у = и1,1,у , V1,0,у = у , ] = 1,...,п

Здесь (.,.) - скалярное произведение следующего вида:

N

"\,ку1,кп1п2 к=1

Когда функция / непрерывна в области П, то можно полагать

А,иу = /(хг, у у), ' =1,..., т, у = 1,..., п.

Решение задачи (3) существует и единственно, потому что А > 0. Строим продолжение задачи (3):

и е М^: Би = /, / е М2N, /2 = 0, (4)

когда векторы

V е М2N: V = V, )', здесь блочная, верхняя треугольная матрица Б размерности 2NX 2N такова, что:

Бц =Л, Би = 0, Б21 =вА, Б22 = Мв,

(ul, = Е ад^А Vul, V1е мN.

Б =

где матрицы

Л 0

вА мв

и

вА = лв+вл, Mq= л2 - q2, q=v / vy - v /v x, a=v/v x + v /v y,

матрицы vx, v v с размерностями N x N определяются так:

(v xul, v0=ЕЁ (-(uu+i, j i=i j=i

"i,i, f )hiivi,i, j )hih2, Ui

,i, j^' V2>

i,m+i, j Vi,m+i, j

= 0, j = i,..., и,

III п

(Vй= ЕЕ(-(ы1,»,у+1 -"1,1,у )h2\1,J ии,«+1 = ги,«+1 = 0, 1 = 1,...,т.

1=1 у =1

Вводятся подпространства из векторов в пространстве М2N :

Х1 = {у = (у1',у2')': + Мвг2 = 0}, ¿2 = {у = (У1',у2')': у = о}.

Лемма 1. Решение задачи (4) и е 21 существует, единственно.

Метод итерационных факторизаций при фиктивном продолжении

Дополнительно определяется блочная матрица С размерности 2N X 2N, когда

С11 = С22 = Мв , С12 = вА , С21 = вА,

тогда

C =

мв

-вл

вА Mq

(5)

Чтобы решать задачу (4), предлагается итерационный процесс:

ик е М2N : С (и к - ик-1) = -тк (Бик-1 - /), к е N, ^ > 0,"й° е ^

В итерационном процессе (5) на каждом шаге получается задача с факторизованным оператором такого вида:

йе CN: (ЬЬ*)2й = Ё, Ёе CN, которая расщепляется на простые задачи

а) Т е ^, ЬТ = Ё, Ё е CN,

б) Р е ^, Ь*Р = Т , Т е CN,

в) б е CN , LQ = Р , Р е СN ,

г) U е CN, L*U = Q, б е CN,

здесь матрицы:

L = vX -iv;, LL = L = vx + ivy, LLL = (v;-ivy')(vx + ivv) = A + ie,

поэтому

или

(LL* )2 = (A + ie)2 = Mq+ ieA ,

(Mq + i в л )(«i + U) = fi + i/2

I Мви1 - вАи2 = /, и1 + ш2 = й,

|0А«1 + М в" 2 = /2, ./1 + /2 = Ё, тогда на каждом шаге процесса (5) возникает следующая задача:

Си = /, и = (и;,"2) ', / = (/',/2')'.

Вводится норма:

Теорема 1. Для итерационного процесса (5) в оценке относительной ошибки:

Т-k -

<Ei

— 0 —

если тк, к е N выбирать на основе метода минимальных поправок, то получается, что:

Л

Л

е1 <cpk, cе (0; +~), ре (0;1).

Вывод. По виду матриц L, замечается, что задача (3) с N неизвестными по теореме 1 с помощью итерационного процесса (5) решается при относительной погрешности е1 за O(N 1пе1-1) арифметических действий. А при довольно малых значениях параметров дискретизации решение требует O(Щ операций, а итерационный процесс будет сходиться в несколько итераций и превратится в почти прямой метод.

Алгоритм численного решения модельной задачи для бигармонического уравнения

Решается задача:

Ли1 = Л,

итерационным процессом (5). При выборе итерационных параметров применяется метод минимальных поправок. Тогда для решаемой задачи в итерационном процессе применяется следующий алгоритм вычислений, где хук = ик - и - ошибка:

1) выбирается начальное приближение:

ми = 0е 2Х\

2) вычисляется невязка:

3) ищется поправка:

-к-1: -к-1

= Ли1к-1 - /1, г2к~1 = 0, к е N;

-к -1

Wк-1 е м2N : к-1 = г к-1, к е N;

4) определяется квадрат нормы ошибки:

2

Ек -1 =

У

к-1

= Г

,к-1 ^к-1

= Г

—к-1 —к-1

^WЛ

к е N;

5) ставится условие остановки итераций:

ЕкЕ1 < Е2, Ее (0; 1), ке N

6) вычисляется дополнительно вектор:

к к-1: —-1 =ЛW1k-1, к2к-1 = 0, к е N;

7) находится дополнительно вектор:

Хк-1 е М2N: сх к-1 =к к-1, к е N

8) определяется итерационный параметр:

тк =

DW к-1, W к -1

С ~lDW к -1, DW к-1

г к-1, W к-1

X-1 ^ к-1

кк -1 wk-1 ¿1 , ^

Хк-1 кк-1 Ь1 , ¿1

к е N

9) вычисляется новое приближение:

гк г к-1

и = и

Ее (0; 1) - задаваемая относительная погрешность.

т^к-1, к е N

Программируемый алгоритм численного решения модельной задачи для бигармонического уравнения

Решается задача:

Ли1 = Л,

которая может записываться в виде:

или

если

И4Ли1 = И4 /

а) И2Аи1 = ^ , б) И2А^1 = И4/1,

И = Ъ1/(т + 0,5), И2 = Ъ2/(п + 0,5), 7 = 1,...,т, у = 1,...,п, т = п, Ъ1 = Ъ2, И = И1 = И2,

Б'С 1Б

итерационным процессом (5). При выборе итерационных параметров применяется метод минимальных поправок. Тогда для решаемой задачи в итерационном процессе программируется приведенный ранее алгоритм вычислений. Решаемая задач получается из разностной схемы, которая записывается в следующем виде:

4и1,7, у - и1,7-1, у - и1,7+1, у - и1,7, у-1 - и1,7, у+1 = 4-,7, у , 7 = 1,..., m, ] = n,

где

«1,т+1,у = 0 «1,0,у = «1,1, у, J = П

«1,7,Я+1 = 0, «1,7,0 = «1,7,1 , 7 = т

44U, У - 4U-1, j - 4U+1, j - 4U, у-1 - 4U, у+1 + й4°«1,7, j = ^ A,,, j , 7 = и. ^ j = l,•••, П

где

4-,т+1,у = 0 4-,0,у = 4-,1,у , у' = 1,...,^

Ч1,7,п+1 = 0, 41,7,0 = 4ш, 7 = 1,...,т.

1. Выбирается начальное приближение:

и0 = Ое 2у.

При программировании элементы массива и1°гу7, у = 0,..., п +1. Полагается, что

^7, у = 0 ] = п .

2. Вычисляется невязка:

Ок-1: Гук-1 = И4Лик-1 - И4Л, Г2к-1 = 0, к е N ,

подробнее так:

а) Чк-1 = И2 Аик-1, б) Ц_к-1 = И2 Ачк-1 - И4/1. При программировании полагается, что

к-1 п к -1 к-1 ■ 1 к-1 п к-1 к-1 ■ л

и = 0, и = и , 7 = 1,...,п, и = 0, и = и , ] = 1,...,п.

-,!>+- ' 1,1,0 1,1,1 ''' ' 1,п+1,у ' 1,0,у 1,1,у ' ' '

При программировании ч'-- 7,у = 0,...,п +1 - элементы массива. Элементы этого массива

вычисляются по формуле:

ч к-1 л к-1 к-1 к-1 к-1 к-1 - л

а) Чу,у = 4«1,7,У -«1,7-1,J -«1,7+1,j -«1,7,У-1 -«1,7,У+1, h j = 1,•••,П ,

а затем полагается, что

4к-1 = 0, 4к-1 = 4к-1,7 = 1,...,п, 4к-1 = 0, 4к-1 = 4к-1, у = 1,...,п.

1,7,п+1 ' -'1,1,0 -"Щ 1,п+1,у ' -'1,0,у -'1,1,у'-' ' '

При программировании гк- 7,у = 0,...,п +1 - элементы массива. Элементы этого массива вычисляются по формуле:

к-1 к-1 к-1 к-1 к-1 к-1 , 4 Л-1 ■ ■ ,

б) Г1,,,у = 44-,7,у - 41,7—1,у - Чин-',у - 41,7,у-1 - 41,7,у +1 - И Л,,у, 7,у = и.n,

3. Ищется поправка:

wWк-1 е М2N : И4CWк-1 = Ок-1, к е N, что равносильно последовательному решению четырех систем линейных алгебраических уравнений с треугольными матрицами:

а) Тк-1 е СN: ИЬТк-1 = Ёк-1, б) Рк-1 е СN: ИЬРк-1 = Тк-1. в) ёк-1 е СN : ИЩк-1 = Рк-1, г) Жк-1 е СN : ИЬ*Жк-1 = ёк-1. При программировании решения системы из а) будут в массиве с элементами Тк~1 7, у = 0,..., п + -. При вычислении элементов этого массива сначала полагается

Тк0-1 = Тк-1 = 0,7,у = -,...,п , а затем вычисляются элементы массива по формуле:

T к-1 _

Т У _

(1 + i)r£") + (1 + 1)77---Л- + (1 - 1)7У

и

2

когда индексы принимают последовательно значения 1, у = 1,..., п.

При программировании решения системы из б) будут в массиве с элементами

Р ,■1 1,у = 0,...,п +1. При вычислении элементов этого массива сначала полагается

Зк-1 ' и у рк-1 = р 1,п+1

к-1 п+1, у

= 0,1,у = 1,...,п, а затем вычисляются элементы массива по формуле:

рк -1 =

', у

(1 - хТ-1+(1 -1) Р+-!,+(1+1) рк-1

'2,

1,у • ^ 1+1,у ■ V- ■ /,у +1 _

когда индексы принимают последовательно значения 1, у = п,...,1.

При программировании решения системы из в) будут в массиве с элементами

Оку1 1, у = 0,..., п +1. При вычислении элементов этого массива сначала полагается Оо"1 = <2о-1 = 0,1,у = 1,...,п , а затем вычисляются элементы массива по формуле:

ак-1=[(1+х) ру+(1+1)6^1,+(1 - 1)ак"-1

/2.

когда индексы принимают последовательно значения 1, у =1,...,п.

При программировании решения системы из г) будут в массиве с элементами

0к,11, у = 0,..., п +1. При вычислении элементов этого массива сначала полагается

1, у

гк-1

0^+1 = у = 0,1, у = 1,..., п, а затем вычисляются элементы массива по формуле:

=

1, у

(1 - 1)ак;1+(1 -1)0+-1+(1+1)0у

2,

когда индексы принимают последовательно значения 1, у = п,...,1. 4. Определяется квадрат нормы ошибки:

-к-1

Ек-1 =

У

= < гк-\ммк-1h-2 = Гк-1,Мк-1h-2, ке N.

При программировании у = 0,..., п +1 - элементы массива. При вычислениях

элементов этого массива полагается:

= Яе^А 1, у = 1,..., п + 1.

И

Ек-1 = (Е Егк>и1)h-2, ке N.

1=1 у=1

5. Ставится условие остановки итераций:

Ек_1/Е0 < Е2, Ее (0; 1), ке N,

где Е - задаваемая величина.

6. Дополнительно вычисляется вектор:

Л

к -1 . =-к-1

: Л1к-1 = h4ЛW1k-1, Л2к-1 = 0, к е N ,

подробнее так:

а) Цк-1 = Ъ2Амк-1, б) Л1к-1 = h2 АЦ^-1

При программировании полагается, что:

мк-1 = мк-1,1 = 1,...,п, мк-1 = мк-1, у = 1,...,п.

1,1,0 1,1,1 ' ' ' ' 1,0,у 1,1,у' -1 ' '

к-1

При программировании Цц ] 1,у = 0,...,п +1 - элементы массива. Элементы этого массива вычисляются по формуле:

а) = 4мк-1 - м-, у - wk-!l,} - - мк-1+1,1, у = 1,..., п ,

а затем полагается, что:

чк-1

1 1,1,п+1

= 0, цк-1 = цк-1,1 = 1,...,п, цк-1 = 0, цк-1 = цк-1, у = 1,...,п.

' -4,1,0 ^ 1,1,1 1,п+1,у ' -'1,0,у -'1,1,у'-' ' '

При программировании Лиу!,] = 0,...,п +1 - элементы массива. Элементы этого массива вычисляются по формуле:

лч «к-1 Л к-1 к-1 к-1 к-1 к-1 - л

б) ки,у = 44у,7,у -4и-1,у -417+1,у -4и,у-1 -4l,l,J+l, 7,у = и.п 7. Дополнительно находится вектор:

хОк-1 е М2N: И4СХк-1 =кк-1, ке N, что равносильно последовательному решению четырех систем линейных алгебраических уравнений с треугольными матрицами:

а) Тк-1 е СN: ИЬТк-1 = Нк-1, б) Рк-У е СN: ИЬ*Рк-1 = Тк-1.

в) ёк-1 е СN : ИЩк-1 = Рк-У, г) Нк-У е ^ : ИЬ*Нк-У = О-1.

При программировании решения системы из а) будут в массиве с элементами Тк~1 7, у = 0,..., п + -. При вычислении элементов этого массива сначала полагается

Тк0-1 = Т0к-1 = 0,7,у = -,...,п , а затем вычисляются элементы массива по формуле:

T к-1 _

7 У =

(1+тк-у + (1+1)77-у;1, + (1 -1)77 у-1

когда индексы принимают последовательно значения 7, у = 1,... , п.

При программировании решения системы из б) будут в массиве с элементами

-1

-1

'2,

Рк ^у- 7,у = 0,...,п + -. При вычислении элементов этого массива сначала полагается

рк-1 _ рк-1 7, n+1 ТП+1, j

= 0,7, j = !,•••,n, а затем вычисляются элементы массива по формуле:

р -1 _

7, J

(1 -1)777 У + (1 -1) р+-у. + (1 +1) рк-1

/2.

7,у ■ V- 7+1,у ■ V- ■ 7,у+1 _

когда индексы принимают последовательно значения 7, у = п,...,1.

При программировании решения системы из в) будут в массиве с элементами

О- 1 7, у = 0,..., п + -. При вычислении элементов этого массива сначала полагается О-1 = = 0,7,у = -,...,п , а затем вычисляются элементы массива по формуле:

Qk^ -1 = Г(1 + i) Pj + (1 + 1)йк-1 у + (1 - 1)йк-1

/2.

>, у /*;, у ■ V- ■ ^7-1, у^ ^7, у-1_

когда индексы принимают последовательно значения 7, у = -,..., п.

При программировании решения системы из г) будут в массиве с элементами

-1

Ну 7,у = 0,...,п + -. При вычислении элементов этого массива сначала полагается Нкп+у = Н^-1 у = 0,7,у = -,...,п, а затем вычисляются элементы массива по формуле:

Хк-1

X7, У

о -1)0-1 + (1 - ОН+Уу + (1+1)нк-+-

когда индексы принимают последовательно значения 7, у = п,...,У. 8. Определяется итерационный параметр:

t =

Шк-1, W*^ -1

C -lDw'i -1, Dw к-1

h-1, W*^-1

X-1 'к-1

-к-1 W-1 h1 , w1

Хк-1 h-1 Ь1 , '/1

, ке N •

При программировании Х- 7-- 7, у = 0,..., п + - - элементы массива. При вычислениях элементов этого массива полагается:

Хк-У = ЯеНк-1,7,у = У,...,п.

При программировании:

t= Ш-1, ке N •

7=1 У=1 / 7=1 у=1

Вычисляется новое приближение:

мк = «7к-1 -t -1, ке N •

Это приближение и определяется по соответствующей формуле:

2

йк = йк-1 -ткЖк -1. При программировании ы^^ у .,у = 0,...,п +1 - элементы массива. Следующие элементы этого массива вычисляются по формуле:

К.,у = «Д -*кЧ5 /,У = 1,...,П .

Они принимаются за приближенное решение, если до этого выполнялось условие остановки итераций.

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

В рассматриваемой задаче:

.2 - Г ГЛ А -I -I п ЭАй. I Эй, I А «1 = в О, А«1 |Г1 = «1 |Г1 = 0, -п|г2 = "^1Г2 = 0,

где

О = (0;61)х(0;Ъ2), Г1 ={Й1]х(0;Ъ2)и(0;х^} Г2 ={0}х(0;Ъ2)и(0;^)х{0}, полагается, что Ъ1 = Ъ2 = 2,5,

/1(х;у) = 3((25 - 4х2)(125 - 4х2) + 12(25 - 4х2)(25 - 4у2) + (25 - 4у2)(125 - 4у2))/8 000,

т.е.

«1(х;у) = (25 - 4х2)(125 - 4х2)(25 - 4у2)(125 - 4у2)/1024 000. При аппроксимации этой задачи выбираются значения т = п = 70, 72, 74,...,90. Для численного решения получаемых систем линейных алгебраических уравнений Лй1 = / применим приведенный ранее алгоритм, полагая «1° = 0, Е = 0,001. Для функции к = к(п) числа итераций в зависимости от п приводится график. При условии остановки в итерационном процессе с выбранной относительной погрешностью Ее (0; 1), можно заметить, что вычисления на ЭВМ экспериментально подтверждают асимптотическое поведение, независимость к - количества итераций от т (п) - количества узлов сетки по координатным направлениям осей 0х (0 у) для получения Е - задаваемой относительной погрешности решения линейных систем алгебраических уравнений (см. рис.).

График зависимости к - количества итераций от т(п) - количества узлов сетки по направлениям координат 0х(0у) при заданной относительной погрешности счета Е = 0,001

Литература

1. Самарский, А. А. Методы решения сеточных уравнений / А. А. Самарский. - М.: Наука, 1978. - 591 с.

2. Оганесян, Л.А. Вариационно-разностные методы решения эллиптических уравнений / Л.А. Оганесян, Л.А. Руховец. - Ереван: Изд-во АН АрмССР, 1979. - 235 с.

3. Обэн, Ж.П. Приближённое решение эллиптических краевых задач / Ж.П. Обен. - М.: Мир, -977. - 383 с.

Поступила в редакцию 27 ноября 2018 г.

Bulletin of the South Ural State University Series "Mathematics. Mechanics. Physics" _2019, vol. 11, no. 1, pp. 34-42

DOI: 10.14529/mmph190105

FAST SOLUTION OF A MODEL PROBLEM FOR THE BIHARMONIC EQUATION

A.L. Ushakov

South Ural State University, Chelyabinsk, Russian Federation E-mail: ushakoval@susu.ru

The biharmonic equation in a domain of rectangular shape when boundary conditions are mixed is being considered. Numerical solution of this boundary value problem uses iterative factorization on fictitious continuation after finite-difference approximation of the problem to be solved. Eventually, everything is reduced to solving the linear systems of algebraic equations, the matrices of which are triangular with three or less nonzero elements in lines. If approximation error of the initial problem is sufficiently small, the demanded relative error of the used iterative process gets obtained in several iterations. In this case, the developed iterative method turns out to be the method that has optimal asymptotics by the number of actions in arithmetic operations. The proposed iterative method essentially uses specificities of the obtained model problem. Such a problem can arise in methods of the type of fictious components, regions and spaces, when boundary value problems with elliptic equations in the regions of sufficiently arbitrary shape are being solved. The algorithm at implementation of the iterative process, when the choice of iterative parameters is made automatically using the method of minimal corrections, is given. The criterion for process termination after achieving the preliminarily determined ratio error is specified. Graphic result of a computational experiment that proves the asymptotic optimality of the iterative method in computational outlay is given. Complex analysis gets essentially used when developing the method.

Keywords: fictitious continuation; iterative factorizations.

References

1. Samarskiy A.A. Metody resheniya setochnykh uravneniy (Methods for solving finite-difference equations). Moscow, Nauka Publ., 1978, 591 p. (in Russ.).

2. Oganesyan L.A., Rukhovets L.A. Variatsionno-raznostnye metody resheniya ellipticheskikh uravneniy (Variation-difference methods for solving elliptic equations). Erevan, Izd-vo AN ArmSSR Publ., 1979, 235 p. (in Russ.).

3. Oben Zh.P. Priblizhyennoe reshenie ellipticheskikh kraevykh zadach (Approximate solution of elliptic boundary value problems). Moscow, Mir Publ., 1977, 383 p. (in Russ.). [Aubin J.-P. Approximation of elliptic boundary-value problems. New York, Wiley-Interscience, 1972, 360 p.]

Received November 27, 2018