1 SUBROUTINE intdy (T, K, YH, NYH, DKY, IFLAG)
5 1 icf, ierpj, iersl, jcur, jstart, kflag, l, meth, miter,
6 2 maxord, maxcor, msbp, mxncf, n, nq, nst, nfe, nje, nqu
7 INTEGER I, IC, J, JB, JB2, JJ, JJ1, JP1
8 DOUBLE PRECISION T, YH, DKY
9 DOUBLE PRECISION ROWNS,
10 1 ccmax, el0, h, hmin, hmxi, hu, rc, tn, uround
11 DOUBLE PRECISION C, R, S, TP
13 COMMON /ls0001/ rowns(209),
14 2 ccmax, el0, h, hmin, hmxi, hu, rc, tn, uround,
15 3 iownd(14), iowns(6),
16 4 icf, ierpj, iersl, jcur, jstart, kflag, l, meth, miter,
17 5 maxord, maxcor, msbp, mxncf, n, nq, nst, nfe, nje, nqu
39 IF (k .LT. 0 .OR. k .GT. nq) go to 80
40 tp = tn - hu - 100.0d0*uround*(tn + hu)
41 IF ((t-tp)*(t-tn) .GT. 0.0d0) go to 90
45 IF (k .EQ. 0) go to 15
52 IF (k .EQ. nq) go to 55
58 IF (k .EQ. 0) go to 35
64 40 dky(i) = c*yh(i,jp1) + s*dky(i)
72 80 CALL
xerrwd(
'INTDY-- K (=I1) ILLEGAL ',
73 1 30, 51, 0, 1, k, 0, 0, 0.0d0, 0.0d0)
76 90 CALL
xerrwd(
'INTDY-- T (=R1) ILLEGAL ',
77 1 30, 52, 0, 0, 0, 0, 1, t, 0.0d0)
79 1
' T NOT IN INTERVAL TCUR - HU (= R1) TO TCUR (=R2) ',
80 1 60, 52, 0, 0, 0, 0, 2, tp, tn)
std::string dimension(void) const
subroutine xerrwd(MSG, NMES, NERR, LEVEL, NI, I1, I2, NR, R1, R2)
subroutine intdy(T, K, YH, NYH, DKY, IFLAG)