1-23 RAR 2491-2244 0000-0001-7324-1823 Варин Виктор Петрович Институт прикладной математики им. М.В.Келдыша РАН varin@keldysh.ru Персональная веб-страница: http://keldysh.ru/persons/258.html Varin V.P. Keldysh Institute of Applied Mathematics RAS Функциональное суммирование рядов Functional summation of series. Рассматривается способ суммирования рядов, который сводится к решению некоторых линейных функциональных уравнений. Частичные суммы любого числового ряда удовлетворяют очевидному разностному уравнению. Это уравнение преобразуется в функциональное уравнение на интервале [0,1] для непрерывного аргумента. Далее это уравнение либо решается явно (с точностью до произвольной константы), либо вычисляется асимптотическое разложение решения в нуле. Сумма исходного ряда определяется однозначно как константа, которая необходима для согласования асимптотического разложения решения с частичными суммами исходного ряда. Понятие предела не участвует в данной вычислительной схеме, что позволяет суммировать также расходящиеся ряды. We consider a summation technique which reduces summation of a series to the solution of some linear functional equations. Partial sums of a series satisfy an obvious difference equation. This equation is transformed to the functional equation on the interval [0,1] for the continuous argument. Then this equation is either solved explicitly (to within an arbitrary constant) or an asymptotic expansion of the solution is computed at the origin. The sum of the original series is determined uniquely as a constant needed for the matching of the asymptotic series with partial sums of the original series. The notion of a limit is not involved in this computational technique, which allows summation of divergent series as well. 1. ВВЕДЕНИЕ Суммирование рядов является одной из старейших задач математического анализа и берет свое начало в античности. Эта задача всегда имела два аспекта: практический, т.е. вычислительный, и теоретический, т.е. концептуальный. Иными словами, прежде чем вычислять сумму бесконечного ряда, следует условиться, что следует понимать под этой суммой. Во времена Эйлера понятие суммы ряда не было четко определено, что давало определенный простор в толковании этого понятия и позволило просуммировать многие расходящиеся в обычном понимании ряды. С появлением строгого понятия предела многие математики отказались от использования расходящихся рядов. Достаточно вспомнить знаменитое высказывание Абеля, который усмотрел в расходящихся рядах происки дьявола. Однако довольно скоро выяснилось, что без расходящихся рядов никак не обойтись, что привело к созданию теории расходящихся рядов и появлению методов их суммирования (см. [1]). В настоящее время существует обширный арсенал методов суммирования рядов, многие из которых могут суммировать также и расходящиеся ряды. Однако все эти методы (за исключением формулы Эйлера Миклорени) сводят так или иначе исходно медленно сходящийся или расходящийся процесс к сходящемуся, предел которого и называется суммой ряда. Формула Эйлера Миклорени. а также значительно менее известная формула Эилери Буле здесь стоят особняком, так как они сводят (медленно) сходящийся ряд к расходящемуся асимптотическому. Т.е. эти формулы являются скорее (весьма эффективными) способами вычисления суммы исходно сходящегося ряда. Таким образом, все существующие способы суммирования рядов, как сходящихся, так и расходящихся, используют понятие предела и определяют сумму ряда как предел некоторого сходящегося процесса, который и является алгоритмом суммирования. В то же время сумму ряда можно определить так, что понятие предела нигде не будет использовано, что дает возможность суммировать ряды, суммы которых ранее не были определены и для которых существующие методы суммирования неприменимы. Приведем пример такого ряда. Рассмотрим формальный ряд и определим его сумму при x 1. При х = —1 сумма этого ряда равна константе Эйлера-Гомпертца 6 = eEi(1,1), где Ei() - это интегральная экспонента (см. разд. 5). Если бы этот ряд сходился при х< 1 a lim^i у(х) существовал, то это был бы классический способ суммирования Абеля. Однако этот ряд всюду расходится. Причем ни один из существующих способов суммирования рядов, как то е—алгоритм Винни. Паде-аппроксимации, метод Левина и его модификации, итерационный метод Айткена (см. обзор в [2] и ссылки там), а также т—метод Ланцоша (см. [3]) и т.д., не могут вычислить его сумму. у( х) ренциальному уравнению х2 у' — (1 — х) у +1 = 0, которое интегрируемо. Функция где C - это константа интегрирования, является общим решением этого урав- у(0) = 1 имеют это начальное значение и одно и то же асимптотическое разложение в нуле (так как exp(—1/х) является плоской функцией). у( х) функции Ei() при х ? [0,1] расположен в точности на разрезе комплексной плоскости (—оо, 0), где функция Ei() претерпевает скачок. Константу C ным образом. Возьмем эту константу чисто мнимой и так, чтобы функция у(х) принимала вещественные значения при х ? R+, т.е. возьмем в качестве решения функцию ,W = _Iexp(_I) (ijr + Ei(l,-i)). Тогда в качестве суммы исходного ряда следует принять у(1) = — exp(—1) (in + Ei (1, — 1)) « 0.69717488323505. Заметим, что интегральный метод суммирования Бореля здесь не применим напрямую. Но если его модифицировать и вычислять интеграл в смысле главного значения, то найденная нами константа правильно вычисляется. Данный метод суммирования, т.е. замена числового ряда на степенной, предъявление дифференциального уравнения для его суммы с последующим интегрированием ОДУ (аналитически либо численно) - является весьма мощным средством суммирования, но имеющим весьма ограниченное применение. Нужное дифференциальное уравнение в общем случае может и не существовать. Но этот пример показывает, что сумму ряда можно, в принципе, определить, не прибегая к понятию предела. 2. СУММА РЯДА КАК КОНСТАНТА ИНТЕГРИРОВАНИЯ По аналогии с ОДУ решение разностных уравнений мы будем называть их интегрированием, а общее решение разностного уравнения будем называть первым интегралом разностного уравнения. Эта аналогия хорошо заметна при построении непрерывных итераций дискретных аналитических отображений (см. [4]), где предъявление первого интеграла отображения эквивалентно построению непрерывных итераций и инвариантных кривых этого отображения (см. [5]). При суммировании числового ряда его частичные суммы n s(n)= ? a(k) (2) удовлетворяют очевидному разностному уравнению s(n) — s(n — 1) = a(n), (3) где а(х) - это некоторая аналитическая функция при х ? R, х> const. Это единственное ограничение, которое мы накладываем. При этом ряд может быть и знакопеременным (см. ниже). Начальный индекс суммирования в (2) не играет никакой роли, но обычно k=1 Таким образом, вместо суммы (2), т.е. «интеграла», предлагается изучать разностное уравнение (3), т.е. «ОДУ». Аналогия здесь очевидна. Отмеченная аналогия разностного уравнения и дифференциального на самом деле является весьма глубокой и хорошо изученной. Например, разностное уравнение (3) можно записать как уравнение (exp(D) — 1) s(n) = a(n + 1), где D - это оператор дифференцирования. Если двигаться по этому пути, то мы придем к методу суммирования Эйлера Миклорени. Но мы сосредоточимся на разностном уравнении (3). Разностное уравнение (3) иногда можно преобразовать в однородное трех-звенное линейное соотношение. Для линейных разностных уравнений с полиномиальными коэффициентами существует весьма развитая теория, которая, в частности, позволяет вычислить асимптотику решения при n — оо. Но в общем случае проще работать с уравнением (3). Предположим, что мы нашли какое-то решение s(n) уравнения (3) или, что более реалистично, его асимптотику при n — ос. Функция 5(х), вообще х (по крайней мере локально) задача о непрерывном (либо аналитическом) n — s(n) Пусть исходный ряд (2) сходится. Тогда можно считать, что 5(х) — 0 при х —о констатны С В любом случае именно такую асимптотику и предлагается искать. Таким образом, имеем два решения уравнения (3): одно дискретное (2) и одно непрерывное s^) = С + ё(х) которое, вообще говоря, может быть определено лишь асимптотически при хС ряда (2) при n — ос) однозначно восстанавливается условием согласования (т.е. равенства) частичных сумм исходного ряда со значениями функции s^) в дискретных точках х = n ? N, т.е. При этом, очевидно, достаточно установить это равенство только в одной (любой) точке n =1, 2,.... Во всех остальных точках это равенство будет выполняться автоматически. Таким образом, сумма ряда (2) при n — оо определена однозначно, но никакого предела при этом мы не вычисляли. Задача свелась к решению разностного уравнения (3) либо к вычислению асимптотики его решения. Разумеется, вычисление асимптотики предполагает, строго говоря, вычисление предела на каждом шаге. Но алгоритмически это чаще всего чисто алгебраическая операция. Возьмем в качестве примера Базельскую проблему, т.е. a(n) = 1/n2. Предположим, что мы ничего не знаем об этой проблеме и не умеем решать разностное уравнение (3). Но мы можем тем не менее вычислить асимптотику его решения s(n) методом неопределенных коэффициентов, что является простой алгебраической задачей. Имеем Заметим теперь, что асимптотические ряды такого рода хорошо изучены. В частности, известно, что величина погрешности при вычислении по таким рядам обычно оценивается последним отброшенным членом (см. [6]). Поэтому они суммируются до наименьшего по модулю члена. В данном случае, если взять n = 10, что эквивалентно суммированию первых 10 членов исходного ряда, а затем учесть найденную асимптотическую поправку, то мы найдем константу С = 7г2/6 с шестью десятичными знаками (на самом деле, с восемью). Исходная сумма (2) при n = 10 дает только два десятичных разряда. То, что мы получили, весьма напоминает метод суммирования Эйлера Маклорена, но только без вычисления интеграла. Другая аналогия, которая также очевидна, - это решение краевых задач методом погранслоя, когда решение получается в результате сращивания двух различных разложений, вычисленных в противоположных концах интервала. В случае, если ряд (2) расходится в обычном смысле, это означает, что последовательность частичных сумм {s(n)} не имеет предел а при n — оо, либо стремится к бесконечности. Так же ведет себя, очевидно, и решение s(n) уравнения (3). Однако s(n) (вообще говоря) имеет некоторую асимптотику, которую и предлагается вычислять. Иногда можно найти и явное s(n) s(n) в уравнении (4) могут сократить друг друга (как разрешение неопределенности оо — оо) и сумма расходящегося в обычном смысле ряда будет однозначно определена. Отметим, что смысл указанной операции не столько в назначении числового значения сумме расходящегося в обычном смысле ряда, сколько в установлении связи частичных сумм ряда с их асимптотическим поведением в бесконечности. a( n)=1/n ряда. Он также не может быть просуммирован обычными средствами. a( n)=1/n мы не будем этим пользоваться, а вычислим асимптотику его решения s(n) методом неопределенных коэффициентов. Имеем sin) = Inn Н т Н т + О —г = Inn - У] -г^1, KJ 2п 12п2 120п4 \пЧ itiknk' где B (к) это числа Вернул ли. Щк) Обозначив константу интегрирования С = 7, получаем из (4) знаменитое асимптотическое тождество Эйлера 7 = #(n)-lnn+?-Ay, n G N, (5) где H(n) - это гармоническое число. Метод суммирования, изложенный в этом разделе, нельзя назвать принципиально новым потому, что подобный подход уже использовался Ри.мину жи ном. Хотя Раманужан исходил из формулы Эйлера Миклорени. но он также пришел к сумме ряда как к «константе Раманужана» (см. [7, pp. ix, 3]). В следующих разделах мы приведем ряд примеров суммирования рядов, которые (как мы надеемся) покажут преимущества нашего метода. Но сперва преобразуем разностное уравнение (3) к функциональному уравнению на интервале [0,1] и дадим мотивировку этого преобразования. 3. ФУНКЦИОНАЛЬНОЕ УРАВНЕНИЕ ДЛЯ СУММЫ РЯДА Разностное уравнение (3) определено на интервале [0, оо], что может иметь иногда некоторые преимущества. Однако, как показывает опыт, функциональное уравнение на конечном интервале часто бывает удобнее. Сделаем замену переменных в уравнении (3), n = 1/x, S(x) = s(1/x), A(x) = a(1/x) («Мап muss immer umkehren.» C.G.J. Jacobi). Тогда получим функциональное уравнение A( x)= x2 ем уравнение, которое лишь знаком в левой части отличается от уравнения (А.7) в [9, р. 241]. Это последнее приводилось в [9] в качестве примера уравнения, решением которого является обыкновенная производящая функция x Уравнение (6) легко преобразуется к функциональному уравнению на интервале [0,1]. Для этого сделаем замену x — x/(1 + x). Тогда получим уравнение зШ-з(х)=АШ- (7) Одно из преимуществ функциональной формы записи разностного уравнения (3) - это возможность использовать развитый аппарат вычисления степенных рядов (и более общих разложений), реализованный в системах компьютерной алгебры (CAS). Другое не столь очевидное преимущество появляется, если уравнение (7) S( x) представить, например, в виде разложения по смещенным полиномам Чебы-шева. Так можно, в принципе, решить численно и с контролируемой точностью уравнение (7), т.е. вычислить сумму исходного ряда, даже если функция A(x) Численное решение уравнения (7) являлось первоначальной мотивировкой данного исследования. Однако это является отдельной задачей, которой мы намерены посвятить отдельную статью. Для вычисления асимптотики решения S(x) при x — +0 (которая может быть весьма сложной) уравнение (7) не имеет никаких преимуществ по сравнению с уравнением (6). Так что это последнее мы и будем далее использовать. В случае знакопеременного ряда запишем его в виде n s(n)= ? a(k)cosnk. (8) k=i Тогда, сделав те же замены, что и ранее, и взяв S(x) = U(x)cos п/x, получим уравнение U(x) + U [j^j = А(х) (9) вместо уравнения (6). Асимптотическое решение уравнения (9) имеет некоторые незначительные отличия от решения уравнения (6). Это, как мы полагаем, выгодно отличает наш подход от классического, который в случае знакопеременных рядов приводит к (весьма малоизвестной) формуле Эйлера Буле (см. [8]). Заметим также, что иногда вместо множителя cos пк =(—1)k следует брать множитель exp(ink), который бывает нужен для аналитического продолжения дискретного отображения n — s(n). Разностные уравнения типа (6) и (9) для достаточно простых правых ча-A( x) добные уравнения вводились в [9] комбинаторным путем для обыкновенных производящих функций чисел Бернулли и их обобщений. Многие из таких уравнений решаются в явном виде (см. [10]), хотя в [9] были даны лишь асимптотические разложения решений. Поскольку эти уравнения линейны, то это дает возможность комбинировать уже известные решения с целью получения новых. Приведем пример суммирования ряда с использованием известных решений уравнения (6). Рассмотрим ряд и найдем его сумму нашим методом и двумя способами. Для этого ряда уравнение (6) имеет вид S(x) - S {-^-А =х + 1п(1 - х). (11) Уравнение (11) лишь знаком левой части отличается от уравнения (А.14) в [9, р. 243], для которого мы нашли явное решение в [10]. Решение уравнения (11) имеет вид S(x) = х + In х + Ф (-\ = i х - х2 + -i- xA + 0{ \ x J 2 12 120 x6) Комбинируя полученное решение и формулы (4) и (10), получаем тождество n 1 n 1 k 7 = +1пп-Ф(п) + ]Г -_1п =Я(п-1)-Ф(п), n G N, (12) n k=2 Vк \к 1// которое является явным решением асимптотического тождества Эйлера (5). Заметим теперь, что правая часть уравнения (11) записывается как х + 1п(1 - х) = — У) Решение уравнения (6) с правой частью xn+\ n Е N0, имеет вид (см. [10]) ЗД = -Ц^ф(тг,1 + -У (13) n! V xj Комбинируя полученные решения и формулу (4), получаем тождество ,-.+?({-ь(^))+?!^±» м Используем интегральное представление функции Ф() в формуле (14): Ф(п, 1 + ж)= / eft, n G N. Тогда получим обобщение известной формулы для константы Эйлера: 7 = Н(п) -\пп+ I1 f-i- + —^—) Г"1 eft, n G N. ' v у -Я) Vln t 1 — t/ Ясно, что если суммировать не весь ряд (14) этим способом, а только его часть, то получим множество аналогов формулы Эйлера Миклорени для гармонического ряда с остаточными членами в интегральной форме. Пользуясь выражением гармонического числа через Ф—функцию, получим (см. [2]) 1 о (_i)k k s (к j) Щп + 1) = - + ln„ + Г(п) Е 1±Л- ?? -jLf, Re(n) > 0, где Si - это числа Стирлинга первого рода. и Из этой формулы можно получить бесконечное множество разложений константы Эйлера в (сколь угодно быстро) сходящиеся ряды рацональных чисел (например, подставить п = 2, затем п = 4 и исключить In2). Если теперь вычесть (12) из (14) и подставить туда п = ж, то получим тождество . ~ (-1)nЪ{п,х + 1) n=o (п + 1)! Заметим, что ряд Тейлора для логарифма от Гамма функции дает похожее тождество In* =ЕЖ ^>0. (16) n=o (п + 1)! x=1 (еще Эйлеру) представления константы Эйлера „ л f C(n + i)-i °° (-ir+1C(n + i) п . n=i п + 1 n=i п + 1 причем первый из этих рядов сходится значительно быстрее. Подстановка x = 2, 3,... в формулы (15) и (16) не дает, как легко проверить, ничего нового, так как все логарифмы ln п, п> 2 там сокращаются. Данный пример показывает, что решения уравнений (6) и (9), пользуясь их линейностью, можно составлять из уже известных решений. В следующем разделе мы приведем сводку формул для некоторых известных решений этих уравнений. Далее мы дадим ряд примеров суммирования рядов с их помощью. 4. РЕШЕНИЕ ФУНКЦИОНАЛЬНЫХ УРАВНЕНИЙ Решение уравнений S(x) - S (^-1 = xs+l (11 x при s = п ? N0 мы нашли в [10]. Это решение (13), что легко проверяется. Однако это решение не продолжается на вещественные либо комплексные s. s с нестандартными правыми частями можно решать с помощью подходящей замены переменных. Для уравнения (18) следует сделать замену S(x) = xs V(x). Тогда для функции V(x) получим уравнение V(x)-(l-x)-sV i-^A =х, которое не имеет логарифмической особенности в нуле. Для него следует использовать биномиальные тождества. Пользуясь известным асимптотическим разложением (см., например, [10 / 1 \ 00 Ф n,l + - = (-1)п+1(п-1)! J2 В{к)С{к + п-1,к)хп+\ п е N, (19) /г=0 где В(n) - это числа Бернулли и С() - это биномиальный коэффициент, для произвольного комплексного s = 0, получим (с помощью аналитического продолжения) асимптотическое решение S(x) = Y.B{n)C{n + s-l,n)xn. (20) s n=0 Для отрицательных целых s функция S(х) в (20) выражается через полиномы Бернулли: (-l)n S(x) = ± J- В{п, -1/ж), п = —s ? N. n Для s = 0 асимптотическое разложение (20) непригодно, и его следует заменить на разложение S(x) = Ф 1 + - = -1пж- Е -^*n- (21) V X/ n=1 n s=0 S(x) = -( (o,s + l,l + i), где Z() ~~ это С—функция Гурвица в обозначении, принятом в Maple, т.е. где первый аргумент - это номер производной этой функции. Последняя формула, впрочем, может быть использована скорее для вычисления функции Гурвица с использованием предыдущих формул, чем наоборот. Перейдем к решению уравнений (9) для знакопеременных рядов. Решение уравнения и(х) + U {jzr^ = %s+1 (22) для s = n ( N0 имеет вид ^=^('(»-1+Э4'(»-1+й)+{Г'»Й <»> Пользуясь разложением (19), получаем асимптотическое решение уравнения (22) для s G N, а также (аналитическим продолжением) и для произвольного комплексного s = 0: U(x) = — Y.B{n)C{n + s-l,n){2n-l)xn. (24) s n=1 Для отрицательных целых s функция U(x) в (24) выражается через полиномы Эйлера: (-1)n U(x) = ^—^-Е(п,-1/х), п = -s - 1 G N0, 2 что соответствует появлению этих полиномов в формуле Эйлера Буле для знакопеременных рядов (см. [7, р. 135] и [8]). Для s = 0 асимптотическое разложение (24) непригодно, и его следует заменить на разложение Щх) = Ш 2 - Ф f 1 + i) + Ф f 1 + -1) = t В{П) (1 " 2"' Л s=0 f/W = c(o,, + i,i + i)-ic(o,s + i,i + ^;). Приведенные формулы позволяют суммировать обычный ряд для Z—функции Римана во всей комплексной плоскости, кроме особой точки. Сращивая частичные суммы с асимптотическим разложением, получим по формуле (4) c(a+i) = f:J- + _LiBin)c(n+s-i,n) e^ N (25) n=l n sm n=0 m Данный ряд всюду расходится, но если его суммировать до наименьшего по модулю члена, то функцию Римана можно вычислять по формуле (25) в любой точке s G C, s = 0 с контролируемой точностью. Причем при выбранном s чем больше m, тем больше членов ряда можно использовать. Формула (25) аналогична той, которая получается по формуле Эйлера Маклорена, и которая ранее была основным вычислительным инструментом для функции Римана (см. [11]). s вестную формулу m 1 C(s + i) = Е —+ + -seN, n=1 n которая дает ((0) = —1/2, а также рациональные значения функции ((s + Г в точках s = —2 k и тривиальные нули в точках s = —2 k — 1, k ? N. И, наконец, при s ? N0 формула (25) принимает вид m 1 (—1)s^(s, 1 + m) ns+1 s! C(s + i)= E —- 1, ^eN0, n=1 где «s) = { tfsj = N',s> 1. ^ Расходящийся асимптотический ряд (25) можно преобразовать в сходящийся факториальный ряд, причем двумя способами (см. [2, 10]). Однако, как мы полагаем, факториальные сходящиеся ряды в данном случае не дают никаких преимуществ. Рассмотрим преобразование степенных рядов, которое (насколько нам известно) не описано в литературе, но может быть полезно для решения разностных уравнений, которые мы рассматриваем. Асимптотический ряд (19) можно записать следующим образом Ф 71,1 + - = -1Г Е п М ~Х > nGN' 27 что можно рассматривать как бесконечную треугольную матрицу, преобразующую вектор <x,x2,... > в вектор функций Ф(). Эта матрица обратима для любого порядка главного минора. Можно проверить, что мы получим сходящиеся ряды 1 оо (_ 1)к+1 / 1\ хп = - — Е A \is,* fc,l + - , n G N, Re(x) > 0, (28) (n — 1)! к=п (k — n +1)! x что, кстати, дает еще одно представление для чисел Бернулли в матричной форме. Формула (28) позволяет получать необычные разложения. Например, 1п(1+ж)-2?(Атф(2'-1л+Э- ™=й^^- Однако использование формулы (28) в (25) (x = 1/m) дает лишь другой асимптотический ряд. 5. ПРИМЕРЫ СУММИРОВАНИЯ ЗНАКОПЕРЕМЕННЫХ РЯДОВ Мы уже привели несколько примеров суммирования для сходящихся и расходящихся знакопостоянных рядов. Дадим примеры суммирования некоторых расходящихся знакопеременных рядов. Рассмотрим ряд оо д = ?(—1)n In n. (29) n=1 Суммирование этого ряда приводилось в [12] как демонстрация эффективности предложенного численного алгоритма. Заметим, что Maple также без труда находит численное значение суммы этого ряда, которое оказывается равным 9 = \ Ь (f) ~ 0-225791352644727- Эвристически угадать это значение можно из представления функции Римана оо ( (l-21-sR(s)= ?- ~—, Re(S)>0. n=1 n Если продифференцировать это тождество по s и подставить туда s = 0 (что, вообще-то, незаконно), то получим как раз нужное значение. Таким образом, ряд (29) «помнит» о своем происхождении. Покажем, что ряд (29) суммируется нашим способом в явном виде (чего Maple сделать не может). Действуя описанным в разд. 3 способом, получим функциональное уравнение U{x) + U (j^-^J = -Ьж. Сделаем в нем подстановку U(x) = V(x) + W(x)lnx. Тогда получим уравнение + v (гУ -w Ш111(1 - *>=- К> + w Ш + х)Ых- Правая часть этого уравнения обнуляется при W(x) = —1/2. Поэтому получаем уравнение Заметим теперь, что последнее уравнение мы получили бы описанным в разд. 3 способом из сходящегося ряда. Поэтому 1 о x ^^?(-i)Mn(i-_). Обычно С AS не может суммировать такие ряды. Поэтому преобразуем его к виду °° (1 + (2 n- 1) x)2 exp(2 V (x)) = П 1(1 + (2п-2)ж)(1 + 2пж) Это произведение вычисляется в явном виде. Получаем ^ = ьгШ-1пГЙ + Й-51п(2^ ж>а Заметим, что если бы целью являлось численное решение задачи (с неограниченной точностью), то мы могли бы вычислить асимптотическое разложение решения методом неопределенных коэффициентов, т.е. т,/ ч 1 1 3 1 5 17 7 31 9 691 11 V (х) = - х т х -\ х х Н х X + . . . v ; 4 24 20 112 36 88 Появление числа 691 в асимптотическом разложении V(x) является явным индикатором того, что в нем присутствуют числа Бернулли и данное разложение расходится. Поэтому другой способ просуммировать такие ряды - это найти достаточно много коэффициентов асимптотического разложения, угадать их явный вид и применить средства CAS. В данном случае v ; n=! 2п(2п-1) Собирая полученные формулы и применяя формулу (4), получим сумму ряда (29) g, как и ожидалось. Однако, как мы отмечали, суммирование расходящегося в обычном смысле ряда имеет вполне прозаический смысл. В данном случае мы нашли частичную сумму ряда (29) в замкнутой форме (чего Maple сделать не может), а также получили аналитическую интерполяцию частичных сумм ряда (29). Эта формула имеет вид На рис. 1 показаны график частичных сумм ряда (29) в дискретных точках и график их интерполяции по данной формуле. В завершии этого примера отметим, что найденная нами функция V(х), которая описывает асимптотическое поведение частичных сумм ряда (29) при n — оо, связана с тангенциальными числами следующей формулой где T(х) - это обыкновенная производящая функция тангенциальных чисел, для которой мы нашли аналитическое представление (см. [10]). Цепная дробь производящей функции тангенциальных чисел дана в [9, р. 257], Нам удалось просуммировать нашим способом все ряды, приведенные в [12] в качестве примеров применения численного алгоритма в [12]. Если принять в качестве постулата, что суммирование асимптотических рядов следует производить до наименьшего по модулю члена (что требует отдельного исследования в каждом конкретном случае), то наш метод суммирует эти ряды с произвольной точностью и, по сути, алгебраическими средствами, программируемыми в С AS. Применим наш метод к классическому ряду (1) при х = —1, т.е. найдем константу Эйлера Гомпертци 6 = e Ei(1,1). Действуя описанным в разд. 3 способом, получим функциональное уравнение Сделаем в нем подстановку U(х) = Г(1+1/х) W(ж). Тогда для W(х) получим уравнение Данное уравнение отличается от стандартных функциональных уравнений, которые мы рассматривали, но алгоритмически (и асимптотически) решается тем же способом. Получаем W (х) = 1 — х + ж2 — ж4 — ж5 + 2 ж6 + 9 ж7 + 9 ж8 — 50 ж9 — 267 ж10 + 0(хп). Последовательность целых коэффициентов этого разложения можно найти в энциклопедии целых последовательностей (см. [13, А000587]). Она связана с числами Белла. С помощью формулы суммирования Бореля получаем 00 W(x) = 1 — x J exp(—t)exp(1 — exp(tx)) dt, x > 0. (31) 0 Фактически мы уже решили задачу. Однако Maple не может вычислить x>0 еще один способ решения функционального уравнения (30), который обладает большей универсальностью. Сделаем в уравнении (30) замены, обратные тем, которые мы применили в начале разд. 3. Тогда получим из (30) полиномиальное разностное уравнение nw(n) + w(n — 1) = n, n G N. (32) Maple умеет решать такие уравнения в явном виде. Получаем _ n-l + eEi(n,l) еЕ1(1,1)-1 + Ц0) 1 ]~ п +[ } Г(п + 1) w(0) Осциллирующая часть этого выражения удовлетворяет однородному уравнению (32). Но осциллирующей части здесь не должно быть, так как она w(0) = 1 — e Ei(1,1). Таким образом, имеем еще одно решение уравнения (30): W(x) = 1 - х + х е Ei ^i, 1 j . Это решение совпадает с предыдущим, т.е. мы вычислили в явном виде пре- w(0) определилась бы также из совпадения этих решений. Собирая полученные формулы и применяя формулу (4), получим интерполяционную формулу для суммы ряда (1) при x = — 1: х^/ ^kl, , Л 1 eEi(n, 1)\ . . y^(-\f k\ = 5 + п\ 1 1 cos(7rn), k=o V n n I где 5 = e Ei(1, 1) - это принятая сумма данного ряда, т.е. константа Эйлера Гомпертца. Таким образом, как мы отмечали, сумма расходящегося ряда - это некоторая однозначно определяемая константа, необходимая для согласования дискретной последовательности частичных сумм ряда с их непрерывной (либо аналитической) интерполяцией в бесконечности. С этой точки зрения никакого различия между сходящимися и расходящимися рядами нет. 6. ОБОБЩЕНИЕ ФОРМУЛЫ СТИРЛИНГА Здесь мы рассмотрим пример суммирования расходящегося знакопостоянного ряда, который приводит к обобщению формулы Стирлинга. Рассмотрим ряд то h = ?ln n. (33) n=1 Суммирование этого ряда с помощью формулы Эйлера Миклорени дает асимптотическую формулу Стирлинга для факториала. Мы получим эту формулу в виде сходящегося ряда. Ряд (33) также «помнит» о своем происхождении, что видно, если продифференцировать ряд для функции Римана ж i C(s)=?- Re(s)>l, n=1 n и подставить туда s = 0 (что, разумеется, незаконно). Тогда получим h = -('(О) = ^1п(2тг) « 0.918938533204, но мы не будем этим пользоваться. Действуя описанным в разд. 3 способом, получим функциональное уравнение S(x) - S {j^j = -Ьж. Сделаем в нем подстановку S(x) = V(x)+W(x)lnx. Тогда получим уравнение щх) -уШ+w Шln(l -х)=- Ьх) -w Ш+х)Ых- Правая часть этого уравнения обнуляется при W(x) = —1/x. Поэтому получаем уравнение _ у ( * \ = (1-х) \п(1-х) = ~ хп+1 У } \l-xj х ^пй(п + 1)(п + 2)' Далее мы можем применить асимптотические формулы (20) и (21) и собрать с их помощью асимптотическое решение для функции V (x). Но мы применим формулу (13) и получим функцию V(x) в виде сходящегося ряда. Комбинируя найденные решения и формулу (4), получаем сумму ряда (зз) ^ k h=±ink + n(i-inn)-f: »€N, k=1 k=0 (k + 2)! откуда, подставив п = 1, 2,..., можно получить (с помощью Maple) численное значение константы Н. Чтобы получить точное значение, запишем эту формулу в виде (к + 2)! Правая часть этого выражения стремится к нулю при п — ОО, а левая часть имеет асимптотическое разложение 1 h In 2 1п7гН —- 2 2 бп 24п2 +... поэтому Н = —('(0), как мы нашли в начале этого раздела. Таким образом, получаем вариант формулы Стирлинга в виде сходящегося ряда ! 1 О (_1)к \U(k х + 1 1пГ(ж + 1) = -1п(2тг)-ж(1-1пж) + Е Г. , Ке(ж)>0. (34^ 2 к=0 (k + 2) 2 Данная формула справедлива и при x = 0, если перейти к пределу. Тогда получим тождество - = -ln(2.)-g(fc + 1)(fc + 2). Этот медленно сходящийся ряд можно ускорить с помощью метода Куммера. Формулу (34) можно преобразовать в обычный асимптотический ряд Стир-линга с помощью преобразования (27), сделав в нем замену х — 1/x. Однако существует более общее преобразование, включающее логарифм и функцию Ф(0,х) = Ф(х). ?1)п B (к — п) Рассмотрим две бесконечные верхне-треугольные матрицы (k — п + 1)! (к — п)! v -1)к 1<п<к V,v 'V' 1<п<к 1)! п 1 1 2 ; ; 9; о; Тогда можно проверить, что P = Q—1 и Q = P—1 формально. Причем это справедливо для любого главного минора. Рассмотрим также два бесконечных вектора X In х п * t Y = (ф(п — 1,х + Тогда можно проверить, что Y = Q • X а X = P • Y формально. Причем ком- YX ряды. См. формулы (15), (28). Применив преобразование Р к ряду (34), получим (после некоторых преобразований) формулу Стирлинга в асимптотической форме 1 1 00 B (к + 1) In Г(* + 1) = - 1п(2 тг) - * (1 - In *) + - Нх) + Е к((к + 1)'хк- Расходящийся ряд в формуле (35) можно преобразовать к сходящимся факториальным рядам (см. [2, 10]). Однако это не дает никаких преимуществ, так как в коэффициентах факториальных рядов не удается выявить никаких закономерностей. Эта же проблема возникает, если преобразовать ряд в (35) в (сходящуюся) цепную дробь, что, вероятно, делал еще Стилтьес. Если дифференцировать формулу (35) по ж, то мы получим известные асимптотические разложения для полигамма функций. Однако и для них, за исключением Ф(1,ж + 1)и Ф(2,ж + 1), не удается получить цепные дроби с «квазипериодической» структурой (см. также [2]). Для факториальных рядов та же картина, за исключением ряда для Y = — который мы нашли в [2] и который имеет хотя и сложную, но предсказуемую структуру. Если же дифференцировать формулу (34) по ж, то мы получим сходящиеся разложения полигамма функций. Поскольку значения этих функций в натуральных аргументах выражаются с помощью С—функции Римана, то мы получаем бесконечный набор тождеств для С—функции Римана. Введем вектор-функцию W = {w(n))t1<n ,w(n) n 1 ne N, 1<i<j (J — i + 1)! G G = (0,1,1,2,6,...,(i — 2)!,...У . G дующие компоненты дадут 0 0 0 ЕС(J) —1 = 1, Е(С(J) — —1) = 1, Е(С(J) — — — 2) = 2, .... j =2 j=3 j =4 Эти формулы, вероятно, известны, но мы не смогли найти их в литературе. 7. ЗАКЛЮЧЕНИЕ Мы представили метод суммирования сходящихся и расходящихся рядов, который является аналогом методов Эйлера Миклорени и Эилери Буле, но проще их в реализации и выводе. Метод ориентирован на использование имеющегося весьма развитого инструментария CAS и позволяет в ряде случаев суммировать ряды чисто алгебраическими средствами. Понятие суммы ряда получило толкование, которое не использует понятие предела, что позволяет рассматривать сходящиеся и расходящиеся (в обычном смысле) ряды с единой позиции. Основной вычислительной задачей при этом является вычисление асимптотических разложений решений некоторых классов функциональных уравнений. Эта задача в общем случае является весьма сложной и малоисследованной, но во многом аналогична задаче вычисления асимптотических разложений решений ОДУ. В частности, помимо обычных асимптотических рядов, которые мы получали в данной статье, могут появляться т.н. г ряды (см. [14]), а также более сложные разложения, известные под названием «трансряды» (см. [15]). Развитие этой теории тесно связано с развитием С AS. Суммирование медленно сходящихся рядов со сложно устроенным общим членом было бы, вероятно, непрактично осуществлять с помощью асимптотических разложений. Однако наш метод суммирования позволяет сделать это численно, путем численного решения соответствующих функциональных уравнений. К этой теме мы намерены вернуться в одной из следующих публикаций. 521.1+531.314 10.20948/prepr-2022-28 Непрерывные и дискретные модели математической физики и асимптотические методы, 0017-2019-0004 суммирование рядов функциональные уравнения расходящиеся ряды summation of series functional equations divergent series Харди Г. Расходящиеся ряды. М: (1951). Варин В.П. Преобразование последовательностей в доказательствах иррациональности некоторых фундаментальных констант // Ж. вычисл. матем. и матем. физ. N. 9. (в печати) (2022). Lanczos C. Applied Analysis. Dover Publications, New-York: (1956). Варин В.П. Об интерполяции некоторых рекуррентных последовательностей // Ж. вычисл. матем. и матем. физ. Т. 61. N 6. с. 913-925. (2021). Варин В.П. Инвариантные кривые некоторых дискретных динамических систем // Ж. вычисл. матем. и матем. физ. Т. 62. N 2. с. 199-216. (2022). Olver F.W.J. Asymptotics and special functions. Acad. Press, New-York: (1974). Candelpergher B. Ramanujan Summation of Divergent Series. Lecture Notes in Math., Springer. (2017). Borwein J.M., Calkin N.J., Manna D. Euler--Boole Summation Revisited // The American Mathematical Monthly. V. 116. N. 5. P. 387-412. (2009). Arakawa T., Ibukiyama T., Kaneko M. Bernoulli Numbers and Zeta Functions. Springer, Japan. (2014). Варин В.П. Факториальное преобразование некоторых классических комбинаторных последовательностей // Ж. вычисл. матем. и матем. физ. Т.59. N 6. С. 1747-1770. (2018). Borwein J.M., Bradley D.M., Crandall R.E. Computational strategies for the Riemann zeta function // Journal of Computational and Applied Mathematics. V. 121. p. 247-296. (2000). Cohen H., Villegas F.R., Zagier D. Convergence Acceleration of Alternating Series // Experimental Mathematics, V. 9. N. 1. p. 3--12. (2000). Sloane online encyclopedia of integer sequences. (http://oeis.org). Hille E. Ordinary differential equations in the complex domain. New-York: John Wiley \& Sons. (1976). Edgar G.A. Transseries for beginners // [arxiv:0801.4877v5], (2009). (http://arxiv.org/abs/0801.4877v5). https://keldysh.ru/papers/2022/prep2022_28.pdf