1 SUBROUTINE sewset (N, ITOL, RTOL, ATOL, YCUR, EWT)
25 REAL RTOL, ATOL, YCUR, EWT
26 dimension rtol(*), atol(*), ycur(n), ewt(n)
29 go
to(10, 20, 30, 40), itol
32 15 ewt(i) = rtol(1)*
abs(ycur(i)) + atol(1)
36 25 ewt(i) = rtol(1)*
abs(ycur(i)) + atol(i)
40 35 ewt(i) = rtol(i)*
abs(ycur(i)) + atol(1)
44 45 ewt(i) = rtol(i)*
abs(ycur(i)) + atol(i)
subroutine sewset(N, ITOL, RTOL, ATOL, YCUR, EWT)
may be zero for pure relative error test tem the relative tolerance must be greater than or equal to
OCTAVE_EXPORT octave_value_list etc The functions then dimension(columns)
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))<