integer function uniform(bottom,top,seed)
      implicit integer(a-z)
      real random
      uniform=int((float(top-bottom)+0.9999999)*random(seed))+bottom
      return
      end

      real function random(seed)
      implicit integer(a-z)
      data a/16807/,m/2147483647/,q/127773/,r/2836/
      seed=mod(seed,q)*a-seed/q*r
      if(seed.le.0)seed=seed+m
      random=float(seed)/m
      return
      end