C     TITLE: DEMONSTRATION OF COMPUTING X=J*H
C
C     TWO METHODS ARE USED TO COMPUTE X=J*H,
C     USING (3.40) AND (3.41) IN THE TEXT.
C     SEE SECTION 3.3 FOR MORE INFORMATION.
C
      DOUBLE PRECISION DX, DH
      DATA N/100/, H/0.1/
C
C     INITIALIZE.
      XSUM = 0.0
      DH = H
      PRINT *, 'FOLLOWING ARE THE VALUES OF J, X=J*H, AND'
      PRINT *, 'THE ERRORS IN THE COMPUTATION OF X USING'
      PRINT *, '(3.40) AND (3.41), RESPECTIVELY.'
      PRINT *, ' '
C
      DO 10 J=1,N
        XSUM = XSUM + H
        XPROD = J*H
        DX = J*DH
        ERRSUM = DX - XSUM
        ERRPRD = DX - XPROD
        PRINT *, J, DX, ERRPRD, ERRSUM
 10     CONTINUE
      STOP
      END