Date of Award
1-2019
Document Type
Dissertation
Degree Name
Doctor of Philosophy (PhD)
Department
Mathematical Sciences
First Advisor
Charles T. Fulton
Second Advisor
Mark Archambault
Third Advisor
Dennis Jackson
Fourth Advisor
Tariel Kiguradze
Abstract
This dissertation is mainly concerned with three new and very different numerical methods for computing eigenvalues of fourth-order Sturm-Liouville differential operators on finite intervals with regular boundary conditions. The first method is a fourth-order shooting method based on Magnus expansions (MG4). This method is similar to the SLEUTH method of Greenberg and Marletta [34] which uses MG2 shooting for the integrator (also known as second-order Pruess method [71]). The second method utilizes two first-order systems of nonlinear Riccati equations, which are integrated alternately to avoid integrating into singularities (conjugate points) when shooting across the Sturm-Liouville interval, and this provides good stabilization for the numerical integration, a fact also noticed by Scott, Shampine and Wing [77] for second-order equations. The third method, presented here for the first time, converts to a linear 7×7 first-order system which computes the quantities whose zeros determine the eigenvalues for a certain special class of boundary conditions; it also exhibits relatively good numerical stability. All three methods can very often achieve high accuracy on the order of machine precision, and some comparisons of their performance against the well known SLEUTH software package [34] are presented.
Recommended Citation
Alalyani, Ahmad Ali, "Numerical Methods for Eigenvalue Computation of Fourth-Order Self-Adjoint Ordinary Differential Operators" (2019). Theses and Dissertations. 985.
https://repository.fit.edu/etd/985