rand-seq.clj

(defn rand-seqn [prog size seed]
  (loop [ i     0
          seqn  [seed]]		
    (if (>= i size)
      seqn
      (recur  (inc i)
              (conj seqn (first 
                (execute prog [(last seqn) 1.0 1.0 1.0])))))))


(defn- rev-map [ma]
  (reduce (fn [rev pair] (conj rev pair)) 
    [] (map reverse (seq ma))))


(defn seqn-distr [size]
  (loop [ data {}
          seqn (rand-seqn (gen-rand 80) 10 1.0) 
          cntr 0 ] 
    (if (> cntr size)
      (sort (fn [a b] (> (first a) (first b))) (rev-map data))
      (if (contains? data seqn)
        (recur (assoc data seqn (inc (data seqn)))
               (rand-seqn (gen-rand 80) 10 1.0) 
               (inc cntr))
        (recur (assoc data seqn 1)
               (rand-seqn (gen-rand 80) 10 1.0) 
               (inc cntr))))))