В чем измеряется пи. Вычисление N-го знака числа Пи без вычисления предыдущих

История числа Пи начинается еще с Древнего Египта и идет параллельно с развитием всей математики. Мы же впервые встречаемся с этой величиной в стенах школы.

Число Пи является, пожалуй, самым загадочным из бесконечного множества других. Ему посвящены стихи, его изображают художники, о нем даже снят фильм. В нашей статье мы рассмотрим историю развития и вычисления, а также области применения константы Пи в нашей жизни.

Число Пи – это математическая константа равная отношению длины окружности к длине ее диаметра. Первоначально оно называлось лудольфово числом, а обозначать его буквой Пи было предложено британским математиком Джонсом в 1706 году. После работ Леонарда Эйлера в 1737 году это обозначение стало общепринятым.

Число Пи является иррациональным, то есть его значение не может быть точно выражено в виде дроби m/n, где m и n - целые числа. Впервые это доказал Иоганн Ламберт в 1761 году.

История развития числа Пи насчитывает уже порядка 4000 лет. Еще древнеегипетским и вавилонским математикам было известно, что отношение длины окружности к диаметру одинаково для любой окружности и значение его равно чуть больше трех.

Архимед предложил математический способ вычисления Пи, в котором он вписывал в окружность и описывал около неё правильные многоугольники. По его расчетам Пи примерно равнялась 22/7 ≈ 3,142857142857143.

Во II веке Чжан Хэн предложил два значения числа Пи: ≈ 3,1724 и ≈ 3,1622.

Индийские математики Ариабхата и Бхаскара нашли приблизительное значение 3,1416.

Самым точным приближением числа Пи на протяжении 900 лет было вычисление китайского математика Цзу Чунчжи, проведенное в 480-х годах. Он вывел, что Пи ≈ 355 / 113 , и показал, что 3,1415926 < Пи < 3,1415927.

До II тысячелетия было вычислено не более 10 цифр числа Пи. Лишь с развитием математического анализа, а особенно с открытием рядов, были осуществлены последующие крупные продвижения в вычислении константы.

В 1400-х годах Мадхава смог вычислить Пи=3,14159265359. Его рекорд удалось побить персидскому математику Аль-Каши в 1424 году. Он в своём труде «Трактат об окружности» привёл 17 цифр числа Пи, 16 из которых оказались верными.

Голландский математик Людольф ван Цейлен дошел в своих вычислениях до 20-ти чисел, отдав на это 10 лет жизни. После его смерти в его записях были обнаружены еще 15 цифр числа Пи. Он завещал, чтобы эти цифры были высечены на его надгробии.

С появлением компьютеров число Пи на сегодняшний день насчитывает несколько триллионов знаков и это не предел. Но, как подмечено в книге «Fractals for the Classroom», при всей важности числа Пи «трудно найти сферы в научных расчетах, где потребовалось бы больше двадцати десятичных знаков».

В нашей жизни число Пи используется во многих научных областях. Физика, электроника, теория вероятностей, химия, строительство, навигация, фармакология - это лишь некоторые из них, которые просто невозможно представить себе без этого загадочного числа.

Хотите знать и уметь, больше и сами?

Мы предлагаем Вам обучение по направлениям: компьютеры, программы, администрирование, сервера, сети, сайтостроение, SEO и другое. Узнайте подробности сейчас!

По материалам сайта Calculator888.ru - Число Пи - значение, история, кто придумал .

Введение

В статье присутствуют математические формулы, поэтому для чтения перейдите на сайт для их корректного отображения. Число \(\pi \) имеет богатую историю. Данная константа обозначает отношение длины окружности к ее диаметру.

В науке число \(\pi \) используют в любых расчетах, где есть окружности. Начиная от объема банки газировки, до орбит спутников. И не только окружности. Ведь в изучении кривых линий число \(\pi \) помогает понять периодические и колебательные системы. Например, электромагнитные волны и даже музыку.

В 1706 году в книге «Новое введение в математику» британского ученого Уильяма Джонса (1675-1749 гг.) для обозначения числа 3,141592… впервые была использована буква греческого алфавита \(\pi \). Это обозначение происходит от начальной буквы греческих слов περιϕερεια — окружность, периферия и περιµετρoς — периметр. Общепринятым обозначение стало после работ Леонарда Эйлера в 1737 году.

Геометрический период

Постоянство отношения длины любой окружности к её диаметру было замечено уже давно. Жители Междуречья применяли довольно грубое приближение числа \(\pi \). Как следует из древних задач, в своих расчетах они используют значение \(\pi ≈ 3 \).

Более точное значение для \(\pi \) использовали древние египтяне. В Лондоне и Нью-Йорке хранятся две части древнеегипетского папируса, который называют «папирус Ринда». Папирус был составлен писцом Армесом примерно между 2000-1700 гг. до н.э.. Армес в своем папирусе написал, что площадь круга с радиусом \(r\) равна площади квадрата со стороной, равной \(\frac{8}{9} \) от диаметра окружности \(\frac{8}{9} \cdot 2r \), то есть \(\frac{256}{81} \cdot r^2 = \pi r^2 \). Отсюда \(\pi = 3,16\).

