home *** CD-ROM | disk | FTP | other *** search
/ NetNews Usenet Archive 1992 #20 / NN_1992_20.iso / spool / comp / sys / hp48 / 4501 < prev    next >
Encoding:
Internet Message Format  |  1992-09-08  |  1.6 KB
  1. Path: sparky!uunet!usc!sol.ctr.columbia.edu!ira.uka.de!math.fu-berlin.de!news.netmbx.de!Germany.EU.net!news.Hamburg.Germany.EU.net!mwhh!dmhh!detlef
  2. From: detlef@dmhh.hanse.de (Detlef Mueller)
  3. Newsgroups: comp.sys.hp48
  4. Subject: Re: HP's RAND algorithm
  5. Distribution: world
  6. Message-ID: <715705766snx@dmhh.hanse.de>
  7. References: <1992Sep4.183731.130@news.uiowa.edu>
  8. Date: Sat, 05 Sep 92 15:09:26 GMT
  9. Organization: Who is Organized ?
  10. Lines: 34
  11.  
  12. jollie@zambini.cs.uiowa.edu writes in article <1992Sep4.183731.130@news.uiowa.edu>:
  13. > 
  14. > I just got my 48SX, and I love it! But, I was wondering what algorithm
  15. > HP used for the random number generator. I'm hoping that it is a
  16. > prime modulus multiplicative linear conguential generator (I'm not
  17. > making this up) with a=16807 and m=2^31 - 1. I've written one of my own,
  18. > but would prefer to use the built in RAND.
  19. > 
  20. > Jeff Ollie
  21. > jollie@zambini.cs.uiowa.edu "The Happy User"
  22. > 
  23.  
  24. Don't worry, RAND is basing on a linear congruential method (see pg. 9 of
  25. 'The Art of ...', Vol.2 by D.E. Knuth):
  26.  
  27.     Xn+1 = (a * Xn + c) mod m
  28.  
  29. where
  30.     a = 2851130928467
  31.     c = 0
  32.     m = 1E15
  33.  
  34. The intial value of X0 after a memory lost is 999500333083533.
  35.  
  36. If you choose an X0 which isn't a multiple of 2 or 5 then you'll get a
  37. period of 5E13 (see pg. 20 of above book - a mod 200 is 67).
  38.  
  39. Bye,
  40.     8-Detlef
  41. --
  42. +------------------------------------+--------------------------------------+
  43. |`What a depressingly stupid machine'|             Detlef Mueller           |
  44. | -- Marvin                          |          detlef@dmhh.hanse.de        |
  45. +------------------------------------+--------------------------------------+
  46.