On entry: the type of error control. At each step in the numerical solution an estimate of the local error,
est, is made. For the current step to be accepted the following condition must be satisfied:
where
τr and
τa are defined by
where
ε is a small machine-dependent number and
ei is an estimate of the local error at
yi, computed internally. If the appropriate condition is not satisfied, the step size is reduced and the solution is recomputed on the current step. If you wish to measure the error in the computed solution in terms of the number of correct decimal places, then
RELABS should be set to 'A' on entry, whereas if the error requirement is in terms of the number of correct significant digits, then
RELABS should be set to 'R'. If you prefer a mixed error test, then
RELABS should be set to 'M', otherwise if you have no preference,
RELABS should be set to the default 'D'. Note that in this case 'D' is taken to be 'M'.