1: $e$ – double Input: On entry: $e$ , the coefficient of $z^{4}$ .

Constraint: $e \neq 0.0$ .
2: $a$ – double Input: On entry: $a$ , the coefficient of $z^{3}$ .
3: $b$ – double Input: On entry: $b$ , the coefficient of $z^{2}$ .
4: $c$ – double Input: On entry: $c$ , the coefficient of $z$ .
5: $d$ – double Input: On entry: $d$ , the constant coefficient.
6: $zeror [4]$ – double Output
7: $zeroi [4]$ – double Output: On exit: $zeror [i - 1]$ and $zeroi [i - 1]$ contain the real and imaginary parts, respectively, of the $i$ th root.
8: $errest [4]$ – double Output: On exit: $errest [i - 1]$ contains an approximate error estimate for the $i$ th root.
9: $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_C02_NOT_CONV: The iterative procedure used to determine the eigenvalues has failed to converge.
NE_C02_OVERFLOW: The companion matrix $H$ cannot be formed without overflow.
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.
NE_REAL: On entry, $e = 0.0$ .
Constraint: $e \neq 0.0$ .

7 Accuracy

fail = NE_NOERROR

on exit, then the

i

th computed root should have approximately

|\log_{10} (errest [i - 1])|

correct significant digits.

c02alc is not threaded in any implementation.

The method used by the function consists of the following steps, which are performed by functions from LAPACK.

(a)Form matrix $H$ .
(b)Apply a diagonal similarity transformation to $H$ (to give $H^{'}$ ).
(c)Calculate the eigenvalues and Schur factorization of $H^{'}$ .
(d)Calculate the left and right eigenvectors of $H^{'}$ .
(e)Estimate reciprocal condition numbers for all the eigenvalues of $H^{'}$ .
(f)Calculate approximate error estimates for all the eigenvalues of $H^{'}$ (using the 1-norm).

To find the roots of the quartic equation

z^{4} + 2 z^{3} + 6 z^{2} - 8 z - 40 = 0 .

Interfaces: FL CL AD

c02al: FL CL AD