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

( Analysis of singularities of gradient of a solution in the boundary value problem for the Laplace equation outside cut in a plane with mixed boundary condition
Preprint, Inst. Appl. Math., the Russian Academy of Science)

Крутицкий П.А., Cгибнев А.И.
(P.A.Krutitskii, A.I.Sgibnev)

ИПМ им. М.В.Келдыша РАН

Москва, 2004
Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований (проект № 02-01-01067)

Аннотация

Рассматривается смешанная задача для уравнения Лапласа на плоскости вне разреза. В качестве граничных условий задается значение искомой функции на одной стороне разреза и значение ее косой производной на другой стороне. Проводится подробное исследование градиента решения на концах разреза, выводятся асимптотические формулы, обсуждается эффект исчезновения особенности.

Abstract

The mixed problem for the Laplace equation outside cut in a plane is considered. The Dirichlet condition is posed on one side of the cut and the skew derivative condition is posed on the other side. The gradient of the solution at the ends of the cut is studied in details. The asymptotic formulas are obtained. The effect of disappearing of singularity is discussed.


E-mail: conf@phys.msu.su



1  Постановка задачи



На плоскости x=(x1,x2)
R2 рассмотрим простую разомкнутую кривую G класса C2,l,  l (0,1]. Пусть кривая G параметризована и в качестве параметра выступает длина дуги s:  G = {xx=x(s)=(x1(s),x2(s)),  s [a,b]} . Отрезок оси Os, отвечающий контуру G, будем также обозначать G. Вектор касательной к G в точке x(s) обозначим tx=(cosa(s),sina(s)), а вектор нормали nx=(sina(s),-cosa(s)), т.е. cosa(s)=x1(s),  sina(s)=x2(s). Пусть плоскость разрезана вдоль контура G. Через G+ обозначим ту сторону контура G, которая остается слева при возрастании параметра s, а через G- - другую. Индексами + и - будем обозначать предельные значения функций на G+ и на G- соответствено.

Будем говорить, что функция u(x) принадлежит классу G, если:

1) u(x) C0([`(R2\G)])C2(R2\G) и u(x) непрерывна на концах разреза G,

2) u(x) C0([`(R2\G)]\X), где X=(x(a)x(b)) - множество концов разреза G,

3) при x x(d) X справедливо неравенство

u| const|x-x(d)|d,

(1)

где константа const > 0, число d > -1 и d=a либо d=b .

