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)
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)
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)