CC BY
34
УДК 517.977
К ВОПРОСУ ОБ ОТЛИЧИЯХ В ПОВЕДЕНИИ РЕШЕНИЙ ЛИНЕЙНОГО И НЕЛИНЕЙНОГО УРАВНЕНИЙ ТЕПЛОПРОВОДНОСТИ1
Л.И. Рубина2, О.Н. Ульянов3
Ранее предложенным геометрическим методом исследования нелинейных уравнений в частных производных исследуются линейное и нелинейное уравнения теплопроводности. Показано, чем обусловлено отличие в поведении решений рассматриваемых уравнений и что в случае нелинейного уравнения приводит к обострению. Выделен некоторый класс решений линейного уравнения, представляющий поверхности уровня рассматриваемого нелинейного уравнения теплопроводности.
Ключевые слова: нелинейные уравнения в частных производных, уравнения теплопроводности, точные решения, поверхности уровня.
Введение
Известно (см., например, серию работ [1-5]), что нелинейные уравнения теплопроводности и некоторые другие нелинейные уравнения в частных производных описывают режимы, значительно отличающиеся от тех режимов, которые наблюдаются, если используется линейная модель процесса. Часто в случае нелинейных моделей наблюдаются так называемые катастрофы, при которых решение неограниченно возрастает за конечный промежуток времени [1].
В данной работе приведены некоторые точные решения для линейного и нелинейного уравнений теплопроводности [1], которые являются хорошей иллюстрацией описанной выше проблемы. Полученные одним и тем же геометрическим методом [6, 7] решения имеют одинаковые поверхности уровня, но их поведение значительно отличается, так как сами решения удовлетворяют разным обыкновенным дифференциальным уравнениям (ОДУ), к которым сведены первоначальные уравнения теплопроводности. Аналогично при рассмотрении характеристик уравнений отличие в представлении их решений связано с отличием ОДУ, которые задают условия совместности для получения решения вдоль характеристик. Вид выписанных явно точных решений делает наглядными их отличия и позволяет легко увидеть, почему в нелинейном случае наблюдается обострение, а в линейном случае обострение отсутствует. Приводятся картины течений нелинейного уравнения теплопроводности в зависимости от параметров задачи, которые показывают, как можно отодвинуть по времени момент обострения в решении или, переключаясь на другое решение при подходе к обострению, избежать катастрофы.
Сведение линейного и нелинейного уравнения теплопроводности к обыкновенным дифференциальным уравнениям
Будем рассматривать два (линейное и нелинейное) уравнения теплопроводности [1]
ut = kuxx + qu, к = const, q = const. (1)
ut = kuxx + qu -au3, a = const. (2)
Здесь и далее нижние индексы указывают на независимую переменную, по которой вычисляется производная.
В уравнениях (1) и (2) сделаем замену [6] u = Qx , тогда получим уравнения
Qxt = kQxxx + qQx, (3)
Qxt = kQxxx + qQx -aQl. (4)
1 Работа выполнена в рамках программы межрегиональных и межведомственных исследований УрО РАН (проект 12-С-1-1001).
2 Рубина Людмила Ильинична - кандидат физико-математических наук, старший научный сотрудник, ИММ УрО РАН.
E-mail: [email protected]
3 Ульянов Олег Николаевич - кандидат физико-математических наук, старший научный сотрудник, ученый секретарь института, ИММ УрО РАН.
E-mail: [email protected]
Рубина Л.И., К вопросу об отличиях в поведении решений
Ульянов О.Н. линейного и нелинейного уравнений теплопроводности
Считаем, что существует такая система координат [6, 7], в которой функция Q зависит от одной независимой переменной: Q = Q(¥(x, t)), тогда ¥(x, t) = const - поверхность уровня решения Q(x, t). В этой системе координат уравнения (3) и (4) соответственно имеют вид
kQ"r3 + Q (-Vx¥t + 3k¥x¥xx) + Q (-¥xt + k¥xxx + q¥x) = ^ (5)
(kQ"-aQ'3)¥l + Q" (-¥x¥t + 3k¥x¥xx) + Q'(-¥xt + k¥xxx + q¥x) = 0. (6)
Здесь и далее штрих (') указывает на производную по переменной ¥ . Сравнивая (5) и (6), замечаем, что для определения функции ¥(x, t) имеем одни и те же уравнения, если положим, что ¥х Ф 0 [6, 7]
~¥ + f¥xx = f1(¥), ~¥xt + к¥Т + q¥x = f2(¥). (7)
¥x ¥x
Здесь f1(¥), f2(¥) - произвольные функции. Будем устанавливать, когда система уравнений (7) совместна. Из первого уравнения системы ¥xx = (¥t + f¥x)/(3k). Продифференцируем это соотношение по переменной x и полученную производную ¥xxx подставим во второе уравнение системы. Разрешив полученное соотношение относительно производной ¥xt , имеем
¥xx = ¥t +f¥x , ¥xt = 3 q¥x + f3¥l + f1¥x¥t, f3 =1 f'+ ^ f1 - 3 f2 . (8)
3k 2 3k 2 3k 2
Соотношения (8) будут задавать производные второго порядка одной и той же функции, если смешанные производные равны (¥xxt = ¥xtx). Это условие выполняется, когда
¥tt = 3 q¥t + -1 ¥ + 3f3¥t¥x2 + f4¥x4, f4 = 2f1f3 + 3kf3. (9)
2 3k
Чтобы все полученные вторые производные были производными ОДНОЙ функции ¥xt , их смешанные производные должны быть равны. Требуем, чтобы ¥xtt =¥ttx. Условие будет выполняться, если
6f4¥t¥x2 + 6f3¥t2 + 3k(f4 + flf^^k¥ = 0 . (10)
Замечаем, что если f3 = 0, соотношение (10) обращается в тождество. Полученные результаты приводят к следующему утверждению:
Утверждение 1. Система (7) совместна, если вторые производные функции ¥(x, t) определяются из соотношении (8), (9) и f2 = f1 /3 + 2f-^ /(9k), f1(¥) - произвольная функция. В этом случае уравнения (5), (6) имеют вид соответственно
kQ" + Q fx + Q [f / 3 + 2f12 /(9k)] = 0; (kQ''' - aQ'3) + Q fx + Q [f / 3 + 2f12 /(9k)] = 0 .
Если f3 Ф 0, то, выписав дифференциальные следствия соотношения (10) и подставив в полученные выражения вторые производные из (8), (9), получим, что ¥x = g1(¥), ¥t = g2(¥). Потребовав равенства смешанных производных, будем иметь g2 (¥) = Cg1 (¥), C = const. Тогда справедливо
Утверждение 2. Если f3 Ф 0, то ¥ = ¥(ax + bt), a = const, b = const. Тогда, так как u = Qx = Q (¥)¥x (¥), mo u = u(ax + bt) и уравнения (1) и (2) имеют вид соответственно
-buz + ka2uzz + qu = 0, -buz + ka2uzz + qu - au3 = 0, z = ax + bt.
Отыскание поверхностей уровня для уравнений теплопроводности
Положим, что произвольное f1 = const и f2 = 2f12 /(9k) = const. Тогда условие (10) выполняется. Из соотношения (9) находим ¥t, считая, что x - параметр. Предварительно полагая, что ¥t Ф 0, запишем (9) в виде (1/¥t)t +1,5q(1/¥t) = f1 /(3k). Решаем это линейное уравнение и получаем, что ¥t = 1,5 qC (x)exp(3qt/2)/[1 - fC (x)exp(3qt)/(3k)]. Затем из второго соотношения (8) определяем ¥x . Получаем, что
/x = Cx exp(3qt / 2) /[1 - f1C(x) exp(3qt / 2) /(3k)].
Требуя тождественного выполнения первого соотношения (8), приходим к уравнению для определения C(x): Cxx = qC/(2k). Отсюда C(x) = Axexp[±xjq/(2k)] + A2exp[+x^Jq/(2k)]. Здесь
A1 = const, A2 = const. Далее положим, что A1 = 1, A2 = 0 и выпишем для этого случая окончательный вид / и /х
3q exp[(3qt /2) ± x,Jq /(2k) ]/2
/t =-
±yj q /(2k) exp[(3qt/2) ± x*Jq /(2k) ] f1exp[(3qt/2) ± x«Jq /(2k) ]/(3k) ¥x 1 - f1 exp[(3qt/2) ± x*Jq /(2k) ]/(3k)
В выражении для /х знак перед числителем выражения совпадает со знаком в показателе степени экспоненты. Далее, находим, что поверхность уровня в этом случае имеет вид
y/ = -3k ln{1 - f1 exp[(3qt /2) ±x^q /(2k) ]/(3k)}/ f1. (11)
Положим f2 = 0. Тогда f/) = 3k/(2/ + /o), /o = const (см. утверждение 1). Требуем далее, чтобы вторые производные функции /(x, t) удовлетворяли условиям (8), (9). Решая уравнение (9) и второе уравнение из (8) и считая при этом, что x - параметр, получаем, что
/t =V2/ + /o [M(x) + 3q^/2/ + /o/2], /x = N(x)^j2/ + /o exp(3qt/2).
Требуя тождественного выполнения первого условия (8), окончательно получаем
V-
±.
exp
— qt ± x. —
2 V 2k
M
'3k
, M = const, /o = const.
(12)
Заметим, что знаки перед экспонентой и в показателе степени экспоненты совпадают.
Из условия f2 = 0 следует, что /1 - k/x - q/ = A = const. Нетрудно проверить, что если
/o = 4M2 /(9q2), то подстановка функции /(x, t) из (2) в уравнение (1) дает A = 0 . Итак, имеем
Следствие 1. Решение u(x,t) = (2k2/q2){±^/q/(2k)exp[3qt/2±x^Jq/(2k)]-M/(3k)}2 -/o/2
линейного дифференциального уравнения ut = kuxx + qu + C, где M = const, /o = const,
C = C (q, M ,/o), является поверхностью уровня уравнения (2).
Следствие 2. Решение u(x,t) = (2k2/q2){±^/q/(2k) exp[3qt/2±x^q/(2k)]-M/(3k)}2 -/o/2 линейного дифференциального уравнения (1) является поверхностью уровня уравнения (2), если /0 = 4M 2 /(9q 2).
Положим, что f2 = aQ 2 (/), тогда согласно утверждению 1 система уравнений (7) совместна, если
f'/3 + 2f11(9k) = f2 =aQ2 . (13)
С другой стороны, уравнение (6) в этом случае будет иметь вид kQ + Q f = 0 . Выражая отсюда f и подставляя полученное значение в (13), приходим к зависимости
Q"Q - 5Q'"2 /3 + 3aQ"2Q'2 /k = 0.
Решая это уравнение, получаем частное решение Q = +V2k/а(3/ + /0), где /0 = const. Далее из (13) определяем, что f = 6k/(3/-/0). Требуем, чтобы выполнялись соотношения (8), (9). В результате окончательно получаем, что
V-
3q
exp
с = const.
(14)
Точные решения уравнений теплопроводности
Определяем Q(/) из уравнений (случай f = const)
kQ" + Q'f + Q [ f' / 3 + 2fx2 /(9k)] = 0 ; (kQ" - aQ'3) + Q"^ + Q' [ f' / 3 + 2 fx2 /(9k)] = 0 . Решая линейное уравнение для Q(/), находим, что
1
Q = Cj exp[-f/ /(3k)] + C2 exp[-2 f//(3k)], Cj = const, C2 = const. Учитывая, что решение уравнения (1) u = Qx = Qy/X, получаем ±yj q /(2k) exp[3qt /2 ± x^q /(2k) ]
u=
1 - f exp[3qt/2 ± x^Jq /(2k) ]/(3k) Подставив сюда выражение для / из (11), имеем
{C1 exp[-f1//(3k)] + C2 exp[-2f1//(3k)]}.
u = ±J—exp
V2k
-t ±J±'
2 \2k
1 —f^exp 3k
-t ±xJ^~
2 V 2k
Решаем нелинейное уравнение для Q(/). Полагаем, что Q = p(/), а затем полагаем, что p = r (p). В результате приходим к уравнению krrp + fxr = ap3 - 2 f2 p /(9k). Находим частное решение данного уравнения вида r = ap2 + bp , где a = const, b = const. Получаем, что a = ±yja/(2k), b = f /(3k). Возвращаясь к первоначальному уравнению, имеем
f /(3k)
Q = -
C = const.
±у]а/(2к) - С ехр[/^/(3к)]'
Учитывая, что решение уравнения (2) и = Qx = Q у/х, выпишем окончательно решение уравнения (2)
u = ±[f /(3k)]y[qTaexp[3qt/2±xjq/(2k)] C =±C 1 - C0 - [fx /(3k)]exp[3qt /2 ± xy/q/(2k)], у 2k
(15)
Здесь знак числителя не зависит от знака в показателе экспоненты, то есть фактически имеется четыре решения
= [ f /(3k )]^/ q /a exp[3qt/2 + x^Jq /(2k) ] = -[ f1 /(3k )^ q /a exp[3qt/2 + x^jq /(2k) ]
1 - C0 - [f_/(3k)]exp[3qt/2 ± x^/q/(2k)], 1 - C0 - [fl/(3k)]exp[3qt/2 ± x^Jq/(2k)],
_ [f /(3kХЦ/q/aexp[3qt/2 - x„Jq/(2k)] = -[f\ /(3k)Ц/q/aexp[3qt/2 - x„Jq/(2k)]
1 - C0 -[f /(3k)]exp[3qt/2±x^Jq/(2k)] 1 - C0 -[fx/(3k)]exp[3qt/2± x„Jq/(2k)]
Для случая, когда f2 = 0 уравнение (5) имеет вид Q + 3Q /(2/ + /o) = 0. Решая это уравнение, определяем, что Q = Cj - C2 / у]2/+ /o , Cj = const, C2 = const. Тогда окончательно получаем, учитывая (12), что
2k
u =—Cx exp
— qt ± x. — 2 \2k
(■
±
exp
— qt ± x. — 2 \2k
M
3k
- C2 exp
Л
— qt ± x. — 2 \2k
Если /2 = 0, то уравнение (6) имеет вид ^ -aQ3 /к) + 3Q /(2^+^0) = 0. Выпишем частное решение этого уравнения: Q = ±у]2к/а[1/(2^ + уо)]. Соответствующее решение уравнения (6) будет иметь вид
q 2k u — exp
k\ a
(■
— qt ± x. — 2 \2k
±./— exp 2k
— qt ± x. —
2 V 2k
M_
2k
(17)
Если поверхность уровня определяется выражением (14), решение уравнения (2) имеет вид
u = ±у/q/aexp[3qt/2±xyjq/(2k)]{exp[3qt/2±xy/q/(2k)]-с}-1, с = const. (18)
Здесь знак перед выражением для u может не совпадать со знаком в показателе экспоненты (имеем, фактически, четыре решения).
Другой подход к изучению поведения решений нелинейного уравнения теплопроводности
Покажем, что характеристиками уравнений (1) и (2) являются линии t = const .
Чтобы доказать это утверждение, в уравнении (2) перейдем от функции u(x, t) к функции t(u, x) [8]. Получим
t2 + k(t2txx - 2txtutxu + t2Juu ) - (qu-auЧ3 = 0 . (19)
В уравнении (19) сделаем замену независимых переменных u - ср(x) = ;, x = Т):
t2 + k[t2 (tnn - 2tn9x + t ;9x -1 ;9xx ) - 2t I; (tn - tl;Vx )(t^n - t ;;9x ) + t ; (tn - t%@x ) ]
-[q( ;+9) -a(;+^)3]t| = 0.
Полагая, что ; = const - характеристика уравнения (19), выпишем выражение перед производной t и приравняем его нулю. Получим, что П = 0. Отсюда следует, что на характеристике ; = const имеем t = const. Тогда, чтобы уравнение (19) имело решение, вдоль характеристики должно выполняться условие совместности 12 - ktgp^ - [q( ; + 9) -a(; + 9)3]t| = 0. Решая это оду, в случае, когда a = 0 (уравнение (1)), и заменяя ; + 9 = u, получаем, что вдоль характеристики t = const должно выполняться соотношение
u (t, x) = { w0 (t) + c0 (t) sin[±( x + c1 (t ))yl q / k ]}/ q,
c0 (t) = const, w0 (t) = const, с (t) = const. (20)
Подробнее остановимся на решении уравнения (2). Выпишем вид уравнения, считая, что x(u, t) - независимая переменная
xtxl - kxuu + (qu -au 3) xl =0. (21)
Когда a* 0, требуя выполнения вдоль характеристики t = const условия совместности, приходим к выражению
x = Cj ±! I 4 d2! . (22)
■sjau /(4k) - qu /(4k) + w0u + c0 Здесь, вообще говоря, можно положить с1 = c1 (t), с2 = с2 (t), w0 = w0 (t) и, подставляя полученное выражение в уравнение (21), получить соотношения для определения неизвестных функций с1 = c1 (t), с2 = с2 (t), w0 = w0 (t). Чтобы представить в этом случае характер изменения интересующей нас функции u(x, t), преобразуем выражение под знаком корня в интеграле (22):
au4 /(4k) - qu2 /(4k) + w0u + c0 = [a/(4k)][u2 + 2a(t)u - (q/a- 2a2(t) -b(t))][u2 - 2a(t)u
-(q / a - 2a2 (t) + b(t))]
и таким образом вычисление интеграла (22) приведем к выражению его через эллиптический интеграл первого рода, обратной функцией которого является функция Вейерштрасса [9]. Следовательно вдоль любой характеристики t = const имеем u(x,t) =p(x) (p(x) - функция Вейерштрасса). Известно также [10], что функция Вейерштрасса имеет полюсы, приближение к которым, очевидно, будет приводить к обострению в решении, что не наблюдается в решении (20) линейного уравнения (1).
Заключение
Из полученных решений уравнения (2) следует, что катастрофа в процессах, которые описываются данным уравнением, возникает тогда, когда знаменатель решения (а он присутствует в решениях нелинейного уравнения (15)-(18) и у функции Вейерштрасса) стремится к нулю.
На рис. 1 показано поведение решения нелинейного уравнения в зависимости от времени (1: t = 0; 2: t = 0,5; 3: t = 1; 4: t = 1,5; 5: t = 2; 6: t = 2,5). На рис. 2 имеем вид решения в зависимости от времени для линейного уравнения (1: t = 0; 2: t = 0,1; 3: t = 0,5; 4: t = 1; 5: t = 1,5).
В выражениях (15), (17), (18) сразу виден управляющий параметр, выбор которого позволяет, по крайней мере, отодвинуть катастрофу в случае нелинейного уравнения - это произвольная постоянная в знаменателе решения нелинейного уравнения (см. рис. 3, вид u(x,t) для разных произвольных постоянных с в знаменателе решения при t = 1 (1: с = 100; 2: с = 1000; 3: с = 2500; 4: с = 4000; 5: с = 7000).
Величина произвольной постоянной в знаменателе связана со значением решения в точке {x = 0,t = 0} . Так, в решении (18) с = \±^Jq/a -u(0,0)]/u(0,0). Если u(0,0) = +yjq/a, то с = 0 и
u(x,t) = const = +yjq/a . Катастрофы удастся избежать также, если при стремлении знаменателя выражения к нулю сменить знак в числителе на противоположный (рис. 4, вид решения u(x,t) в случае задания разных произвольных постоянных с в знаменателе решения, 1: с = 1600; 2: с = 2800; 3: с = 4900). Здесь возрастающие ветви решений получены, когда перед выражением (15) знак плюс, а убывающие участки решений получены, когда перед выражением (15) задавался знак минус (см. (16)). При каждом заданном значении u(0,0) увеличение времени процесса ведет к приближению катастрофы (рис. 1, сравните с поведением решения в зависимости от времени в случае линейного уравнения рис. 2).
Рис. 1. Вид решения нелинейного уравнения Рис. 2. Вид решения линейного уравнения
Рис. 3. Вид решения при различных постоянных Рис. 4. Вид решения при смене знака
Рис. 5. Зависимость решения от параметра д Рис. 6. Зависимость решения от параметра а
На рис. 5 показана зависимость решения от параметра q . Увеличение этого параметра приближает катастрофу (1: q = 1; 2: q = 2; 3: q = 3; 4: q = 4; 5: q = 5). Изменение параметра а слабо влияет на решение, но большие его значения заметно отодвигают катастрофу (см. рис. 6; 1: а = 1; 2: а = 40; 3: а = 200; 4: а = 500; 5: а = 1500).
Литература
1. Курдюмов, С.П. Нестационарные структуры, динамический хаос, клеточные автоматы. Новое в синергетике. Загадки мира неравновесных структур / С.П. Курдюмов, Г.Г. Малинецкий, А.Б. Потапов. - М.: Наука, 1996. - 111 с.
2. Режимы с обострением в задачах для квазилинейных параболических уравнений / Самарский А.А., Галактионов В.А. и др. - М.: Наука, 1987. - 477 с.
3. Vazquez, J.L. A Stability Technique for Evolution Partial Differential Equations. A Dynamical System Approach / J.L. Vazquez, V. - Birkhauser Verlag, 2004. - 377 p.
4. Беркович, Л.М. Некоторые аналитические методы нелинейной динамики / Л.М. Беркович // Вестник СамГУ. Естественнонаучная серия. - 2005. - № 2(36). - С. 32-64.
5. Куркина, Е.С. О режимах с обострением в уравнениях ut = div(uCTgradu) + ue / E.C. Куркина, ИМ. Никольский // Дифференциальные уравнения. Функциональные пространства. Теория приближений: Тезисы докладов международной конференции, посвященной 100-летию со дня рождения Сергея Львовича Соболева. - Новосибирск, 2008. - С. 512.
6. Рубина, Л.И. Один геометрический метод решения нелинейных уравнений в частных производных / Л.И. Рубина, О.Н. Ульянов // Труды Института математики и механики УрО РАН. Екатеринбург: ИММ УрО РАН, 2010. - Т. 16, № 2. - C. 209-225.
7. Рубина, Л.И. О решении уравнения потенциала / Л.И. Рубина, О.Н. Ульянов // Труды Института математики и механики УрО РАН. Екатеринбург: ИММ УрО РАН, 2008. - Т.14, № 1. -C.130-145.
8. Рубина, Л.И. О характеристиках и решениях одномерного нестационарного уравнения фильтрации / Л.И. Рубина // ПММ. - 2005. - Т. 69. - Вып. 5. - с. 829-836.
9. Ломкаци, Ц.Д. Таблицы эллиптической функции Вейерштрасса. Теоретическая часть / Ц.Д. Ломкаци; под ред. В.М. Белякова, К.А. Карпова. - М.: ВЦ АН СССР, 1967. - 88 с.
10. Градштейн, И.С. Таблицы интегралов, сумм, рядов и произведений / И.С. Градштейн, ИМ. Рыжик. - М.: Физматлит, 1962. - 1100 с.
TOWARDS THE DIFFERENCES IN BEHAVIOUR OF SOLUTIONS OF LINEAR AND NON-LINEAR HEAT-CONDUCTION EQUATIONS
L.I. Rubina\ O.N. Ul’yanoV
The linear and non-linear heat-conduction equations are analyzed by the previously initiated geometrical method of analyzing linear and non-linear equations in partial derivatives. The reason of the difference in behavior of solutions of equations under consideration was stated, as well as the reason of aggravation of the non-linear equation. A class of solutions of linear equations that represents the surfaces of the levels of non-linear heat-conduction equations was excluded.
Keywords: non-linear equations in partial derivatives, heat-conduction equations, exact solutions, surfaces of the level.
References
1. Kurdyumov S.P., Malinetskiy G.G., Potapov A.B. Nestatsionarnye struktury, dinamicheskiy khaos, kletochnye avtomaty. Novoe v sinergetike. Zagadki mira neravnovesnykh struktur (Unsteady structures, dynamic chaos, cellular automata. New in synergetics. Mysteries of the world of nonequilibrium structure). Moscow: Nauka, 1996. 111 p. (in Russ.).
2. Samarskiy A.A., Galaktionov V.A., Kurdyumov S.P., Mikhaylov A.P. Rezhimy s obostreniem v zadachakh dlya kvazilineynykh parabolicheskikh uravneniy (Blow-up regimes in problems for quasilin-ear parabolic equations). Moscow: Nauka, 1987. 477 p. (in Russ.).
3. Vazquez J.L., Galaktionov V. A Stability Technique for Evolution Partial Differential Equations. A Dynamical System Approach. Birkhauser Verlag, 2004. 377 p. (ISBN: 0-8176-4146-7)
4. Berkovich L.M. Vestnik SamGU - Estestvennonauchnaya seriya. 2005. no. 2(36). pp. 32-64. (in Russ).
5. Kurkina E.S., Nikol'skiy I.M. O rezhimakh s obostreniem v uravneniyakh ut = div(uCT gradu) + ue (About blow-up regimes in equations ut = div(uCT gradu) + ue). Different-sial'nye uravneniya. Funktsional'nye prostranstva. Teoriya priblizheniy: Tezisy dokladov mezhdunarod-noy konferentsii, posvyashchennoy 100-letiyu so dnya rozhdeniya Sergeya L'vovicha Soboleva. (Abstracts of the International Conference dedicated to the 100th anniversary of the birth of Sobolev “Differential Equations. Functional Space. Theory of Approximation”) Novosibirsk, 2008. p. 512.
6. Rubina L.I., Ul’ianov O.N. Trudy Inst. Mat. i Mekh. UrO RAN. 2010. Vol. 16, no. 2. pp. 209225. (in Russ.).
7. Rubina L.I., Ul’yanov O.N. Proceedings of the Steklov Institute of Mathematics (Supplementary issues). 2008, Vol. 261. Suppl. 1. pp. 183-200.
8. Rubina L.I. Journal of Applied Mathematics and Mechanics. 2005. Vol. 69. Issue 5. pp. 829-836. (in Russ.).
9. Lomkatsi Ts.D. Tablitsy ellipticheskoy funktsii Veyershtrassa. Teoreticheskaya chast' (Weier-strass elliptic function charts. Theoretical part). Moscow: VTs AN SSSR, 1967. - 88 p. (in Russ.).
10. Gradshteyn I.S., Ryzhik I.M. Tablitsy integralov, summ, ryadov i proizvedeniy (Table of integrals, sums, series and compositions). Moscow : Fizmatlit, 1962. 1100 p.
Поступила в редакцию 17 сентября 2012 г.
1 Rubina Liudmila Ilinichna is Cand. Sc. (Physics and Mathematics), Senior Staff Scientist, Institute of Mathematics and Mechanics of the Russian Academy of Sciences (Ural branch).
E-mail: [email protected]
2 Ul’yanov Oleg Nikolaevich is Cand. Sc. (Physics and Mathematics), Senior Staff Scientist, University’s academic secretary, Institute of Mathematics and Mechanics of the Russian Academy of Sciences (Ural branch).
E-mail: [email protected]
