*DECK RKFAB SUBROUTINE RKFAB (NCOMP, XPTS, NXPTS, NFC, IFLAG, Z, MXNON, P, + NTP, IP, YHP, NIV, U, V, W, S, STOWA, G, WORK, IWORK, NFCC) C***BEGIN PROLOGUE RKFAB C***SUBSIDIARY C***PURPOSE Subsidiary to BVSUP C***LIBRARY SLATEC C***TYPE SINGLE PRECISION (RKFAB-S, DRKFAB-D) C***AUTHOR Watts, H. A., (SNLA) C***DESCRIPTION C C ********************************************************************** C C Subroutine RKFAB integrates the initial value equations using C the variable-step RUNGE-KUTTA-FEHLBERG integration scheme or C the variable-order ADAMS method and orthonormalization C determined by a linear dependence test. C C ********************************************************************** C C***SEE ALSO BVSUP C***ROUTINES CALLED BVDER, DEABM, DERKF, REORT, STOR1 C***COMMON BLOCKS ML15TO, ML17BW, ML18JR, ML8SZ C***REVISION HISTORY (YYMMDD) C 750601 DATE WRITTEN C 890831 Modified array declarations. (WRB) C 890921 Realigned order of variables in certain COMMON blocks. C (WRB) C 891214 Prologue converted to Version 4.0 format. (BAB) C 900328 Added TYPE section. (WRB) C 910722 Updated AUTHOR section. (ALS) C***END PROLOGUE RKFAB C DIMENSION P(NTP,*),IP(NFCC,*),U(NCOMP,NFC,*), 1 V(NCOMP,*),W(NFCC,*),Z(*),YHP(NCOMP,*), 2 XPTS(*),S(*),STOWA(*),WORK(*),IWORK(*), 3 G(*) C C ********************************************************************** C COMMON /ML8SZ/ C,XSAV,IGOFX,INHOMO,IVP,NCOMPD,NFCD COMMON /ML15TO/ PX,PWCND,TND,X,XBEG,XEND,XOT,XOP,INFO(15),ISTKOP, 1 KNSWOT,KOP,LOTJP,MNSWOT,NSWOT COMMON /ML18JR/ AE,RE,TOL,NXPTSD,NIC,NOPG,MXNOND,NDISK,NTAPE,NEQ, 1 INDPVT,INTEG,NPS,NTPD,NEQIVP,NUMORT,NFCCD, 2 ICOCO COMMON /ML17BW/ KKKZPW,NEEDW,NEEDIW,K1,K2,K3,K4,K5,K6,K7,K8,K9, 1 K10,K11,L1,L2,KKKINT,LLLINT C EXTERNAL BVDER C C ********************************************************************** C INITIALIZATION OF COUNTERS AND VARIABLES. C