УДК 591.65
ОБ УТОЧНЕНИИ ОДНОЙ АСИМПТОТИЧЕСКОЙ ФОРМУЛЫ ДЛЯ КОНСТАНТЫ ЛЕБЕГА
И. А. Шакиров
Кандидат физико-математических наук, проректор по дополнительному образованию, Набережночелнинский институт социально-педагогических технологий и ресурсов, [email protected]
Для константы Лебега классического полинома Лагранжа, определенного в четном числе узлов интерполяции, получена строгая двусторонняя оценка. На этой основе конкретизирована неопределенная величина 0( 1) в известном асимптотическом равенстве для константы Лебега. Решены две актуальные задачи теории интерполирования, связанные с оптимальным выбором 0(1).
Ключевые слова: интерполяционный полином Лагранжа, верхняя и нижняя оценки константы Лебега, асимптотическое равенство, погрешность интерполяции.
ВВЕДЕНИЕ
Хорошо известны теоретическое и прикладное значения тригонометрических полиномов Лагранжа. Достаточно полные сведения о них, а также об их интерполяционных характеристиках можно найти в монографиях В. Л. Гончарова [1], И. П. Натансона [2], С. Б. Стечкина, Ю. Н. Субботина [3], Н. П. Корнейчука [4], В. К. Дзядыка [5], А. А. Привалова [6], посвященных различным вопросам теории функций.
Для приближенного представления функции х = х(£) £ С2п = С[0,2п] в приложениях часто используют интерполяционный полином Лагранжа:
1 2п~1 1 2п , . ч
ФП(х,0 = - £ х«к те & - о = - £х((* те «к - о [о-п (и) = . 0)
к=0 к = 1
определенный в четном числе N = 2n равномерно распределенных узлов вида tk = nk/n (k = 0, 2n — 1 V k = 1, 2n, n £ N). При этом норма отклонения полинома (1) от функции x(t) оценивается согласно фундаментальному неравенству:
llx — ФПx||с[0,2п] < (1 + а;)e;(x) (ФП : CW ^ Н с CW), (2)
где Е;(x) — наилучшее равномерное приближение функции x(t) тригонометрическими полиномами из множества
í n-1 1
НП = \ a0 + ^ (ak cos kt + bk sin kt) + an cos nt; ak , bk £ M f , k=i
1 2к — 1
АП = (п £ М) (3)
к=1
— константа Лебега полинома (1), АП = ||ФП||с2п(см., например, [5, с. 43] и [7]).
Согласно неравенству (2) поведение последовательности величин АП = А* (п) (п £ М) имеет первостепенное значение при исследовании вопросов сходимости полиномов (1) к интерполируемой функции х(¿). Первоначально константа Лебега оценивалась сверху как АП ^ 0(1п п) либо АП ^ А+В 1п п, например, АП < 4+^п [1], АП < 8+(4/п) 1пп [2]. Затем коэффициенты А, В многократно уточнялись.
Было установлено асимптотическое равенство:
2
АП = — 1п п + 0(1) (В = 2/п, п ^то) (4)
п
с неопределенной (но ограниченной) величиной 0(1). Используя явный вид константы Лебега (3), в работе [8, с. 106] уточнено поведение асимптотического равенства (4) при произвольно выбранных значениях п, т.е. получено неравенство вида АП < 4/п + (2/п) 1п(4п/п) (п £ М). В [9, с. 215] для неопределенной константы приведено соотношение 0(1) = 2(1 — 1/п)6>П (6>П £ (0,1)), согласно которому 0(1) ^ 0 при приближении параметра 9П к левой границе интервала (0,1), что указывает
на не строгое оценивание величины O(1) снизу. В случае интерполирования по узлам Чебышева в работах [10, 11] для неопределенной константы O(1) содержится более точное выражение вида O(1) = (2/n)[Y + ln(1/n)] + o(1), где y = 0.577215 ... — константа Эйлера.
Отметим, что в большинстве работ, посвященных данной проблематике, нижние оценки для неопределенной величины O(1) и константы Лебега в основном отсутствуют либо являются грубыми.
Из вышесказанного следует, что исследование приближенного представления константы Лебега является одним из важнейших задач теории интерполирования функции. Она до сих пор окончательно не решена, т. е. в равенстве (4) не определено конкретное значение константы O(1), которое сохраняло бы хорошую близость его левой и правой частей при произвольно взятых натуральных значениях параметра n. В связи с этим в данной работе ставятся и решаются две актуальные задачи, связанные с указанными проблемами.
Задача 1. Определить наименьшее значение A* неопределенной константы O(1) в (4) так, чтобы для всех натуральных значений параметра n было выполнено неравенство АП < < Л* + (2/п) ln n (n £ N), и оценить допущенную при этом погрешность
е* = sup{A* + (2/п) ln n - АП | n £ N}. (5)
Задача 2. В случае положительного решения задачи 1 найти другую (отличную от Л*) константу Л2 так, чтобы величина абсолютной равномерной погрешности
Г = sup{|A* + (2/п) lnn - АП|| n £ N} (6)
была бы меньше ранее найденного значения погрешности (5).
Ясно, что первая из них является экстремальной задачей, связанной с определением наилучшей равномерной верхней оценки для константы Лебега. Вторая задача по своей постановке не является экстремальной. Она связана с проблемой уменьшения величины погрешности (6) относительно (5), используя при этом свободу выбора неопределенной константы O(1).
1. ВСПОМОГАТЕЛЬНЫЕ РЕЗУЛЬТАТЫ
Введем классы дискретных функций V+ и V-, которые будут существенно использованы при доказательстве необходимых лемм, а также при решении поставленных во введении задач.
Определение 1. Строго монотонная функция ^ = <^(n) (n £ D) дискретного аргумента, имеющая малое изменение 5 области значений R(^), называется функцией с малой монотонной вариацией; множество таких функций обозначим через символ V±, где D = D(^) С N; 5 = 5(<р) = sup{^(n) | n £ D} - inf{^(n) | n £ D} (5 £ (0,0.2)) — вариация функции в рассматриваемой области; знак «+» для класса возрастающих в области D функций, знак «—» для класса убывающих функций.
Заметим, что в процессе доказательства лемм и некоторых теорем будем использовать непрерывные продолжения функций ^ = <^(n) (n £ D) из классов V± на область D = D(^) = = (inf D(^),supD(^)) с R, состоящую из предельных точек (не дискретную область), для которых суть и формулировка определения 1, обозначения классов V+ и V— полностью сохраняются.
Лемма 1. Функция дискретного аргумента
ап = a(n) = (1/n) cos ec(n/2n) (a : D ^ R, D = N \ {1} = {n| n > 2}) (7)
принадлежит классу V -, а для ее образа и вариации верны соотношения:
R(an) = (2/п, V2/2] с (0.636, 0.708), 5 = 5(an) < 0.071, V n > 2.
Доказательство. Функция an = a(n) (n ^ 2) является дискретным аналогом равномерно убывающей в области D = [2, непрерывно дифференцируемой функции y = (1/x)cosec(n/2x). Для ее образа и вариации верны соотношения: R(y) = (2/п, ^2/2] с (0.636,0.708) и 5 = 5(y) = = л/2/2 - 2/п < 0.071. Следовательно, an £ V—.
Аналогичный результат получим, если функцию (7) представим в несколько другом виде an =
2 /sinn/2n\-1 , „ ^
= — -----(n > 2) и при этом используем свойства первого замечательного предела. □
п \ n/2n )
Лемма 2. Следующие далее зависимости ^ = ^ (п) (п £ Б), выраженные через функцию (7),
принадлежат классам V±:
^i(n) = ап • (1 + cos(—/2n)) е V+, ) = [1/2 + ^2/2, 4/—) С (1.207,1.274);
2
(n) = - ln ап е V-, R(^2) =
—
2 —'
- п ln 2, - - ln ^ ) С (-0.288, -0.220);
^s(n) = ((9 - V2)/8)^i е V+, R(^s) = [(7 + 8л/2)/16, (9 - v^)/2—) С (1.144,1.208); <^(n) = ((л/2 - 1)/4)^i е V+, ) = [1/8, (л/2 - 1)/—) С [0.125,0.132).
Доказательство. С целью исследования поведения функции = (n) (n ^ 2) при помощи производной непрерывно продолжим ее на область Б, где она к тому же является гладкой функцией, т. е. <pi (n) е C'[2,
Покажем, что производная функции (n) положительна в области Б, определим ее образ и оценим вариацию 5 = 5):
(n) = (Y = (1 (-/4n)A' = -(1 - (W-) sin(—/2n)) =
V n sin(-/2n) y Vn V ) 4n3 sin2 (—/4n)
— / A2
= 4n (an)
1
sin(-/2n)
—/2n
> 0, V n е Б;
Я(<^1) = [1/2 + ^2/2,4/п), 5 = 5(^1) < 0.067. Следовательно, £ .
Согласно лемме 1 функция ап = а(п) (п ^ 2) принадлежит классу V—, и на основании известного свойства логарифма имеем (п) = (2/п)1п ап £ V—.
Остальные утверждения леммы 2 очевидны. □
Лемма 3. В области Б для функций (п) = 1—<^(п) —(п) £ и <^6(п) = (п)+^5(п) £ V/ верны следующие соотношения:
Я(^5 ) =
Я(^6) =
1 f* + 1 ln 2,1 - 4 + 2ln — ) С (0.013,0.015), 2 — — — 2 /
Ц^2 + 1ln2,1 - + 2ln - С (0.138,0.147).
8 — — — 2 /
(8)
Доказательство. Изучим поведение производной функции <^5 = <^5 (п) в расширенной области
Б = [2,
, = / 2 1 \ ' /1 + ес8(п/2п)\' = 2вт(п/2п) — (п/п) ео8(п/2п) п ( ^5) = ^П 1П ап) ^ п в1п(п/2п) = пп в1п(п/2п) 2^ +
+
sin(—/2n) + sin(—/2n) cos(—/2n) - (—/2n) cos(—/2n) - (—/2n) cos2(—/2n)
n2 sin2(—/2n)
1
[4n2 sin2(—/2n) + 2—n sin(—/2n) - —2(1 + cos(—/2n)] =
2—n3 sin2 (—/2n)
= 2n )2Í
sin(—/2n)^ sin(—/2n)
+--—--1
—/2n
—/2n
- cos(—/2n) > .
Разложение Тейлора (1/90) • (п/2п)4 + о((п/2п)4) функции в фигурных скобках имеет только положительные значения для всех п £ Б, поэтому производная исследуемой функции всюду положительна. Следовательно, = (п) монотонно возрастает в рассматриваемой области,
и п2
R( ч> 5) =
^ + 1 ln2,1 - 4 + 2ln —, 5 = 5Ы = ^ - 4 , -2 — ' — — 2 2 — — 8
- z + zln^T < 0.001, ^ (n) е V+.
Функция ^ 6 = ^ 6 (п)(п £ Б) принадлежит классу как линейная комбинация функций из этого
же класса.
□
Замечание 1. Леммы 1-3 доказаны в области I) = [2, с использованием элементов дифференциального исчисления. Следовательно, все они справедливы в ее дискретной подобласти ) = {п| п ^ 2}, чего достаточно для доказательства теорем следующего пункта.
2. ОСНОВНЫЕ РЕЗУЛЬТАТЫ
В этом параграфе докажем три теоремы, которые позволят решить поставленные во введении задачи.
Теорема 1. Для константы Лебега (3) справедлива двусторонняя оценка:
5 — v/— 2 2
---<£1(п) + <р2(п) + 1п п < АП < <Мп) + ^2(п) + 1п п, V п > 2,
4 п п
где (п), <^2 (п) — введенные в лемме 2 функции.
Доказательство. Оценим константу АП(п > 2) сверху и снизу:
(9)
а; = -£ ctg
tk — n/2n 1
k=i
2
п
n ii — n/2n
-ctg-~-
n 2
+ V - ctgtk - n/2n
k=2
n
2
cos(n/4n) 2 fn/2 7 2 cos2(n/4n) 2, + - / ctg tdt =-. ; L / + - ln
1 n 1 Г t 1 < -ct^ — ^ - / ctg - dt = n 4n n J t, 2
n
n sin(n/4n) n Jn/2n
n sin(n/2n) n n sin(n/2n)
22
= an (1 + cos(n/2n)) +— ln an +— ln n;
П П
л* 13 п n АП = - i 7 ' " ctg — + n П 4 n 4n
1 п п n tk — n/2n
4 * nctg 4n + £ n
k=2
>
> 1 Г3n t n п I 4n 4n
t
п П Л 1 , n 3n\ 1 Г 4nctg 4n' ■ 2 + (4nctg4^J ■ 2+1 ctg2 dt
n/2
3 П 1 П 1 П 4 cos2 (n / 4n) — 3 2 " 4nctg 4n + ctg 4n + ctg 4n ' 3 — 4 sin2 (n/4n) + П Jn/2n
ctg tdt >
7 3 — 2^2 A 1
+
- n 2 2 5 — V2
-ctg ---+ - ln an + - ln n = ---an
n 4n п n 4
88
где при проведении нижней оценки дополнительно использовано неравенство
4С°82(п/24п) — 3 > 3 — 2^2 (п > 2). 3 — 4ип2(п/4п) 1 У
Итак, для константы Лебега имеет место двусторонняя оценка: 5 —
Л п \ 2 2,
1 + cos — +— ln an +— ln n, v 2n/ п n
-an
(^1 + cos —^ +— lnan +— lnn < АП < an (1 + cos —^ +— lnan +— lnn (n ^ 2), V 2n / п n П V 2n/ п n
□
4 v 2n/ п n v 2n/ п n
равносильная неравенству (9) (см. обозначения функций в лемме 2).
Замечание 2. В формуле (3) при значении параметра n =1 для константы Лебега верно равенство А* = 1, используя которое двойное неравенство (9) перепишется в виде (5 — \/2)/4 < А1 < 1 (n = 1). Другими словами, (9) останется непротиворечивым и при n = 1, если в нем строгое верхнее неравенство заменить на нестрогое:
((5 — V2)/4)^i (n) + «(n) + (2/п) ln n < АП < «i(n) + «^(n) + (2/п) ln n, V n £ N.
С учетом сказанного рассмотренный тривиальный случай включим в содержание предыдущих теорем для полноты их формулировок.
Полученные в теореме 1 результаты позволяют константу А; (n > 2) приближенно заменить полусуммой ее верхней и нижней оценок, а затем и функциональной зависимостью дП = M*(n) = = 1 + (2/п) ln n:
5 —
I a„ I I + cos !
2n/ n
А*
2 ■+
4
П \ 2 2
an ( 1 + cos — ) +— ln an +— ln n =
9 — . . . , 2,
---«1 (n) + «2 (n) + - ln n =
8 n
22
= ^2 (n) + ¥з (n) + — ln n « 1 + ln n = мП (n > 2). (10)
——
Если в (10) положим n = 1, то имеем непротиворечивые приближенные равенства:
Ai « ¥2(1) + ¥з(1) + 2 ln 1 « М* 1 « « ^ ,
причем А* = м* .
Уточним поведение приближенной замены АП ~ мП (n е N) и оценим допущенную при этом погрешность.
Теорема 2. Для константы Лебега (3) равномерно относительно значений аргумента n справедливо неравенство:
АП < 1 + (2/—) lnn = мП (n е N), (11)
в котором равенство достигается только при n = 1; для погрешности еП = МП - АП верна равномерная оценка
е* = sup{< | n е N} = 1 - (5 - V2)/— + (2/—) ln(—/2) < 0.147. (12)
Доказательство. Если в неравенстве (11) положим n = 1, то его левая и правая части совпадают, т. е. Ai = м* = 1, следовательно, е* = 0.
Пусть n > 2. Тогда с целью двусторонней оценки функции погрешности еП = е* (n) (n > 2) в
2
двойном неравенстве (9) вычтем всюду мП = 1 +—ln n и после некоторых преобразований получим
—
(см. леммы 2 и 3):
((5 - У2)/4)¥i(n) + ¥2(n) - 1 < АП - мП < ¥i(n)+ ¥2(n) - 1 ^ ^ 1 - ¥i (n) - ¥2(n) < мП - АП < 1 - ¥2(n) - ((5 - V^)/4)¥i(n) ^
^ 1 — ^1(п) — ^2(п) < < < 1 — (п) — ^2(п) + ((л/2 — 1)/4)<Р1 (п) ^
^ ^5(п) << < ^6(п) ( 5(п), ^6(п) £ V/, п > 2). (13)
Используя теперь содержащиеся в лемме 3 сведения о функциях (8), преобразуем неравенство (13):
inf{¥5(n)| n е I} < ¥5(n) < еП < ¥6(n) < sup{¥e(n)| n е I}, V n > 2 ^
^ 1 -^ + 1 ln2 <еП < 1 - 5 - ^ + 2ln — ^ 0.013 <мП - АП < 0.147, V n > 2. 2 — — — 2
Учитывая результаты тривиального случая (n =1 ^ е* = м* - А* =0), уточним последнее неравенство:
0 < мП - АП < 0.147, V n е N ^ 0 < sup^| n е N} < 0.147 ^ 0 < е* < 0.147. Из этих соотношений легко получим справедливость оценок (11) и (12). □
Следствие. Единственным решением экстремальной задачи
AnfR{A* | АП < A* + (2/—) ln n, V n е N} (14)
является константа A* = Ai = 1.
Доказательство. Если в правой части неравенства (11) единицу заменим другим вполне определенным числом A, больше единицы, то согласно теореме 2 в (11) соответственно появится строгое неравенство, которое будет выполнено при любых натуральных значениях n, и при этом значение соответствующей погрешности увеличится по сравнению с (12). Следовательно, A (A > 1) никак не может быть решением задачи (14). Полагая обратное (A есть решение (14)), сразу же обнаружим, что решением рассматриваемой экстремальной задачи является константа A' (1 < A' < A). Вывод: решение задачи (14) не может быть больше единицы.
Если же константу А выбрать таким, что А < 1, то тогда неравенство (11) не будет выполнено для всех натуральных значений параметра п. Например, при значении параметра п = 1 всегда имеем А* < А+ (2/п)1п1(^ 1 < А). Полученное противоречие вида 1 < А < 1 показывает, что А не может быть меньше единицы. Следовательно, единственным решением задачи 1 является константа
А* = А* = 1. □
Теорема 3. Для константы Лебега (3) имеет место приближенное равенство вида
АП « 3/п + (2/п)1пп = ДП (п £ М),
в котором для абсолютной погрешности еП = |дП — АП I равномерно относительно значений аргумента п верна оценка
е* = Бир{еП| п £ М} = Бир{|3/п + (2/п) 1пп — АП| | п £ М} < 0.102.
Доказательство. На этот раз для константы Лебега (3) согласно ранее использованной схеме (10) получим следующие приближенные равенства:
АП ~ Р2 (п) + Рз (п) + (2/п) 1п п « 3/п + (2/п) 1п п = ДП (п £ М).
(15)
Выбранная в приближенной замене (15) константа А* = 3/п (А* < А* ^ 3/п < 1) удовлетворяет основному условию существования решения задачи 2: А2 £ [1,1 — е*], где е* = 1 — (5 — \/2)/п + (2/п) 1п(п/2) « 0.147 — значение погрешности, определенное в предыдущей теореме 2.
Теперь оценим допущенную в приближенных заменах (15) погрешность, используя результаты теорем 1 и 2. Для этого в двойном неравенстве (9) всюду вычтем ДП = 3/п + (2/п) 1п п (п > 2), затем, используя обозначения лемм 2 и 3, последовательно упростим полученные соотношения:
((5 — У2)/4)^1 (п) + р2(п) — 3/п < АП — дП < Р1 (п) + Р2(п) — 3/п ^ ^ 1 — р1 (п) — р2 (п) — (1 — 3/п) < дП — АП < < 1 — (п) — Р2 (п) + ((^2 — 1)/4)р1 (п) — (1 — 3/п) ^ ^ Р5 (п) — (1 — 3/п) <еП <Рб (п) — (1 — 3/п) ^ (п) <еП < рв (п)
(р7 (п) = Р5 (п) — (1 — 3/п), рв (п) = Р6 (п) — (1 — 3/п), п > 2). (16)
Ясно, что вновь введенные функции р7(п), р8(п) принадлежат области , так как согласно лемме 3 таковыми являются функции р5(п), р6 (п) (п ^ 2). Для их образов и вариаций верны соотношения:
Я(Р7) = Я(Р8) =
3 — 2 1п п — I) с (—0.032, —0.030), 5 = 5(р7) < 0.001;
п 2 п п 2 п /
3 — 3+о4л/2 + 11п2, 2 1п п — с (0.093, 0.102), 5 = 5(рв) < 0.008.
п 8 п п 2 п
Используя полученные сведения о функциях р7(п), р8(п) (п > 2), преобразуем неравенство (16):
т£ р7(п) < р7(п) < еП < Р8 (п) ^ яир рв(п), V п > 2.
П^2 П>2
Отсюда
|еП < тах{
т£ р7(п) БИр рв (п) > = тах <
П>2 п>2 ' 1
^ + 11п 2
2п
2 п 2 — л/2
— 1п —--
п 2 п
2 п 2-^/2
= _ 1п---< 0.102, V п > 2.
п 2 п
Учитывая результаты тривиального случая (п = 1 ^ е* =0), уточним последнее неравенство:
|еП I < 0.102, V п £ М Теорема доказана.
•Ф»
яир еП < 0.102
< 0.102.
□
—*
Замечание 3. Выбранная в условиях теоремы 3 константа А 2 = 3/п действительно является решением задачи 2, так как она обеспечивает требуемое уменьшение погрешности (е* < е*). Однако из результатов теорем 2 и 3 следует, что внутри интервала (1,1 — е*) существует «лучшая» константа А*, которая позволяет уменьшить теоретически полученное значение погрешности е* (см. (12)) в два раза, т. е. до значения е*/2 = 1/2 + (1/п) 1п(п/2) — (5 — \/2)/2п « 0.073.
Библиографический список
1. Гончаров В. Л. Теория интерполирования и приближения функций. М. ; Л. : ГТТИ, 1934.
2. Натансон И. П. Конструктивная теория функций. М. ; Л. : Гостехиздат, 1949.
3. Стечкин С. Б., Субботин Ю. Н. Сплайны в вычислительной математике. М. : Наука, 1976.
4. Корнейчук Н. П. Точные константы в теории приближения. М. : Наука, 1987.
5. Дзядык В. К. Аппроксимационные методы решения дифференциальных и интегральных уравнений. Киев : Наук. думка, 1988.
6. Привалов А. А. Теория интерполирования функций : в 2 кн. Саратов : Изд-во Сарат. ун-та, 1990.
7. Шакиров И. А. О фундаментальных характеристиках семейства интерполяционных полиномов Лагранжа // Изв. Сарат. ун-та. Нов. сер. Сер. Ма-
тематика. Механика. Информатика. 2013. Т. 13, вып. 1, ч. 2. С. 99-104.
8. Габдулхаев Б. Г. Оптимальные аппроксимации решений линейных задач. Казань : Изд-во Казан. ун-та, 1980.
9. Бабенко К. И. Основы численного анализа. М. ; Ижевск : НИЦ Регулярная и хаотическая динамика, 2002.
10. Brutman L. Lebesgue functions for polynomial interpolation - a survey // Ann. Numer. Math. 1997. Vol. 4. P. 111-127.
11. Vertesi P. On the Lebesgue function and Lebesgue constant : a tribute to Paul Erdos // Bolyai Society of Mathematical Studies. 2002. Vol. 11. P. 705-728.
On a Refinement of the Asymptotic Formula for the Lebesgue Constants
I. A. Shakirov
Naberezhnye Chelny Institute of Social Pedagogical Technologies and Resources, 28, Nizametdinov st., 423806, Naberezhniye Chelny, Tatarstan, Russia, [email protected]
For the Lebesque constant of the classical Lagrange polynomial defined in the even number of nodes of interpolation, strict two-sided estimation is received. On this basis, an undefined value O( 1) is refined in the well-known asymptotic equality for the Lebesque constant. Two actual problems in the interpolation theory associated with the optimal choice of O(1) are solved.
Key words: Lagrange interpolation polynomial, upper and lower assessment of the Lebesque constant, asymptotic equality, error of interpolation.
References
1. Goncharov V. L. Teoriia interpolirovaniia i priblizheniia funktsii [Interpolation theory and approximations of functions]. Moscow; Leningrad, GTTI, 1934 (in Russian).
2. Natanson I. P. Constructive function theory. Vol. 1-3. New York, F. Ungar Publ. Co., 19641965.
3. Stechkin S. B., Subbotin Yu. N. Splainy v vychislitel'noi matematike [Splins in computational mathematics]. Moscow, Nauka, 1976 (in Russian).
4. Korneichuk N. P. Tochnye konstanty v teorii priblizheniia [Constants in approximation theory]. Moscow, Nauka, 1987 (in Russian).
5. Dzyadyk V. K. Approksimatsionnye metody resheniia differentsial'nykh i integral'nykh uravnenii [Approximation Methods for solving Differential and Integral Equations]. Kiev, Naukova Dumka, 1988 (in Russian).
6. Privalov A. A. Teoriia interpolirovaniia funktsii [Interpolation of functions theory]. Saratov, Saratov Univ. Press, 1990 (in Russian).
7. Shakirov I. A. About the fundamental characteristics of the lagrange interpolation polynomials family. Izv. Saratov Univ. (N. S.), Ser. Math. Mech. Inform., 2013, vol. 13, iss. 1, pt. 2, pp. 99104 (in Russian).
8. Gabdulhaev B. G. Optimal'nye approksimatsii reshenii lineinykh zadach [Optimal approximations of linear problems solutions]. Kazan, Kazan Univer. Press, 1980 (in Russian).
9. Babenko K. I. Osnovy chislennogo analiza [Fundamentals of numerical analysis]. Moscow, Izhevsk, NIC Regular and chaotic dynamics, 2002 (in Russian).
10. Brutman L. Lebesgue functions for polynomial interpolation - a survey. Ann. Numer. Math., 1997, vol. 4, pp. 111-127.
11. Vertesi P. On the Lebesgue function and Lebesgue constant : a tribute to Paul Erdos. Bolyai Society of Math. Studies, 2002, vol. 11, pp. 705-728.