program test_random

 use random_mod
use nr
 implicit none
 integer :: seed = 0. 
 integer i,nmuons/10000/
integer j, bin(10)
 real(rp) harvest,x,dum,y,z
                                                                   
!call ran_seed_put(seed)

bin = 0
do i=1,nmuons

  call ran_uniform(harvest)
j = 10*harvest + 1
bin(j) = bin(j)+1
!  call ran1(dum)  

!  x = log(harvest)     
!  y=  log(dum)
! call expdev(z)

write(11,'(i8, es12.4)') i, harvest !,    dum, x,y,z
end do


do j = 1, 10
   print *, j, bin(j)
enddo
end