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)
std::string dimension(void) const
subroutine dqagpe(F, A, B, NPTS2, POINTS, EPSABS, EPSREL, LIMIT, RESULT, ABSERR, NEVAL, IER, ALIST, BLIST, RLIST, ELIST, PTS, IORD, LEVEL, NDIN, LAST)