2 DOUBLE PRECISION FUNCTION dcsevl (X, CS, N)
38 DOUBLE PRECISION B0, B1, B2, CS(*), ONEPL, TWOX, X, D1MACH
43 IF (first) onepl = 1.0d0 + d1mach(4)
45 IF (n .LT. 1) CALL
xermsg(
'SLATEC',
'DCSEVL',
46 + .LE.
'NUMBER OF TERMS 0', 2, 2)
47 IF (n .GT. 1000) CALL
xermsg(
'SLATEC',
'DCSEVL',
48 + .GT.
'NUMBER OF TERMS 1000', 3, 2)
49 IF (
abs(x) .GT. onepl) CALL
xermsg(
'SLATEC',
'DCSEVL',
50 +
'X OUTSIDE THE INTERVAL (-1,+1)', 1, 1)
59 b0 = twox*b1 - b2 + cs(ni)
62 dcsevl = 0.5d0*(b0-b2)
subroutine xermsg(LIBRAR, SUBROU, MESSG, NERR, LEVEL)
OCTAVE_EXPORT octave_value_list return the value of the option it must match the dimension of the state and the relative tolerance must also be a vector of the same length tem it must match the dimension of the state and the absolute tolerance must also be a vector of the same length The local error test applied at each integration step is xample roup abs(local error in x(i))<