home *** CD-ROM | disk | FTP | other *** search
/ Oakland CPM Archive / oakcpm.iso / sigm / vol266 / tstrnd.for < prev    next >
Encoding:
Text File  |  1986-05-19  |  935 b   |  45 lines
  1. C    [TSTRND.FOR of JUGPDS Vol.8]
  2. C
  3.     program    tstrnd
  4. C
  5. C    This is a test program for DRAND and RANDU.
  6. C
  7.     integer*4    isum
  8.     dimension    irand(500),rand(500)
  9.     data     bel/Z'07'/
  10. C
  11.     write(1,100)
  12.   100      format(1h1,5x,'** Test of Random Number Routine **'/)
  13.    10    write(1,110)
  14.   110       format(1h ,'Enter number of sets to be generated (I4): ')
  15.     read(1,200) iset
  16.   200     format(I4)
  17.     if (iset.le.500) goto 20
  18.     write(1,120) bel
  19.   120      format(1h ,A1,'Too big (ISET<=500)!'/)
  20.                     goto 10
  21.    20    write(1,130)
  22.   130      format(1h ,'Enter the seed (I4): ')
  23.     read(1,200) iseed
  24. C
  25.     isum = 0
  26.     sum = 0.0
  27.     ix = iseed
  28.     do 300 i=1,iset
  29.       call    drand(ix,rnum)
  30.       ix = ix * rnum
  31.       irand(i) = ix
  32.       rand(i) = rnum
  33.       isum = isum + ix
  34.       sum = sum + rnum
  35.   300    continue
  36.     write(1,400) (irand(i),rand(i),i=1,iset)
  37.   400      format(1h ,4(i10,f10.4))
  38.     imean = isum /float(iset)
  39.     xmean = sum / iset
  40.     write(1,500) imean,xmean
  41. 500    format(1h0,'Mean: ',I10,2x,f6.4)
  42.     stop
  43.     end
  44.  
  45.