298 research outputs found
A matrix method for fractional Sturm-Liouville problems on bounded domain
A matrix method for the solution of direct fractional Sturm-Liouville
problems on bounded domain is proposed where the fractional derivative is
defined in the Riesz sense. The scheme is based on the application of the
Galerkin spectral method of orthogonal polynomials. The order of convergence of
the eigenvalue approximations with respect to the matrix size is studied. Some
numerical examples that confirm the theory and prove the competitiveness of the
approach are finally presented
Large-scale computation of pseudospectra using ARPACK and eigs
ARPACK and its MATLAB counterpart, eigs, are software packages that calculate some eigenvalues of a large non-symmetric matrix by Arnoldi iteration with implicit restarts. We show that at a small additional cost, which diminishes relatively as the matrix dimension increases, good estimates of pseudospectra in addition to eigenvalues can be obtained as a by-product. Thus in large-scale eigenvalue calculations it is feasible to obtain routinely not just eigenvalue approximations, but also information as to whether or not the eigenvalues are likely to be physically significant. Examples are presented for matrices with dimension up to 200,000
- …