圓周率

圓周率是数学常数，为圆的周长和其直径的比，近似值约3.14159265，常用符号${\displaystyle \pi }$表示。

圓周率

圓周率
數表—无理数 ${\displaystyle \color {blue}{\sqrt {2}}}$　${\displaystyle \color {blue}\varphi }$　${\displaystyle \color {blue}{\sqrt {3}}}$　${\displaystyle \color {blue}{\sqrt {5}}}$　${\displaystyle \color {blue}\delta _{S}}$　${\displaystyle \color {blue}e}$　${\displaystyle \color {blue}\pi }$
命名
數字	π
名稱	圓周率
識別
種類	無理數 超越數
符號	${\displaystyle \pi }$
位數數列編號	A000796
性質
定義	${\displaystyle \pi ={\frac {C}{d}}.}$，其中${\displaystyle C}$為圓周長、${\displaystyle d}$為直徑 ${\displaystyle \pi =\int _{-1}^{1}{\frac {dx}{\sqrt {1-x^{2}}}}.}$
以此為根的多項式或函數	${\displaystyle e^{ix}+1=0}$
表示方式
值	${\displaystyle \pi \approx }$3.14159265
無窮級數	${\displaystyle \pi =\sum _{k=0}^{\infty }{\frac {4(-1)^{k}}{2k+1}}}$
二进制	11.00100100001111110110…^[1]
十进制	3.14159265358979323846…
十六进制	3.243F6A8885A308D31319…^[2]^:242
六十进制	3;8,29,44,0,47,25,53,7,24,57,36…^[3]^[4]
查 论 编

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

几个文明古国很早就須计算出π的精确值以便于生产的计算。西元5世纪，中國劉宋数学家祖冲之用几何方法将圆周率计算到小数点后7位。大约同时，印度数学家也将圆周率计算到小数点后5位。史上首條π的精确无穷级数公式（即π的莱布尼茨公式）直到约1000年后才由印度数学家发现。^[5][6]微積分出現，π的位數很快計到數百位，足以滿足任何科學工程的計算需求。在20和21世纪，计算机技术快速发展，π的计算精度急速提高。截至2022年6月，π的十进制精度已達1×10¹⁴位。^[7]当前计算π的值主要都为打破记录、测试超级计算机的计算能力和高精度乘法算法，因为几乎所有科学研究对π的精度要求都不超过几百位。^[2]:17[8]

π的定义涉及圆，所以π在三角学和几何学的许多公式，特别是在圆形、球形或椭球形相關公式中广泛应用。^[9]在更近代的數學分析裡，π改由實數系統譜性質中的特征值或週期來定義，不再指涉幾何。所以它也在一些和圓之幾何無甚相關的数学和科學領域中出現，像是數論、統計以及幾乎物理學中所有領域。π的广泛用途使它成为科学界内外最广为人知的常数。几本专门介绍π的书籍经已出版，圆周率日（3月14日）和π值计算突破记录也往往会成为报纸的新闻头条。^[10]此外，背诵π值的世界记录已达10萬位。^[11]

直徑為一的圓的周長是π

基本概念

名称

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

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

定义

 

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

π常用定义为圆的周长${\displaystyle C}$ 与直径${\displaystyle d}$ 的比值：^[2]:8

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

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

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

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

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

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

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

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

并且其中包括独特的正实数${\displaystyle \pi }$ 。^[15][21]

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

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

无理及正规性

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

統計隨機性（英语：statistical randomness）检验，包括正规数检验，可验证${\displaystyle \pi }$ 的位數沒有明顯的固定模式。${\displaystyle \pi }$ 的小数中任意固定长度的序列（如3位數000，001……999）出現機率都相同^[25]。不過有關π是正规数的猜想既無證明，亦無证伪^{[2]:22-23[25]}。

