[updated]C語言和python實現(xiàn)蒙特卡里算法求pi值

高中數(shù)學(xué)課本上曾經(jīng)寫過用計算機算得pi的近似值胳徽,用的是蒙特卡里方法。但是那時的我還不知道編程為何物江耀,今天突然又想到這個問題,嘗試了一下诉植,發(fā)現(xiàn)并不是那么難祥国,用C語言和python各實現(xiàn)了一遍。
不多說上代碼晾腔。


c語言實現(xiàn)

#include <stdlib.h>
#include <stdio.h>
#include <math.h>
#include <time.h>
#define MAXN 1200000//總隨機點數(shù)舌稀,這個數(shù)值越大,pi值計算的越精確灼擂,相應(yīng)的耗費更多的時間

 int main(void)
 {
    int i=0;
    double count=0.0;
    double x=0.0,y=0.0,dist=0.0;
    double my_pi=0;
    srand((unsigned)time(NULL));
    //隨機數(shù)種子被初始化
    clock_t start = clock();//時間起點
    for(i=0;i<MAXN;i++)
     {
          x = ((double)rand()/RAND_MAX);//生成0-1的隨機數(shù)
          y = ((double)rand()/RAND_MAX);
          dist = sqrt(x*x + y*y);//計算(x壁查,y)到原點的距離
          if(dist <= 1.0)
             count++;//記錄落在半圓內(nèi)的隨機點的個數(shù)
     }
    my_pi = 4*(count/MAXN);
    clock_t finish = clock();//函數(shù)運行結(jié)束時間記錄
    printf("%lf\n",my_pi);
    double duration = (double)(finish - start) / CLOCKS_PER_SEC;//計算出函數(shù)所耗費的時間
    printf("It costs %lf s",duration);
    return 0;
 }

python版 python版本為3.5X

# -*- coding:utf-8 -*- 
from random import random
from math import sqrt
from time import clock

MAXN = 120000
count =0
clock()
for i in range (1,MAXN):
    x, y = random(),random()
    dist = sqrt(x**2+y**2)
    if dist <= 1.0:
        count = count + 1
pi = 4*(count/MAXN)
print("PI的值為") % pi
print("程序的運行時間是%-5.5ss") % clock()



pi.PNG

看的出來計算出的pi值離精確的pi值還有一定距離,并且數(shù)值會波動剔应,畢竟這是用概率去模擬的睡腿,樣本點足夠多,才可能精確不是领斥?
第一次使用Markdown寫東西嫉到,將就看吧o(╯□╰)o

