圆周率

圆周率

圆周率
数表—无理数 $\color {blue}{\sqrt {2}}$ - $\color {blue}\varphi$ - $\color {blue}{\sqrt {3}}$ - $\color {blue}{\sqrt {5}}$ - $\color {blue}\delta _{S}$ - $\color {blue}e$ - $\color {blue}\pi$
识别
种类	无理数 超越数
符号	$\pi$
位数数列编号	A000796
性质
定义	$\pi ={\frac {C}{d}}.$，其中$C$为圆周长、$d$为直径 $\pi =\int _{-1}^{1}{\frac {dx}{\sqrt {1-x^{2}}}}.$
连分数	$\pi =3+\textstyle {\frac {1}{7+\textstyle {\frac {1}{15+\textstyle {\frac {1}{1+\textstyle {\frac {1}{292+\textstyle {\frac {1}{1+\textstyle {\frac {1}{1+\textstyle {\frac {1}{1+\ddots }}}}}}}}}}}}}}$
以此为根的多项式或函数	$e^{ix}+1=0$
表示方式
值	$\pi \approx$3.14159265
无穷级数	$\pi =\sum _{k=0}^{\infty }{\frac {4(-1)^{k}}{2k+1}}$
二进制	11.00100100001111110110…
十进制	3.14159265358979323846…
十六进制	3.243F6A8885A308D31319…
六十进制	3;8,29,44,0,47,25,53,7,24,57,36…

圆周率是数学常数，等于任何圆的周长和其直径的比，一个常见的近似值等于3.14159265，常用符号$\pi$表示。

$\pi$是无理数，不能用分数表示出来（即它的小数部分是无限不循环小数），但近似${\textstyle {\frac {22}{7}}}$等有理数。学界认为π的数字序列在统计上是随机分布，但迄今未能证明。此外，π还是超越数，亦即它不是任何有理系数多项式的根；化圆为方的问题不可能用尺规作图解决。

几个文明古国很早就须计算出π的精确值以便于生产的计算。西元5世纪，中国刘宋数学家祖冲之用几何方法将圆周率计算到小数点后7位。大约同时，印度数学家也将圆周率计算到小数点后5位。史上首条π的精确无穷级数公式（即π的莱布尼茨公式）直到约1000年后才由印度数学家发现。微积分出现，π的位数很快计到数百位，足以满足任何科学工程的计算需求。在20和21世纪，计算机技术快速发展，π的计算精度急速提高。截至2024年3月，π的十进制精度已达105万亿位。几乎所有科学研究对π的精度要求都不超过几百位，当前计算π的值主要都为打破记录、测试超级计算机的计算能力和高精度乘法算法。

π的定义涉及圆，在三角学和几何学的许多公式，特别是广泛应用在圆形、球形或椭球形相关公式中。在近代数学分析里，π改由实数系统谱性质中的特征值或周期来定义，其他数学领域如数论、统计以及几乎所有物理学领域均有出现，π的广泛用途使它成为科学界内外最广为人知的数学常数。几本专门介绍π的书籍经已出版，圆周率日（3月14日）和π值计算突破记录也往往会成为报纸的新闻头条。此外，背诵π值的世界记录已达10万位。

直径为一的圆的周长是π（3.14159265...）

基本概念

名称

数学家用小写希腊字母$\pi$表示圆周和其直径之比，有时也将其拼写为“Pi”，来自希腊语“περίμετρος”（周长）的首字母。英语π的发音与英文单词“Pie”（/paɪ/，西式馅饼）相同。π的小写字母（或其无衬线体）在数学要和表示连乘积的大写Π相区分开。

关于选择符号π的原因，请参见引入π符号一节。

定义

圆周长略大于其直径的三倍；精确的比例称为$\pi$

π常用定义为圆的周长$C$与直径$d$的比值：

$\pi ={\frac {C}{d}}$。

无论圆的大小如何，比值${\frac {C}{d}}$为恒值。如果圆的直径变为原先的二倍，周长也变为二倍，比值${\frac {C}{d}}$不变。π目前的定义暗地用了欧几里得几何的一些定理，虽然圆的定义可扩展到任意曲面（即非欧几里得几何），但这些圆不符合定律$\pi ={\frac {C}{d}}$。

这里，圆的周长指其圆周的弧长，弧长这概念可以不依赖几何学，而是用微积分学的极限来定义。例如，若想计算笛卡儿坐标系中单位圆$x^{2}+y^{2}=1$上半部分的弧长，需要用到积分：

$\pi =\int _{-1}^{1}{\frac {dx}{\sqrt {1-x^{2}}}}.$

上述积分是由卡尔·魏尔斯特拉斯于1841年对π的积分定义。

π这些依赖周长、且暗地依赖积分的定义如今在文献中并不常见。雷默特（Remmert (1991)）解释说现代教微积分时，大学一般将微分学课程安排在积分学课程之前，所以不依赖于后者的π的定义就很有必要了。其中一种定义由理查·巴尔策提出，由爱德蒙·兰道推广，其表述如下：π是两倍于能使余弦函数等于零的最小正数。余弦函数可以由独立于几何之外的幂级数定义，或者使用微分方程的解来定义。

在相似的启发下，π可以用关于复变量$z$的复指数函数$\exp(z)$来定义。复指数类似余弦函数，可用多种方式定义。令函数$\exp(z)$值为一的复数集合是如下所示的（虚）等差数列：

$\{\dots ,-2\pi i,0,2\pi i,4\pi i,\dots \}=\{2\pi ki|k\in \mathbb {Z} \}$，

并且其中包括独特的正实数$\pi$。

基于同样想法但更抽象的定义运用了精巧的拓扑学和代数学概念，用以下定理描述：存在唯一的从加法模数整数组成的实数群R/Z到绝对值为1的复数组成的乘法群的连续同态（拓扑学概念，指在拓扑空间之间的一种态射）。数字π定义为此同态派生的模的一半。

周长固定，圆会围成最大面积，π同样表述为等周不等式中出现的常数（乘四分之一）。此外，在很多其他紧密相关的方程中，π作为某些几何或者物理过程的特征值出现；详见下文。

无理及正规性

π是无理数，无法表示成两整数之比的形式（形如${\textstyle {\frac {22}{7}}}$的分数常用来近似表达π，但是没有任何普通分数（指整数的比）可以取到π的精确值）。由于$\pi$是无理数，故可表示为无限不循环小数。有多种方法能证明π是无理数，这些证明也都要用到微积分学和反证法。$\pi$可以用有理数来近似的程度还无法准确得知（称为无理性度量），不过估计其无理性度量比e或ln(2)的要大，但是小于刘维尔数的无理性度量。

统计随机性检验，包括正规数检验，可验证$\pi$的位数没有明显的固定模式。$\pi$的小数中任意固定长度的序列（如3位数000，001……999）出现机率都相同。不过有关π是正规数的猜想既无证明，亦无证伪。

电脑出现后可生成大量π的不同位数，并统计分析之。金田康正详细统计分析了π的十进制数字，并验证了其分布正规：例如，假设检定0到9十个数的出现频率，找不到有特定重复规律的证据。根据无限猴子定理，任何任意长度、由随机内容组成的子序列看起来都有可能像不随机生成。因此，就算π的小数序列通过了随机性统计测试，其中也可能有几位的数字看起来似有规律可循而非随机数，例如π的十进制写法在小数第762位后开始出现了连续六个9。

超越性

由于π是超越数，不能利用尺规作图化圆为方。

$\pi$不仅是无理数，还是超越数，即$\pi$不是任何有理系数多项式的根。（比方说，试图解有限项方程${\textstyle {\frac {x^{5}}{120}}-{\frac {x^{3}}{6}}+x=0}$来求$\pi$的值）

$\pi$的超越性衍生出一些重要的结果：$\pi$不能经有限次四则运算和开平方运算有理数来获得，因此不是规矩数。换言之，尺规作图作不出长度为$\pi$的线段，也就不可能用尺规方法做出与已知圆面积相等的正方形。后者即为有名的化圆为方问题，该问题早在古典时代即已提出，曾困扰人数千年之久。直至今天，依然有民间数学爱好者声称他们解决了这问题。

连分式

$\pi$像所有无理数一样无法表示成分数，但$\pi$等全部无理数都能表示成一系列叫连分数的连续分数形式：

$\pi =3+\textstyle {\frac {1}{7+\textstyle {\frac {1}{15+\textstyle {\frac {1}{1+\textstyle {\frac {1}{292+\textstyle {\frac {1}{1+\textstyle {\frac {1}{1+\textstyle {\frac {1}{1+\ddots }}}}}}}}}}}}}}$

在这连分数的任意一点截断化简，都能得到π的近似值；前四位近似值是3、${\frac {22}{7}}$、${\frac {333}{106}}$、${\frac {355}{113}}$。这些数在历史上是π最广为人知且广为使用的几个近似值。用以上方式得出的$\pi$的近似值要比任何有相同或更小的整数分母的其他整数分数近似值更接近π。π是超越数，据定义来说它不是代数数，又因此不可能是二次无理数；是故π不能表示为循环连分数。尽管$\pi$的简单连分数没有表现出任何其他明显规律，数学家发现了数条广义连分数能表示π，例如：

$\pi =\textstyle {\cfrac {4}{1+\textstyle {\frac {1^{2}}{2+\textstyle {\frac {3^{2}}{2+\textstyle {\frac {5^{2}}{2+\textstyle {\frac {7^{2}}{2+\textstyle {\frac {9^{2}}{2+\ddots }}}}}}}}}}}}=3+\textstyle {\frac {1^{2}}{6+\textstyle {\frac {3^{2}}{6+\textstyle {\frac {5^{2}}{6+\textstyle {\frac {7^{2}}{6+\textstyle {\frac {9^{2}}{6+\ddots }}}}}}}}}}=\textstyle {\cfrac {4}{1+\textstyle {\frac {1^{2}}{3+\textstyle {\frac {2^{2}}{5+\textstyle {\frac {3^{2}}{7+\textstyle {\frac {4^{2}}{9+\ddots }}}}}}}}}}$