Под C0([`(R2\G)]) понимается класс функций, которые имеют предельные значения на разрезе G слева и справа, но эти значения могут быть различны во внутренних точках контура G, т.е. функции могут иметь скачок на контуре G.

Сформулируем смешанную задачу с косой производной для уравнения Лапласа вне разреза на плоскости.

Задача S . Найти функцию u(x) из класса G, удовлетворяющую в R2\G уравнению Лапласа

Du(x)=0,    x R2\G,

(2)

граничным условиям

u(x)|x(s) G+=f+(s),

(3)



ж
и
  ∂ u

nx
+b   ∂ u

tx
ц
ш
к
к


x(s) О G- 
=f-(s),
(4)

и условиям на бесконечности


u(x)=Aln|x|+O(1),      ∂ u

|x|
=  A

|x|
+O(|x|-2),     |x| →∞ .
(5)

Считаем, что f+(s) C1,l(G), f-(s) C0,l(G) - известные функции, A и b - заданные константы. В случае A=0 получим классическое условие ограниченности на бесконечности. В случае b = 0 получим смешанную задачу Дирихле - Неймана [2].

В более общей постановке задача S изучалась в [0], где рассматривался случай нескольких разрезов. Из результатов [0] вытекает, что решение задачи S существует и единственно. Более того, для случая нескольких разрезов в [0] было получено интегральное представление решения этой задачи в виде потенциалов, плотности в которых удовлетворяют однозначно разрешимой системе интегральных уравнений. Как показано ниже, в случае задачи S , т.е. в случае одного разреза, эту систему интегральных уравнений можно упростить, а значит, упростится и интегральное представление для решения. Эти результаты будут использованы в настоящей работе при выводе асимптотических формул, описывающих поведение градиента решения на концах разреза. Будем искать решение задачи S в виде суммы потенциала простого слоя V[m](x), углового потенциала T[n](x) (см. [G, K1] ) и константы c:

u[m,n](x)=V[m](x)+T[n](x)+c,

(6)



V[m](x)=-  1

2p

у
х
G 
m(s)ln|x-y(s)|ds,   T[n](x)=-  1

2p

у
х
G 
n(s)y(x,y(s))ds.

Неизвестные функции m(s), n(s), заданные на G, будем разыскивать в пространстве Cvq(G),  v (0,1], q [0,1).

Будем говорить, что функция F(s), определенная на G, принадлежит банахову пространству Cvq(G), v (0,1], q [0,1), если

F0(s)=F(s)|(s-a)(s-b)|q C0,v(G).

Норма в пространстве Cvq(G) задается соотношением

||F(s)||Cvq(G)=||F0(s)||C0,v(G).

Ядро углового потенциала определяется с точностью до 2pm   (m целое) формулами

cosy(x,y)=

 x1-y1


|x-y|

,   siny(x,y)=

 x2-y2


|x-y|

.

Очевидно, y(x,y) - угол между вектором [( ) || ( yx)] и направлением оси Ox1. С другой стороны, y(x,y) - многозначная гармоническая функция, сопряженная с ln|x-y| относительно соотношений Коши-Римана.

Пусть x - произвольная фиксированная точка плоскости, лежащая вне контура G, и y=y(s) G, тогда под y(x,y(s)) будем понимать произвольную фиксированную ветвь этой функции, непрерывно изменяющуюся по s вдоль G. При таком определении y(x,y(s)) угловой потенциал - многозначная функция. Для его однозначности необходимо потребовать выполнения следующего условия [G , K1] :



у
х
G 
n(sds=0.
(7)

При выполнении условия (7) гармонический угловой потенциал можно записать в виде гармонического потенциала двойного слоя


T[n](x)=  1

2p

у
х
G 
w(s)   ∂ ln|x-y(s)|

ny
 ds,

где


w(s)= s
у
х
a 
n(sds,    s О [a,b].

Из условия (7) вытекает, что w(a)=w(b)=0 и w(s) C0(G).

Если потребовать



у
х
G 
m(sds=-2pA,
(8)

то функция u[m, n](x) из (6) удовлетворит условиям на бесконечности (5).

Если плотности m(s),n(s) пpинадлежат Cvq(G),  v (0,1],  q [0,1), и удовлетвоpяют условиям (7), (8), то функция u[m, n](x) из (6) принадлежит классу G и удовлетворяет всем условиям задачи S , за исключением граничных условий [G, K1, K2]. В частности, неpавенство (1) будет выполнено с показателем d = -q, если q (0,1).

Заменим условие (3) эквивалентными граничными условиями на G:

 

 u


x

ú
ú



x(s) G+ 

=f ′ +(s),    f ′ +(s)

 d f+(s)


ds

C0,l(G),

(9)


u(x(a))=f+(a).

(10)

В условии (10) учтено, что u+(x(a))=u-(x(a))=u(x(a)) в силу непрерывности функции u(x) на концах контура G. Кроме того, при выбранной параметризации /t x=/s в любой точке x(s) G.

Плотности потенциалов будем разыскивать в виде


n(s)=
r1(s)+(
Ц
 

1+b2
 
+b)r2(s)

2
Ц

1+b2
,
m(s)=
r2(s)-(
Ц
 

1+b2
 
+b)r1(s)

2
Ц

1+b2
,
(11)

где r1(s),r2(s) Cvq(G),  v (0,1],  q [0,1). Подставляя функцию u[m,n](x) в граничные условия (4), (9) и учитывая условия (7), (8), (10), получим систему уравнений относительно функций r1(s),  r2(s) и константы c:


r1(s)-(
Ц
 

1+b2
 
+b)  1

p

у
х
G 
 r1(s)

s-s
 ds = F1(s),    s О G,
(12)

r2(s)+(
Ц
 

1+b2
 
-b)  1

p

у
х
G 
 r2(s)

s-s
 ds = F2(s),    s О G.
(13)


у
х
G 
(r1(s)+(
Ц
 

1+b2
 
+b)r2(s)) ds=0,
(14)


у
х
G 
(r2(s)-(
Ц
 

1+b2
 
+b)r1(s)) ds=-4pA
Ц
 

1+b2
 
,
(15)

 1

2
Ц

1+b2
ж
и
V[r2-(
Ц
 

1+b2
 
+b)r1](x(a))+

+T[r1+(
Ц
 

1+b2
 
+b)r2](x(a)) ц
ш
+c=f+(a),
(16)

где функции


F1(s)=2
Ц
 

1+b2
 
fў+(s)+2f-(s)+

+(
Ц
 

1+b2
 
+b) ж
и
-
у
х
G 
r2(s)Y1(s,sds+
у
х
G 
r1(s)Y2(s,sds ц
ш
,      s О G,      

F2(s)=2(
Ц
 

1+b2
 
-b)(
Ц
 

1+b2
 
fў+(s)-f-(s))-

-(
Ц
 

1+b2
 
-b) ж
и

у
х
G 
r1(s)Y1(s,sds+
у
х
G 
r2(s)Y2(s,sds ц
ш
,      s О G,   
(17)

принадлежат C0,l(G). Функции Y1(s,s), Y2(s,s) принадлежат C0,l(G×G) (см.[K1]) и имеют вид:


Y1(s,s)=  1

p
 cosf0(x(s),y(s))

|x(s)-y(s)|
,

Y2(s,s)=  1

p
й
л
 sinf0(x(s),y(s))

|x(s)-y(s)|
-  1

s-s
щ
ы
.

Через f0(x(s),y(s)) обозначен угол между вектором [( ) || ( xy)] и направлением нормали nx в точке x G. Этот угол считается положительным, если он отсчитывается от вектора nx против часовой стрелки, и отрицательным, если по часовой стрелке. Кроме того, при x y угол f0(x,y) полагается непрерывным по x,y на G. Углы f0(x(s),y(s)),  y(x(s),y(s)) с точностью до 2pm (m целое) связаны соотношением f0(x(s),y(s))=y(x(s),y(s))-a(s)-p/2.

Решая уравнение (12) относительно r1(s), а (13) - относительно r2(s) согласно лемме 2 из [0] и считая функции F1(s), F2(s) известными, найдем:


r1(s)=
(
Ц
 

1+b2
 
-b)

2
Ц

1+b2
F1(s)+

+  1

2p
Ц
 

1+b2
 
Q1(s)

у
х
G 
 F1(s)Q1(s)

s-s
ds-  B1cosph

Q1(s)
,     s О G,
(18)

r2(s)=
(
Ц
 

1+b2
 
+b)

2
Ц

1+b2
F2(s)-

-  1

2p
Ц
 

1+b2
 
Q2(s)

у
х
G 
 F2(s)Q2(s)

s-s
ds-  B2sinph

Q2(s)
,     s О G,
(19)

где B1, B2 - произвольные вещественные константы,

Q1(s)=(s-a)1-h(b-s)h,   Q2(s)=(s-a)1/2-h(b-s)1/2+h,   s G,

Q1(s),Q2(s) C0,h0(G),   h0=min{h,1/2-h}, константа h опpеделяется pавенством

h =

 1


2p

arcctg b ∈ (0,1/2)

и удовлетворяет соотношениям


cosph = ж
з
и

Ц
 

1+b2
 
+b

2
Ц

1+b2
ц
ч
ш
1/2


 
,   sinph = ж
з
и

Ц
 

1+b2
 
-b

2
Ц

1+b2
ц
ч
ш
1/2


 
,

tgph = (
Ц
 

1+b2
 
-b),     ctgph = (
Ц
 

1+b2
 
+b).

Подставляя функции r1(s), r2(s) из (18), (19) в уравнения (14), (15) и пользуясь формулами для интегралов из [0, п.6], получим для констант B1, B2 систему линейных алгебраических уравнений:


B1+(
Ц
 

1+b2
 
+b)B2=0,

(
Ц
 

1+b2
 
+b)B1-B2=-4A
Ц
 

1+b2
 
.
Отсюда
B1=-4A
Ц
 

1+b2
 
sinphcosph,   B2=4A
Ц
 

1+b2
 
sin2ph.
(20)

Раскрывая выражения F1(s),F2(s), получим для r1(s),r2(s) систему регуляризованных уравнений:


r1(s)+  1

Q1(s)

у
х
G 
r1(s)K11(s,sds+  1

Q1(s)

у
х
G 
r2(s)K12(s,sds =  F1(s)

Q1(s)
,
(21)

r2(s)+  1

Q2(s)

у
х
G 
r1(s)K21(s,sds+  1

Q2(s)

у
х
G 
r2(s)K22(s,sds =  F2(s)

Q2(s)
,
(22)

где s G;


K11(s,s)=-  Y2(s,s)Q1(s)

2
Ц

1+b2
-
(
Ц
 

1+b2
 
+b)

2p
Ц

1+b2

у
х
G 
 Y2(x,s)

x-s
Q1(xdx,

K12(s,s)=  Y1(s,s)Q1(s)

2
Ц

1+b2
+
(
Ц
 

1+b2
 
+b)

2p
Ц

1+b2

у
х
G 
 Y1(x,s)

x-s
Q1(xdx,

K21(s,s)=  Y1(s,s)Q2(s)

2
Ц

1+b2
-
(
Ц
 

1+b2
 
-b)

2p
Ц

1+b2

у
х
G 
 Y1(x,s)

x-s
Q2(xdx,

K22(s,s)=  Y2(s,s)Q2(s)

2
Ц

1+b2
-
(
Ц
 

1+b2
 
-b)

2p
Ц

1+b2

у
х
G 
 Y2(x,s)

x-s
Q2(xdx,

F1(s) = (
Ц
 

1+b2
 
-b) ж
з
и
fў+(s)+  f-(s)


Ц

1+b2
ц
ч
ш
Q1(s)+

+  1

p

у
х
G 
 Q1(s)

s-s
ж
з
и
fў+(s)+  f-(s)


Ц

1+b2
ц
ч
ш
 ds,
(23)

F2(s) = ж
з
и
fў+(s)-  f-(s)


Ц

1+b2
ц
ч
ш
Q2(s)-

-
(
Ц
 

1+b2
 
-b)

p

у
х
G 
  Q2(s)

s-s
ж
з
и
fў+(s)-  f-(s)


Ц

1+b2
ц
ч
ш
 ds.
(24)

Заметим, что плотности сингулярных интегралов в выражениях для функций Kpj(s,s), Fp(s)   (p=1,2, j=1,2) являются гельдеровыми функциями на G (причем плотности в Kpj(s,s) гельдеровы по обеим переменным). В частности, эти плотности гельдеровы на G по переменной x с показателем w = min{l,h, 1/2-h}   (равномерно по s в случае Kpj(s,s)) и обращаются в ноль, если x - конец G (так как функции Q1(x), Q2(x) принадлежат классу C0,h0(G) и обращаются в ноль на концах G). Из этих рассуждений и из свойств сингулярных интегралов [mus,§ 18] следует, что функции Kpj(s,s)   (p=1,2, j=1,2) являются гельдеровыми на G по обеим переменным. В частности, эти функции гельдеровы на G по переменной s с показателем

w = min{l,h, 1/2-h}

(25)

равномерно по s G. Если f+(s),f-(s) C0,l(G), то функции F1(s),F2(s) являются гельдеровыми на G с показателем w.

Очевидно, что если функции r1(s), r2(s) из пространства Cvq(G) с v (0,1], q [0,1) дают решение интегральных уравнений (21), (22), то эти функции представимы в виде rj(s)=rj*(s)/Qj(s),   j=1,2, где r1*(s),r2*(s) C0,w(G) ,  w беpется из (25). Поэтому будем искать функции r1(s), r2(s) именно в таком виде. Из этого представления, в частности, следует, что r1(s), r2(s) Cvq(G), где


q= max
{1/2+h,1-h},   v = м
п
н
п
о
min
{w,|2h-1/2|},
h 1/4,
w,
h = 1/4.
(26)

Ниже будем считать, что v и q определяются в (26). Отметим, что 0 < w 1/4, 0 < v 1/4, 1/2 < q < 1. Если b > 0, то 0 < h < 1/4 и q=1-h; если b < 0, то 1/4 < h < 1/2 и q=1/2+h; если же b = 0, то h = 1/4, q=3/4.

Далее введем функции rj*(s)=Qj(s)rj(s) C0,w(G) и операторы Kpj:


Kpj[v](s)=
у
х
G 
Kpj(s,s)Q-1j(s)v(sds,     p=1,2,  j=1,2.
(27)

В итоге уравнения (21), (22) сведутся к одному векторному уравнению относительно неизвестного вектор-столбца [`(r)]=(r1*(s),r2*(s))T, принадлежащего банахову пространству C0(GC0(G) с нормой
||[`(r)]||C0(GC0(G) = ||r1*||C0(G)+||r2*||C0(G) :


(I+R)

r
 
=

F
 
,      R= ж
з
з
и
K11
K12
K21
K22
ц
ч
ч
ш
,
(28)

где [`(F)]=(F1(s),F2(s))T C0(GC0(G), функции F1(s), F2(s) заданы в (23), (24), I - единичный оператор, отображающий пространство C0(GC0(G) в себя, операторы Kpj определены в (27). Из свойств гладкости ядер в уравнениях (28) и гладкости F1(s), F2(s) в (23) вытекает, что если r1*(s),r2*(s) C0(G) - решение уравнения (28), то r1*(s),r2*(s) C0,w(G) с w из (25).

Уравнение (28) является фредгольмовым в пространстве C0(GC0(G) и имеет единственное решение в этом пространстве. Доказательство этих утверждений проводится так же, как и для уравнения (32) в  [0].

Тем самым, решение задачи S  существует, единственно и дается формулой (6), в которой функции m(s),n(s) Cvq(G) даются формулами (11), где rj(s)=rj*(s)Q-1j(s),   j=1,2, а функции r1*(s), r2*(s) C0,w(G) определяются как единственное решение уравнения (28) в C0(GC0(G). Константа c в (6) находится из условия (16):


c=f+(a)-  1

2
Ц

1+b2
V й
л
 r2*

Q2
-(
Ц
 

1+b2
 
+b)  r1*

Q1
щ
ы
(x(a))-

-  1

2
Ц

1+b2
T й
л
 r1*

Q1
+(
Ц
 

1+b2
 
+b)  r2*

2Q2
щ
ы
(x(a)).
(29)

В дальнейшем существенным для нас будет то, что функции r1(s),r2(s) представимы в виде rj(s)=rj*(s)/Qj(s),   j=1,2, где r1*(s), r2*(s) C0,w(G) с w из (25). Кроме того, существенно то, что функции r1(s),r2(s) удовлетворяют уравнениям (18), (19), где функции F1(s), F2(s) принадлежат C0,l(G) и даются формулами (17), а B1, B2 берутся из (20).



2 Поведение градиента решения на концах контура



Через u будем обозначать решение задачи 
S , построенное в параграфе 1. Используя введенные выше обозначения, исследуем поведение u вблизи концов контура

G = {x: x=x(s)=(x1(s),x2(s)),  s ∈ [a,b]} .


Пусть x(d) - один из таких концов. Введем в окрестности x(d) полярную систему координат

x1=|x-x(d)|cosj,   x2=|x-x(d)|sinj.


Напомним, что
a(s) - угол между направлением оси Ox1 и вектором касательной tx в точке x(s) G. Будем считать, что j (a(d),a(d)+2p) , если d=a, и j (a(d)-p, a(d)+p), если d=b. (Полагаем по непрерывности a(a)=a(a+0),   a(b)=a(b-0)). Таким образом, угол j меняется непрерывно в окрестности точки x(d), разрезанной вдоль контура G.

Выделим в явном виде степенные особенности в функциях m(s),n(s):

1) на конце a:



m(s)=  1

2
Ц

1+b2
ж
и
 ra2(s)

|s-a|1-h
-(
Ц
 

1+b2
 
+b)  ra1(s)

|s-a|1/2-h
ц
ш
,
n(s)=  1

2
Ц

1+b2
ж
и
(
Ц
 

1+b2
 
+b)  ra2(s)

|s-a|1-h
+  ra1(s)

|s-a|1/2-h
ц
ш
,
(30)

2) на конце b:



m(s)=  1

2
Ц

1+b2
ж
и
 rb2(s)

|s-b|h
-(
Ц
 

1+b2
 
+b)  rb1(s)

|s-b|1/2+h
ц
ш
,
n(s)=  1

2
Ц

1+b2
ж
и
(
Ц
 

1+b2
 
+b)  rb2(s)

|s-b|h
+  rb1(s)

|s-b|1/2+h
ц
ш
,
(31)

Здесь функции

 

ra1(s)=

 r1*(s)


|s-b|1/2+h

=r1(s)|s-a|1/2-h,

ra2(s)=

 r2*(s)


|s-b|h

=r2(s)|s-a|1-h

 

(32)

являются гельдеровыми на G в окрестности a, а функции

 

rb1(s)=

 r1*(s)


|s-a|1/2-h

=r1(s)|s-b|1/2+h,

rb2(s)=

 r2*(s)


|s-a|1-h

=r2(s)|s-b|h

 

(33)

являются гельдеровыми на G в окрестности b. Функции r1*(s), r2*(s) из C0,w(G) - решение уравнения (28)   (w берется из (25)).

Пользуясь результатами из [mus] о поведении интегралов типа Коши на концах контура G, приходим к следующей теореме.

Теорема 1. Пусть x x(d) X, где d=a или d=b; тогда в окрестности точки x(d) для производных решения задачи S справедливы формулы:

1) при d=a:


  ∂ u

x1
к
к


xx(a) 
= -  ra1(a)sinJa(1/2-h)

2|x-x(a)|1/2-h
+

+(
Ц
 

1+b2
 
+b)  ra2(a)sinJa(1-h)

2|x-x(a)|1-h
+O(1),

  ∂ u

x2
к
к


xx(a) 
=  ra1(a)cosJa(1/2-h)

2|x-x(a)|1/2-h
-

-(
Ц
 

1+b2
 
+b)  ra2(a)cosJa(1-h)

2|x-x(a)|1-h
+O(1);

2) при d=b:


  ∂ u

x1
к
к


xx(b) 
=  rb1(b)cosJb(1/2+h)

2|x-x(b)|1/2+h
+(
Ц
 

1+b2
 
+b)  rb2(b)sinJb(h)

2|x-x(b)|h
+O(1),

  ∂ u

x2
к
к


xx(b) 
=  rb1(b)sinJb(1/2+h)

2|x-x(b)|1/2+h
-(
Ц
 

1+b2
 
+b)  rb2(b)cosJb(h)

2|x-x(b)|h
+O(1).

Здесь

Ja(g)=gj+(1-g)a(a),   Jb(g)=gj+(1-g)a(b)-ph,

через O(1) обозначены функции, непрерывные в окрестности точки x(d), разрезанной вдоль контура G. (Здесь и далее подразумевается, что функции O(1) непрерывны и в самой точке x(d).) Функции rd1(s) и rd2(s)   (d=a,b) введены в (32), (33).

Доказательство. Наряду с декартовыми координатами x=(x1,x2) введем комплексную координату [(x)\tilde]=x1+ix2. Используем представление производных углового и логарифмического потенциалов через интеграл типа Коши на комплексной плоскости [K1]:


  ∂ V[m](x)

x1
=-БW[m](
~
x
 
),       ∂ V[m](x)

x2
=-ВW[m](
~
x
 
),
  ∂ T[n](x)

x1
=ВW[n](
~
x
 
),       ∂ T[n](x)

x2
=-БW[n](
~
x
 
),
(34)

где


W[m](
~
x
 
)=  1

2pi

у
х
G 
m(s)e-ia(s)
d
~
y
 
(s)

~
y
 
(s)-
~
x
 
,   
~
y
 
=y1+iy2.

Из (6) и (34) вытекает, что для изучения поведения u достаточно изучить поведение интегралов типа Коши

W[ra1(s)/|s-a|1/2-h](

~

x

 

),   W[ra2(s)/|s-a|1-h](

~

x

 

),


W[rb1(s)/|s-b|1/2+h](

~

x

 

),   W[rb2(s)/|s-b|h](

~

x

 

)

на комплексной плоскости вблизи концов контура G.

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


W й
л
 rb2(s)

|s-b|h
щ
ы
(
~
x
 
) = W[Yb2(s)](
~
x
 
)+rb2(b)W й
л
 1

|s-b|h
щ
ы
(
~
x
 
),
(35)

где Yb2(s)=(rb2(s)-rb2(b))/|s-b|h.

Пусть d=a либо d=b. Преобразуем плотность интеграла W[1/|s-d|g]( [(x)\tilde]) с g (0,1):

 

 e-ia(s)


|s-d|g

=

~

f

 

(

~

y

 

(s))


(

~

y

 

(s)-

~

x

 

(d))g

,

где


~
f
 
(
~
y
 
(s))=f(s)=e-ia(s) й
к
л
ж
з
и
~
y
 
(s)-
~
x
 
(d)

|s-d|
ц
ч
ш
g

 
щ
ъ
ы
=

= e-ia(s) к
к
к
~
y
 
(s)-
~
x
 
(d)

s-d
к
к
к
g

 
exp{igarg(
~
y
 
(s)-
~
x
 
(d))}.
(36)

Пользуясь техникой, развитой в [K1,лемма 1], нетрудно показать, что функция [(g)\tilde](s)=([(y)\tilde](s)-[(x)\tilde](d))/|s-d| принадлежит классу C1(G) и не обращается в ноль ни при каком s G. Заметим, что функция ([(g)\tilde])g гельдерова с показателем 1 на любой гладкой дуге, не содержащей точки [(g)\tilde]=0. Согласно теореме о гельдеровости сложной функции можно утверждать, что функция в квадратных скобках в (36) гельдерова по s на контуре G в окрестности d c показателем 1. Очевидно, e-ia(s)=(y1(s)-iy2(s)) C1(G). Поэтому функция f(s) гельдерова по s на G вблизи конца d с показателем 1, т.е. для любых s1,s2, лежащих на G в окрестности точки d, имеем

|

~

f

 

(

~

y

 

(s2))-

~

f

 

(

~

y

 

(s1))|=|f(s2)-f(s1)| const|s2-s1|.

(37)

Покажем, что функция [(f)\tilde]([(y)\tilde](s))=f(s) гельдерова по [(y)\tilde] на контуре G в окрестности [(x)\tilde](d). В силу [K1,лемма 1] функция |s2-s1|/|[(y)\tilde](s2)-[(y)\tilde](s1)| принадлежит C1(G×G) C0(G×G) по обеим переменным s1, s2 и поэтому равномерно ограничена по s1, s2 G константой c*:

 

 |s2-s1|


|

~

y

 

(s2)-

~

y

 

(s1)|

c*.

(38)

Отсюда: |s2-s1| c*|[(y)\tilde](s2)-[(y)\tilde](s1)|. Используя эту оценку в (37), находим

|

~

f

 

(

~

y

 


2 

)-

~

f

 

(

~

y

 


1 

)| c0|

~

y

 


2 

-

~

y

 


1 

|,

(39)

где [(y)\tilde]2=[(y)\tilde](s2),   [(y)\tilde]1=[(y)\tilde](s1) - произвольные точки, лежащие на контуре G в окрестности конца [(x)\tilde](d), и c0 - некоторая константа. Тем самым, функция [(f)\tilde]([(y)\tilde]) гельдерова на контуре G по [(y)\tilde] вблизи конца [(x)\tilde](d) с показателем 1.

Преобразуем интеграл W[1/|s-d|g]([(x)\tilde]) следующим образом:


W й
л
 1

|s-d|g
щ
ы
(
~
x
 
) =  1

2pi

у
х
G 
~
f
 
(
~
y
 
)

(
~
y
 
-
~
x
 
(d))g
d
~
y
 

~
y
 
-
~
x
 
=

=
~
f
 
(
~
x
 
(d))

2pi

у
х
G 
 1

(
~
y
 
-
~
x
 
(d))g
d
~
y
 

~
y
 
-
~
x
 
+  1

2pi

у
х
G 
~
f
 
(
~
y
 
)-
~
f
 
(
~
x
 
(d))

(
~
y
 
-
~
x
 
(d))g
d
~
y
 

~
y
 
-
~
x
 
,

где функция [(f)\tilde]([(y)\tilde]) введена в (36). Второе слагаемое в последней формуле в силу оценки (39) и результатов [mus,§ 22] представляет функцию O(1), непрерывную в окрестности x(d), разрезанной вдоль контура G. К первому слагаемому применим результаты [mus]. Поведение интеграла типа Коши


I(
~
x
 
)=  1

2pi

у
х
G 
 1

(
~
y
 
-
~
x
 
(d))g
d
~
y
 

~
y
 
-
~
x
 
,

где g (0,1), изучалось в [mus,§ 23]. Для [(x)\tilde] достаточно близких к [(x)\tilde](d), но не расположенных на контуре G, справедливо представление [mus,§ 23]:

I(

~

x

 

)=(-1)j(d)

 exp{(-1)j(d)gpi}


2isingp

 

 1


(

~

x

 

-

~

x

 

(d))g

+O(1),

(40)

где ([(x)\tilde]-[(x)\tilde](d))g обозначает ветвь, голоморфную вблизи [(x)\tilde](d) на разрезанной вдоль контура G плоскости и принимающую значение ([(x)\tilde](s)-[(x)\tilde](d))g на левой стороне разреза; под O(1) понимается обозначение, введенное в теореме 1; j(a)=2,   j(b)=1.

Если s G и s d, то

|[(y)\tilde](s)- [(x)\tilde](d)|/|s-d → 1,       arg([(y)\tilde](s)-[(x)\tilde](d))a (d)+p(2-j(d)).


Учитывая эти соотношения, вычислим предел функции [(
f)\tilde]([(y)\tilde]) из (36) при [(y)\tilde] [(x)\tilde](d) и [(y)\tilde] G. Используя (40), получим


W й
л
 1

|s-d|g
щ
ы
(
~
x
 
) = (-1)j(d)  exp{(-1)j(d)gpi}

2isingp
 1

(
~
x
 
-
~
x
 
(d))g
×

×exp{i ж
и
-a(d)+g(a(d)+p(2-j(d))) ц
ш
}+O(1) =

=  1

2singp
 1

|x-x(d)|g
exp(iqd(g))+O(1),
(41)

где qa(g)=-gj-(1-g)a(a)+(g-1/2)p и qb(g)=-gj-(1-g)a(b)+p/2.

Теперь изучим поведение функции W[Yb2(s)]([(x)\tilde]). Для этого обратимся к исследованию функции Yb2(s). Заметим, что функция

r2(s)=rb2(s)/|s-b|h


удовлетворяет уравнению (
19). Используя это уравнение и результаты [mus,§ 22], можно убедиться в том, что справедлива

Лемма 1. Если F2(s) из (19) - гельдерова функция на G, то Yb2(s) также является гельдеровой функцией переменной s на G вблизи b. Кроме того, Yb2(b)= limsbYb2(s)=0.

Доказательство леммы будет дано в следующем параграфе. Из леммы 8 и оценки (38) вытекает, что функция [(Y)\tilde]b2( [(y)\tilde](s))=Yb2(s), рассматриваемая как функция комплексной переменной [(y)\tilde], гельдерова на контуре G вблизи [(x)\tilde](b). Кроме того, [(Y)\tilde]b2([(x)\tilde](b))=0. Тогда согласно [mus,§ 22], функция W[Yb2(s)]([(x)\tilde]) непрерывно продолжима на конец [(x)\tilde](b). Используя обозначение O(1), введенное в теореме 1, при [(x)\tilde] [(x)\tilde](b) имеем W[Yb2(s)]=O(1).

Возвращаясь к исходному интегралу (35) и учитывая (41), для [(x)\tilde] [(x)\tilde](b) получим


W й
л
 rb2(s)

|s-b|h
щ
ы
(
~
x
 
) =  rb2(b)

2sinph
 1

|x-x(b)|h
eiqb(h)+O(1).
(42)

Проводя аналогичные рассуждения и выкладки, нетрудно показать, что при [(x)\tilde] [(x)\tilde](d)


W й
л
 rdj(s)

|s-d|g
щ
ы
(
~
x
 
) =  rdj(d)

2sinpg
 1

|x-x(d)|g
eiqd(g)+O(1),
(43)

где d=a,   j=2,   g = 1-h, либо d=a,   j=1,   g = 1/2-h, либо d=b,   j=1,   g = 1/2+h.

Используя формулы (30), (42), (43), для [(x)\tilde] [(x)\tilde](d) имеем:

1) при d=a:


W[m](
~
x
 
)=  1

2
Ц

1+b2
ж
и
 ra2(a)

|x-x(a)|1-h
 eiqa(1-h)

2sinph
-

-  ra1(a)

|x-x(a)|1/2-h
 eiqa(1/2-h)

2sinph
ц
ш
+O(1),

W[n](
~
x
 
)=  1

2
Ц

1+b2
ж
и
 ra2(a)

|x-x(a)|1-h
 eiqa(1-h)cosph

2sin2ph
+

+  ra1(a)

|x-x(a)|1/2-h
 eiqa(1/2-h)

2cosph
ц
ш
+O(1);

2) при d=b:


