SUBROUTINE DQDAG(F,A,B,ERRABS,ERRREL,IRULE,RESULT,ERREST)
      IMPLICIT DOUBLE PRECISION(A-H,O-Z)
      PARAMETER (LIMIT=500,LENW=3000)
      DIMENSION WORK(LENW),IWORK(LIMIT) 
      EXTERNAL F
      CALL DQAG(F,A,B,ERRABS,ERRREL,IRULE,RESULT,ABSERR,NEVAL,IER,
     +LIMIT,LENW,LAST,IWORK,WORK)
      RETURN
      END