|
PROGRAM TEST
|
IMPLICIT NONE
|
C
|
INTEGER N
|
PARAMETER (N = 3)
|
C
|
INTEGER ICOL, INCX, IROW
|
DOUBLE PRECISION SCALE, X(N,N)
|
C
|
EXTERNAL DSCAL
|
C
|
C Initialize the array X to store the vector x shown below.
|
C
|
C 5.0 5.5 6.0
|
C x = 2.5 3.0 3.5
|
C 1.0 1.5 2.0
|
C
|
DATA X / 5.0D0, 2.5D0, 1.0D0, 5.5D0, 3.0D0, 1.5D0,
|
$ 6.0D0, 3.5D0, 2.0D0 /
|
C
|
PRINT 1000
|
PRINT 1010, ((X(IROW,ICOL), ICOL = 1, N), IROW = 1, N)
|
INCX = N
|
SCALE = X(1,1) / X(2,1)
|
CALL DSCAL (N, SCALE, X(2,1), INCX)
|
PRINT 1020
|
PRINT 1010, ((X(IROW,ICOL), ICOL = 1, N), IROW = 1, N)
|
|
|
C
|
1000 FORMAT (1X, 'Before scaling row 2:')
|
1010 FORMAT (3(3X, F3.1))
|
1020 FORMAT (/1X, 'After scaling row 2:')
|
C
|
END
|
|