1 SUBROUTINE dqagp(F,A,B,NPTS2,POINTS,EPSABS,EPSREL,RESULT,ABSERR,
2 * neval,ier,leniw,lenw,last,iwork,work)
187 DOUBLE PRECISION A,ABSERR,B,EPSABS,EPSREL,POINTS,RESULT,WORK
188 INTEGER IER,IWORK,LAST,LENIW,LENW,LIMIT,LVL,L1,L2,L3,L4,NEVAL,
191 dimension iwork(leniw),points(npts2),work(lenw)
203 IF(leniw.LT.(3*npts2-2).OR.lenw.LT.(leniw*2-npts2).OR.npts2.LT.2)
208 limit = (leniw-npts2)/2
214 CALL
dqagpe(f,a,b,npts2,points,epsabs,epsrel,limit,result,abserr,
215 * neval,ier,work(1),work(l1),work(l2),work(l3),work(l4),
216 * iwork(1),iwork(l1),iwork(l2),last)
221 10
IF(ier.EQ.6) lvl = 1
222 IF(ier.GT.0) CALL
xerror(
'ABNORMAL RETURN FROM DQAGP',26,ier,lvl)
subroutine xerror(MESSG, NMESSG, NERR, LEVEL)
subroutine dqagp(F, A, B, NPTS2, POINTS, EPSABS, EPSREL, RESULT, ABSERR, NEVAL, IER, LENIW, LENW, LAST, IWORK, WORK)
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 dqagpe(F, A, B, NPTS2, POINTS, EPSABS, EPSREL, LIMIT, RESULT, ABSERR, NEVAL, IER, ALIST, BLIST, RLIST, ELIST, PTS, IORD, LEVEL, NDIN, LAST)