近似值

圆周率近似值包括：

• 整数：3
• 分数（依准确度顺序排列）：13/4、16/5、19/6、22/7、179/57、267/85、333/106、355/113、52163/16604、53228/16943、55358/17621、57843/18412、60328/19203、103993/33102、245850922/78256779（选自 A063674 及 A063673。）
• 小数（整数后首80位）：3.14159265358979323846264338327950288419716939937510582097494459230781640628620899...（另见 A000796）

其他进位制的近似值

• 二进制（整数后首48位）：11.001001000011111101101010100010001000010110100011…
• 十六进制（整数后首20位）：3.243F6A8885A308D31319…
• 六十进制（整数后首20位）：3;8,29,44,0,47,25,53,7,24,57,36,17,43,4,29,7,10,3,41,17…

复数与欧拉恒等式

欧拉公式给出了e的复指数与复平面上以原点为圆心的单位圆上的点之间的关系。

任何复数（以$z$为例）都可以表示为一组实数对：极坐标系用实数$r$表示半径，代表复平面上复数$z$离原点的距离；实数$\varphi$则表示夹角，即这条半径（复平面上复数$z$与原点的连线）与正实轴经顺时针转动的夹角。这样一来，$z$就可写成

$z=r\cdot (\cos \varphi +i\sin \varphi )$，这里$i$代表虚数单位，即$i^{2}$＝-1。

复分析中，欧拉公式将三角函数与复指数函数糅合在一起：

$e^{i\varphi }=\cos \varphi +i\sin \varphi$，这里数学常数e是自然对数的底数。

欧拉公式确立了$e$的复指数与复平面上以原点为圆心的单位圆上的点之间的关系，而且当$\varphi =\pi$时，欧拉公式就能改写为欧拉恒等式的形式：

${{{e}^{{i}\,{\pi }}}+{1}}=0$。此等式亦称“最奇妙的数学公式”（英语：the most remarkable formula in mathematics），全因它将五个最基本的数学常数简洁联系起来。

欧拉等式亦可用于求出方程$z^{n}=1$的$n$个不同复数根（这些根叫做$n$次单位根”），可以根据以下公式求得：

$e^{2\pi ik/n}\qquad (k=0,1,2,\dots ,n-1)$。

谱特征

震荡弦的泛音是二次微分的本征函数，会形成泛音列。对应的本征值会形成由π整数倍组成的等差数列

$\pi$常出现在有关几何的问题中。然而，不少和几何无关的问题也可看到$\pi$的身影。

$\pi$在许多用处中都会以特征值形式出现。例如理想的振动弦问题可以建模为函数$f$在单位区间$[0,1]$的图形，固定边界值为$f(0)=f(1)=0$。弦振动的模态会是微分方程的$f^{n}(x)+\lambda ^{2}f(x)=0$，此处λ是相关的特征值。受施图姆－刘维尔理论限制，$\lambda$只能是一些特定的数值。而$\lambda =\pi$即为一个特征值，因为函数$f(x)=\sin(\pi x)$满足边界条件及微分方程$\lambda =\pi$。

依照第一代开尔文男爵威廉·汤姆森所述的一篇传说，古迦太基城的外形是等周长问题的一项解(Thompson 1894)。这些包围著海的区域由迦太基女王狄多所围，城不靠海的边界须用指定大小的牛皮围住，后来是将牛皮剪成小段

$\pi$是上述方程的最小特征值，也和弦振动的基本模式有关。一种让弦振动的方式是提供弦能量，能量会满足维廷格函数不等式，其中提到若函数$f:[0,1]\rightarrow \mathbb {C}$使得$f(0)=f(1)=0$，且$f$和$f'$都是平方可积函数，则以下的不等式成立：

$\pi ^{2}\int _{0}^{1}|f(x)|^{2}\,dx\leq \int _{0}^{1}|f'(x)|^{2}\,dx,$

此例中等号成立的条件恰好是$f$为$\sin(\pi x)$倍数的时候。因此$\pi$似乎是维尔丁格不等式的最佳常数，也是最小的特征值（根据雷利商数的计算方式）

$\pi$在更高维度的分析也有类似的角色，出现在其他类似问题的特征值中。就如以上所述，$\pi$的一项特点是等周定理中的最佳常数：周长为$P$的平面若尔当曲线，所围面积$A$满足以下的不等式

$4\pi A\leq P^{2}$，

$A=\pi r^{2}$及$P=2\pi r$，故等号成立的条件是曲线为圆形。

圆周率π也和庞加莱不等式的最佳常数有关，$\pi$是一维及二维的狄氏能量特征向量最佳值中最小，会出现在许多经典的物理现象中，例如经典的位势论。其一维的情形即为维廷格不等式。

圆周率π也是傅里叶变换的重要常数，傅里叶变换属于积分变换，将实数线上有复数值、可积分的函数，转换为以下形式：

${\hat {f}}(\xi )=\int _{-\infty }^{\infty }f(x)e^{-2\pi ix\xi }\,dx.$

傅里叶变换有几种不同的写法，但不论怎么写，傅里叶变换及反傅里叶变换中，一定会有某处出现$\pi$。不过上述的定义是最经典的，因为其描述了L²空间中唯一的幺正算符，也是$L^{1}$空间到$L^{\infty }$空间的代数同态。

不确定性原理也用到$\pi$。不确定性原理提出了可以将函数在空间及在频域中局部化程度的下限，用傅立叶转换的方式表示：

$\int _{-\infty }^{\infty }x^{2}|f(x)|^{2}\,dx\ \int _{-\infty }^{\infty }\xi ^{2}|{\hat {f}}(\xi )|^{2}\,d\xi \geq \left({\frac {1}{4\pi }}\int _{-\infty }^{\infty }|f(x)|^{2}\,dx\right)^{2}$。

物理的结果，有关量子力学中同时观测位置及动量的不确定性，见下文。傅立叶分析中出现π是史东－凡纽曼定理的结果，证实了海森伯群的薛定谔表示是唯一。

高斯积分

高斯函数${\textstyle f(x)=e^{-x^{2}}}$的图像，函数下方与X轴围成的阴影部分面积为${\textstyle {\sqrt {\pi }}}$。

高斯积分是对高斯函数$e^{-x^{2}}$在整条实轴上的积分，即函数下方与X轴围成的面积，其结果为${\sqrt {\pi }}$，

$\int _{-\infty }^{\infty }e^{-x^{2}}\,dx={\sqrt {\pi }}$

此积分的计算可以先计算$f(x)=e^{-x^{2}}$对整条实轴的积分的平方，通过转换笛卡尔坐标系为极坐标系从而求得

$\left(\int _{-\infty }^{\infty }e^{-x^{2}}\,dx\right)^{2}=\iint _{\mathbf {R} ^{2}}e^{-(x^{2}+y^{2})}\,dxdy=\int _{0}^{2\pi }\int _{0}^{\infty }e^{-r^{2}}r\,dr\,d\theta =\pi$

其他计算方法可参阅高斯积分。高斯函数更一般的形式为${\textstyle f(x)=a\exp {\frac {-(x-b)^{2}}{2c^{2}}}}$，求一般形式的高斯积分均可通过换元积分法转化为求${\textstyle f(x)=e^{-x^{2}}}$的积分。

另外，当高斯函数为以下形式时，它则是平均数为$\mu$和标准差为$\sigma$的正态分布的机率密度函数：

$f(x)={1 \over \sigma {\sqrt {2\pi }}}\,\exp {\frac {-(x-\mu )^{2}}{2\sigma ^{2}}}$

这函数是概率密度函数，函数下方与X轴围成的面积必须为1，令$\mu =0$和$\sigma =1$即可变换得出${\textstyle \int _{-\infty }^{\infty }e^{-x^{2}}\,dx={\sqrt {\pi }}}$。概率论与统计学领域经常使用正态分布来作为复杂现象的简单模型：例如科学家通常假设大多数试验观测值的随机误差都是服从正态分布。

由一维布朗运动的反正弦定律，可以通过试验正信号相对于负信号领先权过零点的分布反过来推算π

概率论与统计学中的中心极限定理解释了正态分布以及$\pi$的核心作用，这定理本质上是联系着$\pi$的谱特征与海森堡不确定性原理相关的特征值，并且在不确定性原理中有

$\sigma _{x}\sigma _{p}\geq \hbar /2$，

这里的$\sigma _{x}$与$\sigma _{p}$分别为位置与动量的标准差，$\hbar$是约化普朗克常数，而不等式的等号当且仅当粒子的波函数为高斯函数使成立。

同样地，$\pi$作为唯一独特的常数使得高斯函数等于其自身的傅里叶变换，此时的高斯函数形式为$f(x)=e^{-\pi x^{2}}$。根据豪（Howe）的说法，建立傅里叶分析基本定理的“全部工作（whole business）”简化为高斯积分。

历史

主条目：π的近似值

参见：π计算年表

远古时期

圆周率在远古时期（西元前一千纪）已估算至前两位（3.1）。有些埃及学家声称，远至古王国时期时期的古埃及人已经用${\textstyle {\frac {22}{7}}}$作为圆周率的约数，但这说法受到质疑。

最早有记载的对圆周率估值在古埃及和巴比伦出现，两估值都与圆周率的正确数值相差不到百分之一。巴比伦曾出土一块西元前1900至1600年的泥板，泥板上的几何学陈述暗示人们当时把圆周率视同${\textstyle {\frac {25}{8}}}$（等于3.125）。埃及的莱因德数学纸草书（鉴定撰写年份为西元前1650年，但抄自一份西元前1850年的文本）载有用作计算圆面积的公式，该公式中圆周率等于${\textstyle ({\frac {16}{9}})^{2}}$（≈3.1605）。