W[m](
~
x
 
)=  1

2
Ц

1+b2
ж
и
 rb2(b)

|x-x(b)|h
 eiqb(h)

2sinph
-

-  rb1(b)

|x-x(b)|1/2+h
 eiqb(1/2+h)

2sinph
ц
ш
+O(1),

W[n](
~
x
 
)=  1

2
Ц

1+b2
ж
и
 rb2(b)

|x-x(b)|h
 eiqb(h)cosph

2sin2ph
+

+  rb1(b)

|x-x(b)|1/2+h
 eiqb(1/2+h)

2cosph
ц
ш
+O(1).

Используя (6) и (34), для [(x)\tilde] [(x)\tilde](d) получим:

1) при d=a:


  ∂ u

x1
к
к


xx(a) 
=

=  ra1(a)

2
Ц
 

1+b2
 
|x-x(a)|1/2-h
м
н
о
 sinqa(1/2-h)

2sinph
+  cosqa(1/2-h)

2cosph
ь
э
ю
+

+  ra2(a)

2
Ц
 

1+b2
 
|x-x(a)|1-h
м
н
о
-  sinqa(1-h)

2sinph
+  cosphcosqa(1-h)

2sin2ph
ь
э
ю
+O(1),

  ∂ u

x2
к
к


xx(a) 
=

