1 SUBROUTINE intdy (T, K, YH, NYH, DKY, IFLAG)
4 INTEGER ILLIN, INIT, LYH, LEWT, LACOR, LSAVF, LWM, LIWM,
5 1 mxstep, mxhnil, nhnil, ntrep, nslast, cnyh,
6 2 ialth, ipup, lmax, meo, nqnyh, nslp
7 INTEGER ICF, IERPJ, IERSL, JCUR, JSTART, KFLAG, L, METH, MITER,
8 2 maxord, maxcor, msbp, mxncf, n, nq, nst, nfe, nje, nqu
9 INTEGER I, IC, J, JB, JB2, JJ, JJ1, JP1
10 DOUBLE PRECISION T, YH, DKY
11 DOUBLE PRECISION CONIT, CRATE, EL, ELCO, HOLD, RMAX, TESCO,
12 1 ccmax, el0, h, hmin, hmxi, hu, rc, tn, uround
13 DOUBLE PRECISION C, R, S, TP
15 COMMON /ls0001/ conit, crate, el(13), elco(13,12),
16 1 hold, rmax, tesco(3,12),
17 2 ccmax, el0, h, hmin, hmxi, hu, rc, tn, uround,
18 2 illin, init, lyh, lewt, lacor, lsavf, lwm, liwm,
19 3 mxstep, mxhnil, nhnil, ntrep, nslast, cnyh,
20 3 ialth, ipup, lmax, meo, nqnyh, nslp,
21 4 icf, ierpj, iersl, jcur, jstart, kflag, l, meth, miter,
22 5 maxord, maxcor, msbp, mxncf, n, nq, nst, nfe, nje, nqu
44 IF (k .LT. 0 .OR. k .GT. nq) go
to 80
45 tp = tn - hu - 100.0d0*uround*(tn + hu)
46 IF ((t-tp)*(t-tn) .GT. 0.0d0) go
to 90
50 IF (k .EQ. 0) go
to 15
57 IF (k .EQ. nq) go
to 55
63 IF (k .EQ. 0) go
to 35
69 40 dky(i) = c*yh(i,jp1) + s*dky(i)
77 80 CALL
xerrwd(
'INTDY-- K (=I1) ILLEGAL ',
78 1 30, 51, 0, 1, k, 0, 0, 0.0d0, 0.0d0)
81 90 CALL
xerrwd(
'INTDY-- T (=R1) ILLEGAL ',
82 1 30, 52, 0, 0, 0, 0, 1, t, 0.0d0)
84 1
' T NOT IN INTERVAL TCUR - HU (= R1) TO TCUR (=R2) ',
85 1 60, 52, 0, 0, 0, 0, 2, tp, tn)
subroutine xerrwd(MSG, NMES, NERR, LEVEL, NI, I1, I2, NR, R1, R2)
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 intdy(T, K, YH, NYH, DKY, IFLAG)