西元前4世纪的《百道梵书》的天文学运算把${\textstyle {\frac {339}{108}}}$（≈3.139，精确到99.91%）用作圆周率估值。西元前150年前其他印度文献把圆周率视为${\textstyle {\sqrt {10}}}$（≈3.1622）。

割圆时代

π可以透过计算圆的外切多边形及内接多边形周长来估算

第一条有纪录、严谨计算π数值的演算法是用正多边形的几何算法，在西元前250年由希腊数学家阿基米德发明。这算法用了有一千年之久，因而有时π亦称阿基米德常数。阿基米德的算法是在计算圆的外切正六边形及内接正六边形的边长，以此计算$\pi$的上限及下限，之后再将六边形变成十二边形，继续计算边长，一直计到正96边形为止。他根据多边形的边长证明${\textstyle {\frac {223}{71}}<\pi <{\frac {22}{7}}}$（也就是${\textstyle 3.1408<\pi <3.1429}$）。阿基米德得到的上限${\textstyle {\frac {22}{7}}}$也造成常见误解，认为$\pi$就等于${\textstyle {\frac {22}{7}}}$。在西元前150年，希腊罗马的科学家克劳狄乌斯·托勒密在《天文学大成》一书中提到π的数值是3.1416，可能来自阿基米德，也可能来自阿波罗尼奥斯。数学家在1630年利用多边形的方式计算π到第39位小数，一直到1699年，其他数学家才利用无穷级数的方式打破其纪录，计算到第71位小数。

阿基米德发展了用多边形近似π的计算方式

中国历史上，$\pi$的数值有3、3.1547（西元前一世纪）、${\sqrt {10}}$（西元前100年，数值约3.1623）及${\textstyle {\frac {142}{45}}}$（第三世纪，数值约3.1556）。大约在西元265年，曹魏数学家刘徽创立割圆术，用正3072边形计算出π的数值为3.1416。他后来又发明了较快的算法，利用边数差两倍的正多边形，其面积的差值会形成等比数列，其公比为${\textstyle {\frac {1}{4}}}$的原理，配合96边形算出π的值为3.14。祖冲之在西元480年利用割圆术计算12288边形边长，得到π的值在3.1415926和3.1415927之间。他同时提出了π的约率${\textstyle {\frac {22}{7}}}$和密率${\textstyle {\frac {355}{113}}}$。在之后的八百年内，这都是π最准确的估计值。为纪念祖冲之对圆周率发展的贡献，日本数学家三上义夫将这推算值命名为“祖冲之圆周率”，简称“祖率”。

印度天文学家阿耶波多在西元499年的著作《阿里亚哈塔历书》中使用了3.1416的数值。斐波那契在大约1220年用独立于阿基米德多边形法，计算出3.1418。义大利作家但丁·阿利吉耶里用的数值则是${\textstyle 3+{\frac {\sqrt {2}}{10}}\approx 3.14142}$。

波斯天文学家卡西在1424年利用3×2²⁸边的多边形，计算到六十进制的第9位小数，相当十进制的第16位小数。这一突破成为当时的纪录，延续了约180年。法国数学家弗朗索瓦·韦达在1579年用3×2¹⁷边形计算到第9位小数，佛兰芒数学家阿德里安·范·罗门在1593年计算到第15位小数。荷兰数学家鲁道夫·范·科伊伦在1596年计算到第20位小数，他之后又计算到第35位小数（因此在二十世纪初之前，圆周率在德国会称为鲁道夫数）。荷兰科学家威理博·司乃耳在1621年计算到第34位小数，而奥地利天文学家克里斯托夫·格林伯格在1630年用10⁴⁰边形计算到第38位小数，至今这仍是利用多边形算法可以达到最准确的结果。

无穷级数

比较几条曾用来计π的无穷级数的收敛情形。S_n是只取前n项的近似值。每张图都是对应前一张图的阴影部份，然后放大横轴10倍。（点击察看细节）

16及17世纪时，开始改用无穷级数的方式去计π。无穷级数是一组无穷数列的和。无穷级数让数学家可以计算出比阿基米德以及其他用几何方式计算的数学家更准确的结果。虽然詹姆斯·格雷果里及戈特弗里德·莱布尼茨等欧洲数学家利用无穷数列计算π而使得该方法为大家所知，但这种方法最早是由印度科学家在大约1400到1500年之间发现。第一个记载用无穷级数计算π的人是约西元1500年左右时，印度天文学家尼拉卡莎·萨默亚士在他的著作《系统汇编》中用梵语诗所记录。当时没有这数列对应的证明，而证明出现在另一本较晚的印度作品《基本原理》，年代约在西元1530年。尼拉卡莎将该数列归功于更早期的印度数学家桑加马格拉马的马德哈瓦（1350–1425）。相关的无穷级数有许多，包括有关$\sin$、$\tan$及$\cos$的，现在称为马德哈瓦数列或π的莱布尼茨公式。玛达瓦在1400年用无穷级数计算π到第11位小数，但在1430年一位波斯数学家卡西利用多边形算法否定了他算的结果。

艾萨克·牛顿利用无穷级数计算π到第15位，后来写道：“我很羞愧的告诉你我为了这个计算用了多少个数字。”

欧洲发现的第一条无穷项圆周率公式是无穷乘积（和一般用来计算π的无穷级数不同），由法国科学家弗朗索瓦·韦达在1593年发现：

${\frac {2}{\pi }}={\frac {\sqrt {2}}{2}}\cdot {\frac {\sqrt {2+{\sqrt {2}}}}{2}}\cdot {\frac {\sqrt {2+{\sqrt {2+{\sqrt {2}}}}}}{2}}\cdots$

约翰·沃利斯在1655年发现了沃利斯乘积，是欧洲发现的第二条无穷项圆周率公式：

${\frac {\pi }{2}}={\frac {2}{1}}\cdot {\frac {2}{3}}\cdot {\frac {4}{3}}\cdot {\frac {4}{5}}\cdot {\frac {6}{5}}\cdot {\frac {6}{7}}\cdot {\frac {8}{7}}\cdot {\frac {8}{9}}\cdots$

微积分学由英国科学家艾萨克·牛顿及德国数学家戈特弗里德·莱布尼茨在1660年代发明，许多计π的无穷级数出现。牛顿自己就用反正弦（$\arcsin$）数列在1655年或1666年将π近似到第15位小数，后来写到“我很羞愧告诉你我为了计算它用了多少数字，我当时没有做其他事。”

苏格兰数学家詹姆斯·格雷果里在1671年发现了马德哈瓦公式，莱布尼茨也在1674年发现：

$\arctan z=z-{\frac {z^{3}}{3}}+{\frac {z^{5}}{5}}-{\frac {z^{7}}{7}}+\cdots$

这公式即为格雷果里－莱布尼茨公式，在$z=1$时数值为${\frac {\pi }{4}}$。1699年时英国数学家亚伯拉罕·夏普用格雷果里－莱布尼茨公式，在$z={\frac {1}{\sqrt {3}}}$时计算，计算到π的第71位小数，打破由多边形算法得到的第39位小数的记录。格雷果里－莱布尼茨公式在$z=1$时非常简单，但收敛到最终值的速度非常慢，现在不会再用此公式来计π。

约翰·梅钦在1706年用格雷果里－莱布尼茨级数产生了可以快速收敛的公式：

${\frac {\pi }{4}}=4\,\arctan {\frac {1}{5}}-\arctan {\frac {1}{239}}$

梅钦用这公式计到π第100位小数后来其他数学家也发展了一些类似公式，现在称为梅钦类公式，创下了许多计算π位数的纪录。在进入电脑时代时，梅钦类公式仍然是耳熟能详可以计算π的公式，而且在约250年的时间里，很多有关π位数的纪录都是梅钦类公式所得，比如在1946年时由达尼尔·弗格森（Daniel Ferguson）用这类公式计到第620位小数，是没有计算设备辅助的最佳纪录。

1844年，计算天才扎卡里亚斯·达斯在德国数学家卡尔·弗里德里希·高斯的要求下以梅钦类公式心算了π的200位小数，并创下纪录。英国数学家威廉·谢克斯花了15年的时间计算π到小数707位，不过第528位小数出错，后面的小数也都不正确。

收敛速度

有些π的无穷级数收敛的比其他级数要快，数学家一般会选用收敛速度较快的级数，可以在较少的计算量下计算π，且达到需要的准确度。以下是π的莱布尼茨公式：

$\pi ={\frac {4}{1}}-{\frac {4}{3}}+{\frac {4}{5}}-{\frac {4}{7}}+{\frac {4}{9}}-{\frac {4}{11}}+{\frac {4}{13}}-\cdots$

随著一项一项的值加入总和中，只要项次够多，总和最后会慢慢接近π。不过此数列的收敛速度很慢，要到50万项之后，才会精确到π的第五位小数。

尼拉卡莎在15世纪发展了π的另一条无穷级数，收敛速度比格雷果里－莱布尼茨公式快很多：

$\pi =3+{\frac {4}{2\times 3\times 4}}-{\frac {4}{4\times 5\times 6}}+{\frac {4}{6\times 7\times 8}}-{\frac {4}{8\times 9\times 10}}+\cdots$

以下比较两条级数的收敛速率：

