版塊導航: 正在加載中...

客戶端APP下載

論文輔導

調(diào)劑小程序

登錄注冊

應《網(wǎng)絡安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

北京石油化工學院2026年研究生招生接收調(diào)劑公告

返回列表

zyj8119

木蟲 (著名寫手)

模擬EPI: 10
應助: 65 (初中生)
貴賓: 0.003
金幣: 915.1
散金: 1440
紅花: 35
帖子: 2936
在線: 1329.4小時
蟲號: 664177
注冊: 2008-11-29
性別: GG
專業(yè): 理論和計算化學

[交流] 【轉(zhuǎn)帖】用fortran編的投針問題的小程序

蒲豐氏投針問題的模擬過程，隨機數(shù)發(fā)生器也是自編的，以供大家參考和提出建議。謝謝。（seed1和seed2最好選擇3和5，為了使投針次數(shù)達到1000000，CVF進行如下設置Project->settings->link->
output，將stack allocations reserve:設為1000000000）

CODE:

       program getpi  

      implicit none 

      real,parameter::a=5,L=4,pi=3.14159   

      integer::n1,i,counter=0  

      real,allocatable::R1(:),R2(:)

      real::theta,x,pi1 

      write(*,*) 'input the size of the array:'  

      read(*,*) n1 

      allocate(R1(n1)) 

      allocate(R2(n1))   

      call random(n1,R1,R2) 

      do i=1,n1

       x=a*(2*R1(i)-1) 

       theta=pi*R2(i)  

      if(abs(x)
      end do 

      pi1=(1.0*n1)/(counter*1.0)*(2.0*L/a) 

      write(*,*) 'this is PI:'   

      write(*,*) pi 

      write(*,"('this is ratio of counter to total number',F10.6)") 1.0

     &*counter/n1

      stop   

      end program  

      subroutine random(n,R1,R2) 

      implicit none  

      integer n,seed1,seed2,i,little,m 

      real::R1(n),R2(n)  

      integer::temp1(n),temp2(n)

      write(*,*) 'input seed1'  

      read(*,*) seed1  

      write(*,*) 'input seed2'

      read(*,*) seed2  

      m=2**30   

      m=m*2-1   

      temp1(1)=397204094  

      temp2(1)=temp1(1)  

      R1(1)=(1.0*temp1(1))/(1.0*m)   

      R2(1)=(1.0*temp2(1))/(1.0*m) 

      little=0 

      if(R1(1)<0.5) little=little+1  

      do i=1,n-1 

      temp1(i+1)=mod(seed1*temp1(i),m)  

      R1(i+1)=(1.0*temp1(i+1))/(1.0*m) 

      temp2(i+1)=mod(seed2*temp2(i),m)  

      R2(i+1)=(1.0*temp2(i+1))/(1.0*m) 

      if(R1(i+1)<0.5) little=little+1 .  

      end do ; 

      write(*,*) 'ratio of number which is little than 0.5'  

      write(*,*) 1.0*little/n 

      return 

      end subroutine

[ Last edited by zyj8119 on 2010-11-12 at 16:09 ]

回復此樓

» 猜你喜歡

283求調(diào)劑已經(jīng)有5人回復
一志愿哈爾濱工業(yè)大學材料與化工方向336分已經(jīng)有7人回復
數(shù)一英一271專碩（085401）求調(diào)劑，可跨已經(jīng)有4人回復
0856，材料與化工321分求調(diào)劑已經(jīng)有11人回復
0703化學求調(diào)劑已經(jīng)有8人回復
277跪求調(diào)劑已經(jīng)有9人回復
材料277求調(diào)劑已經(jīng)有6人回復
085602 307分求調(diào)劑已經(jīng)有7人回復
一志愿南昌大學324求調(diào)劑已經(jīng)有7人回復
材料與化工272求調(diào)劑已經(jīng)有6人回復

好好學習，天天向上。

1樓 2010-10-17 01:08:42

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

智能機器人

Robot (super robot)

我們都愛小木蟲

找到一些相關的精華帖子，希望有用哦~

在linux下編FORTRAN程序已經(jīng)有14人回復
我編的Simpson積分法fortran程序給不出結(jié)果，大俠們看看哪里出了問題？已經(jīng)有4人回復
求助Fortran編程已經(jīng)有8人回復
求助計算機高人，編寫個Fortran程序已經(jīng)有13人回復
請教一個fortran小程序編譯出錯的問題，謝謝已經(jīng)有9人回復
如何用fortran編寫高斯白噪聲程序已經(jīng)有6人回復
有關fortran的一次而問題，希望大家能幫幫忙，謝謝已經(jīng)有4人回復
寫了一個fortran90的小程序，編譯通不過，請大俠幫忙已經(jīng)有59人回復
【求助】請問一下fortran的可視化編程已經(jīng)有7人回復

點擊這里搜索更多相關資源

科研從小木蟲開始，人人為我，我為人人

相關版塊跳轉(zhuǎn) 我要訂閱樓主 zyj8119 的主題更新

返回列表

普通表情龍兔虎貓高級回復 (可上傳附件)

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 0703化學求調(diào)劑 +7	奶油草莓. 2026-03-22	8/400	2026-03-28 09:59 by 神馬都不懂
[材料工程] 一志愿C9材料與化工專業(yè)總分300求調(diào)劑 +8	曼111 2026-03-24	9/450	2026-03-28 07:58 by YYYYX1234
[考研] 學碩274求調(diào)劑 +5	Li李魚 2026-03-26	5/250	2026-03-27 20:51 by 熱情沙漠
[考研] 化學調(diào)劑 +4	愛吃番茄的旭 2026-03-24	5/250	2026-03-27 17:50 by kiokin
[考研] 一志愿西北大學總分282 英語一62 求調(diào)劑 +7	18419759900 2026-03-25	8/400	2026-03-27 16:38 by 18419759900
[考研] 材料與化工085600，總分304，本科有兩篇sci參與，求調(diào)劑 +10	幸運的醬醬 2026-03-22	12/600	2026-03-27 16:08 by muchong357
[考研] 0856調(diào)劑 +5	求求讓我有書讀?/a> 2026-03-26	6/300	2026-03-27 15:12 by caszguilin
[考研] 化學308分求調(diào)劑 +8	你好明天你好 2026-03-23	9/450	2026-03-27 14:01 by 楊光于青云
[考研] 071000生物學求調(diào)劑，初試成績343 +6	小小甜面團 2026-03-25	6/300	2026-03-26 23:01 by 不吃魚的貓
[考研] 321求調(diào)劑 +6	Ymlll 2026-03-24	6/300	2026-03-26 20:50 by 不吃魚的貓
[考研] 297求調(diào)劑 +6	田洪有 2026-03-26	6/300	2026-03-26 15:55 by 不吃魚的貓
[考研] 291 求調(diào)劑 +7	化工2026屆畢業(yè)?/a> 2026-03-21	8/400	2026-03-26 11:25 by AlenQIN.
[考研] 285求調(diào)劑 +3	AZMK 2026-03-24	3/150	2026-03-25 12:23 by userper
[考研] 340求調(diào)劑 +5	話梅糖111 2026-03-24	5/250	2026-03-25 06:53 by ilovexiaobin
[考研] 284求調(diào)劑 +3	yanzhixue111 2026-03-23	6/300	2026-03-23 22:58 by pswait
[考研] 361求調(diào)劑 +3	Glack 2026-03-22	3/150	2026-03-23 22:03 by fuyu_
[考研] 一志愿東華大學化學070300，求調(diào)劑 +7	2117205181 2026-03-21	8/400	2026-03-22 22:55 by chixmc
[考研] 求調(diào)劑一志愿海大，0703化學學碩304分，有大創(chuàng)項目，四級已過 +6	幸運哩哩 2026-03-22	10/500	2026-03-22 20:10 by edmund7
[考研] 306求調(diào)劑 +5	來好運來來來 2026-03-22	5/250	2026-03-22 16:17 by BruceLiu320
[考研] 275求調(diào)劑 +6	shansx 2026-03-22	8/400	2026-03-22 15:27 by barlinike

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時熱門版塊排行榜

zyj8119

[交流] 【轉(zhuǎn)帖】用fortran編的投針問題的小程序

» 猜你喜歡

智能機器人

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频