naginterfaces.library.zeros.cubic_​complex

naginterfaces.library.zeros.cubic_complex(u, r, s, t)

cubic_complex determines the roots of a cubic equation with complex coefficients.

For full information please refer to the NAG Library document for c02am

https://support.nag.com/numeric/nl/nagdoc_30/flhtml/c02/c02amf.html

Parameters

ucomplex

$u$, the coefficient of $z^3$.

rcomplex

$r$, the coefficient of $z^2$.

scomplex

$s$, the coefficient of $z$.

tcomplex

$t$, the constant coefficient.

Returns

zerocomplex, ndarray, shape $\left(3\right)$

$zero[i-1]$ contains the $i$th root.

errestfloat, ndarray, shape $\left(3\right)$

$errest[i-1]$ contains an approximate error estimate for the $i$th root.

Raises

NagValueError

(errno $1$)

On entry, $u=(0.0,0.0)$.

Constraint: $u\ne (0.0,0.0)$

(errno $2$)

The companion matrix $H$ cannot be formed without overflow.

(errno $3$)

Failure to converge in lapackeig.zhseqr.

Notes

No equivalent traditional C interface for this routine exists in the NAG Library.

cubic_complex attempts to find the roots of the cubic equation

$$u{z}^3+r{z}^2+sz+t=0\text{,}$$

where $u$, $r$, $s$ and $t$ are complex coefficients with $u\ne 0$. The roots are located by finding the eigenvalues of the associated $3\times 3$ (upper Hessenberg) companion matrix $H$ given by

$$\begin{array}{c}H=\left(\begin{array}{ccc}0& 0& -t/u\\ 1& 0& -s/u\\ 0& 1& -r/u\end{array}\right)\text{.}\end{array}$$

The eigenvalues are obtained by a call to lapackeig.zhseqr. Further details can be found in Further Comments.

To obtain the roots of a quadratic equation, quadratic_complex() can be used.

References

Golub, G H and Van Loan, C F, 1996, Matrix Computations, (3rd Edition), Johns Hopkins University Press, Baltimore