電腦出現後可生成大量π的不同位数，并統計分析之。金田康正詳細統計分析了π的十進制數字，并验证了其分布正规：例如，假設檢定0到9十個數的出現頻率，找不到有特定重复规律的證據^{[2]:22, 28–30}。根據無限猴子定理，任何任意長度、由隨機內容組成的子序列看起來都有可能像不隨機生成。因此，就算π的小数序列通過了隨機性統計測試，其中也可能有幾位的數字看起來似有规律可循而非隨機数，例如π的十進制写法在小數第762位后开始出现了連續六個9^[2]:3。

超越性

 

由于π是超越數，不能利用尺规作图化圓為方。

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

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

连分式

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

${\displaystyle \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、${\displaystyle {\frac {22}{7}}}$ 、${\displaystyle {\frac {333}{106}}}$ 、${\displaystyle {\frac {355}{113}}}$ 。这些数在历史上是π最广为人知且广為使用的几个近似值。用以上方式得出的${\displaystyle \pi }$ 的近似值要比任何有相同或更小的整数分母的其他整数分数近似值更接近π。^[30]π是超越数，据定义来说它不是代數數，又因此不可能是二次無理數；是故π不能表示为循环连分数。尽管${\displaystyle \pi }$ 的简单连分数没有表现出任何其他明显规律，^[31]数学家發現了数條广义连分数能表示π，例如：^[32]

${\displaystyle \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^[30]（选自 A063674 及 A063673。）
• 小數（整数后首80位）：3.14159265358979323846264338327950288419716939937510582097494459230781640628620899...^[2]:240（另见 A000796）

其他进位制的近似值

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

复数与欧拉恒等式

 

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

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

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

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

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

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

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

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

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

谱特征

 

震盪弦的泛音是二次微分的本徵函數，會形成泛音列。對應的本徵值會形成由π整數倍組成的等差数列

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

${\displaystyle \pi }$ 在許多用處中都會以特征值形式出現。例如理想的振動弦（英语：vibrating string）問題可以建模為函數${\displaystyle f}$ 在單位區間${\displaystyle [0,1]}$ 的圖形，固定邊界值为${\displaystyle f(0)=f(1)=0}$ 。弦振動的模態會是微分方程的${\displaystyle f^{n}(x)+\lambda ^{2}f(x)=0}$ ，此處λ是相關的特徵值。受施图姆－刘维尔理论限制，${\displaystyle \lambda }$ 只能是一些特定的數值。而${\displaystyle \lambda =\pi }$ 即為一個特征值，因為函數${\displaystyle f(x)=\sin(\pi x)}$ 滿足邊界條件及微分方程${\displaystyle \lambda =\pi }$ ^[37]。

 

依照第一代开尔文男爵威廉·汤姆森所述的一篇傳說，古迦太基城的外形是等周長問題的一項解（Thompson 1894）。這些包圍著海的區域由迦太基女王狄多所圍，城不靠海的邊界須用指定大小的牛皮圍住，後來是將牛皮剪成小段

${\displaystyle \pi }$ 是上述方程的最小特征值，也和弦振動的基本模式（英语：fundamental mode）有關。一種讓弦振動的方式是提供弦能量，能量會滿足維廷格函數不等式^[38]，其中提到若函數${\displaystyle f:[0,1]\rightarrow \mathbb {C} }$ 使得${\displaystyle f(0)=f(1)=0}$ ，且${\displaystyle f}$ 和${\displaystyle f'}$ 都是平方可積函數，則以下的不等式成立：

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

此例中等號成立的條件恰好是${\displaystyle f}$ 為${\displaystyle \sin(\pi x)}$ 倍數的時候。因此${\displaystyle \pi }$ 似乎是維爾丁格不等式的最佳常數，也是最小的特征值（根據雷利商數（英语：Rayleigh quotient）的計算方式）

${\displaystyle \pi }$ 在更高維度的分析也有類似的角色，出現在其他類似問題的特徵值中。就如以上所述，${\displaystyle \pi }$ 的一項特點是等周定理中的最佳常數：周長為${\displaystyle P}$ 的平面若尔当曲线，所圍面積${\displaystyle A}$ 滿足以下的不等式

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

${\displaystyle A=\pi r^{2}}$ 及${\displaystyle P=2\pi r}$ ，故等號成立的條件是曲线為圓形^[39]。

圓周率π也和庞加莱不等式的最佳常數有關^[40]，${\displaystyle \pi }$ 是一維及二維的狄氏能量（英语：Dirichlet energy）特征向量最佳值中最小，會出現在許多經典的物理現象中，例如經典的位势论^[41][42][43]。其一維的情形即為維廷格不等式。

圓周率π也是傅里叶变换的重要常數，傅里叶变换屬於积分变换，將實數線上有複數值、可積分的函數，轉換為以下形式：

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

傅里叶变换有幾種不同的寫法，但不論怎麼寫，傅里叶变换及反傅里叶变换中，一定會有某處出現${\displaystyle \pi }$ 。不過上述的定義是最經典的，因為其描述了L²空間中唯一的幺正算符，也是${\displaystyle L^{1}}$ 空間到${\displaystyle L^{\infty }}$ 空間的代數同態^[44]。

不确定性原理也用到${\displaystyle \pi }$ 。不确定性原理提出了可以將函數在空間及在頻域中局部化程度的下限，用傅立葉轉換的方式表示：

${\displaystyle \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}}$ 。

物理的結果，有關量子力学中同時觀測位置及動量的不確定性，見下文。傅立葉分析中出現π是史東－凡紐曼定理（英语：Stone–von Neumann theorem）的結果，證實了海森伯群的薛定諤表示（英语：Schrödinger representation）是唯一^[45]。

高斯积分

 

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

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

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

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

${\displaystyle \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}}}$ 的积分。

