Простая одержимость. Бернхард Риман и величайшая нерешенная проблема в математике. - Дербишир Джон

Чтобы посмотреть, что же получается в результате этого процесса, и разобраться в том, почему Риман назвал этот вторичный член «периодическими членами», поупражняемся немного в арифметике, используя конкретные значения буквы x. Как и раньше, возьмем x = 20; тем самым мы вычисляем величину J(20) — что, как несложно проверить из исходного определения функции J, равно 97/12 т.е. 9,5833333…. Вот как это получается.

Сначала возводим 20 в степень 1/2 + 14,134725i. В результате получаем точку, которая на рисунке 21.2 помечена как 1 и численно выражается как −0,302303 − 4,46191i. Интегральный логарифм от этого — т.е. функция Li — дает самую западную точку на рисунке 21.3, выражаемую числом −0,105384 + 3,14749i. Теперь разберемся с сопряженным членом из этой пары нулей. Возводим 20 в степень 1/2 − 14,134725i. Результат равен −0,302303 + 4,46191i. Он показан на средней картинке на рисунке 21.4. Это зеркальный образ точки, помеченной на рисунке 21.2 как 1, относительно вещественной оси. Берем интегральный логарифм и получаем ответ −0,105384 − 3,14749i — точку, лежащую глубоко на юге в правой части рисунка 21.4. Складывая два ответа, получаем −0,210768. Мнимые части, разумеется, сократились. Вот и все с первой парой сопряженных нулей.

Повторим все это для второй пары, 1/2 + 21,022040i и 1/2 − 21,022040i. На этот раз окончательный ответ будет равен 0,0215632. Для третьей пары он равен −0,0535991. С тремя парами мы разобрались, но впереди бесконечность!

После 50 таких вычислений получаем (таблицу следует читать по колонкам):

Первое значение представляет собой некоторую аномалию, поскольку самая западная точка на рисунке 21.3 отстоит от вертикальной оси более чем в два раза дальше, чем остальные. Однако затем числа в таблице уменьшаются по мере того, как значения, соответствующие северной половине критической прямой, по спирали приближаются к πi. И взгляните на их знаки — имеется примерно равное число положительных и отрицательных.[199] Это хорошая новость, потому что, хотя ответы и становятся меньше, они делают это не очень быстро, и нам потребуется вся возможная помощь, которую могут нам оказать сокращения между положительными и отрицательными значениями. Не будем забывать, что все это происходит под знаком суммы — эти 50 чисел предстоит еще сложить друг с другом. (Сумма равна −0,343864, что, кстати, составляет не более 8 процентов от полной бесконечной суммы. Не так плохо для всего лишь 50 слагаемых.)

Рисунок 21.5. Первые 50 значений, полученных путем взятия нетривиального нуля и его комплексно сопряженного, вычисления значений функции Li(20z)и их последующего суммирования.

Из рисунка 21.5 видно, почему Риман назвал эти компоненты вторичного члена «периодическими». Они изменяются нерегулярным образом (что означает, если уж быть совсем скрупулезным, что они не строго «периодические», а только «колебательные») вверх и вниз от положительных к отрицательным значениям и обратно.[200] Причина этого совершенно ясна из рисунка 21.3. Колебательная природа вторичных членов связана с тем, что, как видно из рисунка 21.3, функция Li(xρ) скручивает критическую прямую во все более и более плотную спираль. Значения функции, соответствующие нулям дзета-функции, могут при этом оказаться где угодно на этой спирали; определяющая причина состоит в том, что для больших x критическая прямая чрезвычайно сильно растягивается перед закручиванием. Закручивание настолько плотное, что высоко расположенный отрезок критической прямой отображается в нечто очень близкое по форме к окружности. В силу этого получается, что значения функции Li(xρ) в нулях дзета-функции выглядят примерно как точки, раскиданные по окружности. Если вы немного знакомы с тригонометрией, то вам известно, что это приводит нас в мир синусов и косинусов, волновых функций, колебаний, вибраций… музыки. Именно отсюда и взялось введенное сэром Майклом Берри понятие «музыка простых чисел».

По мере прибавления новых членов сами они убывают, а положительные и отрицательные до некоторой степени сокращают друг друга при суммировании, так что мы зарабатываем сходимость. Эта сходимость, правда, страшно медленная. Для получения результата с точностью в три значащие цифры приходится складывать более 7000 членов; в четыре цифры — более 86 000. На графике на рисунке 21.6 показаны первые 1000 результатов (хотя некоторые из самых левых при выбранном масштабе оказались за пределами рисунка); на этот раз не делается никаких попыток соединить точки между собой. Видно, что члены под знаком суммы действительно уменьшаются, хотя и делают это с достаточной ленцой.

Рисунок 21.6. То же, что на рисунке 21.5, но показана 1000 значений (точки не соединены между собой).

Окончательный результат равен −0,370816425…. Это, как мы помним, второй член в выражении (21.1). Первый же член — это в нашем случае Li(20), равный 9,90529997763…. Третий равен ln 2, что составляет 0,69314718055994…. И четвертый член, тот самый надоедливый интеграл, добавляет пустячный результат 0,000364111…. Подставим все это в выражение (21.1) и — хлоп! — J(20) = 9,58333333… (что мы, конечно, и так знали).

Закончим тем, что с использованием формулы Римана проведем полное вычисление π(1000 000) — т.е. числа простых чисел в пределах одного миллиона — не ради веселья, хотя веселье и немалое, а для того, чтобы сделать несколько важных замечаний по поводу остаточного члена.

Как мы помним из главы 19.iv,

Сколько же членов в правой части надо вычислять? До тех пор пока числа в скобках не станут меньше 2, потому что J(x) равна нулю, когда x меньше 2. Корень девятнадцатой степени из 1000 000 равен 2,069138…, а корень двадцатой степени 1,995262… Следовательно, можно остановиться на 19. Поскольку число 19 свободно от квадратов и имеет только один простой делитель — самого себя, — функция Мебиуса μ(19) имеет значение −1. Таким образом, последний член в правой части равен −1/19J(19√1000 000). Всего в правой части будет 13 слагаемых, поскольку между 1 и 19 функция Мебиуса принимает ненулевые значения 13 раз — при аргументах 1, 2, 3, 5, 6, 7, 10, 11, 13, 14, 15, 17, 19. Напомним, что функция Мебиуса равна нулю всякий раз, когда аргумент делится на точный квадрат (например, 4 или 9).