3楼: Originally posted by totoymy at 2011-10-30 12:53:52:
帮我说说问题也好。。。应该一共出100组数。。。结果总是不对啊~

MODULE PARA
IMPLICIT NONE

INTEGER*4, PARAMETER:: Ndim=12   

REAL,DIMENSION(Ndim,5):: pp

INTEGER*4::NOUT

END MODULE PARA

PROGRAM CM
!USE MSIMSL
USE PARA
IMPLICIT NONE

INTEGER:: i,j,k
INTEGER::t
real::x,y,z
REAL*8::r
REAL::s,m

OPEN (11, file='input.txt')

READ(11,*) ((pp(i,j),j=1,4),i=1,Ndim)
close(11)
   call random_seed()

t=1
DO
   call random_number(r)
!  write (*,*) r
   
   k=1
   s = 0.0
   do while(k<=Ndim)
      s=s+pp(k,4)
      m=s-pp(k,4)

!     write(*,*) s, m
      
      if (m           x=x+pp(k,1)
          y=y+pp(k,2)
          z=z+pp(k,3)
          open(12,file="result.txt",access='append')
          write(12,*) x,y,z,t,k,r
          close(12)  
          t = t + 1
      end if
      k=k+1
   enddo
   
   if (t>100) exit
enddo

end program cm

1楼: Originally posted by totoymy at 2011-10-28 10:57:29:
MODULE PARA
IMPLICIT NONE

INTEGER*4, PARAMETER:: Ndim=12   

REAL,DIMENSION(Ndim,5):: pp

INTEGER*4::NOUT

END MODULE PARA

PROGRAM CM
USE MSIMSL
USE PARA
IMPLICIT NONE

INTEGER: ...