c   乱数の生成 
      
      iseed=1
      do 10 j=1,100
        nran=irandf(iseed)
        a=randf(iseed)
        write(6,*) j,nran,a
 10   continue
      stop
      end 

c   0から2147483647の一様乱数を生成 
      function irandf(iseed)
      iseed=iseed*48828125
      if(iseed.le.0) iseed=(iseed+2147483647)+1
      irandf=iseed
      return
      end

c   0から1の一様乱数を生成 
      function randf(iseed)
      iseed=iseed*48828125
      if(iseed.le.0) iseed=(iseed+2147483647)+1
      randf=float(iseed)/2147483647
      return
      end