Древнегреческий математик Архимед (287-212 гг. до н.э.) впервые поставил задачу измерения круга на научную почву. Он получил оценку \(3\frac{10}{71} < \pi < 3\frac{1}{7}\), рассмотрев отношение периметров вписанного и описанного 96-угольника к диаметру окружности. Архимед выразил приближение числа \(\pi \) в виде дроби \(\frac{22}{7}\), которое до сих называется архимедовым числом.

Метод достаточно простой, но при отсутствии готовых таблиц тригонометрических функций потребуется извлечение корней. Кроме этого, приближение сходится к \(\pi \) очень медленно: с каждой итерацией погрешность уменьшается лишь вчетверо.

Аналитический период

Несмотря на это, до середины 17 века все попытки европейских учёных вычислить число \(\pi \) сводились к увеличению сторон многоугольника. Так например, голландский математик Лудольф ван Цейлен (1540-1610 гг.) вычислил приближенное значение числа \(\pi \) с точностью до 20-ти десятичных цифр.

На вычисление ему понадобилось 10 лет. Удваивая по методу Архимеда число сторон вписанных и описанных многоугольников, он дошел до \(60 \cdot 2^{29} \) — угольника с целью вычисления \(\pi \) с 20 десятичными знаками.

После смерти в его рукописях были обнаружены ещё 15 точных цифр числа \(\pi \). Лудольф завещал, чтобы найденные им знаки были высечены на его надгробном камне. В честь него число \(\pi \) иногда называли «лудольфовым числом» или «константой Лудольфа».

Одним из первых, кто представил метод, отличный от метода Архимеда, был Франсуа Виет (1540-1603 гг.). Он пришел к результату , что круг, диаметр которого равен единице, имеет площадь:

\[\frac{1}{2 \sqrt{\frac{1}{2}} \cdot \sqrt{\frac{1}{2} + \frac{1}{2} \sqrt{\frac{1}{2}} } \cdot \sqrt{\frac{1}{2} + \frac{1}{2} \sqrt{\frac{1}{2} + \frac{1}{2} \sqrt{\frac{1}{2} \cdots }}}} \]

С другой стороны, площадь равна \(\frac{\pi}{4} \). Подставив и упростив выражение, можно получить следующую формулу бесконечного произведения для вычисления приближенного значения \(\frac{\pi}{2} \):

\[\frac{\pi}{2} = \frac{2}{\sqrt{2}} \cdot \frac{2}{\sqrt{2 + \sqrt{2}}} \cdot \frac{2}{\sqrt{2+ \sqrt{2 + \sqrt{2}}}} \cdots \]

Полученная формула представляет собой первое точное аналитическое выражение для числа \(\pi \). Кроме этой формулы, Виет, используя метод Архимеда, дал с помощью вписанных и описанных многоугольников, начиная с 6-угольника и заканчивая многоугольником с \(2^{16} \cdot 6 \) сторонами приближение числа \(\pi \) с 9 правильными знаками.

Английский математик Уильям Броункер (1620-1684 гг.), используя цепную дробь , получил следующие результаты вычисления \(\frac{\pi}{4}\):

\[\frac{4}{\pi} = 1 + \frac{1^2}{2 + \frac{3^2}{2 + \frac{5^2}{2 + \frac{7^2}{2 + \frac{9^2}{2 + \frac{11^2}{2 + \cdots }}}}}} \]

Данный метод вычисления приближения числа \(\frac{4}{\pi} \) требует довольно больших вычислений, чтобы получить хотя бы небольшое приближение.

Получаемые в результате подстановки значения то больше, то меньше числа \(\pi \), и каждый раз все ближе к истинному значению, но для получения значения 3,141592 потребуется совершить довольно большие вычисления.

Другой английский математик Джон Мэчин (1686-1751 гг.) в 1706 году для вычисления числа \(\pi \) со 100 десятичными знаками воспользовался формулой, выведенной Лейбницем в 1673 году, и применил её следующим образом:

\[\frac{\pi}{4} = 4 arctg\frac{1}{5} — arctg\frac{1}{239} \]

Ряд быстро сходится и с его помощью можно вычислить число \(\pi \) с большой точностью. Формулы подобного типа использовались для установки нескольких рекордов в эпоху компьютеров.

В XVII в. с началом периода математики переменной величины наступил новый этап в вычислении \(\pi \). Немецкий математик Готфрид Вильгельм Лейбниц (1646-1716 гг.) в 1673 году нашел разложение числа \(\pi \), в общем виде его можно записать следующим бесконечным рядом:

\[ \pi = 1 — 4(\frac{1}{3} + \frac{1}{5} — \frac{1}{7} + \frac{1}{9} — \frac{1}{11} + \cdots) \]

Ряд получается при подстановке x = 1 в \(arctg x = x — \frac{x^3}{3} + \frac{x^5}{5} — \frac{x^7}{7} + \frac{x^9}{9} — \cdots\)

Леонард Эйлер развивает идею Лейбница в своих работах, посвященных использованию рядов для arctg x при вычислении числа \(\pi \). В трактате «De variis modis circuli quadraturam numeris proxime exprimendi» (О различных методах выражения квадратуры круга приближенными числами), написанном в 1738 году, рассматриваются методы усовершенствования вычислений по формуле Лейбница.

