c
c Time-stamp: <2000-10-04 17:24:21 baron>
c
      program try
c     -----------
      implicit none
      real*8 :: x,randu
      integer :: iseed=32038,i,n
      print *,"give n"
      read(*,*)n
      do i=1,n
       x=randu(iseed)
       write(7,*)i,x
      enddo
      end program try
c
c
c
      real*8 function randu(iseed)
c     ----------------------------
      implicit none
      integer :: iseed
      integer,parameter :: bigint=huge(1)
      real*8,parameter :: xbigint=bigint,ix=69069.d0
      real*8,parameter :: rbigint=1.d0/xbigint
      save
c
      iseed = mod(ix*iseed,xbigint)
      randu = rbigint*iseed
c
      end function randu