=  ra1(a)

2
Ц
 

1+b2
 
|x-x(a)|1/2-h
м
н
о
 cosqa(1/2-h)

2sinph
-  sinqa(1/2-h)

2cosph
ь
э
ю
+

+  ra2(a)

2
Ц
 

1+b2
 
|x-x(a)|1-h
м
н
о
-  cosqa(1-h)

2sinph
-  cosphsinqa(1-h)

2sin2ph
ь
э
ю
+O(1);

2) при d=b:


  ∂ u

x1
к
к


xx(b) 
=

=  rb1(b)

2
Ц
 

1+b2
 
|x-x(b)|1/2+h
м
н
о
 sinqb(1/2+h)

2sinph
+  cosqb(1/2+h)

2cosph
ь
э
ю
+

+  rb2(b)

2
Ц
 

1+b2
 
|x-x(b)|h
м
н
о
-  sinqb(h)

2sinph
+  cosqb(h)cosph

2sin2ph
ь
э
ю
+O(1),

  ∂ u

x2
к
к


xx(b) 
=

=  rb1(b)

2
Ц
 

1+b2
 
|x-x(b)|1/2+h
м
н
о
 cosqb(1/2+h)

2sinph
-  sinqb(1/2+h)

2cosph
ь
э
ю
+

+  rb2(b)