Эйлер пишет о том, что ряд для арктангенса будет сходиться быстрее, если аргумент будет стремиться к нулю. Для \(x = 1\) сходимость ряда очень медленная: для вычисления с точностью до 100 цифр необходимо сложить \(10^{50}\) членов ряда. Ускорить вычисления можно, уменьшив значение аргумента. Если принять \(x = \frac{\sqrt{3}}{3}\), то получается ряд

\[ \frac{\pi}{6} = artctg\frac{\sqrt{3}}{3} = \frac{\sqrt{3}}{3}(1 — \frac{1}{3 \cdot 3} + \frac{1}{5 \cdot 3^2} — \frac{1}{7 \cdot 3^3} + \cdots) \]

По утверждению Эйлера, если мы возьмем 210 членов этого ряда, то получим 100 верных знаков числа. Полученный ряд неудобен, потому что необходимо знать достаточно точное значение иррационального числа \(\sqrt{3} \). Также Эйлер в своих вычислениях использовал разложения арктангенсов на сумму арктангенсов меньших аргументов :

\[где x = n + \frac{n^2-1}{m-n}, y = m + p, z = m + \frac{m^2+1}{p} \]

Далеко не все формулы для вычисления \(\pi \), которые использовал Эйлер в своих записных книжках, были опубликованы. В опубликованных работах и записных книжках он рассмотрел 3 различных ряда для вычисления арктангенса, а также привел множество утверждений, касающихся количества суммируемых членов, необходимых для получения приближенного значения \(\pi \) c заданной точностью.

В последующие годы уточнения значения числа \(\pi \) происходили все быстрее и быстрее. Так, например, в 1794 году Георг Вега (1754-1802 гг.) определил уже 140 знаков , из который только 136 оказались верными.

Период компьютерных вычислений

XX век ознаменован совершенно новым этапом в вычислении числа \(\pi \). Индийский математик Сриниваса Рамануджан (1887-1920 гг.) обнаружил множество новых формул для \(\pi \). В 1910 году он получил формулу для вычисления \(\pi \) через разложение арктангенса в ряд Тейлора:

\[\pi = \frac{9801}{2\sqrt{2} \sum\limits_{k=1}^{\infty} \frac{(1103+26390k) \cdot (4k)!}{(4\cdot99)^{4k} (k!)^2}} .\]

При k=100 достигается точность в 600 верных цифр числа \(\pi \).

Появление ЭВМ позволило существенно увеличить точность получаемых значений за более короткие сроки. В 1949 году всего за 70 часов с помощью ENIAC группа ученых под руководством Джона фон Неймана (1903-1957 гг.) получила 2037 знаков после запятой числа \(\pi \) . Давид и Грегорий Чудновские в 1987 году получили формулу, с помощью которой смогли установить несколько рекордов в вычислении \(\pi \):

\[\frac{1}{\pi} = \frac{1}{426880\sqrt{10005}} \sum\limits_{k=1}^{\infty} \frac{(6k)!(13591409+545140134k)}{(3k)!(k!)^3(-640320)^{3k}}.\]

Каждый член ряда дает по 14 цифр. В 1989 году было получено 1 011 196 691 цифр после запятой. Данная формула хорошо подходит для вычисления \(\pi \) на персональных компьютерах. На данный момент братья являются профессорами в политехническом институте Нью-Йоркского университета.

Важным событием недавнего времени стало открытие формулы в 1997 году Саймоном Плаффом . Она позволяет извлечь любую шестнадцатеричную цифру числа \(\pi \) без вычисления предыдущих. Формула носит название «Формула Бэйли — Боруэйна — Плаффа» в честь авторов статьи, где формула была впервые опубликована. Она имеет следующий вид:

\[\pi = \sum\limits_{k=1}^{\infty} \frac{1}{16^k} (\frac{4}{8k+1} — \frac{2}{8k+4} — \frac{1}{8k+5} — \frac{1}{8k+6}) .\]

В 2006 году Саймон, используя PSLQ, получил несколько красивых формул для вычисления \(\pi \). Например,

\[ \frac{\pi}{24} = \sum\limits_{n=1}^{\infty} \frac{1}{n} (\frac{3}{q^n — 1} — \frac{4}{q^{2n} -1} + \frac{1}{q^{4n} -1}), \]

\[ \frac{\pi^3}{180} = \sum\limits_{n=1}^{\infty} \frac{1}{n^3} (\frac{4}{q^{2n} — 1} — \frac{5}{q^{2n} -1} + \frac{1}{q^{4n} -1}), \]

где \(q = e^{\pi}\). В 2009 году японские ученые, используя суперкомпьютер T2K Tsukuba System, получили число \(\pi \) c 2 576 980 377 524 десятичными знаками после запятой. Вычисления заняли 73 часа 36 минут. Компьютер был оснащен 640-ка четырех ядерными процессорами AMD Opteron, что обеспечило производительность в 95 триллионов операций в секунду.

