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

Архимед в сочинении “Измерение круга” вычислил отношение длины окружности к диаметру (число ) и нашел, что оно заключено между 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

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

В числе ПИ очень много загадок. Вернее это даже не загадки, а своего рода какая-то Истина, которую за всю историю человечества никто еще не разгадал…

Что такое число Пи? Число ПИ - математическая «константа», выражающая отношение длины окружности к её диаметру. Сначала по невежеству его (это отношение) считали равным трем, что было грубо приближенно, но им хватало. Но когда времена доисторические сменились временами древними (т.е. уже историческими), то удивлению пытливых умов не было предела: оказалось, что число три весьма неточно выражает это соотношение. С течением времени и развитием наук это число стали полагать равным двадцати двум седьмым.

Английский математик Август де Морган назвал как-то число ПИ “…загадочным числом 3,14159…, которое лезет в дверь, в окно и через крышу”. Неутомимые ученые продолжали и продолжали вычислять десятичные знаки числа Пи, что является на самом деле дико нетривиальной задачей, потому что просто так в столбик его не вычислить: число это не только иррациональное, но и трансцендентное (это вот как раз такие числа, которые не вычисляются путем простых уравнений).

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

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

В 1965-м году американский математик М. Улэм, сидя на одном скучном собрании, от нечего делать начал писать на клетчатой бумаге цифры, входящие в число пи. Поставив в центре 3 и двигаясь по спирали против часовой стрелки, он выписывал 1, 4, 1, 5, 9, 2, 6, 5 и прочие цифры после запятой. Попутно он обводил все простые числа кружками. Каково же было его удивление и ужас, когда кружки стали выстраиваться вдоль прямых!

В десятичном хвосте числа пи можно отыскать любую задуманную последовательность цифр. Любая последовательность цифр в десятичных знаках числа пи рано или поздно найдется. Любая!

Ну и что? – спросите вы. А то. Прикиньте: если там есть ваш телефон (а он есть), то ведь там же есть и телефон той девушки, которая не захотела дать вам свой номер. Более того, там есть и номера кредиток, и даже все значения выигрышных номеров завтрашнего тиража лотереи. Да что там, вообще всех лотерей на много тысячелетий вперед. Вопрос в том, как их там отыскать…

Если зашифровать все буквы цифрами, то в десятичном разложении числа пи можно найти всю мировую литературу и науку, и рецепт изготовления соуса бешамель, и все священные книги всех религий. Это строгий научный факт. Ведь последовательность БЕСКОНЕЧНА и сочетания в числе ПИ не повторяются, следовательно она содержит ВСЕ сочетания цифр, и это уже доказано. А раз все, то ВСЕ. В том числе и такие, которые соответствуют выбранной вами книге.

А это опять-таки означает, что там содержится не только вся мировая литература, которая уже написана (в частности и те книги, которые сгорели и т.д.), но и все книги, которые еще БУДУТ написаны. В том числе и Ваши статьи на сайтах. Получается, что это число (единственное разумное число во Вселенной!) и управляет нашим миром. Надо только рассмотреть побольше знаков, найти нужный участок и расшифровать его. Это чем-то сродни парадоксу со стадом шимпанзе, долбящем по клавиатуре. При достаточно долгом (можно даже оценить это время) эксперименте они напечатают все пьесы Шекспира.

Тут же напрашивается аналогия с периодически появляющимися сообщениями о том, что в Ветхом Завете, якобы, закодированы послания потомкам, поддающиеся прочтению с помощью хитроумных программ. Отметать сходу такую экзотическую особенность Библии не совсем мудро, кабаллисты веками занимаются поиском таких пророчеств, но хотелось бы привести сообщение одного исследователя, который с помощью компьютера нашел в Ветхом завете слова о том, что в Ветхом Завете нет никаких пророчеств. Скорее всего, в очень большом тексте, так же, как и в бесконечных цифрах числа ПИ, можно не только закодировать любую информацию, но и “найти” фразы, изначально не заложенные туда.

Для практики, в пределах Земли достаточно 11 знаков после точки. Тогда, зная, что радиус Земли равен 6400 км или 6,4*1012 миллиметров, получится, что мы, отбросив двенадцатую цифру в числе ПИ после точки при вычислении длины меридиана, ошибемся на несколько миллиметров. А при расчете длины Земной орбиты при вращении вокруг Солнца (как известно, R=150*106 км = 1,5*1014 мм) для такой же точности достаточно использовать число ПИ с четырнадцатью знаками после точки, да что уж там мелочиться - диаметр нашей Галактики около 100.000 световых лет (1 световой год примерно равен 1013 км) или 1018 км или 1030 мм., а еще в XVII веке были получены 34 знака числа ПИ, избыточные для таких расстояний, а их на данный момент вычислено до 12411-триллионного знака !!!