另外，当高斯函数为以下形式时，它则是平均数为${\displaystyle \mu }$ 和標準差为${\displaystyle \sigma }$ 的正态分布的機率密度函數^[46]：

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

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

 

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

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

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

这里的${\displaystyle \sigma _{x}}$ 與${\displaystyle \sigma _{p}}$ 分別為位置與動量的標準差，${\displaystyle \hbar }$ 是約化普朗克常数，而不等式的等号当且仅当粒子的波函数为高斯函数使成立^[48]。

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

历史

主条目：π的近似值

参见：π計算年表

远古时期

圓周率在远古时期（西元前一千纪）已估算至前两位（3.1）。有些埃及學家聲稱，遠至古王國時期時期的古埃及人已經用${\textstyle {\frac {22}{7}}}$ 作為圓周率的約數^{[50][註 2]}，但這說法受到質疑。^{[52][53][54][55]}

最早有記載的对圓周率估值在古埃及和巴比伦出现，兩估值都与圆周率的正确数值相差不到百分之一。巴比伦曾出土一塊西元前1900至1600年的泥板，泥板上的幾何學陳述暗示人们当时把圓周率視同${\textstyle {\frac {25}{8}}}$ （等於3.125）。^[2]:167埃及的莱因德数学纸草书（鉴定撰寫年份為西元前1650年，但抄自一份西元前1850年的文本）載有用作計算圓面積的公式，该公式中圓周率等于${\textstyle ({\frac {16}{9}})^{2}}$ （≈3.1605）。^[2]:167

西元前4世紀的《百道梵書（英语：Shatapatha Brahmana）》的天文學運算把${\textstyle {\frac {339}{108}}}$ （≈3.139，精确到99.91%）用作圓周率估值^[56]。西元前150年前其他印度文獻把圓周率視為${\textstyle {\sqrt {10}}}$ （≈3.1622）^[2]:169。

割圆时代

 

π可以透過計算圓的外切多邊形及內接多邊形周長來估算

