| 7 | 1/1 | 返回列表 |
| 查看: 1122 | 回復(fù): 6 | |||
836449366金蟲 (小有名氣)
|
[求助]
計(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就等于圓周率了(即面積之比),我所編寫的程序如下: 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,yend 求蟲友幫幫忙,謝謝 |
|
感覺你代碼的錯(cuò)誤蠻多,然后我改寫了一下。 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/numdeallocate(x,y) stop end program |
金蟲 (小有名氣)
金蟲 (小有名氣)
金蟲 (小有名氣)
| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 336求調(diào)劑 +5 | kiyy 2026-04-01 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 274求調(diào)劑 +5 | 薛定諤的虎。 2026-04-01 | 5/250 |
|
|
[考研] 理學(xué)07化學(xué) 303求調(diào)劑 +12 | 睿08 2026-03-27 | 12/600 |
|
|
[考研] 安全工程 285 求調(diào)劑 +3 | Xinyu56 2026-04-01 | 4/200 |
|
|
[考研] 266分,一志愿電氣工程,本科材料,求材料專業(yè)調(diào)劑 +8 | 哇呼哼呼哼 2026-04-01 | 8/400 |
|
|
[考研] 0710生物學(xué)考研調(diào)劑 +3 | 李多米lee. 2026-03-27 | 4/200 |
|
|
[考研] 調(diào)劑 +5 | 好好讀書。 2026-03-28 | 7/350 |
|
|
[考研] 環(huán)境工程 085701,267求調(diào)劑 +15 | minht 2026-03-29 | 16/800 |
|
|
[教師之家] 張雪峰戛然而止的飛馳人生 +3 | yexuqing 2026-03-26 | 4/200 |
|
|
[考研] 318求調(diào)劑 +3 | 篤行致遠(yuǎn). 2026-03-31 | 3/150 |
|
|
[考研] 材料求調(diào)劑 一志愿哈工大總分298分,前三科223分 +11 | dongfang59 2026-03-27 | 11/550 |
|
|
[考研] 315求調(diào)劑 +6 | akie... 2026-03-28 | 7/350 |
|
|
[考研] 南京大學(xué)化學(xué)調(diào)劑 +11 | 景隨風(fēng) 2026-03-29 | 16/800 |
|
|
[考研] 085600材料與化工調(diào)劑 +16 | kikiki7 2026-03-30 | 16/800 |
|
|
[考研] 083000學(xué)碩274求調(diào)劑 +12 | Li李魚 2026-03-26 | 12/600 |
|
|
[考研] 085601一志愿西北工業(yè)大學(xué)初試346 +4 | 085601初試346 2026-03-30 | 4/200 |
|
|
[有機(jī)交流] 考研調(diào)劑 +8 | watb 2026-03-26 | 8/400 |
|
|
[考研] 297求調(diào)劑 +17 | 田洪有 2026-03-26 | 18/900 |
|
|
[考研] 086000生物與醫(yī)藥調(diào)劑 +5 | Feisty。 2026-03-28 | 9/450 |
|
|
[考研] 316求調(diào)劑 +7 | 江辭666 2026-03-26 | 7/350 |
|