7楼: Originally posted by 19870617 at 2012-08-27 21:54:47
您好！源码是这样的
        program main
     implicit none
     integer :: n,i
     real(8) ::pi,u1,u2,maxr
         real(8) :: r(1)  
    pi=3.1415926535
     call random_seed ()
     n=1000      
    d ...

program main
   implicit none
   integer :: n,i
   real(8) ::pi,u1,u2,maxr
   real(8) :: r(1000)  
   pi=3.1415926535
   call random_seed ()
   n=1000      
   do i=1,n
      call random_number(u1)
      call random_number(u2)
      r(i)=sqrt(-log(u1))*cos(2*pi*u2)  
      print *, r(i)
   enddo
   maxr=maxval(r)
   print *,"the largest is",maxr
end

2楼: Originally posted by snoopyzhao at 2012-08-27 17:56:59
把 r 定义为一个数组就可以了……

找本入门的书看一下吧

3楼: Originally posted by 19870617 at 2012-08-27 18:30:01
r定义为数组？一维大小为1000嘛？试了一下输出的r有3列，感觉不对啊...