29 #define dls001_1 (mdls001_._1)
30 #define dls001_2 (mdls001_._2)
31 #define dls001_3 (mdls001_._3)
33 #define dlsa01_1 (mdlsa01_._1)
34 #define dlsa01_2 (mdlsa01_._2)
35 #define dlsa01_3 (mdlsa01_._3)
37 double d_sign(
const double & a,
const double & b);
49 #define pow_di(__x, __n) pow(*__x, (double) *__n)
57 C_INT yh_dim1, yh_offset, i__1, i__2;
113 yh_offset = 1 + yh_dim1;
128 if ((*t - tp) * (*t -
dls001_1.tn) > 0.)
144 for (jj = jj1; jj <= i__1; ++jj)
154 for (i__ = 1; i__ <= i__1; ++i__)
157 dky[i__] = c__ * yh[i__ +
dls001_1.l * yh_dim1];
168 for (jb = 1; jb <= i__1; ++jb)
182 for (jj = jj1; jj <= i__2; ++jj)
192 for (i__ = 1; i__ <= i__2; ++i__)
195 dky[i__] = c__ * yh[i__ + jp1 * yh_dim1] + s * dky[i__];
211 for (i__ = 1; i__ <= i__1; ++i__)
214 dky[i__] = r__ * dky[i__];
220 msg =
"DINTDY- K (=I1) illegal ";
226 msg =
"DINTDY- T (=R1) illegal ";
227 mxerrwd(msg, &
c__30, &
c__52, &
c__0, &
c__0, &
c__0, &
c__0, &
c__1, t, &
c_b34,
229 msg =
" T not in interval TCUR - HU (= R1) to TCUR (=R2) ";
C_INT dintdy_(double *t, const C_INT *k, double *yh, C_INT *nyh, double *dky, C_INT *iflag)
double d_sign(const double &a, const double &b)