Следующее достижение в вычислении \(\pi \) принадлежит французскому программисту Фабрису Беллару , который в конце 2009 года на своем персональном компьютере под управлением Fedora 10 установил рекорд, вычислив 2 699 999 990 000 знаков после запятой числа \(\pi \). За последние 14 лет это первый мировой рекорд, который поставлен без использования суперкомпьютера. Для высокой производительности Фабрис использовал формулу братьев Чудновских. В общей сложности вычисление заняло 131 день (103 дня расчеты и 13 дней проверка результата). Достижение Беллара показало, что для таких вычислений не обязательно иметь суперкомпьютер.

Всего через полгода рекорд Франсуа был побит инженерами Александром Йи и Сингеру Кондо. Для установления рекорда в 5 триллионов знаков после запятой числа \(\pi \) был также использован персональный компьютер, но уже с более внушительными характеристиками: два процессора Intel Xeon X5680 по 3,33 ГГц, 96 ГБ оперативной памяти, 38 ТБ дисковой памяти и операционная система Windows Server 2008 R2 Enterprise x64. Для вычислений Александр и Сингеру использовали формулу братьев Чудновских. Процесс вычисления занял 90 дней и 22 ТБ дискового пространства. В 2011 году они установили еще один рекорд , вычислив 10 триллионов десятичных знаков числа \(\pi \). Вычисления происходили на том же компьютере, на котором был поставлен их предыдущий рекорд и занял в общей сложности 371 день. В конце 2013 года Александр и Сингеру улучшили рекорд до 12,1 триллиона цифр числа \(\pi \), вычисление которых заняло у них всего 94 дня. Такое улучшение в производительности достигнуто благодаря оптимизации производительности программного обеспечения, увеличения количества ядер процессора и значительного улучшения отказоустойчивости ПО.

Текущим рекордом является рекорд Александра Йи и Сингеру Кондо, который составляет 12,1 триллиона цифр после запятой числа \(\pi \).

Таким образом, мы рассмотрели методы вычисления числа \(\pi \), используемые в древние времена, аналитические методы, а также рассмотрели современные методы и рекорды по вычислению числа \(\pi \) на компьютерах.

Список источников

  1. Жуков А.В. Вездесущее число Пи – М.:Изд-во ЛКИ, 2007 – 216 с.
  2. Ф.Рудио. О квадратуре круга, с приложением истории вопроса, составленной Ф.Рудио. / Рудио Ф. – М.: ОНТИ НКТП СССР, 1936. – 235c.
  3. Arndt, J. Pi Unleashed / J. Arndt, C. Haenel. – Springer, 2001. – 270p.
  4. Шухман, Е.В. Приближенное вычисление числа Пи с помощью ряда для arctg x в опубликованных и неопубликованных работах Леонарда Эйлера / Е.В. Шухман. — История науки и техники, 2008 – №4. – С. 2-17.
  5. Euler, L. De variis modis circuli quadraturam numeris proxime exprimendi/ Commentarii academiae scientiarum Petropolitanae. 1744 – Vol.9 – 222-236p.
  6. Шумихин, С. Число Пи. История длиною в 4000 лет / С. Шумихин, А. Шумихина. — М.: Эксмо, 2011. — 192с.
  7. Борвейн, Дж.М. Рамануджан и число Пи. / Борвейн, Дж.М., Борвейн П.Б. В мире науки. 1988 – №4. – С. 58-66.
  8. Alex Yee. Number world. Access mode: numberworld.org

Понравилось?

Расскажи

Уже много веков и даже, как ни странно, тысячелетий люди понимают важность и ценность для науки математической постоянной, равной отношению длины окружности к ее же диаметру. число Пи, до сих пор неизвестно, но к нему имели отношение самые лучшие математики на протяжении всей нашей истории. Большинство из них хотели выразить его рациональным числом.

1. Исследователи и истинные поклонники числа Пи организовали клуб, для вступления в который требуется знать наизусть достаточно большое количество его знаков.

2. С 1988 года празднуется «День числа Пи», который приходится на 14 марта. Готовят салаты, торты, печенья, пирожные с его изображением.

3. Число Пи уже переложили на музыку, при этом оно весьма неплохо звучит. Ему даже воздвигли памятник в американском Сиэтле перед зданием городского Музея искусств.

В то далекое время число Пи старались вычислить при помощи геометрии. То, что это число постоянно для самых разных окружностей, знали еще геометры в Древнем Египте, Вавилоне, Индии и Древней Греции, утверждавшие в своих работах, что оно всего лишь немного больше трех.

В одной из священных книг джайнизма (древняя индийская религия, которая возникла в VI в. до н. э.) упоминается, что тогда число Пи считалось равным корню квадратному из десяти, что в итоге дает 3,162... .

Древнегреческие математики проводили измерение окружности методом построения отрезка, а вот для того, чтобы измерить круг, им приходилось строить равновеликий квадрат, то есть фигуру, равную ему по площади.

Когда еще не знали десятичных дробей, великий Архимед нашел значение числа Пи с точностью 99,9%. Он открыл способ, который стал основой многих последующих вычислений, вписывал в окружность и описывал вокруг нее правильные многоугольники. В результате Архимед рассчитал значение числа Пи как отношение 22 / 7 ≈ 3,142857142857143.

