In this paper we present and test a numerical method for computing eigenvalues of 4th order Sturm-Liouville (SL) differential operators on finite intervals with regular boundary conditions. This method is a 4th order ...In this paper we present and test a numerical method for computing eigenvalues of 4th order Sturm-Liouville (SL) differential operators on finite intervals with regular boundary conditions. This method is a 4th order shooting method based on Magnus expansions (MG4) which use MG4 shooting as the integrator. This method is similar to the SLEUTH (Sturm-Liouville Eigenvalues Using Theta Matrices) method of Greenberg and Marletta which uses the 2nd order Pruess method (also known as the MG2 shooting method) for the integrator. This method often achieves near machine precision accuracies, and some comparisons of its performance against the well-known SLEUTH software package are presented.展开更多
Based on some recent results for interlacing eigenvalue intervals from 1-parameter families of se- quences of eigenvalue inequalities, a new method is given to solving the index problem for Sturm-Liouville eigenvalues...Based on some recent results for interlacing eigenvalue intervals from 1-parameter families of se- quences of eigenvalue inequalities, a new method is given to solving the index problem for Sturm-Liouville eigenvalues for coupled self-adjoint boundary conditions in the regular case. The key is a new characteristic principle for indices for Sturm-Liouville eigenvalues. The algorithm corresponding on the characteristic princi- ple are discussed, and numerical examples are presented to illustrate the theoretical results and show that the algorithm is valid.展开更多
文摘In this paper we present and test a numerical method for computing eigenvalues of 4th order Sturm-Liouville (SL) differential operators on finite intervals with regular boundary conditions. This method is a 4th order shooting method based on Magnus expansions (MG4) which use MG4 shooting as the integrator. This method is similar to the SLEUTH (Sturm-Liouville Eigenvalues Using Theta Matrices) method of Greenberg and Marletta which uses the 2nd order Pruess method (also known as the MG2 shooting method) for the integrator. This method often achieves near machine precision accuracies, and some comparisons of its performance against the well-known SLEUTH software package are presented.
基金Supported by the National Natural Science Foundation of China(No.11361039 and 11161030)the Natural Science Foundation of Inner Mongolia Province,China(No.2013MS0116)
文摘Based on some recent results for interlacing eigenvalue intervals from 1-parameter families of se- quences of eigenvalue inequalities, a new method is given to solving the index problem for Sturm-Liouville eigenvalues for coupled self-adjoint boundary conditions in the regular case. The key is a new characteristic principle for indices for Sturm-Liouville eigenvalues. The algorithm corresponding on the characteristic princi- ple are discussed, and numerical examples are presented to illustrate the theoretical results and show that the algorithm is valid.
