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

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

24小時(shí)熱門(mén)版塊排行榜

返回列表

當(dāng)前只顯示滿(mǎn)足指定條件的回帖，點(diǎn)擊這里查看本話(huà)題的所有回帖

836449366

金蟲(chóng) (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 3039.8
散金: 62
紅花: 4
帖子: 270
在線: 356.3小時(shí)
蟲(chóng)號(hào): 1478273
注冊(cè): 2011-11-05
專(zhuān)業(yè): 理論和計(jì)算化學(xué)

[求助] 計(jì)算圓周率的小程序出錯(cuò) 已有1人參與

我的程序是計(jì)算圓周率，原理是：在一個(gè)1*1的正方形及其內(nèi)切圓（圓心為（0.5,0.5））中，隨機(jī)數(shù)x和隨機(jī)數(shù)y組成的點(diǎn)位于正方形內(nèi)，但要判斷這個(gè)點(diǎn)是否也在圓內(nèi)，當(dāng)點(diǎn)數(shù)很多時(shí)，圓內(nèi)的點(diǎn)數(shù)除以總數(shù)再乘以4就等于圓周率了（即面積之比），我所編寫(xiě)的程序如下：
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
求蟲(chóng)友幫幫忙，謝謝

回復(fù)此樓

» 猜你喜歡

274求調(diào)劑已經(jīng)有7人回復(fù)
北科281學(xué)碩材料求調(diào)劑已經(jīng)有11人回復(fù)
085600 材料與化工 329分求調(diào)劑已經(jīng)有3人回復(fù)
302求調(diào)劑已經(jīng)有3人回復(fù)
312求調(diào)劑已經(jīng)有5人回復(fù)
網(wǎng)絡(luò)空間安全0839招調(diào)劑已經(jīng)有4人回復(fù)
329求調(diào)劑已經(jīng)有3人回復(fù)
【2026考研調(diào)劑】制藥工程 284分求相關(guān)專(zhuān)業(yè)調(diào)劑名額已經(jīng)有8人回復(fù)
085602 289分求調(diào)劑已經(jīng)有7人回復(fù)
材料學(xué)碩333求調(diào)劑已經(jīng)有4人回復(fù)

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

請(qǐng)教一個(gè)fortran小程序編譯出錯(cuò)的問(wèn)題，謝謝已經(jīng)有9人回復(fù)

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

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

fish.yfyh

銅蟲(chóng) (小有名氣)

應(yīng)助: 42 (小學(xué)生)
金幣: 607.9
紅花: 4
帖子: 108
在線: 35.9小時(shí)
蟲(chóng)號(hào): 2867094
注冊(cè): 2013-12-12
專(zhuān)業(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后加一個(gè)一個(gè)變量輸出時(shí)，改為
” 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é)計(jì)數(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的回帖

查看全部 7 個(gè)回答

fish.yfyh

銅蟲(chóng) (小有名氣)

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

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

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

感覺(jué)你代碼的錯(cuò)誤蠻多，然后我改寫(xiě)了一下。
   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

金蟲(chóng) (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 3039.8
散金: 62
紅花: 4
帖子: 270
在線: 356.3小時(shí)
蟲(chóng)號(hào): 1478273
注冊(cè): 2011-11-05
專(zhuān)業(yè): 理論和計(jì)算化學(xué)

送紅花一朵

引用回帖:

2樓: Originally posted by fish.yfyh at 2013-12-12 20:30:56
感覺(jué)你代碼的錯(cuò)誤蠻多，然后我改寫(xiě)了一下。
   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后加一個(gè)一個(gè)變量輸出時(shí)，改為
” write(*,"(a,F6.3,E10.2)"

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

贊一下

回復(fù)此樓

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

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

836449366

金蟲(chóng) (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 3039.8
散金: 62
紅花: 4
帖子: 270
在線: 356.3小時(shí)
蟲(chóng)號(hào): 1478273
注冊(cè): 2011-11-05
專(zhuān)業(yè): 理論和計(jì)算化學(xué)

引用回帖:

4樓: Originally posted by fish.yfyh at 2013-12-13 11:21:49
因?yàn)閚um是整數(shù)呀，E10.2是科學(xué)計(jì)數(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的回帖

查看全部 7 個(gè)回答

普通表情龍兔虎貓

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

[考研] 347求調(diào)劑 +4	L when 2026-03-25	4/200	2026-03-25 13:37 by cocolv
[考研] 289求調(diào)劑 +9	懷瑾握瑜l 2026-03-20	9/450	2026-03-25 11:02 by userper
[考研] 306求0703調(diào)劑一志愿華中師范 +10	紙魚(yú)ly 2026-03-21	11/550	2026-03-24 17:22 by qingfeng258
[考研] 材料292調(diào)劑 +8	橘頌思美人 2026-03-23	8/400	2026-03-24 16:33 by laoshidan
[考研] 一志愿211 初試270分求調(diào)劑 +5	谷雨上岸 2026-03-23	6/300	2026-03-24 16:32 by laoshidan
[考研] 【雙一流院校新能源、環(huán)境材料，材料加工與模擬招收大量調(diào)劑】 +4	Higraduate 2026-03-22	7/350	2026-03-24 11:23 by 種大樹(shù)
[考研] 279分求調(diào)劑一志愿211 +18	chaojifeixia 2026-03-19	20/1000	2026-03-24 10:34 by dolphin_ycj
[考研] 284求調(diào)劑 +3	yanzhixue111 2026-03-23	6/300	2026-03-23 22:58 by pswait
[考研] 336化工調(diào)劑 +4	王大坦1 2026-03-23	5/250	2026-03-23 18:32 by allen-yin
[考研] 工科0856求調(diào)劑 +5	沐析汀汀 2026-03-21	5/250	2026-03-23 17:56 by 海瑟薇-
[考研] 一志愿西安交通大學(xué)材料工程專(zhuān)業(yè) 282分求調(diào)劑 +11	楓橋ZL 2026-03-18	13/650	2026-03-22 20:26 by edmund7
[考研] 324求調(diào)劑 +6	lucky呀呀呀鴨 2026-03-20	6/300	2026-03-22 16:01 by ColorlessPI
[考研] 298求調(diào)劑一志愿211 +3	上岸6666@ 2026-03-20	3/150	2026-03-22 15:50 by ColorlessPI
[考研] 285求調(diào)劑 +6	ytter 2026-03-22	6/300	2026-03-22 12:09 by 星空星月
[考研] 材料學(xué)碩301分求調(diào)劑 +7	Liyouyumairs 2026-03-21	7/350	2026-03-21 22:31 by peike
[考研] 279求調(diào)劑 +5	紅衣隱官 2026-03-21	5/250	2026-03-21 14:59 by lature00
[考研] 考研調(diào)劑求學(xué)校推薦 +3	伯樂(lè)29 2026-03-18	5/250	2026-03-20 22:59 by JourneyLucky
[考研] 290求調(diào)劑 +7	^O^乜 2026-03-19	7/350	2026-03-20 21:43 by JourneyLucky
[考研] 求調(diào)劑一志愿南京航空航天大學(xué)289分 +3	@taotao 2026-03-19	3/150	2026-03-20 21:34 by JourneyLucky
[考研] 本科鄭州大學(xué)物理學(xué)院，一志愿華科070200學(xué)碩，346求調(diào)劑 +4	我不是一根蔥 2026-03-18	4/200	2026-03-19 09:11 by 浮云166