π的无穷级数	第1项	前2项	前3项	前4项	前5项	收敛到
$\pi ={\frac {4}{1}}-{\frac {4}{3}}+{\frac {4}{5}}-{\frac {4}{7}}+{\frac {4}{9}}-{\frac {4}{11}}+{\frac {4}{13}}\cdots .$	4.0000	2.6666…	3.4666…	2.8952…	3.3396…	3.1415…
$\pi ={3}+{\frac {4}{2\times 3\times 4}}-{\frac {4}{4\times 5\times 6}}+{\frac {4}{6\times 7\times 8}}\cdots .$	3.0000	3.1666…	3.1333…	3.1452…	3.1396…

计算前五项后，格雷果里－莱布尼茨级数的和跟π的误差为0.2，而尼拉卡莎级数和的误差为0.002。尼拉卡莎级数收敛快很多，也甚为适合用来计π的值。收敛更快的级数有梅钦类公式及楚德诺夫斯基算法，后者每计一项就可以得到14位正确的小数位。

无理与超越性

参见：林德曼-魏尔斯特拉斯定理

并非所有和π有关的研究都旨在提高计算它的准确度。1735年，欧拉解决了巴塞尔问题，建立了所有平方数倒数和与π的关系。之后欧拉发现了欧拉乘积公式，得到了π、素数的重要关联，对日后黎曼ζ函数的研究影响深远。

${\frac {\pi ^{2}}{6}}={\frac {1}{1^{2}}}+{\frac {1}{2^{2}}}+{\frac {1}{3^{2}}}+{\frac {1}{4^{2}}}+\cdots$

1761年，瑞士数学家约翰·海因里希·朗伯用正切函数的无穷连分数表达式证明了π是无理数。1794年，法国数学家阿德里安－马里·勒让德证明了$\pi ^{2}$也是无理数。1882年，德国数学家费迪南德·冯·林德曼证明了对任何非零代数数$\alpha$，$e^{\alpha }$都是超越数，该结论后来由魏尔斯特拉斯推广为林德曼－魏尔斯特拉斯定理。据此定理和欧拉公式，π只能是超越数，进而证实了勒让德和欧拉提出的π超越性猜想。哈代在其著作《数论导引》中则称此证明在提出后，经过希尔伯特、施瓦兹和其他一些人化简过。

引入π符号

莱昂哈德·欧拉在他1736年到1748年的作品中开始用希腊字母π表示圆周率，数学界也开始广为使用

在用π专指“圆周率”之前，希腊字母即已用于几何概念中。威廉·奥特雷德在1647年起在《数学之钥》（Clavis Mathematicae）就已经用$\pi$及$\delta$（对应p和d的希腊字母）来表示圆的周长及直径的比例。

威廉·琼斯在他1706年出版的《新数学导论》（A New Introduction to the Mathematics）提到了$\pi$，是目前已知最早专门用希腊字母$\pi$表示圆周和其直径比例的人。这希腊字母第一次出现是在书中讨论一块半径1的圆时提到“其圆周长一半（$\pi$）”。琼斯选用$\pi$可能因它是希腊文“周边”一词“περιφέρεια”的首字母。不过琼斯提到，他那些有关$\pi$的算式出自“真正聪明的约翰·梅钦先生”，人们推测在琼斯之前，约翰·梅钦就已开始用$\pi$表示圆周率。

琼斯在1706年开始使用此希腊字母，但直到莱昂哈德·欧拉在其1736年出版的《力学》中开始使用之后，其他数学家才纷纷开始用$\pi$指代圆周率。在此之前，数字家可能用像c或p之类的字母代表圆周率。欧拉与欧洲其他数学家间时常互相写信来往，$\pi$的用法迅速传播开来。1748年欧拉在他的《无穷小分析引论》再一次提到了$\pi$，写道：“简洁起见，我们将此数字写为$\pi$，$\pi$等于半径为1的圆周长的一半。”这表示方式之后也推展到整片西方世界。

现代数值近似

计算机时代与迭代算法

约翰·冯·诺伊曼所属的团队是用数位计算机ENIAC来计π的第一队

高斯－勒让德算法：
一开始设定

$\scriptstyle a_{0}=1\quad b_{0}={\frac {1}{\sqrt {2}}}\quad t_{0}={\frac {1}{4}}\quad p_{0}=1$

迭代计算:$\scriptstyle a_{n+1}={\frac {a_{n}+b_{n}}{2}}\quad \quad b_{n+1}={\sqrt {a_{n}b_{n}}}$

$\scriptstyle t_{n+1}=t_{n}-p_{n}(a_{n}-a_{n+1})^{2}\quad \quad p_{n+1}=2p_{n}$

则π的估计值为

$\scriptstyle \pi \approx {\frac {(a_{n}+b_{n})^{2}}{4t_{n}}}$

二十世纪中期计算机技术发展、革新再次引发了计算π位数的热潮。美国数学家约翰·伦奇及李维·史密斯在1949年用桌上型计算机计算到1120位。同年，乔治·韦斯纳（George Reitwiesner）及约翰·冯·诺伊曼带领的团队利用反三角函数（arctan）的无穷级数，用ENIAC计算到了小数后2037位，花了70小时的电脑工作时间。这纪录后来多次由其他透过arctan级数计算出的结果打破（1957年到7480位小数，1958年到第一万位数，1961年到第十万位小数），直到1973年，小数点后第一百万位小数经已算出。

1980年代有两项发明加速计算了π。第一项是发现了新的迭代法去计π的值，计算速度比无穷级数快很多；另一项是发现了可以快速计算大数字乘积的乘法演算法。电脑大部分的工作时间都是在计乘法，这类演算法对现代计π格外重要。这类演算法包括嘉良对马（Karatsuba）算法、谭曲（Toom－Cook）乘法及以傅里叶变换为基础的乘法演算法（傅里叶乘法）。

迭代演算法最早是在1975年至1976年间分别由美国物理学家尤金·萨拉明及奥地利科学家理查·布兰特独立提出。这两条演算法没有依赖无穷级数来计算。迭代会重复特定计算，将前一次的计算结果作为这一次的输入值，使得计算结果渐渐的趋近理想值。此方式的原始版本其实是在160年前由卡尔·弗里德里希·高斯提出，现在称为算术－几何平均数算法（AGM法）或高斯－勒让德算法。萨拉明及布兰特都曾修改之，这算法也称为萨拉明－布兰特演算法。

迭代演算法收敛速度比无穷级数快很多，在1980年代以后广为使用。无穷级数随著项次的增加，一般来说正确的位数也会增加几位，但迭代演算法每计算多一次，正确位数会呈几何级数增长。例如萨拉明－布兰特演算法每计算多一次，正确位数会是之前的二倍。1984年加拿大人乔纳森·波温及彼得·波温提出迭代演算法，每计算多一次，正确位数会是之前的四倍，1987年时有另一条迭代演算法，每计算多一次，正确位数会是之前的五倍。日本数学家金田康正使用的演算法在1955年及2002年间创下了若干项纪录。不过迭代演算法的快速收敛也有其代价，需要的记忆体明显比无穷级数多。

计算π的意义

当数学家发现新的算法、电脑变得普及时，π的已知小数位急剧增加。注意垂直坐标使用了对数坐标。

一般而言，π值并不需要过于精确便能够满足大部分数学运算的需求。按照约·安（Jörg Arndt）及古里斯佗夫·希奴（Christoph Haenel）的计算，39位精确度已可将可观测宇宙圆周的精确度准确至一粒原子大小，足以运算绝大多数宇宙学的计算需求。尽管如此，和π有关的成就往往成为世界各地的新闻头条；部分人出于对破纪录的冲动，依然奋力算出π小数点后上千甚至上百万位。此外也有测试超级计算机、测试数值分析算法（包括高精度乘法算法）等实际好处。纯粹数学这领域也能计算π的位数评定其随机度。

快速收敛级数

斯里尼瓦瑟·拉马努金的肖像，他在印度独立工作时提出了许多计算π的新颖数列。

现代计算π的程序不仅局限于迭代算法。20世纪80与90年代，出现了可用来计算π的新无穷级数，其收敛速度可与迭代算法媲美，而又有着复杂度、内存密集度更低的优势。印度数学家斯里尼瓦瑟·拉马努金是这方面的先驱，他在1914年发表了许多与π相关的公式，这些公式十分新颖，极为优雅而又颇具数学深度，收敛速度也非常快。下式即为一例，其中用到了模方程：

${\frac {1}{\pi }}={\frac {2{\sqrt {2}}}{9801}}\sum _{k=0}^{\infty }{\frac {(4k)!(1103+26390k)}{k!^{4}(396^{4k})}}.$

这无穷级数收敛速度远快于绝大多数反正切数列，包括梅钦公式。第一位使用拉马努金公式计算π并取得进展的是比尔·高斯珀，他在1985年算得了小数点后一千七百万位。拉马努金公式开创了现代数值近似算法的先河，此后波尔文兄弟和楚德诺夫斯基兄弟进一步发展了这类算法。后者于1987年提出了楚德诺夫斯基公式，如下所示：

${\frac {1}{\pi }}={\frac {12}{640320^{3/2}}}\sum _{k=0}^{\infty }{\frac {(6k)!(13591409+545140134k)}{(3k)!(k!)^{3}(-640320)^{3k}}}.$

此公式每计算一项就能得到π的约14位数值，因而用于突破圆周率的数位的计算。利用这公式，楚德诺夫斯基兄弟于1989年算得π小数点后10亿（10⁹）位，法布里斯·贝拉于2009年算得2.7千亿（2.7×10¹²）位，亚历山大·易和近藤滋在2011年算得一万亿（10¹³）位。类似的公式还有拉马努金-佐藤级数。

2006年，加拿大数学家西蒙·普劳夫利用PSLQ整数关系算法按照以下模版生成了几条计算π的新公式：

$\pi ^{k}=\sum _{n=1}^{\infty }{\frac {1}{n^{k}}}\left({\frac {a}{q^{n}-1}}+{\frac {b}{q^{2n}-1}}+{\frac {c}{q^{4n}-1}}\right)$，

