| 7 | 1/1 | 返回列表 |
| 查看: 1108 | 回復(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ā)表 | |
|---|---|---|---|---|
|
[考研] 材料學(xué)碩333求調(diào)劑 +3 | 北道巷 2026-03-24 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 274求調(diào)劑 +5 | 顧九笙要謙虛 2026-03-24 | 5/250 |
|
|
[考研] 一志愿南航材料專317分求調(diào)劑 +5 | 炸呀炸呀炸薯?xiàng)l 2026-03-23 | 5/250 |
|
|
[材料工程] 一志愿C9材料與化工專業(yè)總分300求調(diào)劑 +4 | 曼111 2026-03-24 | 5/250 |
|
|
[考研] 一志愿華東理工大學(xué)081700,初試分?jǐn)?shù)271 +5 | kotoko_ik 2026-03-23 | 6/300 |
|
|
[考研] 341求調(diào)劑(一志愿湖南大學(xué)070300) +5 | 番茄頭--- 2026-03-22 | 6/300 |
|
|
[考研] 材料專碩英一數(shù)二306 +8 | z1z2z3879 2026-03-18 | 8/400 |
|
|
[考研] 336求調(diào)劑 +4 | 收到VS 2026-03-20 | 4/200 |
|
|
[考研] 303求調(diào)劑 +4 | 元夕元 2026-03-20 | 4/200 |
|
|
[考研] 生物學(xué)一志愿985,分?jǐn)?shù)349求調(diào)劑 +6 | zxts12 2026-03-21 | 9/450 |
|
|
[考研] 求老師收我 +3 | zzh16938784 2026-03-23 | 3/150 |
|
|
[考研] 323求調(diào)劑 +6 | 洼小桶 2026-03-18 | 6/300 |
|
|
[考研] 317求調(diào)劑 +12 | 申子申申 2026-03-19 | 18/900 |
|
|
[考研] 289材料與化工(085600)B區(qū)求調(diào)劑 +3 | 這么名字咋樣 2026-03-22 | 4/200 |
|
|
[考研] 考研調(diào)劑 +3 | 呼呼?~+123456 2026-03-21 | 3/150 |
|
|
[考研] 266求調(diào)劑 +3 | 哇呼哼呼哼 2026-03-20 | 3/150 |
|
|
[基金申請]
學(xué)校已經(jīng)提交到NSFC,還能修改嗎?
40+4
|
babangida 2026-03-19 | 9/450 |
|
|
[考研] 求調(diào)劑 +3 | 白QF 2026-03-21 | 3/150 |
|
|
[考研] 330求調(diào)劑0854 +3 | assdll 2026-03-21 | 3/150 |
|
|
[考研]
|
簡木ChuFront 2026-03-19 | 8/400 |
|