Отсутствие периодически повторяющихся цифр, а именно, исходя их формулы Длина окружности=Пи*D окружность не замыкается, так как нет конечного числа. Этот факт также может тесно быть связан с спиральным проявлением в нашей жизни …

Есть еще гипотеза о том, что все (или некоторые) универсальные постоянные (постоянная Планка, число Эйлера, универсальная гравитационная постоянная, заряд электрона и т.д.) со временем меняют свои значения, так как меняется кривизна пространства из-за перераспределения материи или по другим, не известным нам причинам.

Рискуя навлечь гнев просвещенного сообщества, можем предположить, что и рассматриваемое сегодня число ПИ, отражающее свойства Вселенной, может со временем меняться. Во всяком случае, никто не может нам запретить заново найти значение числа ПИ, подтвердив (или не подтвердив) имеющиеся значения.

10 интересных фактов про число ПИ

1. История числа насчитывает не одно тысячелетие, почти столько, сколько существует наука математика. Конечно, точное значение числа рассчитали не сразу. Поначалу отношение длины окружности к диаметру считали равным 3. Но с течением времени, когда начала развиваться архитектура, потребовалось более точное измерение. Кстати, число существовало, а вот буквенное обозначение оно получило только в начале XVIII века (1706 год) и происходит от начальных букв двух греческих слов, означающих «окружность» и «периметр». Буквой «π» число наделил математик Джонс, а прочно вошла в математику она уже в 1737 году.

2. В разные эпохи и у разных народов число Пи имело разное значение. Например, в Древнем Египте оно равнялось 3,1604, у индусов оно приобрело значение 3,162, китайцы пользовались числом, равным 3,1459. С течением времени π рассчитывали все точнее, а когда появилась вычислительная техника, то есть компьютер, оно стало насчитывать более 4 миллиардов знаков.

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

4. Примечательно, что значение числа Пи пытались вводить даже на уровне государства, то есть посредством закона. В 1897 году в штате Индиана подготовили билль. Согласно документу Пи равнялось 3,2. Однако ученые вовремя вмешались и предотвратили таким образом ошибку. В частности, против билля выступил профессор Пердью, присутствовавший на законодательном собрании.

5. Интересно, что свое имя имеют несколько чисел в бесконечной последовательности Пи. Так, шесть девяток числа Пи носят имя американского физика. Как-то Ричард Фейнман читал лекцию и ошарашил публику замечанием. Он сказал, что хотел бы наизусть выучить цифры числа Пи до шести девяток только для того, чтобы под конец рассказа произнести шесть раз «девять», намекая на то, что его значение рационально. Тогда как на самом деле оно иррационально.

6. Математики всего мира не прекращают вести исследования, связанные с числом Пи. Оно буквально окутано некой тайной. Некоторые теоретики даже полагают, что в нем заключена вселенская истина. Чтобы обмениваться знаниями и новой информацией о Пи, организовали Пи-клуб. Вступить в него непросто, нужно иметь незаурядную память. Так, желающих стать членом клуба экзаменуют: человек должен по памяти рассказать как можно больше знаков числа Пи.

7. Придумали даже различные техники для запоминания числа Пи после запятой. Например, придумывают целые тексты. В них слова имеют то же количество букв, что и соответствующая цифра после запятой. Чтобы еще упростить запоминание такого длинного числа, сочиняют стихи по тому же принципу. Члены Пи-клуба частенько развлекаются таким образом, а заодно тренируют память и сообразительность. Например, такое хобби было у Майка Кейта, который восемнадцать лет назад придумал рассказ, каждое слово в котором равнялось почти четырем тысячам (3834) первых знаков числа Пи.

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

9. День числа Пи отмечают больше четверти века, с 1988 года. Однажды физик из научно-популярного музея в Сан-Франциско Ларри Шоу заметил, что 14 марта по написанию совпадает с числом Пи. В дате месяц и число образуют 3.14.

10. Есть любопытное совпадение. 14 марта родился великий ученый Альберт Эйнштейн, создавший, как известно, теорию относительности.

