2
FUNCTION CSEVL (X, CS, N)
3
C***BEGIN PROLOGUE CSEVL
4
C***PURPOSE Evaluate a Chebyshev series.
5
C***LIBRARY SLATEC (FNLIB)
7
C***TYPE SINGLE PRECISION (CSEVL-S, DCSEVL-D)
8
C***KEYWORDS CHEBYSHEV SERIES, FNLIB, SPECIAL FUNCTIONS
9
C***AUTHOR Fullerton, W., (LANL)
12
C Evaluate the N-term Chebyshev series CS at X. Adapted from
13
C a method presented in the paper by Broucke referenced below.
16
C X value at which the series is to be evaluated.
17
C CS array of N terms of a Chebyshev series. In evaluating
18
C CS, only half the first coefficient is summed.
19
C N number of terms in array CS.
21
C***REFERENCES R. Broucke, Ten subroutines for the manipulation of
22
C Chebyshev series, Algorithm 446, Communications of
23
C the A.C.M. 16, (1973) pp. 254-256.
24
C L. Fox and I. B. Parker, Chebyshev Polynomials in
25
C Numerical Analysis, Oxford University Press, 1968,
27
C***ROUTINES CALLED R1MACH, XERMSG
28
C***REVISION HISTORY (YYMMDD)
30
C 890831 Modified array declarations. (WRB)
31
C 890831 REVISION DATE from Version 3.2
32
C 891214 Prologue converted to Version 4.0 format. (BAB)
33
C 900315 CALLs to XERROR changed to CALLs to XERMSG. (THJ)
34
C 900329 Prologued revised extensively and code rewritten to allow
35
C X to be slightly outside interval (-1,+1). (WRB)
36
C 920501 Reformatted the REFERENCES section. (WRB)
37
C***END PROLOGUE CSEVL
38
REAL B0, B1, B2, CS(*), ONEPL, TWOX, X
42
C***FIRST EXECUTABLE STATEMENT CSEVL
43
IF (FIRST) ONEPL = 1.0E0 + R1MACH(4)
45
IF (N .LT. 1) CALL XERMSG ('SLATEC', 'CSEVL',
46
+ 'NUMBER OF TERMS .LE. 0', 2, 2)
47
IF (N .GT. 1000) CALL XERMSG ('SLATEC', 'CSEVL',
48
+ 'NUMBER OF TERMS .GT. 1000', 3, 2)
49
IF (ABS(X) .GT. ONEPL) CALL XERMSG ('SLATEC', 'CSEVL',
50
+ 'X OUTSIDE THE INTERVAL (-1,+1)', 1, 1)
59
B0 = TWOX*B1 - B2 + CS(NI)