Abstract:La Budde's method computes the characteristic polynomial of a real matrix A in two stages: first it applies orthogonal similarity transformations to reduce A to upper Hessenberg form H, and second it computes the characteristic polynomial of H from characteristic polynomials of leading principal submatrices of H. If A is symmetric, then H is symmetric tridiagonal, and La Budde's method simplifies to the Sturm sequence method. If A is diagonal then La Budde's method reduces to the Summation Algorithm, a Horner-like scheme used by the MATLAB function POLY to compute characteristic polynomials from eigenvalues. We present recursions to compute the individual coefficients of the characteristic polynomial in the second stage of La Budde's method, and derive running error bounds for symmetric and nonsymmetric matrices. We also show that La Budde's method can be more accurate than POLY, especially for indefinite and nonsymmetric matrices A. Unlike POLY, La Budde's method is not affected by illconditioning of eigenvalues, requires only real arithmetic, and allows the computation of individual coefficients.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to develop a new method for computing the characteristic polynomial of a real matrix, namely the La Budde method. Specifically, the La Budde method achieves this goal through two stages: 1. **First stage**: Use an orthogonal similarity transformation to reduce a real matrix \(A\) to an upper Hessenberg form \(H\). If \(A\) is symmetric, then \(H\) will be further simplified to a symmetric tridiagonal matrix. 2. **Second stage**: Compute the characteristic polynomial of the upper Hessenberg matrix \(H\) through a recursive formula. Specifically, it gradually determines the characteristic polynomial of the entire matrix \(H\) by computing the characteristic polynomials of the leading principal sub - matrices of \(H\). ### Key problems and advantages - **Numerical stability**: The La Budde method remains numerically stable during the reduction to Hessenberg form and is not affected by the eigenvalue condition number. - **Computational efficiency**: This method requires only approximately \(5n^3\) floating - point operations and only real - number arithmetic operations, without the need for complex - number arithmetic. - **Selective computation**: It can efficiently compute individual coefficients of the characteristic polynomial, which is very important for applications that only require some coefficients (such as quantum physics). - **Accuracy**: Compared with the `poly` function in MATLAB, the La Budde method has higher accuracy for indefinite matrices and non - symmetric matrices. ### Formula representation The characteristic polynomial is defined as: \[p(\lambda)=\det(\lambda I - A)=\lambda^n + c_1\lambda^{n - 1}+\cdots+c_{n - 1}\lambda + c_n\] where \(I\) is the identity matrix, \(c_1 =-\text{trace}(A)\), \(c_n=(- 1)^n\det(A)\). For a symmetric matrix \(A\), the computation of its characteristic polynomial is simplified to the Sturm sequence method, and for a diagonal matrix \(A\), it is simplified to the Summation Algorithm, which is the basis of the `poly` function in MATLAB. ### Summary The paper mainly solves the problem of how to compute the characteristic polynomial of a real matrix more accurately and efficiently, especially when dealing with indefinite matrices and non - symmetric matrices.

La Budde's Method for Computing Characteristic Polynomials

On the application of Laguerre's method to the polynomial eigenvalue problem

Exchangeable arrays and integrable systems for characteristic polynomials of random matrices

On the characteristic polynomial of Cartan matrices and Chebyshev polynomials

On prescribed characteristic polynomials

Characteristic Polynomials Of The Matrices With -Entry

A MATLAB package computing simultaneous Gaussian quadrature rules for Multiple Orthogonal Polynomials

Generalized characteristic polynomials and Gaussian cubature rules

Average characteristic polynomials for multiple orthogonal polynomial ensembles

On Convoluted Forms of Multivariate Legendre-Hermite Polynomials with Algebraic Matrix Based Approach

On the condition of characteristic polynomials

Structured eigenvalue backward errors for rational matrix polynomials with symmetry structures

Finding roots of complex analytic functions via generalized colleague matrices

Generalized Golub-Kahan bidiagonalization for nonsymmetric saddle point systems

On the Hadamard and Vandermonde Determinants and the Bernoulli–Euler–Lagrange–Aitken Method for Calculating the Roots of Polynomials

CHARACTERISTIC POLYNOMIALS OF THE MATRICES WITH (j,k) -ENTRY Q J±k +T

A Diagrammatic Approach for the Coefficients of the Characteristic Polynomial

Fast and backward stable computation of the eigenvalues and eigenvectors of matrix polynomials

Bidiagonal matrix factorisations related to multiple orthogonal polynomials

A Riemann--Hilbert approach to the perturbation theory for orthogonal polynomials: Applications to numerical linear algebra and random matrix theory

Bidiagonal factorization of the recurrence matrix for the Hahn multiple orthogonal polynomials