$q$为e^$\pi$，$k$是奇数，$a,b,c$是普劳夫计算出的有理常数。

统计模拟法

布丰投针问题，多枚长度为ℓ的针随机地抛掷向平面。

随机地往内切四分之一圆的正方形内抛掷大量的点。

蒙特卡洛方法基于随机试验结果计算$\pi$的近似值

统计模拟法是以概率统计理论为指导的一类非常重要的计数方法，经大量重复试验计算事件发生频率，按照大数定律（即当试验次数充分大时，频率充分接近概率）可以求得$\pi$的近似值。 布芬（Buffon）投针问题就是其中一项实例：长度$l$的针随机往画满间距$t\left(l\leq t\right)$的平行线的平面上抛掷$n$次， 如果针与平行直线相交$m$次，$n$充分大就可根据以下公式算出$\pi$的近似值：

$\pi \approx {\frac {2n\ell }{mt}}$

用统计模拟法计$\pi$的另一例子是随机往内切四分之一圆的正方形内抛掷大量点，落在四分之一圆内的点的数量与抛掷点的总量的比值会近似于${\textstyle {\frac {\pi }{4}}}$。

此外还可用随机游走试验，并用统计模拟法计算$\pi$值，如抛掷一枚均匀的硬币$N$次，并记录正面朝上的次数，所得结果中，正面朝上的次数$n_{N}$服从二项分布且

$\Pr(n_{N}=m)={\binom {N}{m}}({\frac {1}{2}})^{m}({\frac {1}{2}})^{N-m}$

因为硬币均匀，所以N次试验中每次试验结果相互独立。由此可定义一系列独立的随机变量$X_{k}\left(k=1,2,\ldots \right)$，当抛掷结果为正面时$X_{k}=1$否则为-1，且$X_{k}=\pm 1$且取何值有相同概率（即，正面朝上和背面朝上的概率相同）。对随机变量$X_{k}\left(k=1,2,\ldots ,N\right)$求和可得

$W_{N}=\sum _{k=1}^{N}X_{k}$

设k为“硬币正面朝上的次数”减去“硬币反面朝上的次数”，即可得到$m-\left(N-m\right)=k$。变换式子，得$m={\frac {N+k}{2}}$，因此

$\Pr(W_{N}=k)={\binom {N}{\frac {N+k}{2}}}{\frac {1}{2^{N}}}$，其中$k=-N,-N+2,-N+4,\ldots ,N-2,N$。

可证明，

$E(W_{N})=0$，$E(W_{N}^{2})=N$，以及${\textstyle E(|W_{N}|)={\binom {N}{\left\lceil {N/2}\right\rceil {\frac {\left\lceil {N/2}\right\rceil }{2^{N-1}}}}}={\begin{cases}{\frac {(N-1)!!}{(N-2)!!}},&{\text{若 }}N{\text{偶，}}\\{\frac {N!!}{(N-1)!!}},&{\mbox{若 }}N{\mbox{奇。}}\end{cases}}}$

并且当N变大时，${\textstyle E\left(\left\vert W_{N}\right\vert \right)}$的值会渐近于${\textstyle {\sqrt {\frac {2N}{\pi }}}}$，因此当N充分大时可根据以下公式算出$\pi$的近似值：

$\pi \approx {\frac {2N}{|W_{N}|^{2}}}$

和其他计算$\pi$值的方法相比，蒙特卡洛方法收敛速度很慢，而且无论实验多少次，都无从得知$\pi$的估值已经精确到第几位。因此，当追求速度或精度时，蒙特卡洛方法不适合用来估计$\pi$。

阀门算法

1995年引入的两条算法开辟了研究$\pi$的新途径。因为每计算出一位数字，该数就会像流过阀门的水一样不会再出现在后续的计算过程中，这种新进算法叫阀门算法。这就与无穷级数及迭代算法形成对比——无穷级数和迭代算法自始至终的每一步计算都会涉及到之前所有步骤计算出的中间值。

1995年，美国数学家斯坦·瓦格纳和斯坦利·拉比诺维茨（Stanley Rabinowitz）发明了一种简单的阀门算法，其运算速度类似arctan演算法，但速度比迭代算法慢。

贝利－波尔温－普劳夫公式（BBP）是另一条阀门算法，属于一种位数萃取演算法。1995年，西蒙·普劳夫等人发现

$\pi =\sum _{k=0}^{\infty }{\frac {1}{16^{k}}}\left({\frac {4}{8k+1}}-{\frac {2}{8k+4}}-{\frac {1}{8k+5}}-{\frac {1}{8k+6}}\right)$

这公式和其他公式不同，可以计算$\pi$的任何十六进小数位，而不用计算前面全部小数位。十六进数位可计算得到特定二进数位；想要得到八进制数位的话，计算一、两位十六进小数即可。目前也已发现一些这种演算法的变体，不过还没有发现针对十进制、可以快速生成特定小数位的位数萃取演算法。位数萃取演算法的一项重要用途是用来确认声称是计算到$\pi$小数位数的新纪录：若有声称是新纪录的计算结果出现，先将十进制的数值转换到十六进制，再用贝利－波尔温－普劳夫公式去确认最后一些位数（用乱数决定），若这些位数都对，就能有一定把握认为此计算结果是对的。

1998年到2000年间，分布式计算计画PiHex用贝拉公式（贝利－波尔温－普劳夫公式的一种变体）计算$\pi$第10¹⁵位，结果是0。2010年9月，有雅虎员工用公司的Apache Hadoop应用程式在上千台电脑计算π在2×10¹⁵位开始往后256位，其第2×10¹⁵位刚好也是0。

利用伽玛函数计算

伽玛函数，$\Gamma (n)=(n-1)!$，可以被用作计算圆周率。

$\pi =\Gamma ({\frac {1}{2}})^{2}$

用途

$\pi$与圆密切相关，出现在许多几何学和三角学的公式中（特别是与圆、椭圆和球体相关的那些）。 此外，$\pi$也出现在其他学科的重要公式中，比如统计学、物理学，傅立叶分析和数论的公式。

几何学与三角学

圆的面积等于$\pi$乘以阴影部分面积。

$\pi$出现在基于圆的几何图形（如椭圆、球、圆锥与环面）的面积、体积公式中。下面是一些用到π的常见公式：

• 半径$r$的圆周长$2\pi r$。
• 半径$r$的圆面积$\pi r^{2}$。
• 半径$r$的球体积${\textstyle {\frac {4}{3}}\pi r^{3}}$。
• 半径$r$的球面面积$4\pi r^{2}$。

上述公式是n维球的体积与其边界（（n−1）维球的球面）的表面积的特殊情况，具体将在后文给出解释。

描述由圆生成的图形的周长、面积或体积的定积分常涉及π。例如，表示半径为1的半圆的面积的积分为

$\int _{-1}^{1}{\sqrt {1-x^{2}}}\,dx={\frac {\pi }{2}}.$

${\textstyle {\sqrt {1-x^{2}}}}$的积分表示上半圆（此处的平方根由勾股定理得出），从-1到1的积分${\textstyle \int _{-1}^{x}}$可用来计算计算半圆与x 轴间的面积。

正弦和余弦函数的重复周期为 2π。

三角函数要用到角，而数学家常用弧度作角度单位。π在弧度制起重要作用，数学家将周角，即360度定义为2π度。由这条定义可得，180度＝π弧度，1度＝${\textstyle {\frac {\pi }{180^{\circ }}}}$弧度。因此，常用的三角函数的周期为$\pi$的倍数；例如，正弦和余弦周期为π，任何角度$\theta$和任何整数$k$都有

${\textstyle \sin \theta =\sin \left(\theta +2\pi k\right)}$及${\textstyle \cos \theta =\cos \left(\theta +2\pi k\right)}$。

拓扑学

克莱因四次曲面的单值化，亏格为3且欧拉特征值为−4的面，作为双曲面与菲诺平面的对称群PSL(2,7)的商。根据高斯－博内定理，基本域的双曲面积为8π.

常数$\pi$出现在将平面微分几何及其拓扑学联系起来的高斯－博内定理中。具体来说，如果紧曲面Σ的高斯曲率为$K$，那么有

$\int _{\Sigma }K\,dA=2\pi \chi (\Sigma )$，

其中$\chi (\Sigma )$是该曲面的欧拉示性数，是整数。例如，曲率为1（也就是说其曲率半径也为1，对于球面而言此时的曲率半径与半径重合）的球面$S$的表面积。球面的欧拉特征数可以通过其同源组计算，其结果为2。于是，便得出

$A(S)=\int _{S}1\,dA=2\pi \cdot 2=4\pi$

即为半径为1的球面的表面积公式。

常数$\pi$还出现在拓扑学的许多其他的积分公式中，特别是那些涉及通过陈－韦伊同态的特征类。

向量分析

向量分析的方法可以通过分解成球谐函数来理解（图示）

向量分析是与向量场的性质有关的微积分的分支，并有许多物理用途，例如用在电磁学中。位于三维笛卡尔坐标系原点的点源$Q$的牛顿位势为

$V(\mathbf {x} )=-{\frac {kQ}{|\mathbf {x} |}}$

表示位于距原点$\left\vert {\boldsymbol {x}}\right\vert$的单位质量（或电荷）的势能，而$k$是维度常数。在这里由$\mathrm {E}$表示的场可以是（牛顿）引力场或（库仑）电场，是位势的负梯度：

$\mathbf {E} =-\nabla V.$

特殊情况有库仑定律和牛顿万有引力定律。高斯定律表明，通过包含原点的任何平滑、简单、封闭、可定向曲面$S$的场的向外通量等于$4\pi kQ$：