В Китае, математик и придворный астроном, Цзу Чунчжи в V веке до н. э. обозначил более точное значение числа Пи, рассчитав его до семи цифр после запятой и определил его значение между числами 3, 1415926 и 3,1415927. Более 900 лет понадобилось ученым, чтобы продолжить дальше этот цифровой ряд.

Средние века

Известный индийский ученый Мадхава, который жил на рубеже XIV - XV веков, ставший основателем Керальской школы астрономии и математики, впервые в истории стал работать над разложением тригонометрических функций в ряды. Правда, сохранились всего лишь два его труда, а на другие известны лишь ссылки и цитаты его учеников. В научном трактате «Махаджьянаяна», который приписывают Мадхаве, указано, что число Пи равно 3,14159265359. А в трактате «Садратнамала» приведено число с еще большим количеством точных знаков после запятой: 3,14159265358979324. В указанных числах последние цифры не соответствуют правильному значению.

В XV веке самаркандский математик и астроном Ал-Каши вычислил число Пи с шестнадцатью знаками после запятой. Его результат считался наиболее точным в течение последующих 250 лет.

У. Джонсон, математик из Англии, одним из первых смог обозначить отношение длины окружности к ее диаметру буквой π. Пи - это первая буква греческого слова «περιφέρεια» - окружность. Но этому обозначению удалось стать общепринятым лишь после того, как им воспользовался в 1736 году более известный ученый Л. Эйлер.

Заключение

Современные ученые продолжают работать над дальнейшими вычислениями значений числа Пи. Для этого уже используют суперкомпьютеры. В 2011 г. ученый из Сигэру Кондо, сотрудничая с американским студентом Александром Йи, произвели правильный расчет последовательности из 10 триллионов цифр. Но до сих пор так и неясно, кто открыл число Пи, кто впервые задумался над этой проблемой и произвел первые расчеты этого, по-настоящему мистического числа.

Отношение длины окружности к ее диаметру одно и то же для всех окружностей. Это отношение принято обозначать греческой буквой (“пи” - начальная буква греческого слова , которое и означало “окружность”).

Архимед в сочинении “Измерение круга” вычислил отношение длины окружности к диаметру (число ) и нашел, что оно заключено между 3 10/71 и 3 1/7.

Долгое время в качестве приближенного значения использовали число 22/7, хотя уже в V веке в Китае было найдено приближение 355/113 = 3,1415929..., которое было открыто вновь в Европе лишь в XVI веке.

В Древней Индии считали равным = 3,1622….

Французский математик Ф. Виет вычислил в 1579 г. с 9 знаками.

Голландский математик Лудольф Ван Цейлен в 1596 г. публикует результат своего десятилетнего труда – число , вычисленное с 32 знаками.

Но все эти уточнения значения числа производились методами, указанными еще Архимедом: окружность заменялась многоугольником со все большим числом сторон. Периметр вписанного многоугольника при этом был меньше длины окружности, а периметр описанного многоугольника – больше. Но при этом оставалась неясным, является ли число рациональным, т. е. отношением двух целых чисел, или иррациональным.

Лишь в 1767 г. немецкий математик И.Г. Ламберт доказал, что число иррационально.

А еще через сто с лишним лет в 1882 г. другой немецкий математик – Ф. Линдеман доказал его трансцендентность, что означало и невозможность построения при помощи циркуля и линейки квадрата, равновеликого данному кругу.

Простейшее измерение

Начертим на плотном картоне окружность диаметра d (=15 см) , вырежем получившийся круг и обмотаем вокруг него тонкую нить. Измерив длину l (=46,5 см) одного полного оборота нити, разделим l на длину диаметра d окружности. Получившееся частное будет приближенным значением числа , т. е. = l / d = 46,5 см / 15 см = 3,1 . Данный довольно грубый способ дает в обычных условиях приближенное значение числа с точностью до 1.

Измерение с помощью взвешивания

На листе картона начертим квадрат. Впишем в него круг. Вырежем квадрат. Определим массу картонного квадрата с помощью школьных весов. Вырежем из квадрата круг. Взвесим и его. Зная массы квадрата m кв (=10 г) и вписанного в него круга m кр (=7,8 г) воспользуемся формулами

где p и h –соответственно плотность и толщина картона, S – площадь фигуры. Рассмотрим равенства:

Естественно, что в данном случае приближенное значение зависит от точности взвешивания. Если взвешиваемые картонные фигуры будут довольно большими, то возможно даже на обычных весах получить такие значения масс, которые обеспечат приближение числа с точностью до 0,1.

Суммирование площадей прямоугольников, вписанных в полукруг

Рисунок 1

Пусть А (a; 0), В (b; 0). Опишем на АВ полуокружность как на диаметре. Разделим отрезок АВ на n равных частей точками x 1 , x 2 , ..., x n-1 и восстановим из них перпендикуляры до пересечения с полуокружностью. Длина каждого такого перпендикуляра – это значение функции f(x)= . Из рисунка 1 ясно, что площадь S полукруга можно вычислить по формуле

