УДК 517-956.6
ЗАДАЧА ГЕЛЛЕРСТЕДТА ДЛЯ НЕЛИНЕЙНОГО ФУНКЦИОНАЛЬНО -ДИФФЕРЕНЦИАЛЬНОГО УРАВНЕНИЯ СМЕШАННОГО ТИПА
GELLERSTEDT'S TASK FOR THE NONLINEAR FUNCTIONAL AND DIFFERENTIAL EQUATION OF THE MIXED TYPE
Е.В. Чаплыгина, А.Н. Зарубин E.V. Chaplygina, A.N. Zarubin
ФГБОУВО «Орловский государственный университет имени И.С. Тургенева», Россия, 302026, г. Орел, ул. Комсомольская, 95
FSBEI "Orel state University named after I. S. Turgenev", Russia, 302026, Orel, Komsomolskaya str., 95 E-mail: [email protected]; [email protected]
Аннотация. Исследуется краевая задача для нелинейного уравнения смешанного типа с оператором Лаврентьева-Бицадзе и функциональным запаздыванием, опережением. Построено общее решение уравнения. Задача однозначно разрешима. Найдены в явном виде интегральные представления решений.
Resume. The regional task for the nonlinear equation of the mixed type with Lavrentyev-Bitsadze's operator and functional delay, an advancing is investigated. The common decision of the equation is constructed. The task is unambiguously solvable. Integrated submissions of decisions are found in an explicit form.
Ключевые слова: уравнение смешанного типа, разностные уравнения, функция Римана, задача Ко-ши, задача Дирихле..
Key words: equation of mixed type, difference equations, function of Riman, Cauchy problem, Dirichlet
problem.
Постановка задачи
Уравнение
Lu(x, y) = Ux (x, y) + sgn(y)Uyy (x, y) = u(ln x, y)u(ex, y), (1)
рассмотрим в области D = D+ ^D ^ J, где D+ = {(x,y) : x0 < x < x2,0 < y < h} = D^^ Д+ ^ J (0 < h = const),
D = D0 ^ D- - эллиптическая и гиперболическая части области D, причем D; = {(x,y):xk < x < xk+v 0 < y < h} (k = -1,0,1,2), Dk = {(x,y): -y <a\(x) < y + xv -x/2 < y < 0} (k = -1,0,1,2), I = {(x,y): x0 < x < x2,y = 0}, J = {(x,y): x = xi,0 < y < h},
a 1 (x) = a (x) = ex, a0 (x) = x, a^ (x) = a (x) = ln(x), af (x) = a (a (x)) = ln(ln(x)) и x_! = -», x0 = 0, x = 1, x2 = e, x3 = ee.
Пусть Dk = D+k ^ D- ^ , где = {(x, y): xk < x < xk+y = 0} (k = -1,0,1,2) . Тип функционального запаздывания и опережения следует из представлений
u(ln x, y) = u(x - (x - ln x), y) = u( x - tx (x), y),
u(ex, y) = u(x + (ex - x), y) = u(x + t2 (x), y),
Tx (x) = x - ln x > 0, t2 (x) = ex - x > 0.
Задача G. Найти в области D функцию u(x, y) е C(D) n Cl(D \ J) n C2(D \(I ^ J)), удовлетворяющую уравнению (1), краевым условиям
u(x,h) = p(x),x0 < x < x2, (2)
u(x,—x) = (x), x0 < x < x /2, (3)
u(x,lnx — x) = ^1(x),ex'/2 < x < x2, (4)
u(x0, y) = u(x2, y) = 0,0 < y < h, (5)
u(x, y) = cos(ex + Уд/ — sgn(y) ) cos(ex — уд/ — sgn(y)), (x, y) е D—i, (6)
u( x, y) = cos(ln(ln x) + Уд/ — sgn( y) ) cos(ln(ln x) — Уд/ — sgn( y) ), (x, y) е D 2, (7) условиям сопряжения
u(x,0—) = u(x,0+) = co(x), x0 < x < x2, (8)
u^ (x,0—) = u (x,0+) = v(x), x0 < x < x2, x ^ , (9)
условиям согласования
P(x0 ) = (P(x2 ) =¥o(xo) = (x2 ) = 0, (10)
где p( x), щк (x) (к = 0,1) - заданные непрерывные достаточно гладкие функции.
Общее решение уравнения (1)
Уравнение (1) в терминах функций
u± (x, y) = u(x, y), (x, y) е Dk (к = 0,1), (11)
с учетом (6) и (7), можно записать в форме системы
Lu ± (x, y) = cos(x + y^ — sgn(y) ) cos(x — Уд/ — sgn( y) )Au ± (x, y), (x, y) е D0±, (12)
где
u ± (x, y) = (u± (x, y),u± (^x, y))T , A = 0 n , (13)
V1 0 J
которая в характеристических переменных
í = x + Уд/ — sgn(y), r = x — Уд/ — sgn(y) (14)
будет иметь вид матричного уравнения
4^ (í,r) = cosícosrAu± (í,r). (15)
Используя известную [1] для уравнения (15) функцию Римана
R(t, q; í, r) = J 0 (i^Aa(t, q;í,r) ),
где
í r
cr(t, q;í,r) = J cos rdrJ cos sds = (siní — sin t)(sinr — sin q), (16)
t q
согласно [2, с.43] можно записать общее решение уравнения (15) в форме
u ± (í, r) = J 0 (^Aa(0,0;í,r))ü + (0,0) +
í __r _
+ J J0 (ijAa(t ,0;í,r))^1± (t )dt + J J0 (iy¡Aa(0, q;í,r) )&± (q)dq, (17)
0 0
где фф+ (?), ф2+ (?) — произвольные дважды непрерывно дифференцируемые вектор-функции; I = л/—1 , ^ (2) — функция Бесселя [3, с.727] первого рода нулевого порядка.
Поскольку матрица А из (13) имеет различные собственные значения \ = 1, ^ = —1, то она приводима к диагональному виду, т.е. существует матрица Т ^ТА | ^ о) такая, что
T— AT a = Л a =
Л 0 ^ (1 — Л 1 ( 1 Л
V 0 Л j
, причем T^ =
V1 1 J
а T— = -A 2
V—1 1j
Значит,
J о (i4A¿) = J 0 (Ít¡TaA AT-/a) = TJ o )T—
0V\ A A A / A 0*
Г 7 /„• ЛГ^Гч А Л
= T
^(iVV) 0 fYi(ff) Y 2(^)1
T A
0 2g ) J A 2
Y 2(tf) Yi(a)
(18)
где
уи(а) = Jо¡(i^)-(-1)" Уо^^^а) (" = 1,2). (19)
Поэтому из равенства (17), в силу (18), (19), (11) и возвращения к старым переменным по
формулам (14), найдем общее решение уравнения в форме
+
2о _
и+ (а 2(х) у) = |ф+ (0у о(1л1 а0А 2+, 2о+ +
о
—+
2+ I-
+ \ф2+ (? у о (Ца(0,X; 2+, 2++ ) )йг, (х, у) е Д+ (к = 0,1) (20)
или
zk _
щ(x y) = \Ф!(t) J 0(i v ^(tA z°, zk) )dt+
^ I-
+ |ф2±(t)J0(4v(0,t;z±±,^))dt, (x,y) e D¿ (k = 0,1), (20')
0
где a° (x) = x, a\ (x) = a2 (x) = ex; ф± (t), фф (t) — произвольные дважды непрерывно дифференцируемые функции; z I = a k (x) + iy ,z+k=a\ (x) — iy, z — =a k (x) + y, z—=af (x) — y (k = 0,1), причем а0 (x) = x, aj (x) = a (x) = ln x и
+ + 3 3
u± (xj — 0, y) = u± (xj + 0, y) = 0, 0 < y < h, — (u± (— 0, y)) ф — (u± (+ 0, y)), 0 < y < h. (21)
3x 3x
В равенствах (20), (20') учтено, что U ± (0,0) из (17), т.е. в старых переменных U ± (0,0) = 0 в силу условий (5) и (21).
На основании (16) равенство (20) можно записать в форме
±
z0 _
u± (a k (x), y) = |ф± (t)J0 (i^J sin z± (sin z± — sin t))dt +
0
—±
z0 _
+ |ф± (t)J0 (i^sin z± (sinz± — sint))dt, (x,y) e D± (k = 0Д)
0
или, после интегрирования по частям и соответствующих замен, в виде
Uk (a2 (x), y) = M0 (sin z0±, sin z0± ) —
ÍM0(s • sinz±,s • sinz± )—J0(iJ(] — s)sinz± sinz± )ds, (x,y) e D± (k = 0Д). (22)
J Яс
0 3s Для уравнения (22) имеет место формула обращения
M0 (sinz0 ,sinz± ) = u± (a2 (x), y)\_muu± (sinz0 ,sinz± )] + (23)
V* ± • —± j ± / • ± • —± \ + sin ^sin z0± J uk(s •sin z0±,s •sin z0±) 7 (^s(j — s)sin z± sin z0± )ds, (x, y) e D0±
2 J Vs(j — s)
z 0 ) alcs in (s *s in z 0 ,
M0 (s • sin z±, s • sin z± ) = Jф± (r )dr + Jф± (r )dr, (24)
0
0
0
и I (5 • э1п г о, ^ • э1п г * ) = а 2 ^ 1 [агсв1п( 5 • э1п г к ) + агсв1п( 5 • э1п г к )]|,
[агсв1п( ^ • э1п г к ) — агсБ1п( 5 • э1п г* )] | (к = 0,1), (25)
+г(1±1)/2
2
причем при ^ = 1 из (24), (25) получим
го го
М0 (э1п г0, э1п ) = |фг (г )—г + |ф° (г )—г,
о о
ик От го, вт го ) = ик (а2 (х) У);
(5) = (я); /0 (5), (5) — модифицированные функции Бесселя [3, с.730] первого рода нуле-—5
вого и первого порядка.
При у = 0 выражения (22), (23) после преобразований представимы равенствами
иО (а 2( х),0) = Мо0 (х,0) —
х 7 , . , . . _
^^ , ,------/------ _0 < х < х
|Мк (?,0)—(г^э1п х(э1п х — э1п I) )—?, 0 < х < хг, (26)
М о0 (х,0) = иО (а 2( х),0) +
Г к / к / ч ^N »т х • соэ Г б т ,. Г-.-—-:-г. , Л „ ^
I ик (ак (?),0)--/0 (г^эт ?(б1п х — э1П ?))—?, 0 < х < хг, (27)
э1п ? • соэ х бх
где
х
Мо (х,0) = М0 (э1п х, э1п х) = I[фк (г) + фк (г)]—г. (28)
Кроме того, из (22), (23) можно получить следующие формулы взаимного обращения
ик (а2(х),0) = М0у(х,0) —
гку У^ —±у± о у
к , ^ соэ х • э1п? б т ,. г--Т^-:-Тч , ^ . .
— IМ0 (?,0)--(г^тх(э1пх — Э1П ?))—?, 0 < х < , (29)
э1п х • соэ ? б?
М 0у (х,0) = ику (а 2( х),0) +
х О
+ Iик (ак (?),0)—/0 (г^/э1п ?(э1п х — э1п ?))—?, 0 < х < х, (30)
о бх
где
М 0у (х,0) = Моу (81П гк ,81П гок )| =о = г (1±1)/2[фк (х) — фк (х)]. (31)
Однозначная разрешимость задачи О Теорема 1. Если ф(х) е С[х0, х2 ] С2(х0, х2), (х) е С[х0, хг /2] С 2(х0, х /2), ^(х) е С[ех^2,х2]^С2(ех^2,х2), абсолютно интегрируемы на своих промежутках, ф(х0) = ф(х2) = ^(х0) = ^(х2) = 0 и (х) при х ^ х0, щ[(хх) при х ^ х2, допускают интегрируемую особенность, то существует единственное решение и( х, у) задачи О. Единственность решения задачи О следует из утверждений.
о
о
о
Лемма 1. Если и(х, у) - решение уравнения (1) в области О— = О0 ^ О— из класса С(О " ) о С2 (О" ), обращающееся в нуль на у = —х, х0 < х < х /2, у = 1п х — х, ех1 2 < х < х2 и в О—\, О 2—, то
х2
Р = |о(x)v(> 0.
Доказательство леммы аналогично [5, с. 128-130].
Лемма 2. Если и(х, у) - решение уравнения (1) в области О+ из класса С(О + ) о С2 (О+ \ 3), обращающееся в нуль при х = хк (0 < у < И) (к = 0,2), у = к (х0 < х < х2 ) и в областях О ^, О2+, то ( < 0 и
Р + Я ^ (х, y) + u2y (х, y)]d*dy = 0.
Доказательство леммы аналогично [6].
Вопрос существования решения задачи О в области О = О0 ^ О ^ 3 связан с построением в О = Ок ^ Ок ^ /д. (к = 0,1) на основании общих решений (20) функций и° (х, у), (х, у) е Ок (к = 0,1), удовлетворяющих условиям (2)-(10), (11), (21) в которых ф(х),щк (х) (к = 0,1) заданы, а б)(х),у(х) подлежат определению. Поскольку условие (21) на х = х (0 < у < И) известно, то достаточно решить задачу О для уравнения (1) в областях О0 и О1 , то есть найти функции щ (х, у), (х, у) е О° и щ (х, у), (х, у) е О° (или
х yl(х y) e D0
u1(a2(x),y),(хy) e Do±)•
Проведем построение решения задачи G для уравнения (1) в области D0 = D0 ^ D0 ^ /0 , то есть найдем функции щ (х, y), (х, y) e Dq при условиях (2)-(io), (11):
u+ (х, h) = р(х), х± < х < х1, (32)
u+ (х±, y) = u± (х1, y) = 0,0 < y < h, (33)
u± (х,-х) = y/0 (х), х0 < х < х1 / 2, (34)
u0 (х,0-) = u+ (х,0+) = о(х), х0 < х < х1, (35)
U0y (х,0-) = u0y (х,0+) = v(х), х0 < х < х1, (36)
о(х0) = ®(х1) = У0(х0) = 0, Р(х0) = Р(х1) = 0. Задача Коши. Найти в области D- решение u- (х, y) уравнения (1) из класса С(D0 ) о С2 (D- ), удовлетворяющее условиям (35), (36), то есть
u0 (х,0-) = о(х), х < х < х, u0у (х,0-) = v(х), х0 < х < х, где о(х), v(х) - непрерывные достаточно гладкие функции, причем о(х0 ) = о(х ) = 0 •
Теорема 2. Если о(х) e С[х0,х]оC2(х0,х), v(x) e С1(х0,х), Ох) = о(х) = 0, то существует единственное решение задачи Коши u0 (х, y) e С(D0 ) о С 2(D0 ) вида
1 Z0 I-
u- (х, y) = — IJ0 (iyj(sin z- - sin t) sin z- )[p'(t) + r(t)]dt +
2 0
+1 0 J0 (i4sin Z0 (sin Z0- - sin t))[p'(t) - r(t)]dt, (х, y) e D0, (37)
2 *
2 0
х
0
+
D
, N / ч ^ / ч sin x • cos t d т ,. !—.——-:—- .
p(x) = g(x) + I G(t)--/0 (i^/sm t(sin x - sin t))dt, (38)
J sin t • cos x dx
r(x) = v(x) + Jv(t) —10 (iylsin t(sin x - sin t))dt. (39)
и дх
Доказательство следует из (20) (к = о). Действительно, для определения произвольных функций ф[ (?), ф2 (?) учтем в общем решении (20) (к = о) (или (22)) уравнения (1) в области О-условия (35), (36) задачи Коши. Тогда
G)(x) = M- (x,0) - JM- (t,0)—J (i^Jsin x(sin x - sin t))dt, 0 = x0 < x < xx,
0
dt
cos x • sin t д
/ N гт- / ЛЧ Ггг- s cosx • sin t д T ,. !—.-—-:-- . -
v(x) = M0 (x,0) - IM0 (t,0)--J (4sin x(sin x - sin t))dt, 0 = x0 < x < xx.
7 l sin x • cos t dt
На основании (27) и (30) имеем M0 (x,0) = p(x), MQy (x,0) = r(x), x0 < x < x, то есть, в силу (28), (31), дифференцируя первое равенство, приходим к системе
\ф- (x) + ф2 (x) = p'(x),
[fo (x) -ф2- (x) = r(x).
Из этой системы уравнений найдем функции ф[ (x), ф2 (x) , подставляя которые в (20) (или (20')) (к = 0), получаем решение u- (x, y) задачи Коши (37) в области D- .
Функциональное соотношение между g(x) и v(x) , принесенное из D0 на линию изменения типа уравнения (1) y = 0, x0 < x < xx, получим из (37), полагая y = -x и учитывая условие(34)задачи G:
2 2 x
(x) = - J[p (t) - r(t)]dt, 0 = x0 < x < xj /2,
2 0
то есть, после замены x на x /2 и дифференцирования,
p'(x) = r(x) + (x/ 2), 0 = x0 < x < xx. (40)
Выражение (40) является искомым функциональным соотношением.
Задача Дирихле. В области D^ найти решение u^ (x, y) е C(D0+ ) n C 2 (Di) уравнения (1), удовлетворяющее условиям (32), (33), (35), то есть
ui (x,0+) = g(x), ui (x, h) = ((x), x0 < x < xx, < (x0, y) = < (^, y) = 0, 0 < y < h, где ((x),g(x) - непрерывные достаточно гладкие функции, причем
g(x0) = g(x-) = ((x0) = ((x-) = 0.
Теорема 3. Если ((x), g(x) е C[x0, x¡ ] n C2 (x0, x¡) и
g(x0) = g(x1) = ((x0) = ((x ) = 0, то существует единственное решение ui (x, y) е
C(D0+ ) n C 2 (Di) задачи Дирихле вида
(x, y) = J J0 (i4(sin zi - sin t)sin Z0i )[p'(t) - ¿ R2lkh (p (t) - R'P (t))]dt i
0 к=0
Z+
z0 __
+ J J0 (4sin Z0+ (sin Z0+ - sin t))[£R2'hk(p'(t) - Rfp\t))]dt, (x, y) е Di, (41)
0 к=0
0
x
x
где R0 — оператор сдвига по х: R0 q(x) = q(x — 0),
P( x) = H x) +
W^^L j uo (s • sinzg,s • sinzo ) ^ (ф(1 — s). sinz J. sinz- + ^
2
Vs(l — s)
(42)
y=h
а p( х) определяется равенством (38).
Доказательство следует из (20) (k = 0). Действительно, для определения произвольных
функций фф (t), ф2 (t) учтем в общем решении (20) (k = 0) (или (22)) уравнения (1) в области D
условия (32), (35) задачи Дирихле.
Используя формулы обращения (23), (27), придем к системе
Гмо+ (х,0) = p(х), 0 = х0 < х < х,
|M0+ (sin Zg, sin Z0+ )| a = Дх), 0 = х0 < x < x,
которая, в силу (24), (28) и дифференцирования, примет вид
ф (х) + ф2+ (х) = p'( X),
ф (х) + Я2'кф+2 (х) = RXhp(х).
Подставляя из первого уравнения
фф (х) = p (х) - ф2 (х), 0 = х0 < х < х, во второе уравнение системы (43), приходим к разностному уравнению
где
решение [7] которого
а, в силу (44),
ф; (х) = R2x'hфJ (х) + r(х), o = х0 < х < Xj,
r(х) = Р( х) — R^X х),
Ф2 (х) = 2 RfV( х),
к=0
ф; (х) = p( х)—2 Rfkhr( х).
к=0
(43)
(44)
(45)
(46)
(47)
(48)
Учитывая выражения (47), (48) в (20) (или (20')) (к = 0), приходим к требуемому представлению (41) решения задачи Дирихле для уравнения (1) в области .
Найдем функциональное соотношение между 0)(х) и у(х), принесенное из О^ на линию изменения типа у = 0, 0 = х0 < х < х •
Подставляя в равенство (20) (или (22)) (к = 0) условие (36), согласно (29) получим урав-
нение
__х_ cos х • sin t D 1-
v(х) = Mjy (х,0) — jMjy (t,0) —-—J(iylsinх^тх — sint))dt, 0
— х0 < х < х^,
u sin х • cos t Dt
обращая которое относительно M0y (х,0) аналогично (30), в силу (31), придем к равенству фф (х) — ф2 (х) = —ir(х), т.е., на основании (47), (48), (46), к выражению
r( х) = ip\х) — 2i2 R2ikhp '(х) + 2i2 Rf(2k+1)A'( х),
представимому в виде
к=0
к=0
(1 — R2h )r(х) = —i(1 + R2h )p' (х) + 2iR[hp' (х), 0 = х0 < х < хг Выражение (49) является искомым функциональным соотношением.
0
Вопрос существования решения задачи О в области О0 = ^ О0 ^ /0 сводится к разрешимости системы функциональных соотношений (40), (49), то есть к разностному уравнению
(1 + 1Я2'к )г(х) = а(х) = -1(1 +1)(1 + я2'кУо (х/2) + (I + 1)^ХАА'(х), о = хо < х < хг. (50) Решение [7] разностного уравнения (50), аналогично (45), можно записать в виде
г (х) = 2 (-) "Я2тка(х). (51)
п=о
a(x) = (a(£),8(£ - x) - S(£ + x)) = ja(£)[S(£ - x) - S(£ + x)]d£,
Интегральное представление (51) можно [8] найти, поскольку [9, с.7] финитная на промежутке [о, х ] непрерывная функция
х1
+ х)) =
о
1 тл: где 5(2) =- V ехр(/Лт 2) - дельта-функция [4,с.711-714] Дирака, а А =-.
2 х ^^ х
Таким образом, учитывая (51) в (39), применяя формулы взаимного обращения (29), (30) получим
v(x) = r(x) - [r(t)--J0 (i^Jsin x(sin x - sin t) )dt, 0 = x0 < x < xx. (52)
J0 sin x • cos t dt
На основании свойств функций p(x),^0(x),входящих в (42), (50), из (52) следует, что v(x) G C:(x0,xj).
Очевидно, интегрируя (40), подставляя p(x), r(x) из (38), (51) и применяя формулы взаимного обращения (26), (27), найдем cc(x) g C[x0, xx ] о C2 (x0, xx).
Подстановка функций c(x) и v(x) в формулы (37), (41) приводит к окончательному виду решения задачи Коши и задачи Дирихле в областях D- и D^ , то есть в области D0 = D„ и D- и I0.
Список литературы
1. Чуриков Ф.С., Кокинасидии П.Д. 1976. О построении функции Римана для компактных уравнений методом промежуточного аргумента. Труды Кубинского университета, Т. 222: 5-12.
Churikov F. S., Kokinasidiya P. D. 1976. About creation of function of Riman for the compact equations by method of an intermediate argument. Works of the Cuban university, T. 222: 5-12.
2. Векуа И.Н. 1948. Новые методы решения эллиптических уравнений, ОГИЗ, М.-Л.: 296.
Vekua I.N. New methods of the solution of the elliptic equations. - OGIZ, M.-L., 1948. 296 pp.
3. Прудников А.Н., Брычков Ю.А., Маричев О.И. 1983. Интегралы и ряды. Специальные функции. Москва: Наука : 750.
Prudnikov A. N., Brychkov Yu. A., Marichev O. I. 1983. Integrals and ranks. Special functions. Moscow: Science: 750.
4. Тер-Крикоров А.М., Шабунин М.И. 1988. Курс математического анализа. Москва: Наука: 816.
Ter-Krikorov A. M., Shabunin M. I. 1988. Kurs of the mathematical analysis. Moscow: Science: 816.
5. Зарубин А.Н. 1999. Уравнения смешанного типа с запаздывающим аргументом. Орел: издательство ОГУ: 225.
Zarubin A.N. 1999. The equations of the mixed type with the late argument. Orel: 225.
6. Зарубин А.Н. 2014. Краевая задача для опережающе-запаздывающего уравнения смешанного типа с негладкой линией вырождения. Дифференциальные уравнения, Т. 50, №10: 1362-1372.
Zarubin A.N. 2014. A regional task for the operezhayushche-late equation of the mixed type with the rough line of degeneration. Differential equations, T. 50, No. 10: 1362-1372.
7. Зарубин А.Н. 2012. Краевая задача для уравнения смешанного типа с опережающе-запаздывающим аргументом. Дифференциальные уравнения, Т. 48, №10: 1401-1411.
Zarubin A.N. 2012. A regional task for the equation of the mixed type with the operezhayushche-late argument. Differential equations, T. 48, No. 10: 1401-1411.
8. Зарубин А.Н. 2015. Задача Трикоми для опережающе-запаздывающего уравнения смешанного типа с замкнутой линией вырождения. Дифференциальные уравнения, Т. 51, №10: 1315-1327.
Zarubin A. N. 2015. Zadacha Trikomi for the operezhayushche-late equation of the mixed type with the closed line of degeneration. Differential equations, T. 51, №10: 1315-1327.
9. Агранович М.С. 2008. Обобщенные функции, Москва: МЦНМО: 128.
Agranovich M. S. 2008. The generalized functions. Moscow: MCNMO: 128.