Нахождение производной многочлена

;; Производная монома
 
(defun m-derive (c p x sign)
  (cond ((zerop c) nil)
        ((= p 0) nil)  
        ((= p 0) (list (if (plusp sign) '+ '-) c))
        ((> p 1) (list (if (plusp sign) '+ '-) (* c p) x '^ (- p 1)))
        (t (list (if (plusp sign) '+ '-) c))))
 
;; Производная полинома
        
(defun poly-derive (poly x &optional (c 0) (p 0) (sign 1) (flg 0) (res nil))
  (cond ((null poly) (append res (m-derive c p x sign))) 
        ((numberp (car poly)) 
         (cond ((zerop flg) (poly-derive (cdr poly) x (+ (* c 10) (car poly)) p sign flg res))
               (t           (poly-derive (cdr poly) x c (+ (* p 10) (car poly)) sign flg res))))
        ((or (eq (car poly) x) (eq (car poly) '^)) (poly-derive (cdr poly) x c p sign (+ 1 flg) res))
        ((eq (car poly) '+) (poly-derive (cdr poly) x 0 0 1 0 (append res (m-derive c p x sign))))
        ((eq (car poly) '-) (poly-derive (cdr poly) x 0 0 -1 0 (append res (m-derive c p x sign))))))
 
;; Для удобства запуска
 
(defun derive (poly arg)
  (implode (poly-derive (explode poly) arg)))
 
(derive '2x^3-7x^2+13x^1-8 'x)
 
==> +6X^2-14X^1+13
 
(derive '4y^4-56y^1+2y^2 'y)
 
==> +16Y^3-56+4Y^1