1 subroutine qagi(f,bound,inf,epsabs,epsrel,result,abserr,neval,
2 * ier,limit,lenw,last,iwork,work)
158 real abserr, epsabs,epsrel,result,work
159 integer ier,iwork, lenw,limit,lvl,l1,l2,l3,neval
173 if(limit.lt.1.or.lenw.lt.limit*4) go
to 10
181 call
qagie(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.ne.0) call
xerror(
'abnormal return from qagi',26,ier,lvl)
subroutine xerror(MESSG, NMESSG, NERR, LEVEL)
subroutine qagi(f, bound, inf, epsabs, epsrel, result, abserr, neval, ier, limit, 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 qagie(f, bound, inf, epsabs, epsrel, limit, result, abserr, neval, ier, alist, blist, rlist, elist, iord, last)