2
Ц
 

1+b2
 
|x-x(b)|h
м
н
о
-  cosqb(h)

2sinph
-  cosphsinqb(h)

2sin2ph
ь
э
ю
+O(1).

После элементарных тригонометрических преобразований приходим к утверждению теоремы.

Полученные формулы показывают, что частные производные решения задачи S обычно имеют на концах разреза степенные особенности с показателем q=max{1/2+h,1-h}.



3  Доказательство леммы 1



Докажем лемму 1, которая была использована в теореме 1.

Напомним, что Yb2(s)=(rb2(s)-rb2(b))/|s-b|h, а rb2(s)=r2(s)|s-b|h . Как было отмечено выше, функция r2(s) удовлетворяет уравнению (19). Поэтому для rb2(s) справедливо тождество:


rb2(s)=F2(s)|s-b|h-  1

|s-a|1-h

у
х
G 
Q2(x)
~
F
 

2 
(x)

x-s
 dx-  B2sinph

|s-a|1-h
,    s О G,
(44)

где F2(s)=({1+b2}+b)F2(s)/(2{1+b2}),   [( F)\tilde]2(s)=F2(s)/(2{1+b2});   F2(s),[(F)\tilde]2(s) C0,l(G). Заметим, что для s, лежащих на G в окрестности точки b, справедливо равенство

 

 1