S = (b – a) ((f(x 0) + f(x 1) + … + f(x n-1)) / n.

В нашем случае b=1, a=-1 . Тогда = 2 S .

Значения будут тем точнее, чем больше точек деления будет на отрезке АВ. Облегчить однообразную вычислительную работу поможет компьютер, для которого ниже приводится программа 1, составленная на Бейсике.

Программа 1

REM "Вычисление пи"
REM "Метод прямоугольников"
INPUT "Введите число прямоугольников", n
dx = 1 / n
FOR i = 0 TO n - 1
f = SQR(1 - x ^ 2)
x = x + dx
a = a + f
NEXT i
p = 4 * dx * a
PRINT "Значение пи равно ", p
END

Программа была набрана и запущена при различных значениях параметра n . Полученные значения числа записаны в таблице:

Метод Монте-Карло

Это фактически метод статистических испытаний. Свое экзотическое название он получил от города Монте-Карло в княжестве Монако, знаменитого своими игорными домами. Дело в том, что метод требует применения случайных чисел, а одним из простейших приборов, генерирующих случайные числа, может служить рулетка. Впрочем, можно получить случайные числа и при помощи …дождя.

Для опыта приготовим кусок картона, нарисуем на нем квадрат и впишем в квадрат четверть круга. Если такой чертеж некоторое время подержать под дождем, то на его поверхности останутся следы капель. Подсчитаем число следов внутри квадрата и внутри четверти круга. Очевидно, что их отношение будет приближенно равно отношению площадей этих фигур, так как попадание капель в различные места чертежа равновероятно. Пусть N кр – число капель в круге, N кв – число капель в квадрате, тогда

4 N кр / N кв.

Рисунок 2

Дождь можно заменить таблицей случайных чисел, которая составляется с помощью компьютера по специальной программе. Каждому следу капли поставим в соответствие два случайных числа, характеризующих его положение вдоль осей Ох и Оу . Случайные числа можно выбрать из таблицы в любом порядке, например, подряд. Пусть первое четырехзначное число в таблице 3265 . Из него можно приготовить пару чисел, каждое из которых больше нуля и меньше единицы: х=0,32, у=0,65 . Эти числа будем считать координатами капли, т. е. капля как будто попала в точку (0,32; 0,65). Аналогично поступаем и со всеми выбранными случайными числами. Если окажется, что для точки (х; у) выполняется неравенство, то, значит, она лежит вне круга. Если х + у = 1 , то точка лежит внутри круга.

Для подсчета значения снова воспользуемся формулой (1). Ошибка вычислений по этому методу, как правило, пропорциональна , где D – некоторая постоянная, а N –число испытаний. В нашем случае N = N кв. Из этой формулы видно: для того чтобы уменьшить ошибку в 10 раз (иначе говоря, чтобы получить в ответе еще один верный десятичный знак), нужно увеличить N, т. е. объем работы, в 100 раз. Ясно, что применение метода Монте-Карло стало возможным только благодаря компьютерам. Программа 2 реализует на компьютере описанный метод.

Программа 2

REM "Вычисление пи"
REM "Метод Монте-Карло "
INPUT "Введите число капель ", n
m = 0
FOR i = 1 TO n
t = INT(RND(1) * 10000)
x = INT(t \ 100)
y = t - x * 100
IF x ^ 2 + y ^ 2 < 10000 THEN m = m + 1
NEXT i
p = 4 * m / n

END

Программа была набрана и запущена при различных значениях параметра n. Полученные значения числа записаны в таблице:

n
n

Метод “падающей иголки”

Возьмем обыкновенную швейную иголку и лист бумаги. На листе проведем несколько параллельных прямых так, чтобы расстояния между ними были равны и превышали длину иголки. Чертеж должен быть достаточно большим, чтобы случайно брошенная игла не упала за его пределами. Введем обозначения: а - расстояние между прямыми, l – длина иглы.

Рисунок 3

Положение случайным образом брошенной на чертеж иглы (см. рис. 3) определяется расстоянием Х от ее середины до ближайшей прямой и углом j , которой игла образует с перпендикуляром, опущенным из середины иглы на ближайшую прямую (см. рис. 4). Ясно, что

Рисунок 4

На рис. 5 изобразим графически функцию y=0,5 cos . Всевозможные расположения иглы характеризуются точками с координатами (; у ) , расположенными на участке ABCD. Закрашенный участок AED – это точки, которые соответствуют случаю пересечения иглы с прямой. Вероятность события a – “игла пересекла прямую” – вычисляется по формуле:

Рисунок 5

Вероятность p(a) можно приблизительно определить многократным бросанием иглы. Пусть иглу бросали на чертеж c раз и p раз она упала, пересекая одну из прямых, тогда при достаточно большом c имеем p(a) = p / c . Отсюда = 2 l с / a k.

Замечание. Изложенный метод представляет собой вариацию метода статистических испытаний. Он интересен с дидактической точки зрения, так как помогает совместить простой опыт с составлением довольно сложной математической модели.

Вычисление с помощью ряда Тейлора

Обратимся к рассмотрению произвольной функции f(х). Предположим, что для нее в точке x 0 существуют производные всех порядков до n -го включительно. Тогда для функции f(х) можно записать ряд Тейлора:

Вычисления с помощью этого ряда будут тем точнее, чем больше членов ряда будет задействовано. Реализовать данный способ, конечно, лучше всего на компьютере, для чего можно воспользоваться программой 3.

Программа 3

REM "Вычисление пи"
REM "Разложение в ряд Тейлора "
INPUT n
a = 1
FOR i = 1 TO n
d = 1 / (i + 2)
f = (-1) ^ i * d
a = a + f
NEXT i
p = 4 * a
PRINT "значение пи равно"; p
END

Программа была набрана и запущена при различных значениях параметра n . Полученные значения числа записаны в таблице:

Есть очень простые мнемонические правила для запоминания значения числа :

ЧИСЛО p – отношение длины окружности к ее диаметру, – величина постоянная и не зависит от размеров окружности. Число, выражающее это отношение, принято обозначать греческой буквой 241 (от «perijereia » – окружность, периферия). Это обозначение стало употребительным после работы Леонарда Эйлера , относящейся к 1736, однако впервые оно было употреблено Уильямом Джонсом (1675–1749) в 1706. Как и всякое иррациональное число, оно представляется бесконечной непериодической десятичной дробью:

p = 3,141592653589793238462643… Нужды практических расчетов, относящихся к окружностям и круглым телам, заставили уже в глубокой древности искать для 241 приближений с помощью рациональных чисел. Сведения о том, что окружность ровно втрое длиннее диаметра, находятся в клинописных табличках Древнего Междуречья. Такое же значение числа p есть и в тексте Библии: «И сделал литое из меди море, – от края до края его десять локтей, – совсем круглое, вышиною в пять локтей, и снурок в тридцать локтей обнимал его кругом» (3 Цар. 7. 23). Так же считали и древние китайцы. Но уже во 2 тыс. до н.э. древние египтяне пользовались более точным значением числа 241, которое получается из формулы для площади круга диаметра d :

Этому правилу из 50-й задачи папируса Райнда соответствует значение 4(8/9) 2 » 3,1605. Папирус Райнда, найденный в 1858, назван так по имени его первого владельца, его переписал писец Ахмес около 1650 до н.э., автор же оригинала неизвестен, установлено только, что текст создавался во второй половине 19 в. до н.э. Хотя каким образом египтяне получили саму формулу, из контекста неясно. В так называемом Московском папирусе, который был переписан неким учеником между 1800 и 1600 до н.э. с более древнего текста, примерно 1900 до н.э., есть еще одна интересная задача о вычислении поверхности корзины «с отверстием 4½». Неизвестно, какой формы была корзина, но все исследователи сходятся во мнении, что и здесь для числа p берется то же самое приближенное значение 4(8/9) 2 .

Чтобы понять, каким образом древние ученые получили тот или иной результат, нужно попытаться решить задачу, используя только знания и приемы вычислений того времени. Именно так поступают исследователи старинных текстов, однако решения, которые им удается найти, вовсе не обязательно «те самые». Очень часто для одной задачи предлагается несколько вариантов решения, каждый может выбрать себе по вкусу, однако никто не может утверждать, что именно им пользовались в древности. Относительно площади круга кажется правдоподобной гипотеза А.Е.Раик, автора многочисленных книг по истории математики: площадь круга диаметра d сравнивается с площадью описанного вокруг него квадрата, из которого по очереди удаляются малые квадраты со сторонами и (рис. 1). В наших обозначениях вычисления будут выглядеть так: в первом приближении площадь круга S равна разности между площадью квадрата со стороной d и суммарной площадью четырех малых квадратов А со стороной d :

В пользу этой гипотезы свидетельствуют аналогичные вычисления в одной из задач Московского папируса, где предлагается сосчитать

С 6 в. до н.э. математика стремительно развивалась в Древней Греции. Именно древнегреческие геометры строго доказали, что длина окружности пропорциональна ее диаметру (l = 2 p R ; R – радиус окружности, l – ее длина), а площадь круга равна половине произведения длины окружности и радиуса:

S = ½ l R = p R 2 .

Эти доказательства приписывают Евдоксу Книдскому иАрхимеду .

В 3 в. до н.э. Архимед в сочинении Об измерении круга вычислил периметры вписанных в окружность и описанных около нее правильных многоугольников (рис. 2) – от 6- до 96-угольника. Таким образом он установил, что число p находится между 3 10/71 и 3 1/7, т.е. 3,14084 < p < 3,14285. Последнее значение до сих пор используется при расчетах, не требующих особой точности. Более точное приближение 3 17/120 (p » 3,14166) нашел знаменитый астроном, создатель тригонометрии Клавдий Птолемей (2 в.), но оно не вошло в употребление.

Индийцы и арабы полагали, что p = . Это значение приводит так же и индийский математик Брахмагупта (598 – ок. 660). В Китае ученые в 3 в. использовали значение 3 7/50, которое хуже приближения Архимеда, но во второй половине 5 в. Цзу Чун Чжи (ок. 430 – ок. 501) получил для p приближение 355/113 (p » 3,1415927). Оно осталось неизвестно европейцам и было вновь найдено нидерландским математиком Адрианом Антонисом только в 1585. Это приближение дает ошибку лишь в седьмом десятичном знаке.

Поиски более точного приближения p продолжались и в дальнейшем. Например, аль-Каши (первая половина 15 в.) в Трактате об окружности (1427) вычислил 17 десятичных знаков p . В Европе такое же значение было найдено в 1597 году. Для этого ему пришлось вычислять сторону правильного 800 335 168-угольника. Нидерландский ученый Лудольф Ван Цейлен (1540–1610) нашел для него 32 правильных десятичных знака (опубликовано посмертно в 1615), это приближение называется лудольфовым числом.

Число p появляется не только при решении геометрических задач. Со времени Ф.Виета (1540–1603) разыскание пределов некоторых арифметических последовательностей, составляемых по простым законам, приводило к тому же числу p . В связи с этим в определении числа p принимали участие почти все известные математики: Ф.Виет, Х.Гюйгенс , Дж.Валлис, Г.В.Лейбниц , Л.Эйлер . Они получали различные выражения для 241 в виде бесконечного произведения, суммы ряда, бесконечной дроби.

Например, в 1593 Ф.Виет (1540–1603) вывел формулу

В 1658 англичанин Уильям Броункер (1620–1684) нашел представление числа p в виде бесконечной непрерывной дроби

однако неизвестно, как он пришел к этому результату.

В 1665 Джон Валлис (1616–1703) доказал, что

Эта формула носит его имя. Для практического нахождения числа 241 она мало пригодна, но полезна в различных теоретических рассуждениях. В историю науки она вошла как один из первых примеров бесконечных произведений.

Готфрид Вильгельм Лейбниц (1646–1716) в 1673 установил следующую формулу:

выражающую число p /4 как сумму ряда. Однако этот ряд сходится очень медленно. Чтобы вычислить p с точностью до десяти знаков, потребовалось бы, как показал Исаак Ньютон, найти сумму 5 млрд чисел и затратить на это около тысячи лет непрерывной работы.

Лондонский математик Джон Мэчин (1680–1751) в 1706, применяя формулу

получил выражение

которая до сих пор считается одной из лучших для приближенного вычисления p . Чтобы найти те же десять точных десятичных знаков, потребуется всего несколько часов ручного счета. Сам Джон Мэчин вычислил p со 100 верными знаками.

C помощью того же ряда для arctg x и формулы

значение числа p было получено на ЭВМ с точностью до ста тысяч десятичных знаков. Такого рода вычисления представляют интерес в связи с понятием случайных и псевдослучайных чисел. Статистическая обработка упорядоченной совокупности указанного количества знаков p показывает, что она обладает многими чертами случайной последовательности.

Есть несколько забавных способов запомнить число p точнее, чем просто 3,14. Например, выучив следующее четверостишие, можно без труда назвать семь десятичных знаков p :

Нужно только постараться

И запомнить все как есть:

Три, четырнадцать, пятнадцать,

Девяносто два и шесть .

(С.Бобров Волшебный двурог )

Подсчет количества букв в каждом слове следующих фраз так же дает значение числа p :

«Что я знаю о кругах?» (p » 3,1416). Эту поговорку предложил Я.И.Перельман.

«Вот и знаю я число, именуемое Пи. – Молодец!» (p » 3,1415927).

«Учи и знай в числе известном за цифрой цифру, как удачу примечать» (p » 3,14159265359).

Учитель одной из московских школ придумал строку: «Это я знаю и помню прекрасно», а его ученица сочинила забавное продолжение: «Пи многие знаки мне лишни, напрасны». Это двустишие позволяет определить 12 цифр.

А так выглядит 101 знак числа p без округления

3,14159 26535 89793 23846 26433 83279 50288 41971 69399 37510 58209 74944 59230 78164 06286 20899 86280 34825 34211 70679.

В наше время с помощью ЭВМ значение числа p вычислено с миллионами правильных знаков, но такая точность не нужна ни в каких вычислениях. А вот возможность аналитического определения числа ,

В последней формуле в числителе стоят все простые числа, а знаменатели отличаются от них на единицу, причем знаменатель больше числителя, если тот имеет вид 4n + 1, и меньше в противном случае.

Хотя еще с конца 16 в., т.е. с тех пор, как сформировались сами понятия рациональных и иррациональных чисел, многие ученые были убеждены в том, что p – число иррациональное, но только в 1766 немецкий математик Иоганн Генрих Ламберт (1728–1777), основываясь на открытой Эйлером зависимости между показательной и тригонометрической функциями, строго доказал это. Число p не может быть представлено в виде простой дроби, как ни были бы велики числитель и знаменатель.

В 1882 профессор Мюнхенского университета Карл Луиз Фердинанд Линдеман (1852–1939) используя результаты, полученные французским математиком Ш.Эрмитом , доказал, что p – число трансцендентное, т.е. оно не является корнем никакого алгебраического уравнения a n x n + a n– 1 x n– 1 + … + a 1 x + a 0 = 0 с целыми коэффициентами. Это доказательство поставило точку в истории древнейшей математической задачи о квадратуре круга. Тысячелетия эта задача не поддавалась усилиям математиков, выражение «квадратура круга» стало синонимом неразрешимой проблемы. А все дело оказалось в трансцендентной природе числа p .

В память об этом открытии в зале перед математической аудиторией Мюнхенского университета был установлен бюст Линдемана. На постаменте под его именем изображен круг, пересеченный квадратом равной площади, внутри которого начертана буква p .

Марина Федосова