Часть IV · Анализ Глава 29 из 60
Число e и экспонента
Одно и то же число 2,71828… замечали в банке, на графике, под гиперболой, в радиоактивном распаде, на окружности, в гардеробе и на смотринах невесты. Запишем все встречи с ним, как натуралист в полевом дневнике, и выясним, одно ли это число или несколько похожих.
Опирается на: 28 · Интеграл
Вы научитесь
- объяснять, почему сложные проценты, площадь под гиперболой и функция, равная своей производной, приводят к одному числу e
- решать задачи о непрерывном росте, распаде и остывании и находить период полураспада
- вычислять e и степени e с любой точностью по ряду из факториалов
- понимать формулу Эйлера и узнавать 1/e в задачах о шляпах, о лотерее и о выборе лучшего
Прошлая глава оставила трещину в самой удобной своей формуле. Площадь под графиком $y = x^n$ даёт первообразная $\frac{x^{n+1}}{n+1}$ — при любом $n$, кроме одного: при $n = -1$ в знаменателе оказывается ноль. А гипербола $y = \frac1x$ никуда не делась, площадь под ней есть, её можно закрасить и измерить. Была и подсказка: прямоугольники Ферма над отрезками $[1;\ q]$, $[q;\ q^2]$, $[q^2;\ q^3]$, … все одной площади, так что площадь под гиперболой растёт, как число шагов геометрической прогрессии, — как логарифм. Какой именно? И второй вопрос, оставленный напоследок: есть ли функция, которая совпадает со своей производной, то есть растёт ровно с той скоростью, какую показывает её собственное значение?
Оба вопроса приводят к одному числу, $2{,}71828\ldots$, которое обозначают буквой $e$. Странно другое: это число появляется там, где его никто не ждёт. В банковском договоре. На графике показательной функции. Под гиперболой. В древней кости, где распадается радиоактивный углерод. На окружности. В гардеробе, где перепутали шляпы. И даже на смотринах невесты. Будем вести записи, как натуралист, который встречает похожую птицу то в лесу, то на болоте, то посреди города. Сначала отметим, где и при каких обстоятельствах её видели, а потом разберёмся, один ли это вид или несколько похожих.
Банк: проценты на проценты
Первое наблюдение записал Якоб Бернулли, старший из математиков знаменитой базельской семьи. В 1683 году он разбирал задачу, которую мог бы задать любой вкладчик. Банк обещает $100\,\%$ годовых. Кладём рубль — через год получаем два. А если проценты начисляют дважды в год, по $50\,\%$ за полгода? Через полгода на счёте $1{,}5$ рубля, и следующие $50\,\%$ начисляются уже на полтора рубля, а не на один. К концу года выходит $1{,}5^2 = 2{,}25$ рубля. Проценты на проценты выгоднее.
Будем дробить год всё мельче. При $n$ начислениях каждое добавляет $\frac1n$ от суммы, которая лежит на счёте, то есть умножает её на $1 + \frac1n$. За год рубль превращается в $\left(1 + \frac1n\right)^n$ рублей.
| Проценты начисляют | $n$ | Через год, рублей |
|---|---|---|
| раз в год | $1$ | $2$ |
| раз в полгода | $2$ | $2{,}25$ |
| раз в квартал | $4$ | $2{,}441\,406\ldots$ |
| раз в месяц | $12$ | $2{,}613\,035\ldots$ |
| раз в день | $365$ | $2{,}714\,567\ldots$ |
| раз в час | $8760$ | $2{,}718\,127\ldots$ |
| раз в минуту | $525\,600$ | $2{,}718\,279\ldots$ |
Сумма растёт, но прибавки тают. Переход с ежемесячного начисления на ежедневное приносит ещё десять копеек, с ежедневного на ежечасное — треть копейки, с ежечасного на ежеминутное — полторы сотых копейки. Рост явно во что-то упирается. Бернулли увидел, что суммы подходят к пределу, лежащему между двумя и тремя, но самого числа не нашёл. Запишем в дневник: $2{,}71828\ldots$
При ставке $r$ картина та же: $\left(1 + \frac rn\right)^n$ растёт с ростом $n$ и подходит к своему пределу. Такое начисление «бесконечно часто» называют непрерывным. Сыграйте в банк: частоту выбирайте в таблице, ставку — ползунком.
В выражении $\left(1 + \frac1n\right)^n$ основание стремится к единице, а единица в любой степени равна единице. Значит, при $n \to \infty$ оно стремится к $1$?
Основание тянет к единице, показатель — к бесконечности, и исход перетягивания зависит от того, кто быстрее. Такую запись называют неопределённостью вида $1^\infty$. Сравните: $\left(1 + \frac1{n^2}\right)^n \to 1$, потому что основание подходит к единице слишком быстро, а $\left(1 + \frac1n\right)^{n^2}$ растёт без предела. В нашем случае силы уравновешены, и получается $e$.
Пока у нас только таблица. Ни того, что предел существует, ни того, что он меньше трёх, мы не доказали. Вернёмся к этому, когда наблюдений станет больше, а пока — второе место.
Склон: функция, равная своей скорости
Второе место — график. Показательная функция $y = a^x$ из главы о степенях за одинаковые промежутки времени растёт в одно и то же число раз. А с какой скоростью? Скорость — это производная, и её можно посчитать по определению. Разностное отношение раскладывается на два множителя:
$$\frac{a^{x+h} - a^x}{h} = a^x \cdot \frac{a^h - 1}{h}.$$Второй множитель от $x$ не зависит. При $h \to 0$ это наклон графика $y = a^x$ в точке $x = 0$ — число, которое зависит только от основания. Обозначим его $k(a)$.
Пусть у графика $y = a^x$ ($a > 0$) в точке $x = 0$ есть касательная с наклоном $k(a)$. Тогда в каждой точке $x$ производная существует и равна $(a^x)' = k(a)\,a^x$: скорость роста пропорциональна значению.
Хитрость в равенстве $a^{x_0 + h} = a^{x_0} \cdot a^h$: график возле любой точки $x_0$ — это график возле нуля, поднятый в $a^{x_0}$ раз.
Производная показательной функции пропорциональна ей самой. Разные основания отличаются только множителем $k(a)$. У двойки он около $0{,}693$: функция $2^x$ растёт медленнее своего значения, её касательная в точке $(0;\ 1)$ идёт полого. У тройки — около $1{,}099$, быстрее значения. Где-то между двойкой и тройкой прячется основание, у которого $k = 1$, и для него производная совпадает с функцией. Поищите его сами.
Нашли? У того, кто подбирал основание аккуратно, получилось $2{,}718\ldots$ В дневнике второе появление того же числа. Но сходство ещё не тождество: число найдено на ощупь, с тремя знаками, и никакой причины, по которой банку и графику понадобилось одно и то же число, у нас нет. Больше того, мы даже не доказали, что наклон $k(a)$ существует. Нужно третье наблюдение.
Гипербола: площадь, которая складывает
Вернёмся к вопросу, с которого началась глава, — к площади под гиперболой. Обозначим через $A(x)$ площадь под графиком $y = \frac1t$ от единицы до $x$:
$$A(x) = \int_1^x \frac{dt}{t}, \qquad x > 0.$$При $x > 1$ это обычная площадь, при $x < 1$ — площадь от $x$ до единицы, взятая со знаком минус (так интеграл устроен по определению, см. главу об интеграле), и $A(1) = 0$. Первообразную по готовой формуле не найти, но численно площадь считается легко: $A(2) \approx 0{,}693$, $A(3) \approx 1{,}099$, $A(4) \approx 1{,}386$, $A(8) \approx 2{,}079$.
Два совпадения сразу. Первые два числа — те самые наклоны $k(2)$ и $k(3)$ из прошлого раздела. А последние два — удвоенное и утроенное $A(2)$: площадь до четырёх вдвое больше площади до двух, до восьми — втрое. Так ведёт себя логарифм: $\log 4 = 2\log 2$, $\log 8 = 3 \log 2$. Совпадение не случайно.
Для любых положительных $a$ и $b$ выполняется $A(ab) = A(a) + A(b)$.
Хитрость в том, чтобы растянуть фигуру по горизонтали и сжать по вертикали в одно и то же число раз: площадь при этом не меняется, а гипербола переходит сама в себя. Нарисуем случай $a > 1$, $b > 1$.
$A\left(\frac1x\right) = -A(x)$, и $A(x^q) = q\,A(x)$ для любого рационального $q$ и любого $x > 0$.
Из теоремы при $b = \frac1a$ получаем $A(a) + A\left(\frac1a\right) = A(1) = 0$, откуда $A\left(\frac1a\right) = -A(a)$. Для натурального $k$ применим теорему $k - 1$ раз к произведению $x \cdot x \cdot \ldots \cdot x$ и получим $A(x^k) = k A(x)$; с предыдущим равенством это верно и для целых отрицательных $k$. Наконец, если $q = \frac mn$, положим $y = x^{m/n}$. Тогда $y^n = x^m$, и по уже доказанному $n A(y) = A(y^n) = A(x^m) = m A(x)$, то есть $A(x^{m/n}) = \frac mn A(x)$.
Площадь под гиперболой — логарифм. Но по какому основанию? Основание логарифма — то число, логарифм которого равен единице. Ищем $x$, при котором $A(x) = 1$. Площадь растёт вместе с $x$ и меняется без скачков; $A(2) \approx 0{,}693$ меньше единицы, а $A(3) \approx 1{,}099$ больше, поэтому такое $x$ есть, лежит между двумя и тремя и единственно (здесь работает теорема о промежуточном значении из главы о пределах). Численно оно равно $2{,}71828\ldots$ Третье наблюдение — и снова то же число.
Одно и то же число
Наблюдений три. Банк: предел $\left(1 + \frac1n\right)^n$. Склон: основание, у которого наклон в нуле равен единице. Гипербола: число, до которого площадь под ней равна единице. Пора доказать, что это одно число, а не три похожих. Всё решает одна узкая полоска под гиперболой.
Для любого $h > -1$, $h \ne 0$, выполняется $\dfrac{h}{1 + h} < A(1 + h) < h$.
Гипербола на отрезке $[1;\ 1 + h]$ убывает, поэтому зажата между своей самой низкой и самой высокой точкой.
Сначала сверим банк с гиперболой. Обозначим пока через $e$ число, для которого $A(e) = 1$.
При каждом натуральном $n$
$$\left(1 + \frac1n\right)^n < e < \left(1 + \frac1n\right)^{n+1},$$и поэтому $\left(1 + \frac1n\right)^n \to e$ при $n \to \infty$.
Подставим в лемму $h = \frac1n$. Поскольку $\frac{1/n}{1 + 1/n} = \frac1{n+1}$, получаем $\frac1{n+1} < A\left(1 + \frac1n\right) < \frac1n$.
Умножим правое неравенство на $n$. По следствию о степенях $n A\left(1 + \frac1n\right) = A\left(\left(1 + \frac1n\right)^n\right)$, так что $A\left(\left(1 + \frac1n\right)^n\right) < 1 = A(e)$. Площадь растёт вместе с правой границей, поэтому меньшей площади соответствует меньшее число: $\left(1 + \frac1n\right)^n < e$. Точно так же, умножив левое неравенство на $n + 1$, получаем $A\left(\left(1 + \frac1n\right)^{n+1}\right) > 1$ и $\left(1 + \frac1n\right)^{n+1} > e$.
Осталось заметить, что вилка сужается. Верхняя граница больше нижней в $1 + \frac1n$ раз, поэтому их разность равна $\frac1n\left(1 + \frac1n\right)^n < \frac en$. Число $e$ меньше четырёх: $A(2) > \frac12$ (прямоугольник шириной $1$ и высотой $\frac12$ помещается под гиперболой на $[1;\ 2]$), значит, $A(4) = 2A(2) > 1 = A(e)$. Итого $0 < e - \left(1 + \frac1n\right)^n < \frac4n$, и эта разность становится меньше любого положительного числа.
Банк и гипербола говорят об одном и том же числе. Теперь у него есть имя и строгое определение.
Число $e$ — предел последовательности $\left(1 + \frac1n\right)^n$: $e = 2{,}718\,281\,828\,459\ldots$ Площадь под гиперболой $y = \frac1t$ от $1$ до $e$ равна единице.
Теперь назовём площадь её настоящим именем. В главе о логарифмах натуральный логарифм $\ln x$ был логарифмом по основанию $e$.
Для каждого $x > 0$ выполняется $\ln x = A(x) = \displaystyle\int_1^x \frac{dt}{t}$.
Для чисел вида $e^q$ с рациональным $q$ равенство следует из следствия о степенях: $A(e^q) = q A(e) = q = \ln e^q$. Возьмём любое $x > 0$ и допустим, что $A(x) < \ln x$. Между любыми двумя разными числами есть рациональное (дроби стоят на прямой всюду плотно, глава 6), возьмём такое $q$: $A(x) < q < \ln x$. Из $q < \ln x$ следует $e^q < x$, потому что логарифм по основанию $e > 1$ возрастает. Площадь тоже возрастает, поэтому $A(e^q) < A(x)$, то есть $q < A(x)$. Противоречие. Случай $A(x) > \ln x$ разбирается так же. Остаётся $A(x) = \ln x$.
Эту теорему можно прочитать и наоборот и сделать формулу $\ln x = \int_1^x \frac{dt}{t}$ определением логарифма. Тогда всё строится заново без оговорок. Функция $\ln x$ возрастает, меняется без скачков и принимает любые значения: $\ln 2^k = k\ln 2$ уходит в плюс бесконечность, $\ln 2^{-k} = -k\ln 2$ — в минус бесконечность. Значит, для каждого действительного $x$ есть ровно одно число $y > 0$ с $\ln y = x$. Его и называют $e^x$, а степень с любым основанием $a > 0$ определяют равенством $a^x = e^{x \ln a}$. В главе о степенях мы обещали аккуратно разобраться с показателями вроде $\sqrt 2$ — вот это определение. Для рациональных $x$ оно даёт то же, что раньше (по следствию о степенях $\ln a^q = q\ln a$), а правила степеней следуют из свойства площади: логарифмы чисел $e^x e^y$ и $e^{x+y}$ оба равны $x + y$, значит, сами числа равны.
Экспонента — показательная функция $y = e^x$ с основанием $e$; её пишут и как $\exp x$. Она обратна натуральному логарифму: $e^{\ln x} = x$ при $x > 0$ и $\ln e^x = x$ при любом $x$.
Теперь производные. Они уже встречались в главе о производной, но здесь обе функции построены заново, из площади под гиперболой, и производные мы выведем, опираясь только на основную теорему анализа. Логарифм — накопленная площадь под гиперболой, а по основной теореме скорость роста накопленной площади равна высоте графика. Поэтому $(\ln x)' = \frac1x$, и трещина из начала главы закрыта: $\int \frac{dx}{x} = \ln|x| + C$ (при $x < 0$ производная $\ln(-x)$ по правилу цепочки тоже равна $\frac{-1}{-x} = \frac1x$). Производную экспоненты даёт зеркало.
$(e^x)' = e^x$ при каждом $x$.
Графики взаимно обратных функций симметричны относительно прямой $y = x$. Отражение меняет местами пробег и подъём, поэтому наклон касательной переворачивается.
Вот и объяснение второго совпадения. По определению $a^x = e^{x \ln a}$, и по правилу цепочки $(a^x)' = \ln a \cdot e^{x\ln a} = \ln a \cdot a^x$. Множитель из раздела о склоне — натуральный логарифм основания: $k(a) = \ln a = A(a)$. Отсюда $k(2) = A(2)$ и $k(3) = A(3)$, а наклон в нуле равен единице ровно при $\ln a = 1$, то есть при $a = e$. Все три наблюдения — одно число.
Число $e$ — это одновременно предел $\left(1 + \frac1n\right)^n$, граница, до которой площадь под гиперболой $y = \frac1t$ равна единице, и единственное основание, при котором показательная функция равна своей производной. Натуральный логарифм натурален потому, что это площадь под самой простой кривой, $y = \frac1t$.
Осталось ответить на вопрос из прошлой главы: какие функции совпадают со своей производной? Экспонента — одна из них, и любое её кратное $Ce^x$ тоже. Других нет.
Пусть функция $y(x)$ определена на промежутке, содержащем $0$, и на нём $y' = k y$, где $k$ — число. Тогда $y(x) = y(0)\,e^{kx}$. В частности, функции, равные своей производной, — это в точности функции $Ce^x$.
Хитрость в том, чтобы разделить $y$ на ожидаемый ответ и увидеть, что частное не меняется. Положим $u(x) = y(x)\,e^{-kx}$. По правилу произведения и правилу цепочки
$$u' = y' e^{-kx} + y \cdot \left(-k e^{-kx}\right) = (y' - ky)\,e^{-kx} = 0.$$Функция с нулевой производной на промежутке постоянна — это лемма из главы об интеграле, она следует из теоремы Лагранжа: $u(x) - u(0) = u'(c)\,x = 0$ для некоторого $c$ между $0$ и $x$. Значит, $u(x) = u(0) = y(0)$, и $y(x) = y(0)\,e^{kx}$.
Эта теорема — рабочая лошадка анализа. Всякий раз, когда величина меняется со скоростью, пропорциональной ей самой, ответ — экспонента, и никаких других вариантов нет. Первое применение — к самому банку.
Для любого действительного $x$ выполняется $\left(1 + \frac xn\right)^n \to e^x$ при $n \to \infty$.
При $x = 0$ доказывать нечего. Пусть $x \ne 0$ и $n > |x|$, тогда $h = \frac xn$ лежит между $-1$ и $1$. По лемме о полоске $\frac{h}{1 + h} < \ln(1 + h) < h$; умножим на $n > 0$: $\frac{x}{1 + x/n} < n\ln\left(1 + \frac xn\right) < x$. Левая часть отличается от $x$ на $\frac{x^2/n}{1 + x/n}$, и это стремится к нулю, поэтому $\ln\left(1 + \frac xn\right)^n \to x$. Экспонента меняется без скачков (у неё есть производная), поэтому и $\left(1 + \frac xn\right)^n = e^{\ln(1 + x/n)^n} \to e^x$.
Рубль под $100\,\%$ при непрерывном начислении вырастает за год в $e$ раз, а при ставке $r$ за $t$ лет — в $e^{rt}$ раз. Это и есть та плавная кривая, к которой прижимаются ступеньки в банке.
Вклад $100\,000$ рублей положили под $8\,\%$ годовых с непрерывным начислением на $10$ лет. Сколько рублей будет на счёте? Ответ округлите до рубля; $e^{0{,}8}$ можно посчитать на калькуляторе.
$S = 100\,000 \cdot e^{0{,}08 \cdot 10} = 100\,000 \cdot e^{0{,}8} \approx 100\,000 \cdot 2{,}225541 \approx 222\,554$ рубля. Для сравнения: при начислении раз в год было бы $100\,000 \cdot 1{,}08^{10} \approx 215\,892$ рубля. Разница — почти семь тысяч.
Факториалы: портрет с точностью до знака
Предел $\left(1 + \frac1n\right)^n$ — плохой способ считать $e$: при $n = 1000$ верны лишь два знака после запятой. Есть способ несравнимо быстрее, и дорогу к нему прокладывает интегрирование по частям из прошлой главы.
Для любого действительного $x$ и натурального $n$
$$e^x = 1 + x + \frac{x^2}{2!} + \frac{x^3}{3!} + \ldots + \frac{x^n}{n!} + R_n(x),$$где остаток равен
$$R_n(x) = \frac1{n!}\int_0^x (x - t)^n e^t\,dt,$$причём $|R_n(x)| \le e^{|x|}\dfrac{|x|^{n+1}}{(n+1)!}$. Поэтому суммы $1 + x + \frac{x^2}{2!} + \ldots + \frac{x^n}{n!}$ при $n \to \infty$ стремятся к $e^x$.
Хитрость в том, чтобы вынимать из интеграла по одному слагаемому. При $n = 0$ формула верна: $R_0(x) = \int_0^x e^t\,dt = e^x - 1$ по формуле Ньютона — Лейбница, и $e^x = 1 + R_0(x)$.
Шаг от $n$ к $n + 1$ — интегрирование по частям. В интеграле $R_n(x)$ возьмём $u = e^t$ и $dv = \frac{(x - t)^n}{n!}\,dt$, тогда $du = e^t\,dt$ и $v = -\frac{(x - t)^{n+1}}{(n+1)!}$ (проверьте дифференцированием: минус из правила цепочки съедает минус впереди). Формула $\int u\,dv = uv - \int v\,du$ даёт
$$R_n(x) = \left[-e^t\,\frac{(x - t)^{n+1}}{(n+1)!}\right]_0^x + \frac1{(n+1)!}\int_0^x (x - t)^{n+1} e^t\,dt.$$Оставшийся интеграл по определению равен $R_{n+1}(x)$. В скобке на верхнем пределе $(x - x)^{n+1} = 0$, на нижнем $e^0 = 1$, так что скобка равна $\frac{x^{n+1}}{(n+1)!}$, и $R_n(x) = \frac{x^{n+1}}{(n+1)!} + R_{n+1}(x)$. Значит, если формула верна для $n$, то, подставив в неё $R_n = \frac{x^{n+1}}{(n+1)!} + R_{n+1}$, мы получим её для $n + 1$. По индукции она верна для всех $n$.
Оценка остатка. На отрезке между $0$ и $x$ множитель $e^t$ не больше $e^{|x|}$, а $\left|\int_0^x (x - t)^n\,dt\right| = \frac{|x|^{n+1}}{n+1}$. Поэтому $|R_n(x)| \le \frac{e^{|x|}}{n!} \cdot \frac{|x|^{n+1}}{n+1} = e^{|x|}\frac{|x|^{n+1}}{(n+1)!}$. Факториал обгоняет любую степень: как только $n + 1 > 2|x|$, каждый следующий множитель $\frac{|x|}{n+2}$ меньше половины, и дробь $\frac{|x|^{n+1}}{(n+1)!}$ с каждым шагом уменьшается больше чем вдвое. Остаток стремится к нулю, суммы — к $e^x$.
Для любого другого показателя ряд работает так же быстро: $e^{0{,}1} \approx 1 + 0{,}1 + 0{,}005 + 0{,}000\,167 + 0{,}000\,004 = 1{,}105\,171$, и все шесть знаков после запятой верны.
Сколько первых слагаемых ряда $1 + 1 + \frac1{2!} + \frac1{3!} + \ldots$ нужно сложить, чтобы ошибка стала меньше одной миллионной?
Если сложить слагаемые с $k = 0$ до $k = n$, ошибка равна $R_n(1) < \frac{e}{(n+1)!} < \frac{3}{(n+1)!}$. При $n = 9$ это $\frac{3}{10!} \approx 8{,}3 \cdot 10^{-7}$, меньше миллионной. При $n = 8$ оценки не хватает, и настоящая ошибка тоже велика: $e - 2{,}718\,278\,77\ldots \approx 3{,}1 \cdot 10^{-6}$. Значит, нужны слагаемые с $k = 0$ по $k = 9$ — всего десять.
Число $e$ иррационально: доказательство в четыре строки
Эйлер доказал иррациональность $e$ в 1737 году, через цепную дробь. Короткое доказательство, которое приписывают Фурье (его опубликовали в 1815 году), держится на оценке остатка. Допустим, $e = \frac pq$ с натуральными $p$ и $q$. Умножим ряд на $q!$: число $q!\,e = p \cdot (q - 1)!$ целое, и $q!\left(1 + 1 + \frac1{2!} + \ldots + \frac1{q!}\right)$ тоже целое, потому что $q!$ делится на каждый $k!$ при $k \le q$. Значит, целое и число $q! \, R_q(1)$. Но $0 < q!\,R_q(1) < q! \cdot \frac{e}{(q+1)!} = \frac{e}{q+1} < 1$ при $q \ge 2$: целое число строго между нулём и единицей. Противоречие. Случай $q = 1$ исключён сразу: $e$ лежит между $2$ и $3$ и не целое. В 1873 году Шарль Эрмит доказал больше: $e$ трансцендентно.
Распад: сколько осталось
Следующее место, где видели $e$, — лаборатория. У радиоактивного атома нет часов. Он не стареет: атом, проживший тысячу лет, распадается в следующую секунду с той же вероятностью, что и только что родившийся, и атомы не влияют друг на друга. Поэтому в большом образце за короткое время распадается одна и та же доля атомов: вдвое больше атомов — вдвое больше распадов. Скорость убывания пропорциональна количеству: $N' = -\lambda N$, где $\lambda > 0$ — постоянная распада. По теореме о единственности экспоненты $N(t) = N_0\,e^{-\lambda t}$.
Период полураспада $T_{1/2}$ — время, за которое распадается половина атомов образца. Из $N_0 e^{-\lambda T} = \frac{N_0}2$ следует $\lambda T = \ln 2$: период полураспада не зависит от того, сколько атомов было вначале.
На этой формуле стоит радиоуглеродный анализ, который в конце 1940-х годов придумал американский химик Уиллард Либби (в 1960 году он получил за него Нобелевскую премию по химии). Пока организм жив, доля углерода-14 в нём такая же, как в атмосфере: он постоянно обменивается с ней углеродом. После смерти обмен прекращается, углерод-14 распадается, и часы начинают идти. Измерив, какая доля углерода-14 осталась в кости или угле от костра, можно узнать, сколько лет назад прекратился обмен.
В древней кости осталась $\frac18$ того углерода-14, что был при жизни. Сколько лет назад жило животное? Период полураспада — $5730$ лет.
$\frac18 = \left(\frac12\right)^3$ — прошло три периода полураспада: $3 \cdot 5730 = 17\,190$ лет. Через формулу: $2^{-t/5730} = 2^{-3}$, откуда $t = 17\,190$.
Проверить закон можно без реактора — с игральными кубиками. Бросим горсть кубиков и уберём те, на которых выпала шестёрка: они «распались». Остальные бросим снова. У кубика тоже нет памяти: сколько бы бросков он ни пережил, шестёрка выпадает с вероятностью $\frac16$. За бросок остаётся в среднем $\frac56$ кубиков, после $t$ бросков — $N_0\left(\frac56\right)^t = N_0 e^{-t\ln(6/5)}$. Это распад с постоянной $\lambda = \ln\frac65 \approx 0{,}182$ и периодом полураспада $\frac{\ln 2}{\ln(6/5)} \approx 3{,}8$ броска.
Период полураспада йода-131 — около $8$ суток. Какая часть йода останется через $16$ суток?
За первые $8$ суток распадается половина, за следующие $8$ — половина оставшегося: $\frac12 \cdot \frac12 = \frac14$. Ноль не наступает никогда: экспонента не обращается в ноль. На практике, конечно, после нескольких десятков периодов не остаётся ни одного атома, но это уже потому, что атомов конечное число.
Тем же законом Ньютон в 1701 году описал остывание. Его статья вышла без подписи в журнале Лондонского королевского общества. Горячее тело отдаёт тепло тем быстрее, чем сильнее оно нагрето по сравнению с окружающим воздухом, поэтому разница температур убывает по экспоненте: $T(t) = T_{\text{возд}} + \left(T_0 - T_{\text{возд}}\right)e^{-kt}$.
Чашку кипятка ($100\,°$C) оставили в комнате при $20\,°$C. Через $10$ минут вода остыла до $60\,°$C. Через сколько минут от начала она остынет до $30\,°$C?
Следим за разницей с комнатой: вначале $80°$, через $10$ минут — $40°$. Разница уменьшилась вдвое за $10$ минут, и по закону Ньютона она будет уменьшаться вдвое каждые $10$ минут. $30\,°$C — это разница $10°$, в восемь раз меньше исходной, то есть три «периода полуостывания»: $30$ минут.
Окружность: $e^{i\pi} = -1$
Следующее наблюдение похоже на ошибку наблюдателя: число $e$ видели на окружности, в компании с $\pi$ и мнимой единицей. Что вообще может означать $e$ в мнимой степени? Первая попытка: $e^x$ — это $e$, умноженное само на себя $x$ раз. С мнимым $x$ эта мысль ничего не даёт. Вторая: $e^x = \lim\left(1 + \frac xn\right)^n$. Здесь нужны только сложение, умножение и предел, а их для комплексных чисел мы умеем делать. Подставим $x = i\varphi$ и посмотрим на плоскость. При умножении комплексных чисел модули перемножаются, а углы складываются.
Для любого действительного $\varphi$
$$\lim_{n \to \infty}\left(1 + \frac{i\varphi}{n}\right)^n = \cos\varphi + i\sin\varphi.$$Хитрость в том, что умножение на $1 + \frac{i\varphi}{n}$ — это маленький поворот с маленьким растяжением. Проследим за углом и за длиной отдельно; нарисуем $0 < \varphi \le \pi$ (для отрицательных $\varphi$ всё зеркально).
Поэтому естественно положить, как Эйлер во «Введении в анализ бесконечно малых» (1748):
При $\varphi = \pi$ получается равенство, которое часто называют самым красивым в математике:
$$e^{i\pi} + 1 = 0.$$В нём встретились пять главных констант: $0$ и $1$ из арифметики, $\pi$ из геометрии, $i$ из алгебры и $e$ из анализа. А смысл у него простой: если поворачивать единицу мелкими шагами, как описано в доказательстве, на общий угол $\pi$, она приедет в $-1$.
У формулы Эйлера есть и «скоростное» прочтение. Производная $e^{it}$ по $t$ — это $ie^{it}$: скорость равна положению, повёрнутому на $90°$. Точка, скорость которой всегда перпендикулярна радиусу и равна ему по длине, бежит по окружности радиуса $1$ с единичной скоростью и за время $t$ пробегает дугу длины $t$ — ровно угол $t$ в радианах. Та же функция, что описывает рост вклада, при мнимом показателе описывает вращение. Поэтому экспонента и управляет всем, что колеблется: маятниками, переменным током, звуком.
Гардероб и лотерея
Два последних места, где встречается $e$, — там, где считают шансы. Первое — гардероб. В главе о комбинаторике $n$ гостей сдавали шляпы, а гардеробщик возвращал их наугад, и мы нашли вероятность того, что никто не получил свою шляпу — что получился беспорядок:
$$p_n = \frac{D_n}{n!} = 1 - \frac1{1!} + \frac1{2!} - \frac1{3!} + \ldots + \frac{(-1)^n}{n!}.$$Там же было сказано, что эти вероятности стремятся к $\frac1e \approx 0{,}368$. Теперь это следует из ряда для экспоненты при $x = -1$.
$\left|p_n - \frac1e\right| \le \frac1{(n+1)!}$. Число беспорядков $D_n$ при $n \ge 1$ — ближайшее целое к $\frac{n!}{e}$.
По теореме о ряде при $x = -1$ имеем $e^{-1} = 1 - 1 + \frac1{2!} - \ldots + \frac{(-1)^n}{n!} + R_n(-1) = p_n + R_n(-1)$. Оценка остатка по доказательству теоремы: при $x < 0$ множитель $e^t$ на отрезке $[x;\ 0]$ не больше единицы, поэтому $|R_n(-1)| \le \frac{1}{(n+1)!}$. Это первое утверждение. Умножим его на $n!$: $\left|D_n - \frac{n!}{e}\right| \le \frac{1}{n + 1}$. При $n \ge 2$ это не больше $\frac13$, меньше половины, и целое $D_n$ — ближайшее к $\frac{n!}{e}$ целое число. При $n = 1$: $D_1 = 0$, а $\frac1e \approx 0{,}37$ ближе к нулю, чем к единице.
Проверим: $\frac{4!}{e} \approx 8{,}83$, и $D_4 = 9$; $\frac{5!}{e} \approx 44{,}15$, и $D_5 = 44$; $\frac{13!}{e} \approx 2\,290\,792\,932{,}07$, и число беспорядков из тринадцати карт игры Монмора равно $2\,290\,792\,932$. Уже при шести гостях вероятность отличается от $\frac1e$ меньше чем на $\frac1{7!} \approx 0{,}0002$.
Второе — лотерея. Пусть шанс выиграть на один билет равен $\frac1n$, и вы покупаете $n$ билетов в $n$ разных тиражах. Вероятность ни разу не выиграть равна $\left(1 - \frac1n\right)^n$, а по теореме о непрерывных процентах при $x = -1$ это стремится к $e^{-1}$. Билет с шансом один на миллион, купленный миллион раз, оставляет вас без выигрыша с вероятностью $0{,}367\,879\ldots$ Отличие от $\frac1e$ — в седьмом знаке.
Шанс выиграть в моментальную лотерею — один из ста. Вы купили сто билетов, по одному в каждом тираже. Какова вероятность выиграть хотя бы раз?
Ни разу не выиграть: $0{,}99^{100} \approx 0{,}366$ — почти $\frac1e$. Выиграть хотя бы раз: $1 - 0{,}366 = 0{,}634$. Ту же долю, $1 - \frac1e \approx 63{,}2\,\%$, мы уже видели: столько выигрывает банкомёт в игре «Тринадцать».
Смотрины: правило 37 процентов
Невеста принимает женихов по одному, в случайном порядке. Всего их $n$, и это число известно заранее. Каждого она может сравнить со всеми, кого уже видела, но про будущих ничего не знает. Решать нужно сразу: отказ окончательный, а согласие прекращает смотрины. Цель — выбрать лучшего из всех. Если соглашаться наугад, шанс равен $\frac1n$: при сотне женихов это один процент. Можно ли надеяться на большее?
Попробуем стратегию: первых $m$ женихов пропустить, только присматриваясь, а затем согласиться на первого, кто лучше всех предыдущих. Если такого не найдётся, остаться с последним.
При $1 \le m \le n - 1$ эта стратегия выбирает лучшего с вероятностью
$$P_n(m) = \frac mn\left(\frac1m + \frac1{m + 1} + \ldots + \frac1{n - 1}\right).$$Если $n \to \infty$ и доля пропущенных $\frac mn$ стремится к $x \in (0;\ 1)$, то $P_n(m) \to -x\ln x$. Наибольшее значение $\frac1e$ предел принимает при $x = \frac1e$.
Пусть лучший жених пришёл $k$-м по счёту. Все места для него равновероятны, так что это случается с вероятностью $\frac1n$. Если $k \le m$, мы его пропустили. Пусть $k > m$. Мы выберем его, если никто из пришедших с $(m+1)$-го по $(k-1)$-й не оказался лучше всех своих предшественников, — иначе мы остановились бы раньше. А это то же самое, что лучший из первых $k - 1$ женихов оказался среди первых $m$: тогда все, кто пришёл после него и до $k$-го, хуже его и рекордами не станут; если же лучший из первых $k - 1$ пришёл позже $m$-го, то он сам рекорд, и мы согласимся на него. Порядок первых $k - 1$ женихов случаен, их лучший с равными шансами стоит на любом из $k - 1$ мест, поэтому вероятность равна $\frac{m}{k-1}$. Складывая по всем $k$ от $m + 1$ до $n$, получаем
$$P_n(m) = \sum_{k = m+1}^{n} \frac1n \cdot \frac{m}{k - 1} = \frac mn \sum_{j = m}^{n-1} \frac1j.$$Теперь предел. Перепишем сумму как $\sum_{j=m}^{n-1} \frac{1}{j/n} \cdot \frac1n$. Это сумма площадей прямоугольников шириной $\frac1n$ под графиком $y = \frac1t$ на отрезке от $\frac mn$ до $1$ — сумма Римана. При $n \to \infty$ она стремится к $\int_x^1 \frac{dt}{t} = \ln 1 - \ln x = -\ln x$, а множитель $\frac mn$ — к $x$. Итого $P \to -x\ln x$. Производная этой функции $-\ln x - 1$ равна нулю при $\ln x = -1$, то есть при $x = \frac1e$, слева от этой точки она положительна, справа отрицательна. Значит, это максимум, и он равен $-\frac1e \ln\frac1e = \frac1e$.
Можно доказать, что никакая стратегия не даёт в этой игре больше: правило «пропусти примерно $\frac ne$, потом бери первого, кто лучше всех» лучшее из возможных. Широкой публике задачу представил Мартин Гарднер в 1960 году в журнале Scientific American, как «игру в гугол»: на бумажках пишут числа, и нужно угадать, когда перед тобой наибольшее. Под именами «задача о секретаре» и «задача о разборчивой невесте» она живёт в учебниках теории вероятностей, а её вариантами описывают выбор квартиры, продажу машины и момент, когда пора перестать искать и начать действовать.
Куда дальше
Дневник можно закрывать. В банке, на склоне, под гиперболой, в распаде, на окружности, в гардеробе и на смотринах встречалось одно и то же число, а самый точный его портрет дал ряд из факториалов: $e = 1 + 1 + \frac1{2!} + \frac1{3!} + \ldots$ Больше того, вся экспонента оказалась «бесконечным многочленом»: $e^x = 1 + x + \frac{x^2}{2!} + \frac{x^3}{3!} + \ldots$
А другие функции так можно записать? Подставим в этот ряд $x = i\varphi$ и разделим слагаемые на действительные и мнимые: по формуле Эйлера должны получиться $\cos\varphi = 1 - \frac{\varphi^2}{2!} + \frac{\varphi^4}{4!} - \ldots$ и $\sin\varphi = \varphi - \frac{\varphi^3}{3!} + \ldots$ — и это правда. Логарифм, синус, корень — у многих функций есть свой бесконечный многочлен. Но здесь начинаются сюрпризы. Функция $\frac{1}{1 + x^2}$ гладкая на всей прямой, а её ряд $1 - x^2 + x^4 - x^6 + \ldots$ при $x = 2$ даёт $1 - 4 + 16 - 64 + \ldots$ — полную бессмыслицу. А ряд $1 - \frac12 + \frac13 - \frac14 + \ldots$ можно переставить так, что он даст любую наперёд заданную сумму. Когда сумма бесконечного числа слагаемых вообще имеет смысл? Это вопрос главы о рядах.