[编程交流] 巴恩斯利'；s Fern&Sierpin公司

(defun c:fern (/ ptlst pt probability)
 ;; Lee Mac  ~  20.02.10
 (repeat 50000    
   (setq probability (rng))
   (Point (setq pt  (cond (  (< probability 0.01)
                             (iterate pt '((0.0 0.0)
                                           (0.0 0.16)) '(0.0 0.0)))
                          (  (<= 0.01 probability 0.86)
                             (iterate pt '(( 0.85 0.04)
                                           (-0.04 0.85)) '(0.0 1.6)))
                          (  (< 0.86 probability 0.93)
                             (iterate pt '((0.20 -0.26)
                                           (0.23  0.22)) '(0.0 1.6)))
                          (t (iterate pt '((-0.15 0.28)
                                           ( 0.26 0.24)) '(0.0 0.44)))))))
 (princ))
(defun c:sierpinski (/ ptlst pt probability)
 ;; Lee Mac  ~  20.02.10
 (repeat 50000    
   (setq probability (rng))
   (Point (setq pt (cond (  (< probability 0.333)
                            (iterate pt '((0.5 0.0)
                                          (0.0 0.5)) '(0.0 0.0)))
                         (  (<= 0.333 probability 0.666)
                            (iterate pt '((0.5 0.0)
                                          (0.0 0.5)) '(0.5 0.0)))
                         
                         (t (iterate pt '((0.5 0.0)
                                          (0.0 0.5)) '(0.25 0.5)))))))
 (princ))
(defun iterate (point matrix vector)
 (mapcar
   (function +)
     (mapcar
       (function
         (lambda (row)
           (apply (function +)
                  (mapcar (function *) row point)))) matrix) vector))
(defun rng (/ modulus multiplier increment random) ;; Stig
 (if (not seed) (setq seed (getvar "DATE")))
 (setq modulus    4294967296.0 multiplier 1664525 increment 1 
       seed       (rem (+ (* multiplier seed) increment) modulus) 
       random     (/ seed modulus)))
(defun Point (pt)
 (entmakex (list (cons 0 "POINT") (cons 10 pt) (cons 62 40))))

f(f(...f(x)))  ~ n times