版塊導(dǎo)航: 正在加載中...

登錄注冊

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

24小時熱門版塊排行榜

返回列表

836449366

金蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 3039.8
散金: 62
紅花: 4
帖子: 270
在線: 356.3小時
蟲號: 1478273
注冊: 2011-11-05
專業(yè): 理論和計算化學(xué)

[求助] 計算圓周率的小程序出錯已有1人參與

我的程序是計算圓周率，原理是：在一個1*1的正方形及其內(nèi)切圓（圓心為（0.5,0.5））中，隨機(jī)數(shù)x和隨機(jī)數(shù)y組成的點(diǎn)位于正方形內(nèi)，但要判斷這個點(diǎn)是否也在圓內(nèi)，當(dāng)點(diǎn)數(shù)很多時，圓內(nèi)的點(diǎn)數(shù)除以總數(shù)再乘以4就等于圓周率了（即面積之比），我所編寫的程序如下：
program main
   implicit none
   integer*4 :: i,j,n,ii
   common i
   real*4 :: r
   call random_number(x)
   call random_number(y)
   real*4,dimension(i):: x,y
   r=sqrt((x-0.5)**2+(y-0.5)**2)
   ii=0
   do n=2,8,1
   i=10**n
   do j=1,i
   if ( r .le. 0.5 ) then
      ii=ii+1
   endif
   enddo
   enddo
   write(*,"(a,3F6.3)"

"pai=",4*ii/i,x,y
end
求蟲友幫幫忙，謝謝

回復(fù)此樓

» 猜你喜歡

一志愿廈門大學(xué)化學(xué)學(xué)碩307求調(diào)劑已經(jīng)有4人回復(fù)
材料學(xué)碩297已過四六級求調(diào)劑推薦已經(jīng)有13人回復(fù)
一志愿西北大學(xué) 總分282 英語一62 求調(diào)劑已經(jīng)有6人回復(fù)
一志愿西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。已經(jīng)有5人回復(fù)
0856求調(diào)劑已經(jīng)有9人回復(fù)
一志愿南京郵電大學(xué) 288分材料考研求調(diào)劑已經(jīng)有3人回復(fù)
材料與化工304求B區(qū)調(diào)劑已經(jīng)有5人回復(fù)
一志愿北京化工大學(xué)材料與化工（085600）296求調(diào)劑已經(jīng)有5人回復(fù)
275求調(diào)劑已經(jīng)有7人回復(fù)
274求調(diào)劑已經(jīng)有19人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

請教一個fortran小程序編譯出錯的問題，謝謝已經(jīng)有9人回復(fù)

1樓 2013-12-12 19:57:21

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

fish.yfyh

銅蟲 (小有名氣)

應(yīng)助: 42 (小學(xué)生)
金幣: 607.9
紅花: 4
帖子: 108
在線: 35.9小時
蟲號: 2867094
注冊: 2013-12-12
專業(yè): 化學(xué)生物學(xué)與生物有機(jī)化學(xué)

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

感覺你代碼的錯誤蠻多，然后我改寫了一下。
   program main
   implicit none
   integer*4 :: num,i,j,ii
   real*4 :: r
   real*4,allocatable :: x(

,y(

   num = 100000
   allocate(x(num))
   allocate(y(num))
   call random_seed()
   call random_number(x)
   call random_number(y)

   ii = 0
   do i = 1,num
      r = sqrt((x(i)-0.5)**2 + (y(i)-0.5)**2)
      if(r .le. 0.5) then
            ii = ii+1
      end if
   end do

   write(*,"(a,F6.3)"

"pai=",4.0*ii/num

   deallocate(x,y)

   stop
   end program

贊一下

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

836449366

2樓2013-12-12 20:30:56

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

836449366

金蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 3039.8
散金: 62
紅花: 4
帖子: 270
在線: 356.3小時
蟲號: 1478273
注冊: 2011-11-05
專業(yè): 理論和計算化學(xué)

送紅花一朵

引用回帖:

2樓: Originally posted by fish.yfyh at 2013-12-12 20:30:56
感覺你代碼的錯誤蠻多，然后我改寫了一下。
   program main
   implicit none
   integer*4 :: num,i,j,ii
   real*4 :: r
   real*4,allocatable :: x(

,y(

num = 100000
...

非常感謝你的回答，運(yùn)行也正確，但當(dāng)我在write后加一個一個變量輸出時，改為
” write(*,"(a,F6.3,E10.2)"

"pai=",4.0*ii/num,num “，
這樣輸出的結(jié)果為
" pai= 3.141 0.14E-39 "
這樣的輸出的結(jié)果是錯誤的，num輸出的值應(yīng)為1.00E+05
問題出在哪？求蟲友再次解惑下？
我的目的是做個循環(huán)，看看輸入的值越大，得到的圓周率的結(jié)果越準(zhǔn)確。

贊一下

回復(fù)此樓

3樓2013-12-13 11:05:12

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

fish.yfyh

銅蟲 (小有名氣)

應(yīng)助: 42 (小學(xué)生)
金幣: 607.9
紅花: 4
帖子: 108
在線: 35.9小時
蟲號: 2867094
注冊: 2013-12-12
專業(yè): 化學(xué)生物學(xué)與生物有機(jī)化學(xué)

【答案】應(yīng)助回帖

★ ★ ★ ★ ★
836449366: 金幣+5, ★★★★★最佳答案, 非常感謝 2013-12-13 15:05:24

引用回帖:

3樓: Originally posted by 836449366 at 2013-12-13 11:05:12
非常感謝你的回答，運(yùn)行也正確，但當(dāng)我在write后加一個一個變量輸出時，改為
” write(*,"(a,F6.3,E10.2)"

"pai=",4.0*ii/num,num “，
這樣輸出的結(jié)果為
" pai= 3.141 0.14E-39 &q ...

因?yàn)閚um是整數(shù)呀，E10.2是科學(xué)計數(shù)法的浮點(diǎn)數(shù)，所以這樣就可以了。
write(*,"(a,F6.3,E10.2)" ) "pai=",4.0*ii/num,num*1.0

贊一下

回復(fù)此樓

4樓2013-12-13 11:21:49

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

836449366

金蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 3039.8
散金: 62
紅花: 4
帖子: 270
在線: 356.3小時
蟲號: 1478273
注冊: 2011-11-05
專業(yè): 理論和計算化學(xué)

引用回帖:

4樓: Originally posted by fish.yfyh at 2013-12-13 11:21:49
因?yàn)閚um是整數(shù)呀，E10.2是科學(xué)計數(shù)法的浮點(diǎn)數(shù)，所以這樣就可以了。
write(*,"(a,F6.3,E10.2)" ) "pai=",4.0*ii/num,num*1.0...

非常感謝你的回答

贊一下

回復(fù)此樓

5樓2013-12-13 15:06:17

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

836449366

金蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 3039.8
散金: 62
紅花: 4
帖子: 270
在線: 356.3小時
蟲號: 1478273
注冊: 2011-11-05
專業(yè): 理論和計算化學(xué)

引用回帖:

我還想問問，我做了個循環(huán)：
…………
do 100 n=9,12,1
num=10**(n)
…………
100 continue
end
為什么能夠做到這樣的結(jié)果：
pai= 3.141556501388550 0.10E+10（正確）
pai= 3.141569137573242 0.14E+10（錯誤），就是不能做到10^(10)這樣級別的數(shù)來，要想做到更大的數(shù)，不如10^（25）這樣的數(shù)，怎樣修改程序？？？

贊一下

回復(fù)此樓

6樓2013-12-13 21:01:02

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

fish.yfyh

銅蟲 (小有名氣)

應(yīng)助: 42 (小學(xué)生)
金幣: 607.9
紅花: 4
帖子: 108
在線: 35.9小時
蟲號: 2867094
注冊: 2013-12-12
專業(yè): 化學(xué)生物學(xué)與生物有機(jī)化學(xué)

引用回帖:

6樓: Originally posted by 836449366 at 2013-12-13 21:01:02
我還想問問，我做了個循環(huán)：
…………
do 100 n=9,12,1
num=10**(n)
…………
100 continue
end
為什么能夠做到這樣的結(jié)果：
pai= 3.141556501388550 0.10E+10（正確）
pai= 3.14156 ...

因?yàn)閒ortran里面整型默認(rèn)是短整型的,即integer也就是integer(kind=4)
試試定義num為長整型數(shù) integer(kind=8) :: num

贊一下

回復(fù)此樓

7樓2013-12-13 22:10:55

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

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

返回列表

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

[考研] 材料與化工304求B區(qū)調(diào)劑 +3	邱gl 2026-03-26	5/250	2026-03-26 13:33 by 邱gl
[考研] 284求調(diào)劑 +4	junqihahaha 2026-03-26	4/200	2026-03-26 12:36 by 張凱十八號
[考研] 打過很多競賽，085406控制工程300分，求調(diào)劑 +3	askeladz 2026-03-26	3/150	2026-03-26 09:08 by 給你你注意休息
[考研] 321求調(diào)劑 +3	璞玉~~ 2026-03-25	3/150	2026-03-25 19:07 by Zhanglab-TJU
[考研] 材料與化工 322求調(diào)劑 +6	然11 2026-03-19	6/300	2026-03-25 18:37 by haxia
[考研] 網(wǎng)絡(luò)空間安全0839招調(diào)劑 +4	w320357296 2026-03-25	6/300	2026-03-25 17:59 by 255671
[考研] 347求調(diào)劑 +4	L when 2026-03-25	4/200	2026-03-25 13:37 by cocolv
[考研] 303求調(diào)劑 +6	元夕元 2026-03-20	7/350	2026-03-25 12:00 by edmund7
[考研] B區(qū)考研調(diào)劑 +4	yqdszhdap－ 2026-03-22	5/250	2026-03-25 08:51 by baoball
[有機(jī)交流] 有機(jī)合成求助 20+3	FENGSHUJEI 2026-03-23	5/250	2026-03-24 19:31 by 88817753
[考研] 335求調(diào)劑 +4	yuyu宇 2026-03-23	5/250	2026-03-23 23:49 by Txy@872106
[考研] 一志愿國科過程所081700，274求調(diào)劑 +3	三水研0水立方 2026-03-23	3/150	2026-03-23 23:11 by MajorWen
[考研] 材料/農(nóng)業(yè)專業(yè)，07/08開頭均可，過線就行 +3	呵唔哦豁 2026-03-23	4/200	2026-03-23 22:30 by 汪��？！
[考研] 化學(xué)308分求調(diào)劑 +3	你好明天你好 2026-03-23	3/150	2026-03-23 20:11 by macy2011
[考研] 350求調(diào)劑 +6	weudhdk 2026-03-19	6/300	2026-03-23 15:47 by tangyuan0840221
[考研] 求老師收我 +3	zzh16938784 2026-03-23	3/150	2026-03-23 12:56 by ztnimte
[考研] 311求調(diào)劑 +3	26研0 2026-03-20	3/150	2026-03-22 14:46 by ColorlessPI
[考研] 303求調(diào)劑 +5	安憶靈 2026-03-22	6/300	2026-03-22 12:46 by 素顏傾城1988
[考研] 332求調(diào)劑 +3	鳳凰院丁真 2026-03-20	3/150	2026-03-21 10:27 by luoyongfeng
[考研] 材料學(xué)碩318求調(diào)劑 +5	February_Feb 2026-03-19	5/250	2026-03-19 23:51 by 23Postgrad

24小時熱門版塊排行榜

[求助] 計算圓周率的小程序出錯 已有1人參與

» 猜你喜歡

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

【答案】應(yīng)助回帖

» 本帖已獲得的紅花（最新10朵）

【答案】應(yīng)助回帖

[求助] 計算圓周率的小程序出錯已有1人參與

» 本主題相關(guān)價值貼推薦，對您同樣有幫助: