home *** CD-ROM | disk | FTP | other *** search
/ Frozen Fish 1: Amiga / FrozenFish-Apr94.iso / bbs / alib / d3xx / d385 / xlispstat.lha / XLispStat / lisp.lzh / XLisp-Stat / Book / lik-gamma.lsp < prev    next >
Encoding:
Text File  |  1990-10-11  |  285 b   |  14 lines
  1. ; book p.108
  2.  
  3. (require "functions/newton.lsp")
  4.  
  5. (def x (gamma-rand 30 4.5))
  6.  
  7. (def f (let ((s (sum (log x)))
  8.              (n (length x)))
  9.          #'(lambda (a)
  10.              (- (* (- a 1) s) (* n (log-gamma a))))))
  11.  
  12. (def df (make-derivative f .001))
  13. (def ddf (make-derivative df .001))
  14.