1 SUBROUTINE dqagi(F,BOUND,INF,EPSABS,EPSREL,RESULT,ABSERR,NEVAL,
2 * ier,
limit,lenw,last,iwork,work)
158 DOUBLE PRECISION ABSERR,BOUND,EPSABS,EPSREL,RESULT,WORK
159 INTEGER IER,INF,IWORK,LAST,LENW,LIMIT,LVL,L1,L2,L3,NEVAL
173 IF(limit.LT.1.OR.lenw.LT.limit*4) go
to 10
181 CALL
dqagie(f,bound,inf,epsabs,epsrel,limit,result,abserr,
182 * neval,ier,work(1),work(l1),work(l2),work(l3),iwork,last)
187 10
IF(ier.EQ.6) lvl = 1
188 IF(ier.GT.0) CALL
xerror(
'ABNORMAL RETURN FROM DQAGI',26,ier,lvl)
subroutine xerror(MESSG, NMESSG, NERR, LEVEL)
may be zero for pure relative error test tem the relative tolerance must be greater than or equal to
OCTAVE_EXPORT octave_value_list etc The functions then dimension(columns)
subroutine dqagie(F, BOUND, INF, EPSABS, EPSREL, LIMIT, RESULT, ABSERR, NEVAL, IER, ALIST, BLIST, RLIST, ELIST, IORD, LAST)
subroutine dqagi(F, BOUND, INF, EPSABS, EPSREL, RESULT, ABSERR, NEVAL, IER, LIMIT, LENW, LAST, IWORK, WORK)