home *** CD-ROM | disk | FTP | other *** search
/ Programmer 7500 / MAX_PROGRAMMERS.iso / INFO / FORTRAN / LLSQ.ZIP / GEN.FOR < prev    next >
Encoding:
Text File  |  1984-02-23  |  2.2 KB  |  29 lines
  1.       FUNCTION   GEN(ANOISE)                                            GEN00100
  2. C     C.L.LAWSON AND R.J.HANSON, JET PROPULSION LABORATORY, 1972 DEC 15 GEN00200
  3. C     TO APPEAR IN 'SOLVING LEAST SQUARES PROBLEMS', PRENTICE-HALL, 1974GEN00300
  4. C          GENERATE NUMBERS FOR CONSTRUCTION OF TEST CASES.             GEN00400
  5.       IF (ANOISE) 10,30,20                                              GEN00500
  6.    10 MI=891                                                            GEN00600
  7.       MJ=457                                                            GEN00700
  8.       I=5                                                               GEN00800
  9.       J=7                                                               GEN00900
  10.       AJ=0.                                                             GEN01000
  11.       GEN=0.                                                            GEN01100
  12.       RETURN                                                            GEN01200
  13. C                                                                       GEN01300
  14. C     THE SEQUENCE OF VALUES OF J  IS BOUNDED BETWEEN 1 AND 996         GEN01400
  15. C     IF INITIAL J = 1,2,3,4,5,6,7,8, OR 9, THE PERIOD IS 332           GEN01500
  16.    20 J=J*MJ                                                            GEN01600
  17.       J=J-997*(J/997)                                                   GEN01700
  18.       AJ=J-498                                                          GEN01800
  19. C     THE SEQUENCE OF VALUES OF I  IS BOUNDED BETWEEN 1 AND 999         GEN01900
  20. C     IF INITIAL I = 1,2,3,6,7, OR 9,  THE PERIOD WILL BE 50            GEN02000
  21. C     IF INITIAL I = 4 OR 8   THE PERIOD WILL BE 25                     GEN02100
  22. C     IF INITIAL I = 5        THE PERIOD WILL BE 10                     GEN02200
  23.    30 I=I*MI                                                            GEN02300
  24.       I=I-1000*(I/1000)                                                 GEN02400
  25.       AI=I-500                                                          GEN02500
  26.       GEN=AI+AJ*ANOISE                                                  GEN02600
  27.       RETURN                                                            GEN02700
  28.       END                                                               GEN02800
  29.