Недавно на Хабре в одной статье упомянули про вопрос «Что было бы с миром, если бы число Пи равнялось 4?» Я решил слегка поразмышлять на эту тему, используя некоторые (пусть и не самые обширные) знания в соответствующих областях математики. Кому интересно – прошу под кат.

Чтобы представить такой мир, нужно математически реализовать пространство с иным соотношением длины окружности к ее диаметру. Это я и попытался сделать.

Попытка №1.

Оговорим сразу, что рассматривать я буду только двумерные пространства. Почему? Потому что окружность, собственно, определена в двумерном пространстве (если рассмотреть размерность n>2, то отношение меры (n-1)-мерной окружности к ее радиусу даже не будет константой).
Так что для начала я попытался придумать хоть какое-то пространство, где Пи не равно 3.1415… Для этого я взял метрическое пространство с метрикой, в которой расстояние между двумя точками равно максимуму среди модулей разности координат (т.е. расстояние Чебышева).

Какой же вид будет иметь единичная окружность в этом пространстве? Возьмем точку с координатами (0,0) за центр этой окружности. Тогда множество точек, расстояние (в смысле заданной метрики) от которых до центра равно 1, есть 4 отрезка, параллельных осям координат, образующих квадрат со стороной 2 и с центром в нуле.

Да, в некоторой метрике это - окружность!

Посчитаем здесь Пи. Радиус равен 1, тогда диаметр, соответственно, равен 2. Можно также рассмотреть определение диаметра как наибольшего расстояния между двумя точками, но даже так оно равно 2. Осталось найти длину нашей «окружности» в данной метрике. Это сумма длин всех четырех отрезков, которые в данной метрике имеют длину max(0,2)=2. Значит, длина окружности равна 4*2=8. Ну а тогда Пи здесь равно 8/2=4. Получилось! Но нужно ли сильно радоваться? Результат этот практически бесполезен, ведь рассматриваемое пространство абсолютно абстрактно, в нем даже не определены углы и повороты. Вы можете представить себе мир, где по факту не определен поворот, и где окружностью является квадрат? Я пытался, честно, но у меня не хватило воображения.

Радиус равен 1, а вот с нахождением длины этой «окружности» есть некоторые сложности. После некоторых поисков информации в интернете, я пришел к выводу, что в псевдоевклидовом пространстве такое понятие как «число Пи» вообще не может быть определено, что, безусловно, плохо.

Если кто-нибудь в комментариях расскажет мне, как формально считать длину кривой в псевдоевклидовом пространстве, я буду очень рад, ибо моих познаний в дифференциальной геометрии, топологии (а также усердного гугления) для этого не хватило.

Выводы:

Не знаю, можно ли писать о выводах после таких не сильно продолжительных исследований, но кое-что сказать можно. Во-первых, попытавшись представить пространство с иным числом Пи, я понял, что оно будет слишком абстрактно, чтобы быть моделью реального мира. Во-вторых, когда если попытаться придумать более удачную модель (похожую на наш, реальный мир), выходит, что число Пи останется неизменным. Если принять за данность возможность отрицательного квадрата расстояния (что для обычного человека - просто абсурд), то Пи не будет определено вовсе! Все это и наводит на мысль, что, возможно, мира с другим числом Пи и вовсе быть не могло? Ведь не зря же Вселенная именно такая, какая она есть. А может быть, это и реально, только обычной математики, физики и человеческого воображения для этого недостаточно. А вы как считаете?

Upd. Узнал точно. Длина кривой в псевдоевклидовом пространстве может быть определена только на каком-либо его евклидовом подпространстве. То есть, в частности, для получившейся в попытке N3 «окружности» вовсе не определено такое понятие как «длина». Соответственно, Пи там тоже посчитать нельзя.

С недавних пор существует элегантная формула для вычисления числа Пи, которую в 1995 году впервые опубликовали Дэвид Бэйли, Питер Борвайн и Саймон Плафф:

Казалось бы: что в ней особенного — формул для вычисления Пи великое множество: от школьного метода Монте-Карло до труднопостижимого интеграла Пуассона и формулы Франсуа Виета из позднего Средневековья. Но именно на эту формулу стоит обратить особое внимание — она позволяет вычислить n-й знак числа пи без нахождения предыдущих. За информацией о том, как это работает, а также за готовым кодом на языке C, вычисляющим 1 000 000-й знак, прошу под хабракат.