┴┬┴┬/ ̄\_/ ̄\
┬┴┬┴▏  ▏▔▔▔▔\
┴┬┴/\ /      ﹨
┬┴∕      ≡侣濉/  『味瘛)
┴┬▏        ●  ▏
┬┴▏           ▔█ 
┴◢██◣     \___/
┬█████◣      〗狼/  
┴█████████████◣
◢██████████████▆▄
█◤◢██◣◥█████████◤\
◥◢████ ████████◤   \
┴█████ ██████◤      ﹨
┬│   │█████◤        ▏
┴│   │              ▏
┬ ∕    ∕   ∠覆恪/▔▔▔\     ∕
┴/___/﹨   ∕     ﹨  /\
┬┴┬┴┬┴\    \      ﹨/   ﹨
┴┬┴┬┴┬┴ \___\     ﹨/▔\﹨ ▔\
▲△▲▲╓╥╥╥╥╥╥╥╥\   ∕  /▔﹨/▔﹨
』=А**╠╬╬╬╬╬╬╬╬*﹨ ∫呤辍/  //

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
  • 序言:七十年代末碎节,一起剝皮案震驚了整個濱河市捧搞,隨后出現(xiàn)的幾起案子,更是在濱河造成了極大的恐慌狮荔,老刑警劉巖胎撇,帶你破解...
    沈念sama閱讀 222,946評論 6 518
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件,死亡現(xiàn)場離奇詭異殖氏,居然都是意外死亡晚树,警方通過查閱死者的電腦和手機,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 95,336評論 3 399
  • 文/潘曉璐 我一進店門雅采,熙熙樓的掌柜王于貴愁眉苦臉地迎上來爵憎,“玉大人慨亲,你說我怎么就攤上這事”模” “怎么了刑棵?”我有些...
    開封第一講書人閱讀 169,716評論 0 364
  • 文/不壞的土叔 我叫張陵,是天一觀的道長愚铡。 經(jīng)常有香客問我铐望,道長,這世上最難降的妖魔是什么茂附? 我笑而不...
    開封第一講書人閱讀 60,222評論 1 300
  • 正文 為了忘掉前任,我火速辦了婚禮督弓,結(jié)果婚禮上营曼,老公的妹妹穿的比我還像新娘。我一直安慰自己愚隧,他們只是感情好蒂阱,可當(dāng)我...
    茶點故事閱讀 69,223評論 6 398
  • 文/花漫 我一把揭開白布。 她就那樣靜靜地躺著狂塘,像睡著了一般录煤。 火紅的嫁衣襯著肌膚如雪。 梳的紋絲不亂的頭發(fā)上荞胡,一...
    開封第一講書人閱讀 52,807評論 1 314
  • 那天妈踊,我揣著相機與錄音,去河邊找鬼泪漂。 笑死廊营,一個胖子當(dāng)著我的面吹牛,可吹牛的內(nèi)容都是我干的萝勤。 我是一名探鬼主播露筒,決...
    沈念sama閱讀 41,235評論 3 424
  • 文/蒼蘭香墨 我猛地睜開眼,長吁一口氣:“原來是場噩夢啊……” “哼敌卓!你這毒婦竟也來了慎式?” 一聲冷哼從身側(cè)響起,我...
    開封第一講書人閱讀 40,189評論 0 277
  • 序言:老撾萬榮一對情侶失蹤趟径,失蹤者是張志新(化名)和其女友劉穎瘪吏,沒想到半個月后,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體舵抹,經(jīng)...
    沈念sama閱讀 46,712評論 1 320
  • 正文 獨居荒郊野嶺守林人離奇死亡肪虎,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點故事閱讀 38,775評論 3 343
  • 正文 我和宋清朗相戀三年,在試婚紗的時候發(fā)現(xiàn)自己被綠了惧蛹。 大學(xué)時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片扇救。...
    茶點故事閱讀 40,926評論 1 353
  • 序言:一個原本活蹦亂跳的男人離奇死亡刑枝,死狀恐怖,靈堂內(nèi)的尸體忽然破棺而出迅腔,到底是詐尸還是另有隱情装畅,我是刑警寧澤,帶...
    沈念sama閱讀 36,580評論 5 351
  • 正文 年R本政府宣布沧烈,位于F島的核電站掠兄,受9級特大地震影響,放射性物質(zhì)發(fā)生泄漏锌雀。R本人自食惡果不足惜蚂夕,卻給世界環(huán)境...
    茶點故事閱讀 42,259評論 3 336
  • 文/蒙蒙 一、第九天 我趴在偏房一處隱蔽的房頂上張望腋逆。 院中可真熱鬧婿牍,春花似錦、人聲如沸惩歉。這莊子的主人今日做“春日...
    開封第一講書人閱讀 32,750評論 0 25
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽撑蚌。三九已至上遥,卻和暖如春,著一層夾襖步出監(jiān)牢的瞬間争涌,已是汗流浹背粉楚。 一陣腳步聲響...
    開封第一講書人閱讀 33,867評論 1 274
  • 我被黑心中介騙來泰國打工, 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留第煮,地道東北人解幼。 一個月前我還...
    沈念sama閱讀 49,368評論 3 379
  • 正文 我出身青樓,卻偏偏與公主長得像包警,于是被迫代替她去往敵國和親撵摆。 傳聞我的和親對象是個殘疾皇子,可洞房花燭夜當(dāng)晚...
    茶點故事閱讀 45,930評論 2 361

推薦閱讀更多精彩內(nèi)容