$4\pi kQ=$$\oiint$${\scriptstyle S}$$\mathbf {E} \cdot d\mathbf {A}$

标准形式会将$4\pi$的这因子吸收到常数$k$中，但这种说法表明了它必须出现在“某处”。此外，$4\pi$是单位球面的表面积，但并没有假设$S$是球面。然而，作为散度定理的结果，由于远离原点的区域是真空（无源的），只有$\mathrm {R} ^{3}\setminus \left\{0\right\}$中的表面$S$的同调类与计算积分有关，因此可以由相同同调类中的任何方便的表面代替，特别是球形，因为球面坐标可以用于计算积分。

高斯定律的结果之一是位势$V$的负拉普拉斯算子等于狄拉克δ函数的$4\pi kQ$倍：

$\Delta V(\mathbf {x} )=-4\pi kQ\delta (\mathbf {x} ).$

通过卷积就能得到物质（或电荷）的更一般分布，给出泊松方程

$\Delta V(\mathbf {x} )=-4\pi k\rho (\mathbf {x} )$

其中$\rho$是分布函数。

爱因斯坦方程表明，时空的曲率是由其中的物质能量生成

常数$\pi$在与爱因斯坦场方程中的四维势起类似的作用，爱因斯坦方程是形成广义相对论基础的一条基本公式，并且把引力的基本相互作用描述为物质和能量引起的时空弯曲的结果：

$R_{\mu \nu }-{\frac {1}{2}}Rg_{\mu \nu }+\Lambda g_{\mu \nu }={\frac {8\pi G}{c^{4}}}T_{\mu \nu },$

$R_{\mu v}$是里奇曲率张量，$R$是数量曲率，$g_{\mu v}$是度量张量，$\Lambda$是宇宙学常数，$G$是万有引力常数，$c$是真空中的光速，而$T_{\mu v}$是应力－能量张量。爱因斯坦方程的左边是度量张量的拉普拉斯算子的非线性模拟，并化简（reduce）至在弱域的极限，而右边是分布函数的模拟乘以$8\pi$。

柯西积分公式

复杂的解析函数可以以一系列的流线和等电位线（许多以直角相交的曲线）视觉化，图中是伽玛函数的复数对数。

在复分析中，沿复平面若尔当曲线的围道积分是研究解析函数的重要手段之一。简化版的柯西积分公式表明，对任何若尔当曲线$\gamma$内任一点$z_{0}$，以下围道积分给出$2\pi i$：

$\oint _{\gamma }{\frac {dz}{z-z_{0}}}=2\pi i.$

该命题是柯西积分定理的直接推论，后者表明上述围道积分在围道的同伦变换下保持不变，因而沿任一曲线的积分和沿以$z_{0}$为圆心的圆周积分的结果相同。更为一般地，该公式对不通过$z_{0}$点的任意可求长曲线都成立，但等式右边要乘以曲线关于该点的卷绕数。

一般形式的柯西积分公式建立了全纯函数$f(z)$在若尔当曲线$\gamma$上的值与曲线内任意点$z_{0}$处值的关系：

$\oint _{\gamma }{f(z) \over z-z_{0}}\,dz=2\pi if(z_{0})$

柯西积分定理是留数定理的一项特例。根据留数定理，在区域内除去有限个点解析的亚纯函数$g(z)$在边界上的围道积分与函数在这些点的留数之和满足：

$\oint _{\gamma }g(z)\,dz=2\pi i\sum \operatorname {Res} (g,a_{k})$

Γ函数与斯特灵公式

以维拉瑞索圆将三块球面霍普夫纤维化，下方是富比尼－施图迪度量的黎曼球面与其富比尼－施图迪度量（如图所示的三块平行曲面）。恒等式S₃(1)/S₂(1) = π/2可以确定一条数列

阶乘函数$n!$的值等于所有小于等于$n$的正整数之积，它的定义域只包含非负整数。Γ函数则是阶乘的推广。它在复平面的右半平面定义为：

$\Gamma (z)=\int _{0}^{\infty }{\frac {t^{z-1}}{\mathrm {e} ^{t}}}\,{\rm {d}}t$

再利用解析延拓可以将它的定义域扩展到除去非正整数的整块复数域。当自变量$z=n$取正整数时，$\Gamma$函数给出阶乘$\left(n-1\right)!$；当自变量取半整数时，计算结果含有$\pi$。例如${\textstyle \Gamma (1/2)={\sqrt {\pi }}}$，${\textstyle \Gamma (5/2)={\frac {3{\sqrt {\pi }}}{4}}}$。

根据魏尔施特拉斯分解定理，$\Gamma$函数可分解为如下的无穷乘积：

$\Gamma (z)=e^{-\gamma z}\prod _{n=1}^{\infty }{\frac {e^{z/n}}{1+z/n}}$

$\gamma$是欧拉－马斯刻若尼常数。利用该分解公式和$\Gamma$函数在${\textstyle z={\frac {1}{2}}}$的值${\textstyle \Gamma \left({\frac {1}{2}}\right)^{2}=\pi }$，亦可以证明沃利斯乘积式。$\Gamma$函数和黎曼ζ函数、函数行列式的恒等式存在关联，其中$\pi$扮演着重要的角色。

$\Gamma$函数常用于计算$n$维欧氏空间中n 维球的体积和n 维球面的表面积。对$n$维欧氏空间中半径为$r$$n$维球，其体积$V_{n}(r)$和表面积$S_{n-1}(r)$满足：

$V_{n}(r)={\frac {\pi ^{n/2}}{\Gamma ({\frac {n}{2}}+1)}}r^{n}$

$S_{n-1}(r)={\frac {n\pi ^{n/2}}{\Gamma ({\frac {n}{2}}+1)}}r^{n-1}$

两者还满足如下的关系式：

$2\pi r={\frac {S_{n+1}(r)}{V_{n}(r)}}.$

当$n$很大，用$\Gamma$函数可得到阶乘$n!$的近似公式$n!\sim {\sqrt {2\pi n}}\left({\frac {n}{e}}\right)^{n}$，称斯特灵公式，等价于：

$\pi =\lim _{n\to \infty }{\frac {e^{2n}n!^{2}}{2n^{2n+1}}}.$

斯特灵近似的几何应用之一是埃尔哈特体积猜想。将$n$维欧几里得空间的单纯形记作$\Delta _{n}$，$\left(n+1\right)\Delta _{n}$则表示该单纯形的所有面扩大$n+1$。于是

$\operatorname {Vol} ((n+1)\Delta _{n})={\frac {(n+1)^{n}}{n!}}\sim {\frac {e^{n+1}}{\sqrt {2\pi n}}}.$

这是仅含一点晶格点之凸体体积的（最佳）上界。

数论与黎曼ζ函数

全部质数都有其关联的普鲁法群，即圆的算数定域。分析数论的L函数也定域在每个质数p上

基于韦伊猜想的巴塞尔问题的解：$\zeta (2)$ 的数值是模群中一个基本域的双曲面积的2${\textstyle \pi }$倍。

黎曼ζ函数$\zeta (s)$ 在数学的许多领域均有应用。当自变量$s=2$ ，可写作

$\zeta (2)={\frac {1}{1^{2}}}+{\frac {1}{2^{2}}}+{\frac {1}{3^{2}}}+\cdots$

找到这无穷级数的解析解是数学界著名的“巴塞尔问题”。1735年，欧拉解决了这问题，他得到该无穷级数等于${\frac {\pi ^{2}}{6}}$。欧拉的结论可推导出数论中一项结果，即两随机整数互质（无公因数）的概率为${\frac {6}{\pi ^{2}}}$ 。整数可由质数$p$整除的概率为${\frac {1}{p^{2}}}$（例如，连续7个正整数只有一个可以7整除），任取两随机整数都能以质数$p$整除的概率为${\frac {1}{p^{2}}}$，至少有一数不能整除的概率则为$1-{\frac {1}{p^{2}}}$。又，一随机整数能否以两不同质数整除是相互独立事件，两随机整数互质的概率可以表示成关于所有质数$p$的无穷乘积

${\begin{aligned}\prod _{p}^{\infty }\left(1-{\frac {1}{p^{2}}}\right)&=\left(\prod _{p}^{\infty }{\frac {1}{1-p^{-2}}}\right)^{-1}\\&={\frac {1}{1+{\frac {1}{2^{2}}}+{\frac {1}{3^{2}}}+\cdots }}\\&={\frac {1}{\zeta (2)}}={\frac {6}{\pi ^{2}}}\approx 61\%.\end{aligned}}$

这结论可结合随机数生成器，用统计模拟法计$\pi$的近似值。

巴塞尔问题的结论意味着几何导出量$\pi$的数值与质数的分布有着深刻的关联。巴塞尔问题是谷山－志村定理的一种特殊情况，是安德烈·韦伊对玉河数的猜想的一项特例，即猜想一个这种形式的算术量关于所有质数$p$的无穷乘积能够等于一个几何量——某局部对称空间体积的倒易。巴塞尔问题中，这空间是双曲3-流形SL₂(R)/SL₂(Z)。

$\zeta$函数同样满足黎曼方程的公式，其中用到了$\pi$和伽玛公式：

$\zeta (s)=2^{s}\pi ^{s-1}\ \sin \left({\frac {\pi s}{2}}\right)\ \Gamma (1-s)\ \zeta (1-s)\!.$

除此之外， $\zeta$函数导数也满足