|s-a|1-h

=

 1


|b-a|1-h

+O(|s-b|),

где через O(|s-b|) обозначена функция, гельдеровая на G в окрестности конца b с показателем 1. Преобразуем выражение (44) следующим образом:


rb2(s)=F2(s)|s-b|h-  1

(b-a)1-h

у
х
G 
Q2(x)
~
F
 

2 
(x)

x-s
 dx-

-  B2sinph

(b-a)1-h
+O(|s-b|)X(s),
(45)

где функция X(s) является гельдеровой на G в окрестности точки b с показателем min{h,l}. Теперь положим в (44) s=b и найдем из (44), (45) функцию [rb2(s)-rb2(b)], приводя к общему знаменателю подынтегральные выражения:


rb2(s)-rb2(b) = F2(s)|s-b|h+  b-s

(b-a)1-h

у
х
G 
(x-a)1-h
~
F
 

2 
(x)

(b-x)1-h(x-s)
 dx+

+O(|s-b|)X(s),

Тогда функция Yb2(s) имеет вид:


Yb2(s)=F2(s)-  (b-s)1-h

(b-a)1-h

у
х
G 
(x-a)1-h
~
F
 

2 
(x)

(b-x)1-h(x-s)
 dx+O(|s-b|1-h)X(s).

Используя результаты [mus,§ 22, п. 4] и учитывая гельдеровость функций F2(s), [(F)\tilde]2(s) на G, можно утверждать, что функция Yb2(s) является гельдеровой на G вблизи b. Чтобы найти Yb2(b)=lims bYb2(s), учтем [ mus ,§ 22, п. 2, п. 4], что для функции f(x), гельдеровой на G в окрестности конца d, справедлива формула (0 < g < 1):