第一條有紀錄、嚴謹計算π數值的演算法是用正多邊形的幾何算法，在西元前250年由希臘數學家阿基米德發明。^[2]:170這算法用了有一千年之久，因而有時π亦稱阿基米德常數。^{[2]:175、205}阿基米德的算法是在計算圓的外切正六邊形及內接正六邊形的邊長，以此計算${\displaystyle \pi }$ 的上限及下限，之後再將六邊形變成十二邊形，繼續計算邊長，一直計到正96邊形為止。他根據多邊形的邊長證明${\textstyle {\frac {223}{71}}<\pi <{\frac {22}{7}}}$ （也就是${\textstyle 3.1408<\pi <3.1429}$ ）^[57]。阿基米德得到的上限${\textstyle {\frac {22}{7}}}$ 也造成常見誤解，認為${\displaystyle \pi }$ 就等於${\textstyle {\frac {22}{7}}}$ ^[2]:171。在西元前150年，希臘羅馬的科學家克劳狄乌斯·托勒密在《天文学大成》一書中提到π的數值是3.1416，可能來自阿基米德，也可能來自阿波罗尼奥斯。^[2]:176[58]數學家在1630年利用多邊形的方式計算π到第39位小數，一直到1699年，其他數學家才利用無窮級數的方式打破其紀錄，計算到第71位小數^[59]。

 

阿基米德發展了用多邊形近似π的計算方式

中国历史上，${\displaystyle \pi }$ 的數值有3^[60]、3.1547（西元前一世紀）、${\displaystyle {\sqrt {10}}}$ （西元前100年，數值約3.1623）及${\textstyle {\frac {142}{45}}}$ （第三世紀，數值約3.1556）^{[2]:176–177}。大約在西元265年，曹魏數學家刘徽創立割圆术，用正3072邊形計算出π的數值為3.1416。^[61][2]:177他後來又發明了較快的算法，利用邊數差兩倍的正多邊形，其面積的差值會形成等比數列，其公比為${\textstyle {\frac {1}{4}}}$ 的原理，配合96邊形算出π的值為3.14。^[61]祖冲之在西元480年利用割圆术計算12288邊形邊長，得到π的值在3.1415926和3.1415927之间。他同时提出了π的约率${\textstyle {\frac {22}{7}}}$ 和密率${\textstyle {\frac {355}{113}}}$ 。在之後的八百年內，這都是π最準確的估計值。^[2]:178為紀念祖沖之對圓周率發展的貢獻，日本數學家三上義夫將這推算值命名為“祖沖之圓周率”，簡稱“祖率”。^[62]

印度天文學家阿耶波多在西元499年的著作《阿里亞哈塔曆書》中使用了3.1416的數值。^[2]:179斐波那契在大約1220年用獨立於阿基米德多邊形法，計算出3.1418^[2]:180。義大利作家但丁·阿利吉耶里用的數值則是${\textstyle 3+{\frac {\sqrt {2}}{10}}\approx 3.14142}$ 。^[2]:180

波斯天文學家卡西在1424年利用3×2²⁸邊的多邊形，計算到六十進制的第9位小數，相當十進制的第16位小數。^[63][64]這一突破成為當時的紀錄，延續了約180年。^[65]法國數學家弗朗索瓦·韦达在1579年用3×2¹⁷邊形計算到第9位小數^[65]，佛蘭芒數學家阿德里安·范·羅門在1593年計算到第15位小數^[65]。荷蘭數學家鲁道夫·范·科伊伦在1596年計算到第20位小數，他之後又計算到第35位小數（因此在二十世紀初之前，圓周率在德國會稱為鲁道夫數）。^{[2]:182–183}荷蘭科學家威理博·司乃耳在1621年計算到第34位小數^[2]:183，而奧地利天文學家克里斯托夫·格林伯格（英语：Christoph Grienberger）在1630年用10⁴⁰邊形計算到第38位小數^[66]，至今這仍是利用多邊形算法可以達到最準確的結果^[2]:183。

无穷级数

 

比較幾條曾用來計π的無窮級數的收斂情形。S_n是只取前n項的近似值。每張圖都是對應前一張圖的陰影部份，然後放大橫軸10倍。（點擊察看細節）

