5 SUBROUTINE ddasic (X, Y, YPRIME, NEQ, ICOPT, ID, RES, JAC, PSOL,
6 * H, WT, NIC, IDID, RPAR, IPAR, PHI, SAVR, DELTA, E, YIC, YPIC,
7 * PWK, WM, IWM, HMIN, UROUND, EPLI, SQRTN, RSQRTN, EPCONI,
8 * STPTOL, JFLG, ICNFLG, ICNSTR, NLSIC)
83 IMPLICIT DOUBLE PRECISION(a-h,o-z)
84 dimension y(*),yprime(*),id(*),wt(*),phi(neq,*)
85 dimension savr(*),delta(*),e(*),yic(*),ypic(*),pwk(*)
86 dimension wm(*),iwm(*), rpar(*),ipar(*), icnstr(*)
87 EXTERNAL res, jac, psol, nlsic
98 DATA rhcut/0.1d0/, ratemx/0.8d0/
113 IF (nic .EQ. 2) jskip = 1
114 CALL dcopy (neq, y, 1, phi(1,1), 1)
115 CALL dcopy (neq, yprime, 1, phi(1,2), 1)
117 IF (icopt .EQ. 2)
THEN 130 CALL nlsic(x,y,yprime,neq,icopt,id,res,jac,psol,h,wt,jskip,
131 * rpar,ipar,savr,delta,e,yic,ypic,pwk,wm,iwm,cj,uround,
132 * epli,sqrtn,rsqrtn,epconi,ratemx,stptol,jflg,icnflg,icnstr,
135 IF (iernls .EQ. 0)
RETURN 148 iwm(lcfn) = iwm(lcfn) + 1
151 IF (iernls .EQ. -1)
GO TO 350
152 IF (icopt .EQ. 2)
GO TO 350
153 IF (nh .EQ. mxnh)
GO TO 350
159 IF (iernls .EQ. 1)
GO TO 200
161 CALL dcopy (neq, phi(1,1), 1, y, 1)
162 CALL dcopy (neq, phi(1,2), 1, yprime, 1)
OCTAVE_EXPORT octave_value_list etc The functions then dimension(columns)