home *** CD-ROM | disk | FTP | other *** search
/ Programmer 7500 / MAX_PROGRAMMERS.iso / INFO / FORTRAN / SUPERT87.ZIP / RAND.FOR < prev    next >
Encoding:
Text File  |  1986-12-15  |  1.0 KB  |  19 lines
  1. c      SUBROUTINE RAND(IX,IY,YFL)                                        3484
  2. C     ...RAND. NO. ROUTINE FROM IBM 360 SSP                             3485
  3. C     ...IX - ANY ODD INTEGER <= 9 DIGITS ON FIRST CALL.                3486
  4. C     ...AFTERWARDS, IX=PREVIOUS IY.                                    3487
  5. c      IY=IX*65539                                                       3488
  6. c      IF(IY)5,6,6                                                       3489
  7. c 5    IY=IY+2147483647+1                                                3490
  8. c 6    YFL=IY                                                            3491
  9. c      YFL=YFL*.4656613E-9                                               3492
  10. c      IX=IY                                                             3493
  11. c      RETURN                                                            3494
  12. c      END                                                               3495
  13.  
  14.       function rand(xy)
  15.       n=xy*1048576.
  16.       rand=float(mod(1061*n+221589,1048576))/1048576.
  17.       return
  18.       end
  19.