home *** CD-ROM | disk | FTP | other *** search
/ NetNews Usenet Archive 1993 #3 / NN_1993_3.iso / spool / sci / math / symbolic / 3517 < prev    next >
Encoding:
Internet Message Format  |  1993-01-25  |  1.5 KB
  1. Path: sparky!uunet!pmafire!news.dell.com!natinst.com!cs.utexas.edu!zaphod.mps.ohio-state.edu!howland.reston.ans.net!spool.mu.edu!agate!doc.ic.ac.uk!uknet!edcastle!hwcs!pjbk
  2. From: pjbk@cs.hw.ac.uk (Peter JB King)
  3. Newsgroups: sci.math.symbolic
  4. Subject: Mathematica bug?
  5. Message-ID: <1993Jan25.124536.27674@cs.hw.ac.uk>
  6. Date: 25 Jan 93 12:45:36 GMT
  7. Sender: news@cs.hw.ac.uk (News Administrator)
  8. Organization: Dept of Computing & Electrical Engineering, Heriot-Watt University, Scotland
  9. Lines: 33
  10.  
  11. I have found what appears to be some form of bug in Mathematica (
  12. Version 2.1 for SPARC).  The following program includes as part of its
  13. output
  14.  
  15.  Simplify`Gamma`Private`SimpGamma1 
  16.  
  17. which appears to be some sort of internal definition.  I can find no
  18. mention of SimpGamma1 in the documentation.  Intuitively, it looks like
  19. some sort of rule to simplify Gamma fundtions, but it does not appear to
  20. have any arguments.  Is the answer I get correct if I replace this by
  21. '1'?
  22.  
  23. Here is the program that demonstrates the bug.
  24.  
  25. rstar[s_] := tau / (tau + s )
  26. rbar = 1/tau
  27.  
  28. b[t_] := mu E^(-mu t )
  29. bstar[s_] := mu / (mu + s )
  30. bbar = 1/mu
  31.  
  32.  
  33. pstar[s_] := ( s + gamma ) Integrate[ E^(-(s+gamma)t) b[t] /
  34.         (s + gamma(1 - rstar[s](1-E^(-(s+gamma)t)))), { t, 0,
  35.                                 Infinity}] 
  36. Print[pstar[s]]
  37.  
  38.  
  39. -- 
  40. Peter King, Computing & Electrical Eng.    JANET:    pjbk@uk.ac.hw.cee
  41.   Heriot-Watt University        ARPA:    pjbk@cee.hw.ac.uk
  42.   Riccarton, Edinburgh EH14 4AS     UUCP:    ..!uknet!cee.hw.ac.uk!pjbk
  43. Phone: (+44) 31 451 3433                Fax: (+44) 31 451 3431
  44.