1: $e$ – complex scalar: $e$ , the coefficient of $z^{4}$ .

Constraint: $e \neq (0.0, 0.0)$ .
2: $a$ – complex scalar: $a$ , the coefficient of $z^{3}$ .
3: $b$ – complex scalar: $b$ , the coefficient of $z^{2}$ .
4: $c$ – complex scalar: $c$ , the coefficient of $z$ .
5: $d$ – complex scalar: $d$ , the constant coefficient.

Optional Input Parameters

None.

Output Parameters

1: $zeror (4)$ – double array
2: $zeroi (4)$ – double array: $zeror (i)$ and $zeroi (i)$ contain the real and imaginary parts, respectively, of the $i$ th root.
3: $errest (4)$ – double array: $errest (i)$ contains an approximate error estimate for the $i$ th root.
4: $ifail$ – int64int32nag_int scalar: $ifail = 0$ unless the function detects an error (see Error Indicators and Warnings).

Error Indicators and Warnings

Errors or warnings detected by the function:

$ifail = 1$

On entry,

e = (0.0, 0.0)

$ifail = 2$: The companion matrix $H$ cannot be formed without overflow.

$ifail = 3$: The iterative procedure used to determine the eigenvalues has failed to converge.

$ifail = - 99$: An unexpected error has been triggered by this routine. Please contact NAG.

$ifail = - 399$: Your licence key may have expired or may not have been installed correctly.

$ifail = - 999$: Dynamic memory allocation failed.

Accuracy

ifail = 0

on exit, then the

i

th computed root should have approximately

|\log_{10} (errest (i))|

correct significant digits.

Further Comments

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

(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).

Example

This example finds the roots of the quartic equation

z^{4} + 16 i z^{2} - (8 - 8 i) z - 65 = 0 .

Open in the MATLAB editor: c02an_example

function c02an_example


fprintf('c02an example results\n\n');

e =  complex(1);
a =   -10 + 2i;
b =    48 - 10i;
c =  -100 + 28i;
d =  complex(96);
[zr, zi, errest, ifail] = c02an(e, a, b, c, d);

fprintf('  Roots of quartic    error estimates\n');
for j = 1:4
   if (zi(j)<0)
     fprintf('%8.4f - %7.4fi     %8.2e\n',zr(j),abs(zi(j)),errest(j));
   else
     fprintf('%8.4f - %7.4fi     %8.2e\n',zr(j),abs(zi(j)),errest(j));
   end
end

c02an example results

  Roots of quartic    error estimates
  3.0000 -  3.0000i     1.55e-14
  1.0000 -  1.0000i     1.18e-14
  2.0000 -  2.0000i     2.14e-14
  4.0000 -  4.0000i     2.02e-14

PDF version (NAG web site, 64-bit version, 64-bit version)

Chapter Contents

Chapter Introduction

NAG Toolbox