Как же работает алгоритм вычисления N-го знака Пи?
К примеру, если нам нужен 1000-й шестнадцатеричный знак числа Пи, мы домножаем всю формулу на 16^1000, тем самым обращая множитель, стоящий перед скобками, в 16^(1000-k). При возведении в степень мы используем двоичный алгоритм возведения в степень или, как будет показано в примере ниже, возведение в степень по модулю . После этого вычисляем сумму нескольких членов ряда. Причём необязательно вычислять много: по мере возрастания k 16^(N-k) быстро убывает, так что, последующие члены не будут оказывать влияния на значение искомых цифр). Вот и вся магия — гениальная и простая.

Формула Бэйли-Борвайна-Плаффа была найдена Саймоном Плаффом при помощи алгоритма PSLQ , который был в 2000 году включён в список Top 10 Algorithms of the Century . Сам же алгоритм PSLQ был в свою очередь разработан Бэйли. Вот такой мексиканский сериал про математиков.
Кстати, время работы алгоритма — O(N), использование памяти — O(log N), где N — порядковый номер искомого знака.

Думаю, уместно будет привести код на языке Си, написанный непосредственно автором алгоритма, Дэвидом Бэйли:

/* This program implements the BBP algorithm to generate a few hexadecimal digits beginning immediately after a given position id, or in other words beginning at position id + 1. On most systems using IEEE 64-bit floating- point arithmetic, this code works correctly so long as d is less than approximately 1.18 x 10^7. If 80-bit arithmetic can be employed, this limit is significantly higher. Whatever arithmetic is used, results for a given position id can be checked by repeating with id-1 or id+1, and verifying that the hex digits perfectly overlap with an offset of one, except possibly for a few trailing digits. The resulting fractions are typically accurate to at least 11 decimal digits, and to at least 9 hex digits. */ /* David H. Bailey 2006-09-08 */ #include #include int main() { double pid, s1, s2, s3, s4; double series (int m, int n); void ihex (double x, int m, char c); int id = 1000000; #define NHX 16 char chx; /* id is the digit position. Digits generated follow immediately after id. */ s1 = series (1, id); s2 = series (4, id); s3 = series (5, id); s4 = series (6, id); pid = 4. * s1 - 2. * s2 - s3 - s4; pid = pid - (int) pid + 1.; ihex (pid, NHX, chx); printf (" position = %i\n fraction = %.15f \n hex digits = %10.10s\n", id, pid, chx); } void ihex (double x, int nhx, char chx) /* This returns, in chx, the first nhx hex digits of the fraction of x. */ { int i; double y; char hx = "0123456789ABCDEF"; y = fabs (x); for (i = 0; i < nhx; i++){ y = 16. * (y - floor (y)); chx[i] = hx[(int) y]; } } double series (int m, int id) /* This routine evaluates the series sum_k 16^(id-k)/(8*k+m) using the modular exponentiation technique. */ { int k; double ak, eps, p, s, t; double expm (double x, double y); #define eps 1e-17 s = 0.; /* Sum the series up to id. */ for (k = 0; k < id; k++){ ak = 8 * k + m; p = id - k; t = expm (p, ak); s = s + t / ak; s = s - (int) s; } /* Compute a few terms where k >= id. */ for (k = id; k <= id + 100; k++){ ak = 8 * k + m; t = pow (16., (double) (id - k)) / ak; if (t < eps) break; s = s + t; s = s - (int) s; } return s; } double expm (double p, double ak) /* expm = 16^p mod ak. This routine uses the left-to-right binary exponentiation scheme. */ { int i, j; double p1, pt, r; #define ntp 25 static double tp; static int tp1 = 0; /* If this is the first call to expm, fill the power of two table tp. */ if (tp1 == 0) { tp1 = 1; tp = 1.; for (i = 1; i < ntp; i++) tp[i] = 2. * tp; } if (ak == 1.) return 0.; /* Find the greatest power of two less than or equal to p. */ for (i = 0; i < ntp; i++) if (tp[i] > p) break; pt = tp; p1 = p; r = 1.; /* Perform binary exponentiation algorithm modulo ak. */ for (j = 1; j <= i; j++){ if (p1 >= pt){ r = 16. * r; r = r - (int) (r / ak) * ak; p1 = p1 - pt; } pt = 0.5 * pt; if (pt >= 1.){ r = r * r; r = r - (int) (r / ak) * ak; } } return r; }
Какие возможности это даёт? Например: мы можем создать систему распределённых вычислений, рассчитывающую число Пи и поставить всем Хабром новый рекорд по точности вычисления (который сейчас, к слову, составляет 10 триллионов знаков после запятой). Согласно эмпирическим данным, дробная часть числа Пи представляет собой нормальную числовую последовательность (хотя доказать это достоверно ещё не удалось), а значит, последовательности цифр из него можно использовать в генерации паролей и просто случайных чисел, или в криптографических алгоритмах (например, в хэшировании). Способов применения можно найти великое множество - надо только включить фантазию.

