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)
OCTAVE_EXPORT octave_value_list etc The functions then dimension(columns)