1 SUBROUTINE ewset (N, ITOL, RTOL, ATOL, YCUR, EWT)
11 DOUBLE PRECISION RTOL, ATOL, YCUR, EWT
12 dimension rtol(*), atol(*), ycur(n), ewt(n)
14 go
to(10, 20, 30, 40), itol
17 15 ewt(i) = rtol(1)*dabs(ycur(i)) + atol(1)
21 25 ewt(i) = rtol(1)*dabs(ycur(i)) + atol(i)
25 35 ewt(i) = rtol(i)*dabs(ycur(i)) + atol(1)
29 45 ewt(i) = rtol(i)*dabs(ycur(i)) + atol(i)
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)
subroutine ewset(N, ITOL, RTOL, ATOL, YCUR, EWT)