lim
sd 
ж
и
 (s-d)g

2p

у
х
G 
 f(x)

(x-d)g(x-s)
 dx ц
ш
=  (-1)j(d)ctgpg

2
f(d).

Здесь j(a)=2,  j(b)=1. Таким образом,


Yb2(b)=F2(b)+  2p

2(b-a)1-h
ctgp(1-h)(b-a)1-h
~
F
 

2 
(b) =
=
(
Ц
 

1+b2
 
+b)

2
Ц

1+b2
F2(b)+  -ctgph

2
Ц

1+b2
F2(b)=0

(использована формула ctgph = ({1+b2}+b)). Лемма доказана.



4  Исчезновение особенности
u



Как следует из формул теоремы 1, чтобы
u(x) не имел особенностей в точке x(d)   (d=a или d=b), необходимо и достаточно, чтобы rdj(d)=0 для j=1,2. Особенностей не будет на обоих концах разреза G, если выполнены условия

raj(a)=0,    rbj(b)=0,     j=1,2,

(46)

которые составляют четыре интегральных требования на функции f+(s), f-(s).

Лемма 2. 1) Если rd1(d)=rd2(d)=0, где d=a или d=b, то градиент решения задачи S непрерывен на конце x(d) разреза G.

2) Если rd1(d)=rd2(d)=0 для d=a и d=b, то градиент решения задачи S непрерывен на обоих концах разреза G.

Представляется интересным изучить следующий вопрос: при каких функциях f+(s), f-(s) из (3), (4) особенности градиента решения задачи S на концах разреза G исчезают? Аналитическое исследование этого вопроса для разреза произвольной формы затруднительно. Мы будем исследовать один частный случай.

Пусть G представляет собой прямолинейный отрезок, расположенный под углом q к оси Ox1: G = {x:x(s)=(scosq,ssinq),  s [a,b]}. Тогда Y1(s,s)=Y2(s,s) 0 для s,s G, и уравнения (12) - (16) имеют решения


r1(s)=
(
Ц
 

1+b2
 
-b)

2
Ц

1+b2
F1(s)+

+  1

2p
Ц
 

1+b2
 
Q1(s)

у
х
G 
 F1(x)Q1(x)

x-s
dx-  B1cosph

Q1(s)
,

r2(s)=
(
Ц
 

1+b2
 
+b)

2
Ц

1+b2
F2(s)-

-  1

2p
Ц
 

1+b2
 
Q2(s)

у
х
G 
 F2(x)Q2(x)

x-s
dx-  B2sinph

Q2(s)
,

где s G,


Q1(x)=(x-a)1/2-h(b-x)1/2+h,   Q2(x)=(x-a)1-h(b-x)h,     x О G,

F1(s)=2(
Ц
 

1+b2
 
fў+(s)+f-(s)),

F2(s)=2(
Ц
 

1+b2
 
-b)(
Ц
 

1+b2
 
fў+(s)-f-(s)),

константы B1, B2 берутся из (20). Находя функции rdj(s) по формулам (32), (33), имеем:


ra1(a)=  1

2p
Ц
 

1+b2
 
(b-a)1/2+h
b
у
х
a 
F1(x) ж
и
 b-x

x-a
ц
ш
1/2+h

 
dx-  B1cosph

(b-a)1/2+h
,

rb1(b)=-  1

2p
Ц
 

1+b2
 
(b-a)1/2-h
b
у
х
a 
F1(x) ж
и
 x-a

b-x
ц
ш
1/2-h

 
dx-  B1cosph

(b-a)1/2-h
,

ra2(a)=-  1

2p
Ц
 

1+b2
 
(b-a)h
b
у
х
a 
F2(x) ж
и
 b-x

x-a
ц
ш
h

 
dx-  B2sinph

(b-a)h
,

rb2(b)=  1

2p
Ц
 

1+b2
 
(b-a)1-h
b
у
х
a 
F2(x) ж
и
 x-a

b-x
ц
ш
1-h

 
dx-  B2sinph

(b-a)1-h
.

Покажем, что требования (46) можно удовлетворить все одновременно, так что u(x) будет непрерывен в точках x(a) и x(b). Будем искать F1(x), F2(x) в виде

F1(x)=(x-a)1/2+h(b-x)1/2-hG1(x) C0,l(G),


F2(x)=(x-a)h(b-x)1-hG2(x) C0,l(G),     l ∈ (0,1].

(47)

Тогда условия (46) для функций F1(x), F2(x) превратятся в следующие условия для функций G1(x), G2(x):


