C Typical problem setup for Q1DA
C
         A = 0.0
         B = 1.0          
C Set interval endpoints to [0,1]
         EPS = 0.001       
C Set accuracy request for three digits
         CALL Q1DA (A,B,EPS,R,E,KF,IFLAG)
         WRITE(*,*)'Q1DA RESULTS: A, B, EPS, R, E, KF, IFLAG'
         WRITE(*,'(3F7.4,2E16.8,2I4)') A,B,EPS,R,E,KF,IFLAG
         WRITE(*,*)
         WRITE(*,*)'REFERENCE RESULTS ON IBM PC/AT'
         WRITE(*,*)'0.0000 1.0000 0.0010  0.41406750E-01',
     *'  0.55339400E-06  30   0'
         END
C
         FUNCTION F(X)
C Define integrand F
         F = SIN(2.*X)-SQRT(X)     
         RETURN 
         END