16及17世紀時，開始改用無窮级数的方式去計π。無窮级数是一組無窮數列的和^{[2]:185–191}。無窮级数讓數學家可以計算出比阿基米德以及其他用幾何方式計算的數學家更準確的結果。^{[2]:185–191}雖然詹姆斯·格雷果里及戈特弗里德·莱布尼茨等歐洲數學家利用無窮數列計算π而使得该方法为大家所知，但这种方法最早是由印度科學家在大約1400到1500年之間發現。^{[2]:185-186[67]}第一個记载用無窮级数計算π的人是约西元1500年左右时，印度天文學家尼拉卡莎·薩默亞士（英语：Nilakantha Somayaji）在他的著作《系統匯編（英语：Tantrasamgraha）》中用梵語詩所記錄。^[68]當時沒有這數列對應的證明，而證明出現在另一本較晚的印度作品《基本原理》，年代約在西元1530年。尼拉卡莎將該數列歸功於更早期的印度數學家桑加馬格拉馬的馬德哈瓦（英语：Madhava of Sangamagrama）（1350–1425）。^[68]相關的無窮级数有許多，包括有關${\displaystyle \sin }$ 、${\displaystyle \tan }$ 及${\displaystyle \cos }$ 的，現在稱為馬德哈瓦數列（英语：Madhava series）或π的莱布尼茨公式^[68]。瑪達瓦在1400年用無窮级数計算π到第11位小數，但在1430年一位波斯數學家卡西利用多邊形算法否定了他算的的結果^[69]。

 

艾萨克·牛顿利用無窮级数計算π到第15位，後來寫道：「我很羞愧的告訴你我為了這個計算用了多少個數字。」^[70]

歐洲發現的第一條無窮項圓周率公式是無窮乘積（和一般用來計算π的無窮級數不同），由法國科學家弗朗索瓦·韦达在1593年發現^[2]:187[71]：

