derkfs.f
SUBROUTINE DERKFS (F, NEQ, T, Y, TOUT, INFO, RTOL, ATOL, IDID, H,
+ TOLFAC, YP, F1, F2, F3, F4, F5, YS, TOLD, DTSIGN, U26, RER,
+ INIT, KSTEPS, KOP, IQUIT, STIFF, NONSTF, NTSTEP, NSTIFS, RPAR,
+ IPAR)
C***BEGIN PROLOGUE DERKFS
....
....
Warning: this routine is not intended to be user-callable.
....
....
C***SUBSIDIARY
C***PURPOSE Subsidiary to DERKF
C***LIBRARY SLATEC
C***TYPE SINGLE PRECISION (DERKFS-S, DRKFS-D)
C***AUTHOR Watts, H. A., (SNLA)
C***DESCRIPTION
C
C Fehlberg Fourth-Fifth order Runge-Kutta Method
C **********************************************************************
C
C DERKFS integrates a system of first order ordinary differential
C equations as described in the comments for DERKF .
C
C The arrays YP,F1,F2,F3,F4,F5,and YS (of length at least NEQ)
C appear in the call list for variable dimensioning purposes.
C
C The variables H,TOLFAC,TOLD,DTSIGN,U26,RER,INIT,KSTEPS,KOP,IQUIT,
C STIFF,NONSTF,NTSTEP, and NSTIFS are used internally by the code
C and appear in the call list to eliminate local retention of
C variables between calls. Accordingly, these variables and the
C array YP should not be altered.
C Items of possible interest are
C H - An appropriate step size to be used for the next step
C TOLFAC - Factor of change in the tolerances
C YP - Derivative of solution vector at T
C KSTEPS - Counter on the number of steps attempted
C
C **********************************************************************
C
C***SEE ALSO DERKF
C***ROUTINES CALLED DEFEHL, HSTART, HVNRM, R1MACH, XERMSG
C***REVISION HISTORY (YYMMDD)
C 800501 DATE WRITTEN
C 890531 Changed all specific intrinsics to generic. (WRB)
C 890831 Modified array declarations. (WRB)
C 891024 Changed references from VNORM to HVNRM. (WRB)
C 891024 REVISION DATE from Version 3.2
C 891214 Prologue converted to Version 4.0 format. (BAB)
C 900328 Added TYPE section. (WRB)
C 900510 Convert XERRWV calls to XERMSG calls, replace GOTOs with
C IF-THEN-ELSEs. (RWC)
C 910722 Updated AUTHOR section. (ALS)
C***END PROLOGUE DERKFS