1 SUBROUTINE sintdy (T, K, YH, NYH, DKY, IFLAG)
46 1 icf, ierpj, iersl, jcur, jstart, kflag, l,
47 2 lyh, lewt, lacor, lsavf, lwm, liwm, meth, miter,
48 3 maxord, maxcor, msbp, mxncf, n, nq, nst, nfe, nje, nqu
50 1 ccmax, el0, h, hmin, hmxi, hu, rc, tn, uround
51 COMMON /sls001/ rowns(209),
52 1 ccmax, el0, h, hmin, hmxi, hu, rc, tn, uround,
54 3 icf, ierpj, iersl, jcur, jstart, kflag, l,
55 4 lyh, lewt, lacor, lsavf, lwm, liwm, meth, miter,
56 5 maxord, maxcor, msbp, mxncf, n, nq, nst, nfe, nje, nqu
57 INTEGER I, IC, J, JB, JB2, JJ, JJ1, JP1
63 IF (k .LT. 0 .OR. k .GT. nq) go to 80
64 tp = tn - hu - 100.0e0*uround*sign(
abs(tn) +
abs(hu), hu)
65 IF ((t-tp)*(t-tn) .GT. 0.0e0) go to 90
69 IF (k .EQ. 0) go to 15
76 IF (k .EQ. nq) go to 55
82 IF (k .EQ. 0) go to 35
88 40 dky(i) = c*yh(i,jp1) + s*dky(i)
96 80 CALL
xerrwd(
'SINTDY- K (=I1) illegal ',
97 1 30, 51, 0, 1, k, 0, 0, 0.0e0, 0.0e0)
100 90 CALL
xerrwd(
'SINTDY- T (=R1) illegal ',
101 1 30, 52, 0, 0, 0, 0, 1, t, 0.0e0)
103 1
' T not in interval TCUR - HU (= R1) to TCUR (=R2) ',
104 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 sintdy(T, K, YH, NYH, DKY, IFLAG)