${\displaystyle {\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年發現了沃利斯乘积，是歐洲發現的第二條無窮項圓周率公式^[2]:187：

${\displaystyle {\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年代發明，許多計π的無窮級數出現。牛頓自己就用反正弦（${\displaystyle \arcsin }$ ）數列在1655年或1666年將π近似到第15位小數，後來寫到「我很羞愧告訴你我為了計算它用了多少數字，我當時沒有做其他事。」^[70]

蘇格蘭數學家詹姆斯·格雷果里在1671年發現了馬德哈瓦公式，莱布尼茨也在1674年發現：^{[2]:188–189[72]}

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

這公式即為格雷果里－莱布尼茨公式，在${\displaystyle z=1}$ 時數值為${\displaystyle {\frac {\pi }{4}}}$ 。^[72]1699年時英國數學家亚伯拉罕·夏普用格雷果里－莱布尼茨公式，在${\displaystyle z={\frac {1}{\sqrt {3}}}}$ 時計算，計算到π的第71位小數，打破由多邊形算法得到的第39位小數的记录。^[2]:189格雷果里－莱布尼茨公式在${\displaystyle z=1}$ 時非常簡單，但收斂到最終值的速度非常慢，現在不会再用此公式來計π。^[2]:156

約翰·梅欽在1706年用格雷果里－莱布尼茨級數產生了可以快速收斂的公式：^{[2]:192–193}

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

梅欽用這公式計到π第100位小數^[2]:72–74後來其他數學家也發展了一些類似公式，現在稱為梅欽類公式，創下了許多計算π位數的紀錄。^[2]:72–74在進入電腦時代時，梅欽類公式仍然是耳熟能详可以計算π的公式，而且在约250年的时间里，很多有關π位數的紀錄都是梅欽類公式所得，比如在1946年時由達尼爾·弗格森（Daniel Ferguson）用這類公式計到第620位小數，是沒有計算設備輔助的最佳紀錄。^{[2]:192–196, 205}

1844年，計算天才扎卡里亞斯·達斯（英语：Zacharias Dase）在德國數學家卡爾·弗里德里希·高斯的要求下以梅欽類公式心算了π的200位小數，並創下紀錄。^[2]:194-196英國數學家威廉·謝克斯（英语：William Shanks）花了15年的時間計算π到小數707位，不過第528位小數出錯，後面的小數也都不正確。^{[2]:194–196}

收敛速度

有些π的無窮級數收斂的比其他級數要快，數學家一般會選用收斂速度較快的級數，可以在較少的計算量下計算π，且達到需要的準確度^{[73][2]:15–17, 70–72, 104, 156, 192–197, 201–202}。以下是π的莱布尼茨公式：^[2]:69–72

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

隨著一項一項的值加入總和中，只要項次夠多，總和最後會慢慢接近π。不過此數列的收斂速度很慢，要到50萬項之後，才會精確到π的第五位小數^[74]。

尼拉卡莎在15世紀發展了π的另一條無窮級數，收斂速度比格雷果里－萊布尼茨公式快很多：^[75]

${\displaystyle \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項	收斂到
${\displaystyle \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…
${\displaystyle \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位正確的小數位^[73]。

无理与超越性

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

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

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

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

引入π符号

 

萊昂哈德·歐拉在他1736年到1748年的作品中開始用希臘字母π表示圓周率，數學界也開始廣為使用

在用π专指“圆周率”之前，希腊字母即已用於幾何概念中^[2]:166。威廉·奥特雷德在1647年起在《數學之鑰》（Clavis Mathematicae）就已經用${\displaystyle \pi }$ 及${\displaystyle \delta }$ （對應p和d的希臘字母）來表示圓的周長及直徑的比例。

威廉·琼斯在他1706年出版的《新數學導論》（A New Introduction to the Mathematics）提到了${\displaystyle \pi }$ ，是目前已知最早专门用希臘字母${\displaystyle \pi }$ 表示圓周和其直徑比例的人^[80]。這希臘字母第一次出现是在书中討論一塊半徑1的圓時提到「其圓周長一半（${\displaystyle \pi }$ ）」。琼斯選用${\displaystyle \pi }$ 可能因它是希臘文“周边”一词“περιφέρεια”的首字母^[81]。不過琼斯提到，他那些有關${\displaystyle \pi }$ 的算式出自「真正聰明的約翰·梅欽先生」，人们推測在瓊斯之前，約翰·梅欽就已开始用${\displaystyle \pi }$ 表示圓周率^[2]:166。

瓊斯在1706年開始使用此希臘字母，但直到萊昂哈德·歐拉在其1736年出版的《力學（英语：Mechanica）》中開始使用之后，其他数学家才纷纷开始用${\displaystyle \pi }$ 指代圆周率。在此之前，數字家可能用像c或p之類的字母代表圓周率^[2]:166。歐拉與歐洲其他數學家间时常互相写信来往，${\displaystyle \pi }$ 的用法迅速傳播开来^[2]:166。1748年歐拉在他的《无穷小分析引论》再一次提到了${\displaystyle \pi }$ ，写道：「簡潔起見，我們將此數字寫為${\displaystyle \pi }$ ，${\displaystyle \pi }$ 等於半徑為1的圓周長的一半。」这表示方式之後也推展到整片西方世界^[2]:166。

现代数值近似

计算机时代与迭代算法

 

约翰·冯·诺伊曼所屬的團隊是用數位計算機ENIAC來計π的第一隊

高斯－勒让德算法：
一開始設定

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

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

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

則π的估計值為

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

二十世紀中期计算机技术发展、革新再次引发了計算π位數的熱潮。美國數學家约翰·伦奇及李維·史密斯在1949年用桌上型計算機計算到1120位^[2]:205。同年，喬治·韋斯納（George Reitwiesner）及约翰·冯·诺伊曼帶領的團隊利用反三角函数（arctan）的無窮級數，用ENIAC計算到了小數後2037位，花了70小時的電腦工作時間^[82]。這紀錄後來多次由其他透過arctan級數计算出的結果打破（1957年到7480位小數，1958年到第一萬位數，1961年到第十萬位小數），直到1973年，小数点后第一百萬位小數經已算出^[2]:197。

1980年代有两项發明加速計算了π。第一项是發现了新的迭代法去计π的值，計算速度比無窮級數快很多；另一项是發现了可以快速計算大數字乘積的乘法演算法^[2]:15–17。電腦大部分的工作時間都是在計乘法，這類演算法對現代計π格外重要^[2]:131。這類演算法包括嘉良對馬（Karatsuba）算法、譚曲（Toom－Cook）乘法及以傅里叶变换為基礎的乘法演算法（傅里叶乘法）^{[2]:132, 140}。

迭代演算法最早是在1975年至1976年间分别由美國物理學家尤金·薩拉明（英语：Eugene Salamin (mathematician)）及奧地利科學家理查·布蘭特（英语：Richard Brent (scientist)）独立提出^[2]:87。這两條演算法没有依赖無窮級數來計算。迭代會重覆特定計算，将前一次的計算結果作为這一次的輸入值，使得計算結果漸漸的趨近理想值。此方式的原始版本其實是在160年前由卡爾·弗里德里希·高斯提出，現在稱為算术－几何平均数算法（AGM法）或高斯－勒让德算法^[2]:87。薩拉明及布蘭特都曾修改之，这算法也稱為薩拉明－布蘭特演算法。

迭代演算法收斂速度比無窮級數快很多，在1980年代以後廣為使用。無窮級數隨著項次的增加，一般來說正確的位數也會增加幾位，但迭代演算法每計算多一次，正確位數會呈几何级数增长。例如薩拉明－布蘭特演算法每計算多一次，正確位數會是之前的二倍。1984年加拿大人喬納森·波温（英语：Jonathan Borwein）及彼得·波温（英语：Peter Borwein）提出迭代演算法，每計算多一次，正確位數會是之前的四倍，1987年時有另一條迭代演算法，每計算多一次，正確位數會是之前的五倍^[83]。日本數學家金田康正使用的演算法在1955年及2002年間創下了若干項紀錄^[84]。不過迭代演算法的快速收斂也有其代價，需要的記憶體明顯比無窮級數多^[84]。

计算π的意义

 

當數學家發現新的算法、電腦變得普及时，π的已知小數位急剧增加。注意垂直坐标使用了对数坐标。

一般而言，π值并不需要过于精确便能够满足大部分数学运算的需求。按照約·安（Jörg Arndt）及古里斯佗夫·希奴（Christoph Haenel）的计算，39位精確度已可将可觀測宇宙圆周的精确度準確至一粒原子大小，足以運算絕大多數宇宙学的计算需求^[85]。尽管如此，和π有關的成就往往成為世界各地的新聞頭條；部分人出于對破紀錄的冲动，依然奋力算出π小数点后上千甚至上百萬位^{[2]:17–19[86][87]}。此外也有測試超级计算机、測試数值分析算法（包括高精度乘法算法（英语：Multiplication algorithm#Fast multiplication algorithms for large inputs））等實際好處。純粹數學這领域也能计算π的位数评定其隨機度^[2]:18。

快速收敛级数

 

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

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

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

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

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

此公式每计算一项就能得到π的约14位数值^[88]，因而用於突破圆周率的数位的计算。利用这公式，楚德诺夫斯基兄弟于1989年算得π小数点后10亿（10⁹）位，法布里斯·贝拉于2009年算得2.7千亿（2.7×10¹²）位，亚历山大·易和近藤滋在2011年算得一万亿（10¹³）位。^{[2]:110–111, 206[89][90]}类似的公式还有拉马努金-佐藤级数（英语：Ramanujan–Sato series）。

2006年，加拿大数学家西蒙·普勞夫利用PSLQ整数关系算法（英语：integer relation algorithm）^[91]按照以下模版生成了几條计算π的新公式：

${\displaystyle \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)}$ ，

${\displaystyle q}$ 为e^{${\displaystyle \pi }$} ，${\displaystyle k}$ 是奇数，${\displaystyle a,b,c}$ 是普勞夫计算出的有理常数。^[92]

統計模擬法

 

布豐投針問題，多枚长度为ℓ的针随机地抛掷向平面。

 

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

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

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

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

用統計模擬法计${\displaystyle \pi }$ 的另一例子是随机往内切四分之一圆的正方形内抛掷大量点，落在四分之一圆内的点的数量与抛掷点的总量的比值会近似于${\textstyle {\frac {\pi }{4}}}$ 。^{[2]:39–40[95]}

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

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

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

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

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

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

可证明^[96]，

${\displaystyle E(W_{N})=0}$ ，${\displaystyle 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充分大时可根据以下公式算出${\displaystyle \pi }$ 的近似值：^[97]

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

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

阀门算法

1995年引入的兩條算法开辟了研究${\displaystyle \pi }$ 的新途径。因为每计算出一位数字，該數就会像流过阀门的水一样不会再出现在后续的计算过程中，这种新進算法叫阀门算法（英语：spigot algorithm）。^{[2]:77–84[99]}这就与无穷级数及迭代算法形成对比——无穷级数和迭代算法自始至终的每一步计算都会涉及到之前所有步骤计算出的中间值。^[2]:77–84

1995年，美國數學家斯坦·瓦格纳（英语：Stan Wagon）和斯坦利·拉比諾維茨（Stanley Rabinowitz）发明了一种簡單的阀门算法^{[99][2]:77[100]}，其運算速度類似arctan演算法，但速度比迭代算法慢^[2]:77。

贝利－波尔温－普劳夫公式（BBP）是另一條阀门算法，屬於一种位數萃取演算法（英语：digit extraction algorithm）。1995年，西蒙·普勞夫等人發現^{[2]:117, 126–128[101]}

${\displaystyle \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)}$ 

這公式和其他公式不同，可以計算${\displaystyle \pi }$ 的任何十六进小數位，而不用計算前面全部小數位^{[2]:117, 126–128}。十六进数位可计算得到特定二进数位；想要得到八进制数位的话，计算一、两位十六进小數即可。目前也已發現一些這種演算法的變體，不過還沒有发现針對十進制、可以快速生成特定小數位的位數萃取演算法^[102]。位數萃取演算法的一項重要用途是用來確認聲稱是計算到${\displaystyle \pi }$ 小數位數的新紀錄：若有聲稱是新紀錄的計算結果出現，先將十進制的數值轉換到十六進制，再用贝利－波尔温－普劳夫公式去確認最後一些位數（用亂數決定），若這些位數都對，就能有一定把握认为此計算結果是对的^[90]。

1998年到2000年間，分布式计算計畫PiHex（英语：PiHex）用貝拉公式（贝利－波尔温－普劳夫公式的一種變體）計算${\displaystyle \pi }$ 第10¹⁵位，結果是0^[2]:20[103]。2010年9月，有雅虎員工用公司的Apache Hadoop應用程式在上千台電腦計算π在2×10¹⁵位开始往后256位，其第2×10¹⁵位剛好也是0^[104]。

用途

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

几何学与三角学

 

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

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

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

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

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

${\displaystyle \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π度。^[106]由这条定义可得，180度＝π弧度，1度＝${\textstyle {\frac {\pi }{180^{\circ }}}}$ 弧度。^[106]因此，常用的三角函数的周期为${\displaystyle \pi }$ 的倍数；例如，正弦和余弦周期为π，^[107]任何角度${\displaystyle \theta }$ 和任何整数${\displaystyle k}$ 都有

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

拓扑学

 

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

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

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

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

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

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

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

向量分析

 

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

向量分析是与向量場的性质有关的微积分的分支，并有许多物理用途，例如用在电磁学中。位于三维笛卡尔坐标系原点的点源${\displaystyle Q}$ 的牛顿位势（英语：Newtonian potential）为^[110]

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

表示位于距原点${\displaystyle \left\vert {\boldsymbol {x}}\right\vert }$ 的单位质量（或电荷）的势能，而${\displaystyle k}$ 是维度常数。在这里由${\displaystyle \mathrm {E} }$

行進 18, 2023