1 subroutine qagp(f,a,b,npts2,points,epsabs,epsrel,result,abserr,
2 * neval,ier,leniw,lenw,last,iwork,work)
187 real a,abserr,b,epsabs,epsrel,points,result,work
188 integer ier,iwork,leniw,lenw,limit,lvl,l1,l2,l3,neval,npts2
190 dimension iwork(leniw),points(npts2),work(lenw)
202 if(leniw.lt.(3*npts2-2).or.lenw.lt.(leniw*2-npts2).or.npts2.lt.2)
207 limit = (leniw-npts2)/2
213 call
qagpe(f,a,b,npts2,points,epsabs,epsrel,limit,result,abserr,
214 * neval,ier,work(1),work(l1),work(l2),work(l3),work(l4),
215 * iwork(1),iwork(l1),iwork(l2),last)
220 10
if(ier.eq.6) lvl = 1
221 if(ier.ne.0) call
xerror(
'abnormal return from qagp',26,ier,lvl)
subroutine xerror(MESSG, NMESSG, NERR, LEVEL)
std::string dimension(void) const
subroutine qagp(f, a, b, npts2, points, epsabs, epsrel, result, abserr, neval, ier, leniw, lenw, last, iwork, work)
subroutine qagpe(f, a, b, npts2, points, epsabs, epsrel, limit, result, abserr, neval, ier, alist, blist, rlist, elist, pts, iord, level, ndin, last)