DDCST(3)

SYNOPSIS
SUBROUTINE DDCST (MAXORD, MINT, ISWFLG, EL, TQ)

NAME
ddcst - DDCST sets coefficients used by the core integrator DDSTP.

LIBRARY
SLATEC (SDRIVE)

TYPE
DOUBLE PRECISION (SDCST-S, DDCST-D, CDCST-C)

AUTHOR
Kahaner, D. K., (NIST) National Institute of Standards and Technology Gaithersburg, MD 20899 Sutherland, C. D., (LANL) Mail Stop D466 Los Alamos National Laboratory Los Alamos, NM 87545

DESCRIPTION
DDCST is called by DDNTL. The array EL determines the basic method. The array TQ is involved in adjusting the step size in relation to truncation error. EL and TQ depend upon MINT, and are calculated for orders 1 to MAXORD(.LE. 12). For each order NQ, the coefficients EL are calculated from the generating polynomial: L(T) = EL(1,NQ) + EL(2,NQ)*T + ... + EL(NQ+1,NQ)*T**NQ. For the implicit Adams methods, L(T) is given by dL/dT = (1+T)*(2+T)* ... *(NQ-1+T)/K, L(-1) = 0, where K = factorial(NQ-1). For the Gear methods, L(T) = (1+T)*(2+T)* ... *(NQ+T)/K, where K = factorial(NQ)*(1 + 1/2 + ... + 1/NQ). For each order NQ, there are three components of TQ.

ROUTINES CALLED
(NONE)

REVISION HISTORY (YYMMDD)
790601 DATE WRITTEN
900329 Initial submission to SLATEC.