Больше информации по теме вы можете найти в статье самого Дэвида Бэйли, где он подробно рассказывает про алгоритм и его имплементацию (pdf);

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

Таблица значений тригонометрических функций

Примечание . В данной таблице значений тригонометрических функций используется знак √ для обозначения квадратного корня. Для обозначения дроби - символ "/".

См. также полезные материалы:

Для определения значения тригонометрической функции , найдите его на пересечении строки с указанием тригонометрической функции. Например, синус 30 градусов - ищем колонку с заголовком sin (синус) и находим пересечение этой колонки таблицы со строкой "30 градусов", на их пересечении считываем результат - одна вторая. Аналогично находим косинус 60 градусов, синус 60 градусов (еще раз, в пересечении колонки sin (синус) и строки 60 градусов находим значение sin 60 = √3/2) и т.д. Точно так же находятся значения синусов, косинусов и тангенсов других "популярных" углов.

Синус пи, косинус пи, тангенс пи и других углов в радианах

Приведенная ниже таблица косинусов, синусов и тангенсов также подходит для нахождения значения тригонометрических функций, аргумент которых задан в радианах . Для этого воспользуйтесь второй колонкой значений угла. Благодаря этому можно перевести значение популярных углов из градусов в радианы. Например, найдем угол 60 градусов в первой строке и под ним прочитаем его значение в радианах. 60 градусов равно π/3 радиан.

Число пи однозначно выражает зависимость длины окружности от градусной меры угла. Таким образом, пи радиан равны 180 градусам.

Любое число, выраженное через пи (радиан) можно легко перевести в градусную меру, заменив число пи (π) на 180 .

Примеры :
1. Синус пи .
sin π = sin 180 = 0
таким образом, синус пи - это тоже самое, что синус 180 градусов и он равен нулю.

2. Косинус пи .
cos π = cos 180 = -1
таким образом, косинус пи - это тоже самое, что косинус 180 градусов и он равен минус единице.

3. Тангенс пи
tg π = tg 180 = 0
таким образом, тангенс пи - это тоже самое, что тангенс 180 градусов и он равен нулю.

Таблица значений синуса, косинуса, тангенса для углов 0 - 360 градусов (часто встречающиеся значения)

значение угла α (градусов)	значение угла α в радианах (через число пи)	sin (синус)	cos (косинус)	tg (тангенс)	ctg (котангенс)	sec (секанс)	cosec (косеканс)
0	0	0	1	0	-	1	-
15	π/12			2 - √3	2 + √3
30	π/6	1/2	√3/2	1/√3	√3	2/√3	2
45	π/4	√2/2	√2/2	1	1	√2	√2
60	π/3	√3/2	1/2	√3	1/√3	2	2/√3
75	5π/12			2 + √3	2 - √3
90	π/2	1	0	-	0	-	1
105	7π/12		-	- 2 - √3	√3 - 2
120	2π/3	√3/2	-1/2	-√3	-√3/3
135	3π/4	√2/2	-√2/2	-1	-1	-√2	√2
150	5π/6	1/2	-√3/2	-√3/3	-√3
180	π	0	-1	0	-	-1	-
210	7π/6	-1/2	-√3/2	√3/3	√3
240	4π/3	-√3/2	-1/2	√3	√3/3
270	3π/2	-1	0	-	0	-	-1
360	2π	0	1	0	-	1	-

Если в таблице значений тригонометрических функций вместо значения функции указан прочерк (тангенс (tg) 90 градусов, котангенс (ctg) 180 градусов) значит при данном значении градусной меры угла функция не имеет определенного значения. Если же прочерка нет - клетка пуста, значит мы еще не внесли нужное значение. Мы интересуемся, по каким запросам к нам приходят пользователи и дополняем таблицу новыми значениями, несмотря на то, что текущих данных о значениях косинусов, синусов и тангенсов самых часто встречающихся значений углов вполне достаточно для решения большинства задач.