古希臘歐幾裏得《幾何原本》（約公元前3世紀初）中提到圓周率是常數，中國古算書《周髀算經》（ 約公元前2世紀）中有“徑壹而周三”的記載，也認為圓周率是常數。歷史上曾采用過圓周率的多種近似值，早期大都是通過實驗而得到的結果，如古埃及紙草書（約公元前1700）中取π=（4/3）^4≈3.1604 。第壹個用科學方法尋求圓周率數值的人是阿基米德，他在《圓的度量》（公元前3世紀）中用圓內接和外切正多邊形的周長確定圓周長的上下界，從正六邊形開始，逐次加倍計算到正96邊形，得到(3+(10/71)) < π < (3+(1/7)) ，開創了圓周率計算的幾何方法（亦稱古典方法，或阿基米德方法），得出精確到小數點後兩位的π值。
中國數學家劉徽在註釋《九章算術》（263年）時只用圓內接正多邊形就求得π的近似值，也得出精確到兩位小數的π值，他的方法被後人稱為割圓術。他用割圓術壹直算到圓內接正192邊形。
南北朝時代數學家祖沖之進壹步得出精確到小數點後7位的π值（約5世紀下半葉），給出不足近似值3.1415926和過剩近似值3.1415927，還得到兩個近似分數值，密率355/113和約率22/7。其中的密率在西方直到1573才由德國人奧托得到，1625年發表於荷蘭工程師安托尼斯的著作中，歐洲稱之為安托尼斯率。
阿拉伯數學家卡西在15世紀初求得圓周率17位精確小數值，打破祖沖之保持近千年的紀錄。
德國數學家柯倫於1596年將π值算到20位小數值，後投入畢生精力，於1610年算到小數後35位數，該數值被用他的名字稱為魯道夫數。
1579年法國數學家韋達給出π的第壹個解析表達式。
此後，無窮乘積式、無窮連分數、無窮級數等各種π值表達式紛紛出現，π值計算精度也迅速增加。1706年英國數學家梅欽計算π值突破100位小數大關。1873 年另壹位英國數學家尚可斯將π值計算到小數點後707位，可惜他的結果從528位起是錯的。到1948年英國的弗格森和美國的倫奇***同發表了π的808位小數值，成為人工計算圓周率值的最高紀錄。
電子計算機的出現使π值計算有了突飛猛進的發展。1949年美國馬裏蘭州阿伯丁的軍隊彈道研究實驗室首次用計算機（ENIAC）計算π值，壹下子就算到2037位小數，突破了千位數。1989年美國哥倫比亞大學研究人員用克雷－2型和IBM－VF型巨型電子計算機計算出π值小數點後4.8億位數，後又繼續算到小數點後10.1億位數，創下新的紀錄。
除π的數值計算外，它的性質探討也吸引了眾多數學家。1761年瑞士數學家蘭伯特第壹個證明π是無理數。1794年法國數學家勒讓德又證明了π2也是無理數。到1882年德國數學家林德曼首次證明了π是超越數，由此否定了困惑人們兩千多年的「化圓為方」尺規作圖問題。還有人對π的特征及與其它數字的聯系進行研究。如1929年蘇聯數學家格爾豐德證明了eπ 是超越數等等。
圓周率的計算
古今中外，許多人致力於圓周率的研究與計算。為了計算出圓周率的越來越好的近似值，壹代代的數學家為這個神秘的數貢獻了無數的時間與心血。
十九世紀前，圓周率的計算進展相當緩慢，十九世紀後，計算圓周率的世界紀錄頻頻創新。整個十九世紀，可以說是圓周率的手工計算量最大的世紀。
進入二十世紀，隨著計算機的發明，圓周率的計算有了突飛猛進。借助於超級計算機，人們已經得到了圓周率的2061億位精度。
歷史上最馬拉松式的計算，其壹是德國的Ludolph Van Ceulen，他幾乎耗盡了壹生的時間，計算到圓的內接正262邊形，於1609年得到了圓周率的35位精度值，以至於圓周率在德國被稱為Ludolph數；其二是英國的William Shanks，他耗費了15年的光陰，在1874年算出了圓周率的小數點後707位。可惜，後人發現，他從第528位開始就算錯了。
把圓周率的數值算得這麽精確，實際意義並不大。現代科技領域使用的圓周率值，有十幾位已經足夠了。如果用Ludolph Van Ceulen算出的35位精度的圓周率值，來計算壹個能把太陽系包起來的壹個圓的周長，誤差還不到質子直徑的百萬分之壹。以前的人計算圓周率，是要探究圓周率是否循環小數。自從1761年Lambert證明了圓周率是無理數，1882年Lindemann證明了圓周率是超越數後，圓周率的神秘面紗就被揭開了。
現在的人計算圓周率, 多數是為了驗證計算機的計算能力，還有，就是為了興趣。
圓周率的計算方法
古人計算圓周率，壹般是用割圓法。即用圓的內接或外切正多邊形來逼近圓的周長。Archimedes用正96邊形得到圓周率小數點後3位的精度；劉徽用正3072邊形得到5位精度；Ludolph Van Ceulen用正262邊形得到了35位精度。這種基於幾何的算法計算量大，速度慢，吃力不討好。隨著數學的發展，數學家們在進行數學研究時有意無意地發現了許多計算圓周率的公式。下面挑選壹些經典的常用公式加以介紹。除了這些經典公式外，還有很多其它公式和由這些經典公式衍生出來的公式，就不壹壹列舉了。
1、 Machin公式
[這個公式由英國天文學教授John Machin於1706年發現。他利用這個公式計算到了100位的圓周率。Machin公式每計算壹項可以得到1.4位的十進制精度。因為它的計算過程中被乘數和被除數都不大於長整數，所以可以很容易地在計算機上編程實現。
Machin.c 源程序
還有很多類似於Machin公式的反正切公式。在所有這些公式中，Machin公式似乎是最快的了。雖然如此，如果要計算更多的位數，比如幾千萬位，Machin公式就力不從心了。下面介紹的算法，在PC機上計算大約壹天時間，就可以得到圓周率的過億位的精度。這些算法用程序實現起來比較復雜。因為計算過程中涉及兩個大數的乘除運算，要用FFT(Fast Fourier Transform)算法。FFT可以將兩個大數的乘除運算時間由O(n2)縮短為O(nlog(n))。
2、 Ramanujan公式
1914年，印度數學家Srinivasa Ramanujan在他的論文裏發表了壹系列***14條圓周率的計算公式，這是其中之壹。這個公式每計算壹項可以得到8位的十進制精度。1985年Gosper用這個公式計算到了圓周率的17,500,000位。
1989年，David & Gregory Chudnovsky兄弟將Ramanujan公式改良成為：
這個公式被稱為Chudnovsky公式，每計算壹項可以得到15位的十進制精度。1994年Chudnovsky兄弟利用這個公式計算到了4,044,000,000位。Chudnovsky公式的另壹個更方便於計算機編程的形式是：
3、AGM(Arithmetic-Geometric Mean)算法
Gauss-Legendre公式：
這個公式每叠代壹次將得到雙倍的十進制精度，比如要計算100萬位，叠代20次就夠了。1999年9月Takahashi和Kanada用這個算法計算到了圓周率的206,158,430,000位，創出新的世界紀錄。
4、Borwein四次叠代式：
這個公式由Jonathan Borwein和Peter Borwein於1985年發表，它四次收斂於圓周率。
這個公式簡稱BBP公式，由David Bailey, Peter Borwein和Simon Plouffe於1995年***同發表。它打破了傳統的圓周率的算法，可以計算圓周率的任意第n位，而不用計算前面的n-1位。這為圓周率的分布式計算提供了可行性。1997年，Fabrice Bellard找到了壹個比BBP快40％的公式：