应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 计算圆周率的小程序出错
71/1返回列表
查看: 1070  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

836449366

金虫 (小有名气)

[求助] 计算圆周率的小程序出错 已有1人参与

我的程序是计算圆周率，原理是：在一个1*1的正方形及其内切圆（圆心为（0.5,0.5））中，随机数x和随机数y组成的点位于正方形内，但要判断这个点是否也在圆内，当点数很多时，圆内的点数除以总数再乘以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
求虫友帮帮忙，谢谢
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-12-12 19:57:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fish.yfyh

铜虫 (小有名气)

【答案】应助回帖

感谢参与，应助指数 +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
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

836449366
2楼2013-12-12 20:30:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

836449366

金虫 (小有名气)
送红花一朵
引用回帖:
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
...

非常感谢你的回答，运行也正确，但当我在write后加一个一个变量输出时，改为
” write(*,"(a,F6.3,E10.2)" "pai=",4.0*ii/num,num “，
这样输出的结果为
" pai= 3.141  0.14E-39 "
这样的输出的结果是错误的，num输出的值应为1.00E+05
问题出在哪？求虫友再次解惑下？
我的目的是做个循环，看看输入的值越大，得到的圆周率的结果越准确。
赞一下 回复此楼
3楼2013-12-13 11:05:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fish.yfyh

铜虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★
836449366: 金币+5, ★★★★★最佳答案, 非常感谢 2013-12-13 15:05:24
引用回帖:
3楼: Originally posted by 836449366 at 2013-12-13 11:05:12
非常感谢你的回答，运行也正确，但当我在write后加一个一个变量输出时，改为
” write(*,"(a,F6.3,E10.2)" "pai=",4.0*ii/num,num “，
这样输出的结果为
" pai= 3.141  0.14E-39 &q ...

因为num是整数呀，E10.2是科学计数法的浮点数，所以这样就可以了。
write(*,"(a,F6.3,E10.2)" ) "pai=",4.0*ii/num,num*1.0
赞一下 回复此楼
4楼2013-12-13 11:21:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

836449366

金虫 (小有名气)
引用回帖:
4楼: Originally posted by fish.yfyh at 2013-12-13 11:21:49
因为num是整数呀，E10.2是科学计数法的浮点数，所以这样就可以了。
write(*,"(a,F6.3,E10.2)" ) "pai=",4.0*ii/num,num*1.0...

非常感谢你的回答
赞一下 回复此楼
5楼2013-12-13 15:06:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

836449366

金虫 (小有名气)
引用回帖:
4楼: Originally posted by fish.yfyh at 2013-12-13 11:21:49
因为num是整数呀，E10.2是科学计数法的浮点数，所以这样就可以了。
write(*,"(a,F6.3,E10.2)" ) "pai=",4.0*ii/num,num*1.0...

我还想问问，我做了个循环：
…………
do 100 n=9,12,1
       num=10**(n)
…………
100   continue
       end
为什么能够做到这样的结果：
pai=   3.141556501388550 0.10E+10（正确）
pai=   3.141569137573242 0.14E+10（错误），就是不能做到10^(10)这样级别的数来，要想做到更大的数，不如10^（25）这样的数，怎样修改程序？？？
赞一下 回复此楼
6楼2013-12-13 21:01:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fish.yfyh

铜虫 (小有名气)
引用回帖:
6楼: Originally posted by 836449366 at 2013-12-13 21:01:02
我还想问问，我做了个循环：
…………
do 100 n=9,12,1
       num=10**(n)
…………
100   continue
       end
为什么能够做到这样的结果：
pai=   3.141556501388550 0.10E+10（正确）
pai=   3.14156 ...

因为fortran里面整型默认是短整型的,即integer也就是integer(kind=4)
试试定义num为长整型数 integer(kind=8) :: num
赞一下 回复此楼
7楼2013-12-13 22:10:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 836449366 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区文章，我:8 O5 51O 54,科目齐全 +3 5lbyq5wrhb 2026-02-07 4/200 2026-02-08 08:47 by vs90ilomwc
[考博] 售SCI一区文章，我:8 O5 51O 54,科目齐全 +3 5lbyq5wrhb 2026-02-07 4/200 2026-02-08 08:46 by vs90ilomwc
[论文投稿] 售SCI一区文章，我:8 O5 51O 54,科目齐全 +3 3rkserf6qr 2026-02-07 5/250 2026-02-08 08:32 by vs90ilomwc
[硕博家园] 售SCI一区文章，我:8 O5 51O 54,科目齐全 +3 3rkserf6qr 2026-02-07 4/200 2026-02-08 08:27 by vs90ilomwc
[硕博家园] 售SCI一区文章，我:8 O5 51O 54,科目齐全 +5 2h7du0nuhk 2026-02-07 6/300 2026-02-08 08:26 by vs90ilomwc
[考博] 售SCI一区文章，我:8 O5 51O 54,科目齐全 +5 2h7du0nuhk 2026-02-07 6/300 2026-02-08 08:12 by vs90ilomwc
[考博] 售SCI一区文章，我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 7/350 2026-02-08 08:06 by vs90ilomwc
[教师之家] 售SCI一区文章，我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 7/350 2026-02-08 07:52 by vs90ilomwc
[找工作] 售SCI一区文章，我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 7/350 2026-02-08 07:46 by vs90ilomwc
[公派出国] 售SCI一区文章，我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 8/400 2026-02-08 07:32 by vs90ilomwc
[教师之家] 售SCI一区文章，我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 8/400 2026-02-08 07:26 by vs90ilomwc
[硕博家园] 售SCI一区文章，我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 8/400 2026-02-08 07:07 by vs90ilomwc
[硕博家园] 博士延得我，科研能力直往上蹿 +8 偏振片 2026-02-02 8/400 2026-02-08 06:52 by liyeqik
[教师之家] 有院领导为了换新车，用横向课题经费买了俩车 +7 瞬息宇宙 2026-02-04 7/350 2026-02-07 21:47 by tfang
[有机交流] 酰胺脱乙酰基 10+5 chibby 2026-02-03 12/600 2026-02-07 19:29 by 江东闲人
[基金申请] 同年申请2项不同项目，第1个项目里不写第2个项目的信息，可以吗 +4 hitsdu 2026-02-06 4/200 2026-02-07 13:07 by jurkat.1640
[基金申请] 有时候真觉得大城市人没有县城人甚至个体户幸福 +9 苏东坡二世 2026-02-04 10/500 2026-02-07 12:37 by 小毛球
[考博] 天津大学招2026.09的博士生，欢迎大家推荐交流（博导是本人） +4 a793625982 2026-02-05 5/250 2026-02-07 10:57 by a793625982
[公派出国] CSC & MSCA 博洛尼亚大学能源材料课题组博士/博士后招生｜MSCA经费充足、排名优 +4 雨念 2026-02-01 6/300 2026-02-06 23:32 by MelissaPon
[教师之家] 遇见不省心的家人很难过 +18 otani 2026-02-03 22/1100 2026-02-04 11:06 by tangmnt
信息提示
关闭
请填处理意见
关闭