subroutine test_rn() integer rndist(100) do i=1,1000 do 10 j=1,100000 call RANDOM_NUMBER(XX) if(j==1)then print *,"xx=",xx ! on the FIRST one , pauses if(i==1)read(*,*) endif ! random integer 1 to 100 ixx=xx*100.0+1.d0 rndist(ixx)=rndist(ixx)+1 10 enddo enddo ! i loop print 101,rndist 101 format("rndist="/(10I10)) read(*,*) stop end