home *** CD-ROM | disk | FTP | other *** search
/ Programmer 7500 / MAX_PROGRAMMERS.iso / INFO / FORTRAN / LLSQ.ZIP / PROG4.FOR < prev    next >
Encoding:
Text File  |  1984-02-26  |  1.9 KB  |  26 lines
  1. C     PROG4                                                             PR400100
  2. C     C.L.LAWSON AND R.J.HANSON, JET PROPULSION LABORATORY, 1973 JUN 12 PR400200
  3. C     TO APPEAR IN 'SOLVING LEAST SQUARES PROBLEMS', PRENTICE-HALL, 1974PR400300
  4. C          DEMONSTRATE SINGULAR VALUE ANALYSIS.                         PR400400
  5. C                                                                       PR400500
  6.       DIMENSION A(15,5),B(15),SING(15)                                  PR400600
  7.       CHARACTER*4 NAMES
  8.       DATA NAMES/'    '/                                                PR400700
  9. C                                                                       PR400800
  10. C Set math to single precision format (non-IEEE)
  11.       CALL MPBRQQ
  12.       READ (5,10) ((A(I,J),J=1,5),B(I),I=1,15)                          PR400900
  13.       WRITE (6,20)                                                      PR401000
  14.       WRITE (6,30) ((A(I,J),J=1,5),B(I),I=1,15)                         PR401100
  15.       WRITE (6,40)                                                      PR401200
  16. C                                                                       PR401300
  17.       CALL SVA (A,15,15,5,15,B,SING,NAMES,1,D)                          PR401400
  18. C                                                                       PR401500
  19.       STOP                                                              PR401600
  20.    10 FORMAT (6F12.0)                                                   PR401700
  21.    20 FORMAT (46H1PROG4.    DEMONSTRATE SINGULAR VALUE ANALYSIS/53H LISTPR401800
  22.      1ING OF INPUT MATRIX, A, AND VECTOR, B, FOLLOWS..)                 PR401900
  23.    30 FORMAT (1H /(5F12.8,F20.4))                                       PR402000
  24.    40 FORMAT (1H1)                                                      PR402100
  25.       END                                                               PR402200
  26.