$\exp(-\zeta '(0))={\sqrt {2\pi }}.$

最终的结果是$\pi$可以从谐振子泛函行列式中求得。这泛函行列式可以无穷乘积展开式计算，而且这种方法等价于沃利斯乘积公式。这种方法可用于量子力学，尤其是玻尔模型中的变分。

傅里叶级数

π出现在P进数中的表示（如图），它们是普鲁法群的元素。泰特的论文很大程度地利用了这系统。

周期函数的傅里叶级数很自然出现了$\pi$。周期函数即实数的小数部分所构成群${\textstyle T={\frac {R}{Z}}}$上的函数。傅里叶分解指出，$T$上的复值函数$f$可表示为无穷多个$T$的酉特征的线性叠加之和。也就是说，$T$到圆群$U(1)$（模为1的复数组成的乘法群）的映射是连续群同态。$T$的特征都有$e_{n}(x)=e^{2\pi inx}$的形式，是一条定理。

$T$有唯一的特征值，直到复共轭，那是一群同态。在圆群用哈尔测度，常数$\pi$是这特征值的拉东－尼科迪姆导数值的一半。其他的特征值的导数值为$2\pi$的正整数倍。因此，常数$\pi$是独特的数字，以至于配备了其哈尔测度的群$T$，有对于$2\pi$整数倍的点阵的庞特里亚金对偶性。这是泊松和公式的一维版本。

模形式与𝜃函数

常数$\pi$与模形式和Θ函数密切相关——比如，椭圆曲线中的j变量就很大程度涉及楚德诺夫斯基算法（一种快速计算π的方法）。

模形式是以在上半平面的全纯函数的在模群$SL_{2}(\mathbb {Z} )$（或其子群，$SL_{2}(\mathbb {Z} )$是$SL_{2}(\mathbb {R} )$的一格）下的变换特性归纳。Θ函数便是一例：

$\theta (z,\tau )=\sum _{n=-\infty }^{\infty }e^{2\pi inz+i\pi n^{2}\tau }$

它是一种名为雅可比形式的模形式，有时以诺姆$q=e^{\pi i\tau }$表达。

常数$\pi$是特殊常数，它会使雅可比$\Theta$函数形成自守形式，即该函数会以特定方式变换。有若干恒等式在所有自守形式下成立。，例如：

$\theta (z+\tau ,\tau )=e^{-\pi i\tau -2\pi iz}\theta (z,\tau )$

它使得$\theta$必然在离散海森伯群下以表示（representation）变换。一般模形式和其他$\Theta$函数也包含$\pi$，这也是根据史东–冯纽曼定理。

柯西分布与位势论

箕舌线，英文名来自玛利亚·阿涅西（1718–1799），柯西分布的一幅几何构筑图

积分${\textstyle \int _{-\infty }^{\infty }{\frac {1}{x^{2}+1}}\,dx=\pi }$，柯西分布${\textstyle g(x)={\frac {1}{\pi }}\cdot {\frac {1}{x^{2}+1}}}$是概率密度函数，总概率等于1。

柯西分布的香农熵等于$\log(4\pi )$，也含$\pi$。

柯西分布控制做布朗运动的粒子通过膜的通道

柯西分布在位势论中扮演着重要的角色因为它是最简单的福斯坦堡测度和与在半平面上做布朗运动相关联的经典泊松核。共轭谐波函数以及希尔伯特变换与泊松核的渐近线有关。希尔伯特变换$H$是由奇异积分的柯西主值给出的积分变换

$Hf(t)={\frac {1}{\pi }}\int _{-\infty }^{\infty }{\frac {f(x)}{x-t}}$。

常数$\pi$是唯一的（正）归一化因子因此$H$定义了一个在实数轴上的平方可积分实值函数的希尔伯特空间上的线性复结构。 和傅里叶变换一样，希尔伯特变换就其在希尔伯特空间$L^{2}(R)$的变换特性而言可以完全特征化。直到归一化，它是唯一的与正膨胀对易且与实数轴的所有反射反对易有界线性算子。常数$\pi$是唯一能使这变换幺正的归一化因子。

复变动态系统

可以从曼德博集合中计算π，计算方式和计算从(−0.75, ε)点开始，一直到发散之前的次数有关

大卫·波（David Boll）在1991年发现在曼德博集合分形也有π出现。他检查在曼德博集合在$(-0.75,0)$位置的特性。若考虑坐标在“颈部”$(-0.75,\varepsilon )$的点，而$\varepsilon$趋近零，在发散之前迭代的次数和$\varepsilon$相乘，会趋近$\pi$。若是在右侧尖点处附近的点$(0.25,\varepsilon )$也会有类似的特性：在发散之前迭代的次数和$\varepsilon$的平方根相乘，也会趋近$\pi$。

数学之外的π

描述物理现象

$\pi$与圆以及球坐标系关系密切，即使$\pi$不是物理常数，也常出现在描述宇宙的基本原则方程中。比方说，经典力学领域的简单公式给出长L的单摆小幅摆动的近似周期${\textstyle T\approx 2\pi {\sqrt {\frac {L}{g}}}}$，$g$为地球引力加速度常数。

海森堡不确定性原理是量子力学的基本公式，表明测量粒子时，其位置不确定度（$\Delta x$）与动量不确定度（$\Delta p$）不可能同时达到任意小（$h$为普朗克常数）：${\textstyle \Delta x\,\Delta p\geq {\frac {h}{4\pi }}}$。

$\pi$近似三这特性，和电子偶素的半衰期相对较长有密切的联系。其半衰期的倒数和精细结构常数$\alpha$的关系为${\textstyle {\frac {1}{\tau }}=2{\frac {\pi ^{2}-9}{9\pi }}m\alpha ^{6}}$，$m$为电子质量。

许多结构工程的公式也有$\pi$，例如欧拉推导的挫曲公式说明了长度为$L$、截面二次轴矩为I的细长形物体，在不挫曲的条件下可以承受的最大轴向负载$F$：

$F={\frac {\pi ^{2}EI}{L^{2}}}.$

流体动力学的斯托克斯定律中也有$\pi$。斯托克斯定律是半径约为$R$的小球体在黏度$\eta$的流体中以速度$v$运动时会受到的阻力满足：

$F=6\,\pi \,\eta \,R\,v.$

在理想状态下，河的曲折程度（河道本身的长度与源头到入海口的比值）随着时间的推移逐渐趋向于$\pi$。河流外边缘的快速水流弯曲会使河流内边缘加倍侵蚀，河道变得更弯曲，整条河弯折更厉害。然而，这股弯折劲儿最终会导致河流折回一开始弯折的地方，导致“短路”，并形成河迹湖。这两种相反因素使河道长度与源头到入海口的比值的平均值为π。

π的记忆技巧

主条目：π文字学

π文字学（或译作圆周率的语言学）是指记住$\pi$的大量位值，并将其世界纪录载于健力士世界纪录大全的做法。维尔·美拿（Rajveer Meena）于2015年3月21日在印度于9小时27分钟内背诵了7万位的π，创下健力士世界纪录大全认证的世界纪录。2006年，日本退休工程师原口证自称已经背诵了十万位小数，但他未获健力士世界纪录大全认证。

常用于记忆π的一项技巧是背诵以单词长度代表$\pi$数值的故事或诗歌：第一单词有三字母，第二单词有一字母，第三单词有四字母，第四单词有一字母，第五单词有五字母，如此类推。早期例子是英国科学家詹姆士·金斯设计的诗歌：“How I want a drink, alcoholic of course, after the heavy lectures involving quantum mechanics.”这类诗歌有时在英文中称为“piem”。除了英文，用于记忆π的诗歌亦有不同语言的版本。但是，记忆$\pi$的人一般并不以诗歌记忆$\pi$来创下纪录，而是用如记忆数字规律或轨迹法的方法。

有好几位作家仿照上述记忆技巧，用$\pi$的数值创作了新型的约束写作方式，当中单词长度须符合$\pi$的数值。《The Cadaeic Cadenza》以上述技巧包含了$\pi$前3835位的值，一本标准长度的书《Not a Wake》有一万单词，其中各单词亦代表了$\pi$的一位。

大众文化

π派；圆形西式馅饼是常见的π双关语（英语圆周率和派同音）

也许因为$\pi$的公式很简短而且四处可见， $\pi$比其他数学常数在流行文化中更常见。

在2008年由英国公开大学及英国广播公司联合制作的记录片《数学的故事》于2008年十月由英国广播公司第四台播放。影片讲述了英国数学家马库斯·杜·索托伊在到访印度研究当地三角学的贡献时，展示出历史上π最精确的计算公式的信息图形。

巴黎的科学博物馆发现宫有间圆形房间叫“$\pi$房”，墙上刻有$\pi$的707位数，数字贴在圆顶状的天花板上，由大型的木制字符组成。数值是1853年由英国数学家威廉·尚克思计算出来，但是该结果于第528位后开始出现谬误，在1946年发现，1949年修正。

卡尔·萨根的小说《接触未来》中则暗示说，宇宙的创造者在π的数字中暗藏了一则信息。π的数字也用在凯特·布希所出的专辑Aerial中的《Pi》的歌词里。

美国人在3月14日庆祝圆周率日，此节日在学生中很流行。一些自称“数学极客”的人常常用$\pi$与其数位来创作一些数学或技术圈内人士才能领会到的笑话。麻省理工学院则有几句包含“3.14159”的大学欢呼口号。2015年的圆周率日格外重要，按照美式写法，当天的日期时间3/14/15 9:26:53较其他圆周率日包含更多位数的$\pi$。

在北电网络于2011年举行的技术专利拍卖会上，谷歌用了一些包含$\pi$在内的数学或科学常数来竞价。

在1958年，阿尔伯特·伊格尔提议将$\pi$换成τ（tau）以便简化公式。$\tau$在此定义为$\pi$的两倍。然而，没有任何其他作者曾这样使用过$\tau$。有些人使用不同的值，$\tau =2\pi =6.283185\ldots$。这些人称$\tau$不论是作为弧度制下圆周长的1转还是作为弧长与半径的比值（而不是与直径的比值）都比$\pi$自然，也能因此简化许多公式。有媒体报道称，因为$\tau$的值大小约为6.28，现已有人在6月28日庆祝“$\tau$节”，并吃“两个派”；然而，主流数学界还并未使用$\tau$。

1897年，有业余美国数学家试图藉印第安纳州议会来通过后世所谓印第安纳圆周率法案的法案。这法案试图以法律命令强制规定数学常数之值而臭名远播。该法案描述化圆为方的方法，并间接提到了$\pi$的错误值，例如3.2。该法案通过了印第安纳州众议院的表决，但参议院否决之。

参考资料

书籍

• Arndt, Jörg; Haenel, Christoph. Pi Unleashed. Springer-Verlag. 2006 [2013-06-05]. ISBN 978-3-540-66572-4.  English translation by Catriona and David Lischka.
• Ayers, Frank. Calculus. McGraw-Hill. 1964. ISBN 978-0-070-02653-7. 
• Berggren, Lennart; Borwein, Jonathan; Borwein, Peter. Pi: a Source Book. Springer-Verlag. 1997. ISBN 978-0-387-20571-7. 
• Beckmann, Peter. History of Pi. St. Martin's Press. 1989 [1974]. ISBN 978-0-88029-418-8. 
• Borwein, Jonathan; Borwein, Peter. Pi and the AGM: a Study in Analytic Number Theory and Computational Complexity. Wiley. 1987. ISBN 978-0-471-31515-5. 
• Boyer, Carl B.; Merzbach, Uta C. A History of Mathematics 2. Wiley. 1991. ISBN 978-0-471-54397-8. 
• Bronshteĭn, Ilia; Semendiaev, K. A. A Guide Book to Mathematics. Verlag Harri Deutsch. 1971. ISBN 978-3-871-44095-3. 
• Eymard, Pierre; Lafon, Jean Pierre. The Number Pi. American Mathematical Society. 1999. ISBN 978-0-8218-3246-2. , English translation by Stephen Wilson.
• Howe, Roger, On the role of the Heisenberg group in harmonic analysis, Bulletin of the American Mathematical Society, 1980, 3 (2) .
• Joseph, George Gheverghese. The Crest of the Peacock: Non-European Roots of Mathematics. Princeton University Press. 1991 [2013-06-05]. ISBN 978-0-691-13526-7. 
• Posamentier, Alfred S.; Lehmann, Ingmar. Pi: A Biography of the World's Most Mysterious Number. Prometheus Books. 2004. ISBN 978-1-59102-200-8. 
• Reitwiesner, George. An ENIAC Determination of pi and e to 2000 Decimal Places. Mathematical Tables and Other Aids to Computation. 1950, 4 (29): 11–15. doi:10.2307/2002695. 
• Roy, Ranjan. The Discovery of the Series Formula for pi by Leibniz, Gregory, and Nilakantha. Mathematics Magazine. 1990, 63 (5): 291–306. doi:10.2307/2690896. 
• Schepler, H. C. The Chronology of Pi. Mathematics Magazine (Mathematical Association of America). 1950, 23 (3): 165–170 (Jan/Feb), 216–228 (Mar/Apr), and 279–283 (May/Jun). doi:10.2307/3029284. . issue 3 Jan/Feb, issue 4 Mar/Apr, issue 5 May/Jun （页面存档备份，存于互联网档案馆）
• Thompson, William, Isoperimetrical problems, Nature series: Popular lectures and addresses (MacMillan and Company), 1894, II: 571–592 

延伸阅读

• Blatner, David. The Joy of Pi. Walker & Company. 1999. ISBN 978-0-8027-7562-7. 
• Borwein, Jonathan; Borwein, Peter. The Arithmetic-Geometric Mean and Fast Computation of Elementary Functions. SIAM Review. 1984, 26 (3): 351–365. doi:10.1137/1026073. 
• Borwein, Jonathan; Borwein, Peter; Bailey, David H. Ramanujan, Modular Equations, and Approximations to Pi or How to Compute One Billion Digits of Pi. The American Mathematical Monthly. 1989, 96 (3): 201–219. JSTOR 2325206. doi:10.2307/2325206. 
• Chudnovsky, David V. and Chudnovsky, Gregory V., "Approximations and Complex Multiplication According to Ramanujan", in Ramanujan Revisited (G.E. Andrews et al. Eds), Academic Press, 1988, pp 375–396, 468–472
• Cox, David A., "The Arithmetic-Geometric Mean of Gauss", L' Ensignement Mathematique, 30(1984) 275–330
• Delahaye, Jean-Paul, "Le Fascinant Nombre Pi", Paris: Bibliothèque Pour la Science (1997) ISBN 2902918259
• Engels, Hermann. Quadrature of the Circle in Ancient Egypt. Historia Mathematica. 1977, 4 (2): 137–140. doi:10.1016/0315-0860(77)90104-5. 
• Euler, Leonhard, "On the Use of the Discovered Fractions to Sum Infinite Series", in Introduction to Analysis of the Infinite. Book I, translated from the Latin by J. D. Blanton, Springer-Verlag, 1964, pp 137–153
• Hardy, G. H. and Wright E. M., An Introduction to the Theory of Numbers first published 1938, fifth Edition 1979 with additions 2000, Clarendon Press, Oxford U.K.
• Heath, T. L., The Works of Archimedes, Cambridge, 1897; reprinted in The Works of Archimedes with The Method of Archimedes, Dover, 1953, pp 91–98
• Huygens, Christiaan, "De Circuli Magnitudine Inventa", Christiani Hugenii Opera Varia I, Leiden 1724, pp 384–388
• Lay-Yong, Lam; Tian-Se, Ang. Circle Measurements in Ancient China. Historia Mathematica. 1986, 13 (4): 325–340. doi:10.1016/0315-0860(86)90055-8. 
• Lindemann, Ferdinand. Ueber die Zahl pi. Mathematische Annalen. 1882, 20 (2): 213–225 [2016-09-13]. doi:10.1007/bf01446522. （原始内容存档于2015-01-22）. 
• Matar, K. Mukunda; Rajagonal, C. On the Hindu Quadrature of the Circle" (Appendix by K. Balagangadharan). Journal of the Bombay Branch of the Royal Asiatic Society. 1944, 20: 77–82. 
• Niven, Ivan, "A Simple Proof that pi Is Irrational", Bulletin of the American Mathematical Society, 53:7 (July 1947), 507
• Ramanujan, Srinivasa, "Modular Equations and Approximations to π", Quarterly Journal of Pure and Applied Mathematics, XLV, 1914, 350–372. Reprinted in G.H. Hardy, P.V. Seshu Aiyar, and B. M. Wilson (eds), Srinivasa Ramanujan: Collected Papers, 1927 (reprinted 2000), pp 23–29
• Shanks, William, Contributions to Mathematics Comprising Chiefly of the Rectification of the Circle to 607 Places of Decimals, 1853, pp. i–xvi, 10
• Shanks, Daniel; Wrench, John William. Calculation of pi to 100,000 Decimals. Mathematics of Computation. 1962, 16 (77): 76–99. doi:10.1090/s0025-5718-1962-0136051-9. 
• Tropfke, Johannes, Geschichte Der Elementar-Mathematik in Systematischer Darstellung (The history of elementary mathematics), BiblioBazaar, 2009 (reprint), ISBN 978-1-113-08573-3
• Viete, Francois, Variorum de Rebus Mathematicis Reponsorum Liber VII. F. Viete, Opera Mathematica (reprint), Georg Olms Verlag, 1970, pp 398–401, 436–446
• Wagon, Stan, "Is Pi Normal?", The Mathematical Intelligencer, 7:3(1985) 65–67
• Wallis, John, Arithmetica Infinitorum, sive Nova Methodus Inquirendi in Curvilineorum Quadratum, aliaque difficiliora Matheseos Problemata, Oxford 1655–6. Reprinted in vol. 1 (pp 357–478) of Opera Mathematica, Oxford 1693
• Zebrowski, Ernest, A History of the Circle: Mathematical Reasoning and the Physical Universe, Rutgers University Press, 1999, ISBN 978-0-8135-2898-4

外部链接

• 开放目录项目中的“Digits of Pi”
• Wolfram Mathworld上的“Pi”（页面存档备份，存于互联网档案馆）
• Wolfram Alpha上的Representations of Pi（页面存档备份，存于互联网档案馆）
• Pi搜寻引擎 （页面存档备份，存于互联网档案馆）：可搜寻π、√2和e的20亿位
• Eaves, Laurence. π – Pi. Sixty Symbols. Brady Haran for the University of Nottingham. 2009 [2016-09-13]. （原始内容存档于2020-07-06）. 
• Grime, Dr. James. Pi is Beautiful – Numberphile. Numberphile. Brady Haran. 2014 [2016-09-13]. （原始内容存档于2020-12-12）. 
• Demonstration by Lambert (1761) of irrationality of π, online（页面存档备份，存于互联网档案馆） and analyzed BibNum（页面存档备份，存于互联网档案馆） (PDF).
• 寻找π值的计画 （页面存档备份，存于互联网档案馆）
• 前1.2百万位中的部分资料
• 将π前一万位化作音乐旋律 （页面存档备份，存于互联网档案馆）
• SuperPI（页面存档备份，存于互联网档案馆） 计算π值的软体，电脑硬体玩家常用来测试电脑运算速度（日文）
• 计算圆周率
• PiFast （页面存档备份，存于互联网档案馆） 个人电脑上最快的计算π值软体，是个人电脑计算π值纪录保持软体。
• 用大炮求圆周率 （页面存档备份，存于互联网档案馆）：以蒙特卡罗演算法，利用圆的面积求圆周率。
• 展示将π值计算到10万位的网站 （页面存档备份，存于互联网档案馆）