b
у
х
a 
Gj(x)(x-adx = Alj,    b
у
х
a 
Gj(x)(b-xdx = -Alj,     j=1,2,

или


b
у
х
a 
Gj(xdx = 0,    b
у
х
a 
xGj(xdx = Alj,    j=1,2,
(48)

где l1=8p(1+b2)sinphcos2ph,   l2=8p(1+b2)sin3ph.

Заметим (см. (47)), что для выполнения условия F1(s),F2(s) C0,l(G) с некоторым (не заданным заранее) l (0,1] достаточно, чтобы G1(s), G2(s) Cw0h0(G) с некоторым w0 (0,1] и h0=min{h,1/2-h}.

Лемма 3. Пусть G = {x:x=x(s)=(scosq,ssinq),  s [a,b]}, функции F1(s), F2(s) определены в (47), а G1(s), G2(s) удовлетворяют условиям (48). Тогда градиент решения задачи S непрерывен на концах разреза G.

Приведем несколько способов построения функций G1(s), G2(s). Не ограничивая общности, будем считать, что отрезок [a, b] линейной заменой переменных приведен к симметричному относительно точки s=0 виду.

(I) Рассмотрим случай a=-b < 0.

1) Пусть G1(s), G2(s) - четные функции и A=0. Тогда для выполнения однородных условий (48) достаточно потребовать, чтобы 0b Gj(xdx = 0 для j=1,2. Например, можно выбрать G1(x)=G2(x)=C(b-2|x|), где C - произвольная константа.

2) Пусть G1(s), G2(s) - нечетные функции. Тогда, если A=0, для выполнения (48) достаточно потребовать, чтобы 0b xGj(xdx = 0 для j=1,2. В общем случае, когда A - любое число, достаточно потребовать, чтобы


b
у
х
0 
x ж
и
Gj(x)-  (2k+1)lj A

2b2k+1
x2k-1 ц
ш
 dx = 0,

где j=1,2, и k - любое фиксированное натуральное число. В частности, если A 0, то можно выбрать

Gj(s) = (2k+1)2-1lj A x2k-1/b2k+1,    j=1,2.

Аналогично рассматриваются случаи, когда одна из функций G1(s), G2(s) четная, а другая - нечетная.

(II) Пусть далее a=-1,   b=1. Тогда требования (48) удобно записать в терминах полиномов Лежандра [gr2, gr3]. Напомним, что полиномы Лежандра Pn(s),   n=0,1,, образуют замкнутую ортогональную систему в пространстве L2(-1,1). При этом ||Pn||2=2/(2n+1),   P0(s)=1,   P1(s)=s. Полагая


Gj(s)=
е
n=0 
Anj Pn(s),   j=1,2,

и считая, что ряд сходится в среднем, сразу же определим из (48): A0j=0,   A1j=3A lj/2,   j=1,2. Остальные коэффициенты Ajn остаются неопределенными (j=1,2;  n=2,3,).

Итак, пусть Anj - произвольные коэффициенты (j=1,2;  n=2,3,), для которых ряды n=2An1 Pn(s) и n=2An2 Pn(s) сходятся в среднем. Тогда функции


Gj(s)=  3Alj

2
s+
е
n=2 
Anj Pn(s),    j=1,2,
(49)

удовлетворяют условиям (48).

Требования гладкости, накладываемые на функции G1(s), G2(s) в (47), таковы, что эти функции могут иметь в точках s=1 интегрируемые степенные особенности. Но мы будем для упрощения рассуждений рассматривать сходимость рядов в (49) к гельдеровым функциям.

Используя достаточный признак равномерной сходимости Вейерштрасса и факт равномерной ограниченности Pn(s) на отрезке [-1,1] [gr3,гл.4, § 2], можно доказать следующее достаточное условие равномерной сходимости рядов в (49).

Лемма 4. Пусть Anj=O(1/n2+e) при n →∞ , где e > 0,   j=1,2. Тогда ряды в (49) сходятся равномерно, причем функции G1(s), G2(s) принадлежат классу C1(G) и удовлетворяют условиям (48).

Подытоживая пункт (II), сформулируем следующую теорему.

Теорема 2. Пусть

1) G = {x:x=x(s)=(scosq,ssinq),  s [-1,1]},

2) функции F1(s), F2(s) определяются формулами (47), (49),

3) Ajn - произвольные коэффициенты (n=2,3,), такие что ряды в (49) для j=1 и j=2 сходятся в среднем к гельдеровым функциям.

Тогда градиент решения задачи S непрерывен на концах разреза G.

References

[0]

П.А. Крутицкий, А.И. Cгибнев. Смешанная задача для уравнения Лапласа вне разрезов на плоскости с заданием условия Дирихле и условия с косой производной на разных сторонах разрезов. // Препринт Инст. прикл. матем. им. М.В. Келдыша РАН, 2004, N 8, 28 с.

[2]

Крутицкий П.А., Сгибнев А.И. Метод интегральных уравнений в смешанной задаче Дирихле-Неймана для уравнения Лапласа вне разрезов на плоскости. // Дифф. уравнения. 2001. Т. 37. N 10. С. 1299-1310.

[3]

Гостева А.С., Крутицкая Н.Ч., Крутицкий П.А. Смешанная задача в замагниченной полупроводниковой пленке с разрезами вдоль прямой. // Фунд. и приклад. матем. 2000. Т. 6. Вып. 4. С. 1061-1073.

[G]

Габов С.А. Угловой потенциал и его некоторые приложения. // Матем. сб. 1977. Т. 103 (145). N 4. C. 490-504.

[K1]

Крутицкий П.А. Задача Дирихле для уравнения Гельмгольца вне разрезов на плоскости. // Журнал вычислительной математики и математической физики. 1994. Т. 34. N 8-9. C. 1237-1258.

[K2]

Крутицкий П.А. Задача Неймана для уравнения Гельмгольца вне разрезов на плоскости. // ЖВМ и МФ. 1994. Т. 34. N 11. С. 1652-1665.

[mus]

Мусхелишвили Н.И. Сингулярные интегральные уравнения. - М.: Наука, 1968.

[kol]

Колмогоров А.Н., Фомин С.В. Элементы теории функций и функционального анализа. М.: Наука, 1972.

[cr]

Крейн С.Г. (ред.). Функциональный анализ. М.: Наука, 1964.

[lif]

Лифанов И.К. Метод сингулярных интегральных уравнений и численный эксперимент. - М.: ТОО "Янус", 1995.

[s10]

Крутицкий П.А. Смешанная задача для уравнения Лапласа вне разрезов на плоскости. // Дифф. уравнения. 1997. Т. 33. N 9. С. 1181-1190.

[gr2]

Никифоров А.Ф., Уваров В.Б. Специальные функции мат. физики. - М.: Наука, 1984.

[gr3]

Суетин П.К. Классические ортогональные многочлены. - М.: Наука, 1979.

 


File translated from TEX by TTH, version 3.40.
On 21 Dec 2004, 20:02.