5 SUBROUTINE dnsid(X,Y,YPRIME,NEQ,ICOPT,ID,RES,WT,RPAR,IPAR,
6 * DELTA,R,YIC,YPIC,WM,IWM,CJ,EPCON,RATEMX,MAXIT,STPTOL,
7 * ICNFLG,ICNSTR,IERNEW)
77 IMPLICIT DOUBLE PRECISION(a-h,o-z)
84 parameter(lnni=19, llsoff=35)
96 CALL dslvd (neq, delta, wm, iwm)
100 delnrm = ddwnrm(neq,delta,wt,rpar,ipar)
102 IF (fnrm .LE. epcon)
RETURN 107 iwm(lnni) = iwm(lnni) + 1
113 CALL dlinsd (neq, y, x, yprime, cj, delta, delnrm, wt, lsoff,
114 * stptol, iret, res, ires, wm, iwm, fnrm, icopt, id,
115 * r, yic, ypic, icnflg, icnstr, rlx, rpar, ipar)
120 IF (iret .NE. 0)
GO TO 390
124 IF (fnrm .LE. epcon)
RETURN 130 IF (m .GE. maxit)
GO TO 380
135 CALL dcopy (neq, r, 1, delta, 1)
141 380
IF (rate .LE. ratemx)
THEN 148 390
IF (ires .LE. -2)
THEN OCTAVE_EXPORT octave_value_list etc The functions then dimension(columns)