38 REAL B0, B1, B2, CS(*), ONEPL, TWOX, X
43 IF (first) onepl = 1.0e0 + r1mach(4)
45 IF (n .LT. 1) CALL
xermsg(
'SLATEC',
'CSEVL',
46 + .LE.
'NUMBER OF TERMS 0', 2, 2)
47 IF (n .GT. 1000) CALL
xermsg(
'SLATEC',
'CSEVL',
48 + .GT.
'NUMBER OF TERMS 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)
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))<