1: $m$ – Integer Input: On entry: $m$ , the number of terms in the continued fraction.

Constraint: $m \geq 1$ .
2: $a [m]$ – const double Input: On entry: $a [j - 1]$ must be set to the value of the parameter $a_{j}$ in the continued fraction, for $j = 1, 2, \dots, m$ .
3: $u [m]$ – const double Input: On entry: $u [j - 1]$ must be set to the value of the parameter $u_{j}$ in the continued fraction, for $j = 1, 2, \dots, m - 1$ . (The element $u [m - 1]$ is not used).
4: $x$ – double Input: On entry: the value of $x$ at which the continued fraction is to be evaluated.
5: $f$ – double * Output: On exit: the value of the continued fraction corresponding to the value of $x$ .
6: $fail$ – NagError * Input/Output: The NAG error argument (see Section 7 in the Introduction to the NAG Library CL Interface).

6 Error Indicators and Warnings

NE_ALLOC_FAIL: Dynamic memory allocation failed.
See Section 3.1.2 in the Introduction to the NAG Library CL Interface for further information.
NE_BAD_PARAM: On entry, argument $⟨ value ⟩$ had an illegal value.
NE_INTERNAL_ERROR: An internal error has occurred in this function. Check the function call and any array sizes. If the call is correct then please contact NAG for assistance.
See Section 7.5 in the Introduction to the NAG Library CL Interface for further information.
NE_NO_LICENCE: Your licence key may have expired or may not have been installed correctly.
See Section 8 in the Introduction to the NAG Library CL Interface for further information.
NE_POLE_PRESENT: x corresponds to a pole of $R (x)$ , or is very close. $x = ⟨ value ⟩$ .

7 Accuracy

See Section 7 in e01rac.

e01rbc is not threaded in any implementation.

The time taken by e01rbc is approximately proportional to

m

This example reads in the arguments

a_{j}

and

u_{j}

of a continued fraction (as determined by the example for e01rac) and evaluates the continued fraction at a point

x

e01rb: FL CL CPP AD