Peter de Jong Strange Attractor

dejong.clj

(defn compute-dejong
  "Computes a single DeJong 2d point vector for given params and XY pos"
  [a b c d x y]
  (v/vec2
   (+ (Math/sin (* a y)) (Math/cos (* (* b x) x)))
   (+ (Math/sin (* (* c x) x)) (Math/cos (* d y)))))

(->> (range 1e6)
     (reduce
      (fn [[points [x y]] _]
        (let [pos (compute-dejong a b c d x y)]
          [(conj points pos) pos]))
      [[] [(m/random width) (m/random width)]])
     (first))