home *** CD-ROM | disk | FTP | other *** search
/ Dos/V Magazine 2002 July 1 / VMAG130101.iso / ONLINE / monthly / calc / CLPCA511.LZH / ClipCalc / EXTFUNC / loggamma.cef < prev    next >
Encoding:
Text File  |  2002-08-18  |  724 b   |  33 lines
  1. #!/usr/local/bin/clip
  2. #âKâôâ}è╓Éöé╠æ╬Éö logarithm of gamma function
  3. #\!- <val>
  4. :float \16
  5.  
  6. @N = 8
  7.  
  8. # Bernoulli Éö
  9. @@B  0 = 1.0
  10. @@B  1 = \-1.0    / 2.0
  11. @@B  2 = 1.0      / 6.0
  12. @@B  4 = \-1.0    / 30.0
  13. @@B  6 = 1.0      / 42.0
  14. @@B  8 = \-1.0    / 30.0
  15. @@B 10 = 5.0      / 66.0
  16. @@B 12 = \-691.0  / 2730.0
  17. @@B 14 = 7.0      / 6.0
  18. @@B 16 = \-3617.0 / 510.0
  19.  
  20. @v = 1
  21. while @0 < @N; @v *= @0; @0++; endwhile
  22. @w = 1 / (@0 * @0)
  23.  
  24. @@B 16 / (16 * 15)
  25.     * @w + @@B 14 / (14 * 13)
  26.     * @w + @@B 12 / (12 * 11)
  27.     * @w + @@B 10 / (10 *  9)
  28.     * @w + @@B  8 / ( 8 *  7)
  29.     * @w + @@B  6 / ( 6 *  5)
  30.     * @w + @@B  4 / ( 4 *  3)
  31.     * @w + @@B  2 / ( 2 *  1)
  32.     / @0 + 0.5 * log (2 * !pi) - log @v - @0 + (@0 - 0.5) * log @0
  33.