Analytic vs. numerical solutions to a Sturm-Liouville transmission eigenproblem




Sturm-Liouville eigenproblem, discontinuous coefficient, transmission condition, spectral collocation, FEM, Finite Element Method, accuracy
Abstract views: 380


An elliptic one-dimensional second order boundary value problem involving discontinuous coefficients, with or without transmission conditions, is considered. For the former case by a direct sum spaces method we show that the eigenvalues are real, geometrically simple and the eigenfunctions are orthogonal.

Then the eigenpairs are computed numerically by a local linear finite element method (FEM) and by some global spectral collocation methods.

The spectral collocation is based on Chebyshev polynomials (ChC) for problems on bounded intervals respectively on Fourier system (FsC) for periodic problems.

The numerical stability in computing eigenvalues is investigated by estimating their (relative) drift with respect to the order of approximation. The accuracy in computing the eigenvectors is addressed by estimating their departure from orthogonality as well as by the asymptotic order of convergence. The discontinuity of coefficients in the problems at hand reduces the exponential order of convergence, usual for any well designed spectral algorithm, to an algebraic one.

As expected, the accuracy of ChC outcomes overpasses by far that of FEM outcomes.


Download data is not yet available.


K. Aydemir and O.S. Mukhtarov, Qualitative analysis of eigenvalues and eigenfunctions of one boundary value-transmission problem, Bound. Value Probl. 82 (2016). DOI:

I. Babuska and J.E. Osborn, Numerical Treatment of Eigenvalue Problems for Differential Equations with Discontinuous Coefficients, Math. Comput. 32 (1978), pp. 991–1023, DOI:

P.B. Bailey, M.K. Gordon and L.F. Shampine, Automatic Solution of the Sturm-Liouville Problem, ACM. T. Math. Software. 4 (1978), pp. 193–208, DOI:

P.B. Bailey, W.N. Everitt and A. Zettl, Algorithm 810: The SLEIGN2 Sturm-Liouville Code, ACM T. Math. Software. 27 (2001) pp. 143–192, DOI:

J.P. Boyd, Chebyshev and Fourier spectral methods, 2nd rev. ed. Mineola, NY: Dover Publications. 2001.

P.A.M. Boomkamp, B.J. Boersma, R.H.M. Miesen and G.V. Beijnon, A Chebyshev Collocation Method for Solving Two-Phase Flow Stability Problems, J. Comput. Phys. 132 (1997), pp. 191–200, DOI:

S.C. Brenner, L.R. Scott, The mathematical theory of finite element methods, 3rd ed. Texts in Applied Mathematics 15. New-York, NY: Springer 2008, DOI:

P.G. Ciarlet, The Finite Element Method for Elliptic Problems, North-Holland, Amsterdam, New-York, Oxford 1978. DOI:

W.N. Everitt, A. Zettl, Sturm–Liouville differential operators in direct sum spaces, Rocky Mountain J. Math. 16 (1986), pp. 497–516, DOI:

F. Gesztesyt, C. Macdeo and L. Streit, An exactly solvable periodic Schrodinger operator J. Phys. A: Math. Gen. 18 (1985), pp. L503–L507, DOI:

C. I. Gheorghiu, J. Rommes, Application of the Jacobi-Davidson method to accurate analysis of singular linear hydrodynamic stability problems, Int. J. Numer. Meth. Fl. 71 (2012), pp. 358–369, DOI:

C.I. Gheorghiu, Spectral Methods for Non-Standard Eigenvalue Problems. Fluid and Structural Mechanics and Beyond, Springer Cham Heidelberg New York Dondrecht London 2014, DOI:

Y. He, D.P. Nicholls and J. Shen, An efficient and stable spectral method for electromagnetic scattering from a layered periodic sructure, J. Comput. Phys. 231 (2012), pp. 3007–3022, DOI:

C.O. Horgan, J.P. Spence, A. N. Andry, Lower bounds for eigenvalues of Sturm-Liouville problems with discontinuous coefficients: integral equation methods, Q. Appl. Math. 39 (1982), pp. 455–465, DOI:

Multi-precision Computing Toolbox for MATLAB. Yokohama: Advanpix LLC.; 2008-2017.

M. Marletta, J.D. Pryce, Automatic solution of Sturm-Liouville problems using the Pruess method, J. Comput. Appl. Math. 39 (1992), pp. 57–78, DOI:

O.S. Mukhtarov, M. Kadakal and F.S. Muhtarov, On discontinuous Sturm-Liouville problems with transmission conditions, J. Math. Kyoto Univ. (JMKYAZ) 44 (2004), pp. 779–798, DOI:

J. Necas, Direct Methods in the Theory of Elliptic Equations, Springer Monographs in Mathematics. Berlin 2012, DOI:

S. Pruess, C.T. Fulton, Mathematical Software for Sturm-Liouville Problems, ACM T. Math. Software. 19 (1993), pp. 360–376, DOI:

J. D. Pryce, Numerical solutions of Sturm-Liouville problem, Oxford University Press,Oxford, U. K., 1993.

D.G. Shepelsky, The inverse problem of reconstruction of the medium’s conductivity in a class of discontinuous and increasing functions, Adv. Soviet Math. 19 (1994), pp.209–231. DOI:

I. Titeux, Ya. Yakubov, Completeness of root functions for thermal conduction in a strip with piecewise continuous coefficients Math. Models Methods Appl. Sc., 7 (1997), pp. 1035–1050, DOI:

J.A.C. Weideman, S. C. Reddy, A MATLAB Differentiation Matrix Suite, ACM T. Math. Software. 26 (2000), pp. 465–519, DOI:

A. Zettl, Adjoint and Self-Adjoint Boundary Value Problems with Interface Conditions, SIAM J. Appl. Math. 16 (1968), pp. 851–859, DOI:

A. Zettl, Sturm—Liouville Theory, Math. Surveys Monogr., vol. 121, Amer. Math. Soc., Providence, RI 2005.

H. Yserentant, A Short Theory of the Rayleigh-Ritz Method, C.M.A.M. 13 (2013) pp. 496–502, DOI:




How to Cite

Gheorghiu, C.-I., & Zinsou, B. (2019). Analytic vs. numerical solutions to a Sturm-Liouville transmission eigenproblem. J. Numer. Anal. Approx. Theory, 48(2), 159–174.