QUANC8

!     B ПPИMEPE BЫЧИCЛЯETCЯ ИHTEГPAЛ OT ФУHKЦИИ
!     SIN(X)/X ПO X OT 0 ДO 2

!      REAL FUNCTION FUN(X)
!      REAL X
!      IF(X .EQ. 0.0) FUN=1.
!      IF(X .NE. 0.0) FUN=SIN(X)/X
!      RETURN
!      END

!      EXTERNAL FUN
!      INTEGER NOFUN
!      REAL FUN,A,B,RELERR,ABSERR,RESULT,ERREST,FLAG
!      DATA A/0./,B/2.0/,RELERR/1.E-06/,ABSERR/0.0/
!
!      CALL QUANC8(FUN,A,B,ABSERR,RELERR,RESULT, &
!                 ERREST,NOFUN,FLAG)
!
!      PRINT 1,RESULT,ERREST,NOFUN,FLAG
!      STOP
!    1 FORMAT(10X,'RESULT=',E14.7,3X,'ERREST=',E12.5/ &
!         11X,'NOFUN=',I8,11X,'FLAG=',F10.3)
!      END

!     B ПPИMEPE BЫЧИCЛЯETCЯ ИHTEГPAЛ OT ФУHKЦИИ
!     SIN(X)/X ПO X OT 0 ДO 2

      REAL FUNCTION FUN(X)
        REAL X
        FUN=1/(exp(1.)**(1.9*X**2) + X_H)
        RETURN
      END

      EXTERNAL FUN
        INTEGER NOFUN
        REAL FUN,A,B,RELERR,ABSERR,RESULT,ERREST,FLAG,  X_H

        DATA A/0./,B/1.0/,RELERR/1.E-06/,ABSERR/0.0/
        X_H = 0.

        DO I = 1,6
            X_H = X_H + 0.5
            CALL QUANC8(FUN,A,B,ABSERR,RELERR,RESULT, ERREST,NOFUN,FLAG)
            PRINT 1,RESULT,ERREST,NOFUN,FLAG

            DATA A/0./,B/1.0/,RELERR/1.E-06/,ABSERR/0.0/
            NOFUN = 0
            FLAG = 0
        END DO

        STOP
      1 FORMAT(10X,'RESULT=',E14.7,3X,'ERREST=',E12.5/ &
            11X,'NOFUN=',I8,11X,'FLAG=',F10.3)
      END