Next Article in Journal
Use of Correlated Data for Nonparametric Prediction of a Spatial Target Variable
Next Article in Special Issue
Alikhanov Legendre—Galerkin Spectral Method for the Coupled Nonlinear Time-Space Fractional Ginzburg–Landau Complex System
Previous Article in Journal
Latent Class Regression Utilizing Fuzzy Clusterwise Generalized Structured Component Analysis
Previous Article in Special Issue
Existence and Multiplicity of Solutions to a Class of Fractional p-Laplacian Equations of Schrödinger-Type with Concave-Convex Nonlinearities in ℝN
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Solutions of Sturm-Liouville Problems

by
Upeksha Perera
1,*,† and
Christine Böckmann
2
1
Department of Mathematics, University of Kelaniya, Kelaniya 11600, Sri Lanka
2
Institut für Mathematik, Universität Potsdam, 14476 Potsdam, Germany
*
Author to whom correspondence should be addressed.
Current address: Institut für Mathematik, Universität Potsdam, 14476 Potsdam, Germany.
Mathematics 2020, 8(11), 2074; https://doi.org/10.3390/math8112074
Submission received: 27 October 2020 / Revised: 17 November 2020 / Accepted: 18 November 2020 / Published: 20 November 2020
(This article belongs to the Special Issue Mathematical Analysis and Boundary Value Problems)

Abstract

:
This paper further improves the Lie group method with Magnus expansion proposed in a previous paper by the authors, to solve some types of direct singular Sturm–Liouville problems. Next, a concrete implementation to the inverse Sturm–Liouville problem algorithm proposed by Barcilon (1974) is provided. Furthermore, computational feasibility and applicability of this algorithm to solve inverse Sturm–Liouville problems of higher order (for n = 2 , 4 ) are verified successfully. It is observed that the method is successful even in the presence of significant noise, provided that the assumptions of the algorithm are satisfied. In conclusion, this work provides a method that can be adapted successfully for solving a direct (regular/singular) or inverse Sturm–Liouville problem (SLP) of an arbitrary order with arbitrary boundary conditions.

1. Introduction

The inverse Sturm–Liouville theory was originated in 1929 by Ambarzumian [1] and further developed in [2,3,4,5,6].
The 2 m th order, nonsingular, self-adjoint eigenvalue problem (EVP) or Sturm–Liouville problem (SLP) is given by:
( 1 ) m ( p 0 ( x ) y ( m ) ) ( m ) + ( 1 ) m 1 ( p 1 ( x ) y ( m 1 ) ) ( m 1 ) + + ( p m 2 ( x ) y ) ( p m 1 ( x ) y ) + p m ( x ) y = λ w ( x ) y , a < x < b
together with some boundary conditions at a and b, the functions p k , ( 0 k m ) , and  w ( x ) being continuous on the finite closed interval [ a , b ] , and  p 0 having a continuous derivative. In the inverse SLP, the coefficient functions p k , ( 0 k m ) need to be reconstructed, given suitable valid spectral data.
For a discussion of analytical methods and numerical methods for inverse SLP, see [7,8], respectively. Iterative methods [9,10], Rayleigh–Ritz method [11], finite difference approximation [12], Quasi-Newton method [13], shooting method [14], interval Newton’s method [15], finite-difference method [16], boundary value methods [17,18,19,20], Numerov’s method  [21,22,23], least-squares functional  [24], generalized Rundell–Sacks algorithm [25,26], spectral mappings [27], Lie-group estimation method [28], Broyden method [29,30], decent flow methods [31], modified Numerov’s method [32], Newton-type method [33], Fourier–Legendre series [34], and  Chebyshev polynomials [35] are of particular importance among the existing methods to solve inverse SLP.
Numerical algorithms to solve the inverse fourth-order Sturm–Liouville problem (FSLP) are proposed in  [36,37,38].
In [39], authors implemented Barcilon [9]’s iterative algorithm to solve inverse SLPs successfully. The present paper extends the results for the case m = 2 by considering the general-order inverse SLP algorithm as described in [36].

2. Materials and Methods

Inverse SLP Algorithm for the General Order

Setting p 0 ( x ) , w ( x ) 1 in (1) and taking the domain [ 0 , 1 ] , the EVP reads as:
( y ( m ) ) ( m ) ( p 1 y ( m 1 ) ) ( m 1 ) + + p m ( x ) y = λ y , 0 < x < 1 .
According to [40], m + 1 distinct spectra are required in order to determine p 1 ( x ) , p 2 ( x ) , , p m ( x ) . In other words, m + 1 different EVPs with distinct boundary conditions (BC) are needed. In addition, the spectra need to be interlaced, i.e.,
< λ n 1 ( m + 1 ) < λ n ( 1 ) < λ n ( 2 ) < λ n ( 3 ) < < λ n ( m + 1 ) < λ n + 1 ( 1 ) <
Assume that m boundary conditions are common to all the given m + 1 eigenvalue problems, say, conditions at the right end point x = 1 .
  • Combine the m + 1 eigenvalue problems into an equivalent linear differential equation:
    1.1.
    Define a vector function ϕ ( x , λ ) :
    ϕ T ( x , λ ) = [ ϖ χ , ϖ χ , , ϖ ( m 1 ) χ ( m 1 ) ] = [ ϕ ( 1 ) , ϕ ( 2 ) , , ϕ ( m ) ]
    where ϖ ( x , λ ) and χ ( x , λ ) are two solutions of Equation (2).
    1.2.
    Differentiate each ϕ ( k ) , 2 m times.
    1.3.
    Obtain m ( 2 m + 1 ) linear equations for ϖ ( i ) χ ( j ) + ϖ ( j ) χ ( i ) , where i j , j = 0 , 1 , , 2 m 1 .
    1.4.
    Express ϖ ( i ) χ ( j ) + ϖ ( j ) χ ( i ) , i j , j = 0 , 1 , , 2 m 1 , in terms of ϕ ( 1 ) , ϕ ( 2 ) , , ϕ ( m ) and their first 2 m derivatives.
    1.5.
    Differentiate ϕ ( i ) ( 2 m ) , i = 1 , 2 , , m again, to obtain m linear coupled differential equations of order 2 m + 1 , which is in matrix form:
    M ϕ = λ N ϕ
  • Obtain the corresponding boundary conditions for the linear differential equation using the boundary conditions of the m + 1 eigenvalue problems:
    2.1.
    Assume that ϖ ( x , λ ) and χ ( x , λ ) satisfy the m boundary conditions common to all the given m + 1 eigenvalue problems.
    2.2.
    Then from 1.3, find the m ( 3 m + 1 ) / 2 boundary conditions at x = 1 , and the m ( m + 1 ) / 2 boundary conditions at x = 0 for ϕ ( x , λ ) .
  • Solve Equation (4) using above boundary conditions, to obtain the solutions ϕ n .
  • Solve the adjoint system of equations to the above system, and denote the solutions by η n .
  • Find the bi-orthogonal set of functions { y n } 1 to { ϕ n } 1 , using the relation y n ( x ) = N T η n ( x ) .
Some remarks on Algorithm 1 follow:
Line   (Data:)
Infinite eigenvalue (EV) sequence is replaced with the first N EVs. The truncated EV sequences { λ n ( i ) } n = 1 N , i = 1 , 2 , , m + 1 are obtained, using the Magnus method [39] to solve Equation (2).
Line   (Result:)
Output p : m × n matrix of function values, n: the number of subdivisions in the x-axis and m: the number of potentials.
Line   (1)
We take initial guess p ( 0 ) as 0 m × n .
Line (3)
while loop exit condition is replaced with a maximum number of iterations: k M A X and in each iteration the condition n = 1 ( m + 1 ) N | σ ^ n σ n ( k 1 ) | < t o l ( t o l : a fixed tolerance) is checked to exit the loop.
Line (4)
· denotes the integration w.r.t. x from 0 to 1 and e T = [ 1 , 0 ] . y n and ϕ n are kept fixed at y 0 and ϕ 0 , which are the solutions by setting p i = 0 , i . Solutions y 0 and ϕ 0 are calculated using Mathematica [41]. Matlab [42] built-in functions: trapz() and griddedInterpolant() with pchip (piecewise cubic Hermite interpolating polynomial) option are used to approximate integrals and to have p as a function, respectively. The latter is required since, p as a function is input to the Magnus method. Again, the eigenvalue sequence σ ( k ) is calculated using p ( k ) and the iteration repeats.
You may refer Section 2.6 of [39] for the specific computational details of the Algorithm 1 for the case m = 1 and next section for the case m = 2 .
Algorithm 1: Solving inverse Sturm–Liouville problem (SLP).
Mathematics 08 02074 i001

3. Inverse FSLP Algorithm

Setting m = 2 in Equation (2) we obtain FSLP:
L u u ( 4 ) ( p u ) + q u = λ u .
Three spectra { λ n } , { μ n } and { ν n } required to reconstruct the unknown coefficients: p and q will be obtained from the following EVPs:
L u n = λ n u n , u n ( 0 ) = u n ( 0 ) = u n ( 1 ) = u n ( 1 ) = 0 ,
L v n = μ n v n , v n ( 0 ) = v n ( 0 ) = v n ( 1 ) = v n ( 1 ) = 0 ,
L w n = ν n w n , w n ( 0 ) = w n ( 0 ) = w n ( 1 ) = w n ( 1 ) = 0 .
Although, Equations (9) and (10) are self-adjoint, Equation (11) is not. Its adjoint equation is:
L ω n = ν n ω n , ω n ( 0 ) = ω n ( 0 ) p ( 0 ) ω n ( 0 ) = ω n ( 1 ) = ω n ( 1 ) = 0 .
The eigenvalues λ n , μ n and ν n (for p = q = 0 ) are given by the solutions of the following equations:
sin ( s ) = 0 , s 4 = λ
tan ( r ) tanh ( r ) = 0 , r 4 = μ
tan ( t ) + tanh ( t ) = 0 , t 4 = ν
Furthermore,
t n 1 4 π , s = n π , r n + 1 4 π , n N
so that the spectra { λ n } , { μ n } and { ν n } are interlaced, i.e.,
0 < ν 1 < λ 1 < μ 1 < ν 2 < λ 2 < μ 2 <
Since the eigenvalues are interlaced, they can be coalesced into a single one by defining
σ k , ϕ k ψ k 1 = ν n , w n ω n w n ω n 1 for k = 3 n 2 , λ n , u n 2 u n 2 1 for k = 3 n 1 , μ n , v n 2 v n 2 1 for k = 3 n , n = 1 , 2 , .
Denote:
q = q p , ϕ n = ϕ n ψ n .
Let χ ( x , λ ) and ϖ ( x , λ ) be two solutions of (8) satisfying the common boundary conditions of Equations (9)–(12), namely, ( u ( 1 ) = u ( 1 ) = 0 , etc.)
χ ( 1 , λ ) = ϖ ( 1 , λ ) = χ ( 1 , λ ) = ϖ ( 1 , λ ) = 0 .
Guided by (16), define
ϕ ( x , λ ) = ϖ ( x , λ ) χ ( x , λ ) ψ ( x , λ ) = ϖ ( x , λ ) χ ( x , λ ) .
By differentiating ϕ and ψ five times each, and eliminating ϖ , and χ , 9th order linear differential equation for ϕ (or a coupled system of linear differential equations of order 5 for ϕ and ψ ) can be obtained. When p = q = 0 , this reduces to
ϕ ( 9 ) ( x ) 12 λ ϕ ( 5 ) ( x ) = 64 λ 2 ϕ ( 1 ) ( x )
or
3 D 5 10 D 3 0 D 5 ϕ ψ = λ 8 D 0 10 D 3 24 D ϕ ψ
in the form of the Equation (4), and the boundary conditions read as:
ϕ = ϕ = ψ = 0 at x = 0 ,
ϕ = ϕ = ϕ 2 ψ = ϕ = ϕ 4 ψ = ψ = ψ = 0 at x = 1 .
The solution of the system of Equations (21)–(23) using Mathematica is:
ϕ ( x ) = csc ( μ ) 8 ( 2 cosh ( μ ) + cosh ( 2 μ ) ( 3 cos ( μ ) + cos ( 3 μ ) cos ( ( 2 + i ) μ ) + cos ( ( 1 + 2 i ) μ ) + cosh ( ( 2 + i ) μ ) cosh ( ( 1 + 2 i ) μ ) ) ) × ( i cos ( ( 1 + i ) ( ( 1 + i ) x + 1 ) μ ) i cos ( ( 1 + i ) ( ( 1 + i ) x μ + μ ) ) + 16 sin ( 2 μ ) 4 sin ( 4 μ ) ( 10 + 10 i ) sin ( ( 1 + i ) μ ) + ( 2 2 i ) sin ( ( 3 + i ) μ ) ( 4 + 4 i ) sin ( ( 2 + 2 i ) μ ) ( 2 2 i ) sin ( ( 1 + 3 i ) μ ) + ( 2 + 2 i ) sin ( ( 3 + 3 i ) μ ) 2 i sin ( ( 1 + i ) ( ( 3 + i ) x ) μ ) ( 6 + 6 i ) sin ( ( 1 + i ) ( x 2 ) μ ) + ( 2 + 2 i ) sin ( ( 1 + i ) ( x + 2 ) μ ) + 6 i sin ( ( 1 + i ) ( ( 1 i ) + x ) μ ) + 2 sin ( ( 1 + i ) ( ( 3 + i ) + x ) μ ) + 4 sin ( ( 1 + i ) ( 2 i + x ) μ ) + 4 i sin ( ( 1 + i ) ( 2 i + x ) μ ) ( 2 + 4 i ) sin ( ( 1 + i ) ( ( 1 + i ) + x ) μ ) 2 i sin ( ( 1 + i ) ( ( 1 + 3 i ) + x ) μ ) + ( 2 2 i ) sin ( ( 1 + i ) ( ( 2 2 i ) + x ) μ ) 2 sin ( ( 1 + i ) ( ( 1 3 i ) + x ) μ ) + 2 i sin ( ( ( 1 + i ) + 2 x ) μ ) 2 i sin ( ( ( 1 + 3 i ) + 2 x ) μ ) 6 i sin ( ( 2 i + ( 1 i ) x ) μ ) + ( 2 + 2 i ) sin ( ( 4 i + ( 1 i ) x ) μ ) 2 sin ( ( ( 2 + 4 i ) + ( 1 i ) x ) μ ) + 2 sin ( 2 ( i x + 1 ) μ ) + 4 i sin ( ( 1 + i ) ( i x + 2 ) μ ) + 8 sin ( ( ( 1 i ) + 2 i x ) μ ) 2 sin ( ( ( 3 + i ) + 2 i x ) μ ) + 2 sin ( ( ( 1 3 i ) + 2 i x ) μ ) 2 sin ( ( ( 3 3 i ) + 2 i x ) μ ) + sin ( ( 1 + i ) ( ( 1 + i ) x + 1 ) μ ) ( 4 + 2 i ) sin ( ( ( 1 + i ) x + 2 ) μ ) + 8 i sin ( ( ( 1 + i ) 2 x ) μ ) + 2 i sin ( ( ( 3 + i ) 2 x ) μ ) 2 i sin ( ( ( 3 + 3 i ) 2 x ) μ ) 6 sin ( ( 2 ( 1 i ) x ) μ ) + ( 2 + 2 i ) sin ( ( 4 ( 1 i ) x ) μ ) + 2 sin ( 2 ( 1 i x ) μ ) + 4 sin ( ( 1 + i ) ( 2 i x ) μ ) 8 sin ( 2 ( ( 1 + i ) i x ) μ ) + 2 sin ( 2 ( ( 2 + i ) i x ) μ ) + 2 sin ( 2 ( ( 1 + 2 i ) i x ) μ ) + 2 sin ( ( ( 1 i ) 2 i x ) μ ) 2 sin ( ( ( 3 i ) 2 i x ) μ ) + 8 sin ( ( ( 1 + i ) 2 i x ) μ ) + 2 sin ( ( ( 1 + 3 i ) 2 i x ) μ ) 2 sin ( ( ( 3 + 3 i ) 2 i x ) μ ) 6 sin ( ( 2 ( 1 + i ) x ) μ ) + ( 2 2 i ) sin ( ( 4 ( 1 + i ) x ) μ ) ( 4 2 i ) sin ( ( 1 i ) x μ + 2 μ ) + sin ( ( 1 + i ) ( ( 1 + i ) x μ + μ ) ) 8 i sin ( ( 2 + 2 i ) μ 2 x μ ) + 2 i sin ( ( 4 + 2 i ) μ 2 x μ ) + 2 i sin ( ( 2 + 4 i ) μ 2 x μ ) 16 sinh ( 2 μ ) + 4 sinh ( 4 μ ) + ( 10 + 10 i ) sinh ( ( 1 + i ) μ ) ( 2 2 i ) sinh ( ( 3 + i ) μ ) + ( 4 + 4 i ) sinh ( ( 2 + 2 i ) μ ) + ( 2 2 i ) sinh ( ( 1 + 3 i ) μ ) ( 2 + 2 i ) sinh ( ( 3 + 3 i ) μ ) + ( 6 + 6 i ) sinh ( ( 1 + i ) ( x 2 ) μ ) ( 2 + 2 i ) sinh ( ( 1 + i ) ( x + 2 ) μ ) 2 i sinh ( ( 1 + i ) ( ( 3 i ) + x ) μ ) + 2 i sinh ( 2 ( ( 2 i ) + x ) μ ) 8 i sinh ( 2 ( ( 1 i ) + x ) μ ) 2 sinh ( ( 1 + i ) ( ( 3 + i ) + x ) μ ) + 2 i sinh ( 2 ( ( 1 2 i ) + x ) μ ) + 2 sinh ( ( 1 + i ) ( ( 1 3 i ) + x ) μ ) 2 sinh ( 2 ( i x + 1 ) μ ) + 2 sinh ( ( ( 3 + i ) + 2 i x ) μ ) 2 sinh ( ( 1 + i ) ( ( 1 + i ) x + 1 ) μ ) + ( 4 + 2 i ) sinh ( ( ( 1 + i ) x + 2 ) μ ) 2 sinh ( 2 ( 1 i x ) μ ) + 8 sinh ( 2 ( ( 1 + i ) i x ) μ ) 2 sinh ( 2 ( ( 2 + i ) i x ) μ ) 2 sinh ( 2 ( ( 1 + 2 i ) i x ) μ ) 8 sinh ( ( ( 1 + i ) 2 i x ) μ ) 2 sinh ( ( ( 1 + 3 i ) 2 i x ) μ ) + 2 sinh ( ( ( 3 + 3 i ) 2 i x ) μ ) )
ψ ( x ) = μ 2 csc ( μ ) 8 ( 2 cosh ( μ ) + cosh ( 2 μ ) ( 3 cos ( μ ) + cos ( 3 μ ) cos ( ( 2 + i ) μ ) + cos ( ( 1 + 2 i ) μ ) + cosh ( ( 2 + i ) μ ) cosh ( ( 1 + 2 i ) μ ) ) ) × ( i cos ( ( 1 + i ) ( ( 1 + i ) x + 1 ) μ ) i cos ( ( 1 + i ) ( ( 1 + i ) x μ + μ ) ) 16 sin ( 2 μ ) + 4 sin ( 4 μ ) + ( 10 10 i ) sin ( ( 1 + i ) μ ) ( 2 + 2 i ) sin ( ( 3 + i ) μ ) + ( 4 4 i ) sin ( ( 2 + 2 i ) μ ) + ( 2 + 2 i ) sin ( ( 1 + 3 i ) μ ) ( 2 2 i ) sin ( ( 3 + 3 i ) μ ) 2 sin ( ( 1 + i ) ( ( 3 + i ) x ) μ ) ( 6 6 i ) sin ( ( 1 + i ) ( x 2 ) μ ) + ( 2 2 i ) sin ( ( 1 + i ) ( x + 2 ) μ ) + 6 sin ( ( 1 + i ) ( ( 1 i ) + x ) μ ) 2 i sin ( ( 1 + i ) ( ( 3 + i ) + x ) μ ) 4 i sin ( ( 1 + i ) ( 2 i + x ) μ ) + 4 sin ( ( 1 + i ) ( 2 i + x ) μ ) ( 4 2 i ) sin ( ( 1 + i ) ( ( 1 + i ) + x ) μ ) 2 sin ( ( 1 + i ) ( ( 1 + 3 i ) + x ) μ ) ( 2 + 2 i ) sin ( ( 1 + i ) ( ( 2 2 i ) + x ) μ ) + 2 i sin ( ( 1 + i ) ( ( 1 3 i ) + x ) μ ) 2 i sin ( ( ( 1 + i ) + 2 x ) μ ) + 2 i sin ( ( ( 1 + 3 i ) + 2 x ) μ ) + 6 sin ( ( 2 i + ( 1 i ) x ) μ ) ( 2 2 i ) sin ( ( 4 i + ( 1 i ) x ) μ ) 2 i sin ( ( ( 2 + 4 i ) + ( 1 i ) x ) μ ) + 2 sin ( 2 ( i x + 1 ) μ ) 4 sin ( ( 1 + i ) ( i x + 2 ) μ ) + 8 sin ( ( ( 1 i ) + 2 i x ) μ ) 2 sin ( ( ( 3 + i ) + 2 i x ) μ ) + 2 sin ( ( ( 1 3 i ) + 2 i x ) μ ) 2 sin ( ( ( 3 3 i ) + 2 i x ) μ ) + sin ( ( 1 + i ) ( ( 1 + i ) x + 1 ) μ ) ( 2 4 i ) sin ( ( ( 1 + i ) x + 2 ) μ ) 8 i sin ( ( ( 1 + i ) 2 x ) μ ) 2 i sin ( ( ( 3 + i ) 2 x ) μ ) + 2 i sin ( ( ( 3 + 3 i ) 2 x ) μ ) 6 i sin ( ( 2 ( 1 i ) x ) μ ) ( 2 2 i ) sin ( ( 4 ( 1 i ) x ) μ ) + 2 sin ( 2 ( 1 i x ) μ ) + 4 i sin ( ( 1 + i ) ( 2 i x ) μ ) 8 sin ( 2 ( ( 1 + i ) i x ) μ ) + 2 sin ( 2 ( ( 2 + i ) i x ) μ ) + 2 sin ( 2 ( ( 1 + 2 i ) i x ) μ ) + 2 sin ( ( ( 1 i ) 2 i x ) μ ) 2 sin ( ( ( 3 i ) 2 i x ) μ ) + 8 sin ( ( ( 1 + i ) 2 i x ) μ ) + 2 sin ( ( ( 1 + 3 i ) 2 i x ) μ ) 2 sin ( ( ( 3 + 3 i ) 2 i x ) μ ) + 6 i sin ( ( 2 ( 1 + i ) x ) μ ) ( 2 + 2 i ) sin ( ( 4 ( 1 + i ) x ) μ ) ( 2 + 4 i ) sin ( ( 1 i ) x μ + 2 μ ) + sin ( ( 1 + i ) ( ( 1 + i ) x μ + μ ) ) + 8 i sin ( ( 2 + 2 i ) μ 2 x μ ) 2 i sin ( ( 4 + 2 i ) μ 2 x μ ) 2 i sin ( ( 2 + 4 i ) μ 2 x μ ) 16 sinh ( 2 μ ) + 4 sinh ( 4 μ ) + ( 10 10 i ) sinh ( ( 1 + i ) μ ) ( 2 + 2 i ) sinh ( ( 3 + i ) μ ) + ( 4 4 i ) sinh ( ( 2 + 2 i ) μ ) + ( 2 + 2 i ) sinh ( ( 1 + 3 i ) μ ) ( 2 2 i ) sinh ( ( 3 + 3 i ) μ ) ( 6 6 i ) sinh ( ( 1 + i ) ( x 2 ) μ ) + ( 2 2 i ) sinh ( ( 1 + i ) ( x + 2 ) μ ) + 2 sinh ( ( 1 + i ) ( ( 3 i ) + x ) μ ) + 2 i sinh ( 2 ( ( 2 i ) + x ) μ ) 8 i sinh ( 2 ( ( 1 i ) + x ) μ ) 2 i sinh ( ( 1 + i ) ( ( 3 + i ) + x ) μ ) + 2 i sinh ( 2 ( ( 1 2 i ) + x ) μ ) + 2 i sinh ( ( 1 + i ) ( ( 1 3 i ) + x ) μ ) + 2 sinh ( 2 ( i x + 1 ) μ ) 2 sinh ( ( ( 3 + i ) + 2 i x ) μ ) + 2 sinh ( ( 1 + i ) ( ( 1 + i ) x + 1 ) μ ) ( 2 4 i ) sinh ( ( ( 1 + i ) x + 2 ) μ ) + 2 sinh ( 2 ( 1 i x ) μ ) 8 sinh ( 2 ( ( 1 + i ) i x ) μ ) + 2 sinh ( 2 ( ( 2 + i ) i x ) μ ) + 2 sinh ( 2 ( ( 1 + 2 i ) i x ) μ ) + 8 sinh ( ( ( 1 + i ) 2 i x ) μ ) + 2 sinh ( ( ( 1 + 3 i ) 2 i x ) μ ) 2 sinh ( ( ( 3 + 3 i ) 2 i x ) μ ) )
with μ 4 = λ . The adjoint of system of Equations (21)–(23) is
3 D 5 0 10 D 3 D 5 η ζ = λ 8 D 10 D 3 0 24 D η ζ
with
η = η = η = ζ = ζ = ζ = ζ = 0 at x = 0 ,
ζ = ζ + 2 η = ζ 4 η = 0 at x = 1 .
The solutions are
η ( x ) = cos ( 2 μ x ) + 2 cos ( ( 1 + i ) μ x ) + cosh ( 2 μ x ) + 2 cosh ( ( 1 + i ) μ x ) + ( 2 2 i ) ( sin ( μ ) sinh ( μ ) ) ( sin ( μ ) sinh ( 2 μ ) sin ( 2 μ ) sinh ( μ ) ) ( sin ( 2 μ ) + sinh ( 2 μ ) ) ( sin ( μ ) cosh ( μ ) cos ( μ ) sinh ( μ ) ) ( sin ( ( 1 + i ) μ x ) + sinh ( ( 1 + i ) μ x ) ) ( cos ( 2 μ ) + cosh ( 2 μ ) 2 ) sin ( 2 μ ) + sinh ( 2 μ ) ( sin ( 2 μ x ) + sinh ( 2 μ x ) ) 6
ζ ( x ) = 1 μ 2 ( cos ( 2 μ x ) 2 i cos ( ( 1 + i ) μ x ) + cosh ( 2 μ x ) + 2 i cosh ( ( 1 + i ) μ x ) + ( 2 + 2 i ) ( sin ( μ ) sinh ( μ ) ) ( sin ( μ ) sinh ( 2 μ ) sin ( 2 μ ) sinh ( μ ) ) ( sin ( 2 μ ) + sinh ( 2 μ ) ) ( sin ( μ ) cosh ( μ ) cos ( μ ) sinh ( μ ) ) ( sinh ( ( 1 + i ) μ x ) sin ( ( 1 + i ) μ x ) ) + ( cos ( 2 μ ) + cosh ( 2 μ ) 2 ) sin ( 2 μ ) + sinh ( 2 μ ) ( sin ( 2 μ x ) sinh ( 2 μ x ) ) .
The bi-orthogonal set for { ϕ } given by y = N T η is:
y ( x ) = 24 ( 1 + i ) μ ( sin ( 2 μ ) + sinh ( 2 μ ) ) ( sin ( μ ) cosh ( μ ) cos ( μ ) sinh ( μ ) ) × ( 2 2 i ) sin 2 ( μ ) sinh ( 2 μ ) cos ( μ x ) cosh ( μ x ) + cos ( μ ) sinh ( μ ) sin ( 2 μ ) ( sin ( ( 1 + i ) μ x ) ( 2 2 i ) sin ( 2 μ x ) ) + ( 2 2 i ) ( cos ( 2 μ ) 2 ) cos ( 2 μ x ) + sin ( 2 μ ) ( sinh ( ( 1 + i ) μ x ) ( 2 2 i ) sinh ( 2 μ x ) ) + sinh ( 2 μ ) ( ( 2 2 i ) sin ( 2 μ x ) sin ( ( 1 + i ) μ x ) ( 2 2 i ) sinh ( 2 μ x ) + sinh ( ( 1 + i ) μ x ) ) + ( 2 2 i ) sin 2 ( μ ) cos ( ( 1 + i ) μ x ) + ( 2 2 i ) ( cos ( 2 μ ) 2 ) cosh ( 2 μ x ) + ( 2 2 i ) cosh ( 2 μ ) ( cos ( 2 μ x ) + cosh ( 2 μ x ) ) + ( 2 2 i ) sin 2 ( μ ) cosh ( ( 1 + i ) μ x ) + sin ( μ ) cosh ( μ ) sin ( 2 μ ) ( sin ( ( 1 + i ) μ x ) ( 2 2 i ) sin ( 2 μ x ) ) + ( 2 + 2 i ) ( cos ( 2 μ ) 2 ) cos ( 2 μ x ) + sinh ( 2 μ ) ( ( 2 2 i ) sin ( 2 μ x ) + sin ( ( 1 + i ) μ x ) + ( 2 2 i ) sinh ( 2 μ x ) sinh ( ( 1 + i ) μ x ) ) sin ( 2 μ ) ( sinh ( ( 1 + i ) μ x ) ( 2 2 i ) sinh ( 2 μ x ) ) ( 2 2 i ) ( cos ( 2 μ ) 2 ) cosh ( 2 μ x ) ( 2 2 i ) cosh ( 2 μ ) ( cos ( 2 μ x ) + cosh ( 2 μ x ) ) + ( 4 4 i ) sinh 2 ( μ ) cos ( μ x ) cosh ( μ x ) + ( 2 + 2 i ) sin ( 2 μ ) sinh 2 ( μ ) cos ( μ x ) cosh ( μ x )
z ( x ) = 48 μ ( sin ( 2 μ x ) + ( 1 + i ) sin ( ( 1 + i ) μ x ) + sinh ( 2 μ x ) ( 1 i ) sinh ( ( 1 + i ) μ x ) + ( cos ( 2 μ ) + cosh ( 2 μ ) 2 ) sin ( 2 μ ) + sinh ( 2 μ ) ( cos ( 2 μ x ) cosh ( 2 μ x ) ) + 2 i ( sin ( μ ) sinh ( μ ) ) ( sin ( μ ) sinh ( 2 μ ) sin ( 2 μ ) sinh ( μ ) ) ( sin ( 2 μ ) + sinh ( 2 μ ) ) ( sin ( μ ) cosh ( μ ) cos ( μ ) sinh ( μ ) ) ( cosh ( ( 1 + i ) μ x ) cos ( ( 1 + i ) μ x ) ) .

4. Results for Direct Singular and Inverse SLPs

Here, we present some numerical examples of direct and inverse SLPs of orders 2 and 4.
Algorithm 1 was implemented using MATLAB (2014) [42]. The reference solutions of the EVPs are computed using Wolfram Mathematica 11 [41]. Direct problems were solved using the Magnus expansion method [39].
For the numerical calculations, n is the number of subdivisions in the interval [ a , b ] , m is the number of subdivisions in the interval [ λ 0 , λ * ] ; λ * being the maximum eigenvalue searching, L is the number of multisection steps used to calculate each eigenvalue in the characteristic function and M is the number of inverse algorithm steps.
Error values are stated using max-norm which is defined as | | x | | : = max n ( | x n | ) .

4.1. Singular SLPs

Example 1
(Singular SLP: infinite domain). This example extends the Magnus method, as described in [39] to the infinite domain. Consider the harmonic oscillator problem
u + ( λ x 2 ) u = 0 , x ( , )
together with Dirichlet and Neumann BCs, having exact eigenvalues λ k = 2 k + 1 , n = 0 , 1 , [43].
The domain is truncted to [ l , l ] , and using the parameters m = 10 , n = 100 , L = 10 the first three eigenvalues are computed (see Table 1). It can be deduced that the Magnus method has comparable performance to other methods. Moreover, fixed parameters can be changed to obtain better accuracy.
Example 2
(Singular SLP with asymmetric potential). Consider the Morse potential problem
u ( e α x 1 ) 2 u = λ u , x ( , )
together with Dirichlet and Neumann BCs. When λ ( 0 , 1 ) this has a finite number of EV: λ k = β k ( 2 β k ) , with β k = 1 2 α ( 2 k + 1 ) , k = 0 , 1 , , K , where K = 1 α 1 2 [43].
The domain is truncated to [ l , l ] , and using the parameters m = 10 , n = 100 , L = 10 and α = 0.02 , the first three eigenvalues are computed (see Table 2). Magnus method is more accurate when the truncated interval is small. Furthermore, it can be observed that initial eigenvalues are more accurate than the later ones.
Example 3
(Singular SLP in half axis). Consider the hydrogen atom equation
u 2 x 2 1 x u = λ u , x ( 0 , )
together with Dirichlet and Neumann BCs. This has a finite number of eigenvalues lying between 1 < λ < 0 , given by λ k = 1 4 ( k + 2 ) 2 , with k = 0 , 1 , , [44].
The domain is truncated to [ 0 , 100 ] , and using the parameters m = 100 , n = 200 , L = 10 the first few eigenvalues are computed (see Table 3).

4.2. Inverse SLPs

Example 4
(Inverse SLP with a different domain). This example extends Barcilon’s algorithm in [36] by changing the domain [ 0 , 1 ] to [ 0 , π ] :
u + ( λ cos ( x ) ) u = 0 , x ( 0 , π )
Here, p ( x ) is symmetric and normalized. In the inverse Algorithm 1 of [39], w n ( 0 ) ( x ) and ω n ( 0 ) ( x ) need to be changed into { w n ( 0 ) } 1 = cos ( n x ) cos ( n π ) 1 , { ω n ( 0 ) } 1 = { sin ( n x ) } 1 which are the basic solutions in the new domain. The truncated eigenvalue sequences at N = 6 , for Dirichlet boundary conditions (DDBCs) and Dirichlet– Neumann boundary conditions (DNBCs) are calculated using the Magnus method (with m = 20 , n = 20 , L = 15 ) and M = 10 inverse algorithm steps are used.
Figure 1 shows the reconstructed and exact potential, and the log absolute errors in the reconstructed potential, respectively. From Figure 1a, it is obvious that the potential is converging towards the exact one. The max-norm of difference between the reconstructed p and the actual p is ≈ 3.35 × 10 3 and the max-norm of difference between the reconstructed eigenvalues and the actual ones is ≈ 1.62 × 10 3 .
Figure 2 shows reconstructed p starting with perturbed eigenvalues
σ k δ = σ k + δ σ k · r a n d ( s i z e ( σ k ) , k = 1 , , 2 N
where δ is the noise level and N = 5 , m = 20 , n = 100 , L = 5 , M = 6 .
It is obvious that reconstruction is possible even in the presence of a significant noise.
Example 5
(Inverse SLP with a non-smooth potential). This example extends Barcilon’s algorithm in [36] by reconstructing a non-smooth potential:
u + ( λ p ( x ) ) u = 0 , x ( 0 , 2 )
with p ( x ) = | 1 x | 0.5 : symmetric and normalized, but non-smooth. The truncated eigenvalue sequences at N = 4 , for DDBCs and DNBCs are calculated using the Magnus method (with m = 10 , n = 100 , L = 5 ), and used in the reconstruction of the potential.
Figure 3 shows the reconstructed and exact potential, and the log absolute errors in the reconstructed potential, respectively, using one inverse algorithm step. From Figure 3a, it is obvious that the potential is converging towards the exact one. The max-norm of difference between the reconstructed p and the actual p is ≈ 1.13 × 10 1 and the max-norm of difference between the reconstructed eigenvalues and the actual ones is ≈ 3.15 × 10 1 . As anticipated, the error is maximum at the point of non-differentiability and at the boundaries.
Example 6
(FSLP). Consider
y ( 4 ) ( p 1 y ) + p 2 y = λ y , 0 < x < π
with p 1 0 and p 2 ( x ) = π / 2 x . Using the three spectra corresponding to the sets of BCs: { y ( 0 ) = y ( 0 ) = y ( 1 ) = y ( 1 ) = 0 } , { y ( 0 ) = y ( 0 ) = y ( 1 ) = y ( 1 ) = 0 } and { y ( 0 ) = y ( 0 ) = y ( 1 ) = y ( 1 ) = 0 } two potential functions p 1 , p 2 are reconstructed starting with zero initial guesses and N = 4 , m = 20 , n = 100 , L = 5 , M = 6 (Figure 4). The max-norm of difference between the reconstructed p 2 and actual p 2 is ≈ 3.5 × 10 1 and the max-norm of difference between the reconstructed eigenvalues and the actual ones is ≈ 2.1 × 10 1 .
Example 7
(FSLP). Consider
y ( 4 ) ( p 1 y ) + p 2 y = λ y , 0 < x < π
with p 1 ( x ) = | x π / 2 | π / 4 and p 2 0 . Using the three spectra corresponding to the sets of BCs: y ( 0 ) = y ( 0 ) = y ( 1 ) = y ( 1 ) = 0 , y ( 0 ) = y ( 0 ) = y ( 1 ) = y ( 1 ) = 0 and y ( 0 ) = y ( 0 ) = y ( 1 ) = y ( 1 ) = 0 two potential functions p 1 , p 2 are reconstructed starting with zero initial guesses and N = 4 , m = 20 , n = 100 , L = 15 , M = 6 (Figure 5). The max-norm of difference between the reconstructed p 1 and the actual p 1 is ≈ 1.4 × 10 1 and the max-norm of difference between the reconstructed eigenvalues and the actual ones is ≈ 7.5 × 10 1 .

5. Discussion

In the first few examples, the Magnus method [39] is further extended to solve some direct singular SLPs on infinite domain, half-axis and for an asymmetric potential. Later, Barcilon’s algorithm is extended to solve inverse SLPs in a non–unit domain and for a non–smooth potential. The first two examples, which extend the Magnus method in [39] to infinite interval, illustrates it’s applicability to singular SLPs effectively. Furthermore, the present method does not require transforming the singular SLP to a regular one or re-scaling the eigenvalues. It effectively finds the required eigenvalues, just by truncating the interval to a finite one.
Next, the inverse SLP is solved with smooth and non-smooth potential, even in the presence of noise. It is observed that the method is successful even in the presence of significant noise, provided that the assumptions of the algorithm are satisfied.
The last few examples solve the inverse FSLP using Barcilon’s algorithm with the initial knowledge of three spectra. A simplified FSLP is solved keeping one of the unknown potential functions zero.

6. Conclusions

More than 40 years after Barcilon’s paper [36], this paper gives a concrete implementation of the inverse SLP algorithm proposed therein, to our knowledge for the first time. Furthermore, computational feasibility and applicability of this algorithm for solving inverse SLPs of higher order is verified successfully in this paper (for n = 2 and n = 4 ).
In future work, it is possible to extend the results to solve even higher order SLPs.
In conclusion, this work provides a method that can be adapted successfully for solving a direct (regular/singular) or inverse SLP of an arbitrary order with arbitrary BCs.

Author Contributions

All authors contributed equally. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partially funded by UNIVERSITY OF KELANIYA, SRI LANKA.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Ambarzumian, V. Über eine frage der eigenwerttheorie. Z. Phys. A Hadron. Nucl. 1929, 53, 690–695. [Google Scholar]
  2. Tikhonov, A. On the uniqueness of the solution of the problem of electromagnetic sounding. Dokl. Akad. Nauk SSSR 1949, 69, 797–800. [Google Scholar]
  3. Levinson, N. The Inverse Sturm-Liouville Problem. Mat. Tidsskr. B 1949, 1949, 25–30. [Google Scholar]
  4. Marčenko, V. Sturm-Liouville Operators and Their Applications; Birkhauser Verlag: Basel, Switzerland, 1986. [Google Scholar]
  5. Levitan, B. Inverse Sturm-Liouville Problems; VNU Science Press: Utrecht, The Netherlands, 1987. [Google Scholar]
  6. Kravchenko, V.V. Direct and Inverse Sturm-Liouville Problems: A Method of Solution; Springer Nature: Berlin/Heidelberg, Germany, 2020. [Google Scholar]
  7. McLaughlin, J.R. Analytical methods for recovering coefficients in differential equations from spectral data. SIAM Rev. 1986, 28, 53–72. [Google Scholar] [CrossRef]
  8. Freiling, G.; Yurko, V.A. Inverse Sturm-Liouville Problems and Their Applications; NOVA Science Publishers: New York, NY, USA, 2001. [Google Scholar]
  9. Barcilon, V. Iterative solution of the inverse Sturm-Liouville problem. J. Math. Phys. 1974, 15, 429–436. [Google Scholar] [CrossRef]
  10. Sacks, P.E. An iterative method for the inverse Dirichlet problem. Inverse Probl. 1988, 4, 1055–1069. [Google Scholar] [CrossRef]
  11. Hald, O.H. The inverse Sturm-Liouville problem and the Rayleigh-Ritz method. Math. Comput. 1978, 32, 687–705. [Google Scholar] [CrossRef]
  12. Paine, J. A numerical method for the inverse Sturm-Liouville problem. SIAM J. Sci. Stat. Comput. 1984, 5, 149–156. [Google Scholar] [CrossRef]
  13. Rundell, W.; Sacks, P.E. Reconstruction techniques for classical inverse Sturm-Liouville problems. Math. Comput. 1992, 58, 161–183. [Google Scholar] [CrossRef]
  14. Lowe, B.D.; Pilant, M.; Rundell, W. The recovery of potentials from finite spectral data. SIAM J. Math. Anal. 1992, 23, 482–504. [Google Scholar] [CrossRef]
  15. Neher, M. Enclosing solutions of an inverse Sturm-Liouville problem with finite data. Computing 1994, 53, 379–395. [Google Scholar] [CrossRef] [Green Version]
  16. Fabiano, R.H.; Knobel, R.; Lowe, B.D. A finite-difference algorithm for an inverse Sturm-Liouville problem. IMA J. Numer. Anal. 1995, 15, 75–88. [Google Scholar] [CrossRef]
  17. Ghelardoni, P. Approximations of Sturm-Liouville eigenvalues using boundary value methods. Appl. Numer. Math. 1997, 23, 311–325. [Google Scholar] [CrossRef]
  18. Aceto, L.; Ghelardoni, P.; Magherini, C. BVMs for Sturm-Liouville eigenvalue estimates with general boundary conditions. J. Numer. Anal. Ind. Appl. Math 2009, 4, 113–127. [Google Scholar]
  19. Ghelardoni, P.; Magherini, C. BVMs for computing Sturm-Liouville symmetric potentials. Appl. Math. Comput. 2010, 217, 3032–3045. [Google Scholar] [CrossRef]
  20. Aceto, L.; Ghelardoni, P.; Magherini, C. Boundary value methods for the reconstruction of Sturm-Liouville potentials. Appl. Math. Comput. 2012, 219, 2960–2974. [Google Scholar] [CrossRef]
  21. Andrew, A.L. Numerov’s method for inverse Sturm-Liouville problems. Inverse Probl. 2005, 21, 223–238. [Google Scholar] [CrossRef]
  22. Andrew, A.L. Numerical solution of inverse Sturm-Liouville problems. ANZIAM J. 2004, 45, 326–337. [Google Scholar] [CrossRef]
  23. Andrew, A.L. Computing Sturm-Liouville potentials from two spectra. Inverse Probl. 2006, 22, 2069–2081. [Google Scholar] [CrossRef]
  24. Röhrl, N. A least-squares functional for solving inverse Sturm–Liouville problems. Inverse Probl. 2005, 21, 2009–2017. [Google Scholar] [CrossRef] [Green Version]
  25. Rafler, M.; Böckmann, C. Reconstruction method for inverse Sturm–Liouville problems with discontinuous potentials. Inverse Probl. 2007, 23, 933–946. [Google Scholar] [CrossRef]
  26. Efremova, L.; Freiling, G. Numerical solution of inverse spectral problems for Sturm-Liouville operators with discontinuous potentials. Open Math. 2013, 11, 2044–2051. [Google Scholar] [CrossRef]
  27. Ignatiev, M.; Yurko, V. Numerical Methods for Solving Inverse Sturm-Liouville Problems. Results Math. 2008, 52, 63–74. [Google Scholar] [CrossRef]
  28. Liu, C.S. Solving an inverse Sturm-Liouville problem by a Lie-group method. Bound. Value Probl. 2008, 2008, 749865. [Google Scholar] [CrossRef] [Green Version]
  29. Böckmann, C.; Kammanee, A. Broyden method for inverse non-symmetric Sturm-Liouville problems. BIT Numer. Math. 2011, 51, 513–528. [Google Scholar] [CrossRef]
  30. Böckmann, C.; Rattana, A. An inverse fourth order sturm-liouville problem. AIP Conf. Proc. 2015, 1648, 850030. [Google Scholar] [CrossRef]
  31. Gao, Q. Decent flow methods for inverse Sturm–Liouville problem. Appl. Math. Model. 2012, 36, 4452–4465. [Google Scholar] [CrossRef]
  32. Gao, Q.; Cheng, X.; Huang, Z. Modified Numerov’s method for inverse Sturm–Liouville problems. J. Comput. Appl. Math. 2013, 253, 181–199. [Google Scholar] [CrossRef]
  33. Drignei, M.C. A Newton-type method for solving an inverse Sturm-Liouville problem. Inverse Probl. Sci. Eng. 2014, 23, 851–883. [Google Scholar] [CrossRef]
  34. Kravchenko, V.V. On a method for solving the inverse Sturm-Liouville problem. J. Inverse Ill Posed Probl. 2018. [Google Scholar] [CrossRef]
  35. Neamaty, A.; Akbarpoor, S.; Yilmaz, E. Solving Symmetric Inverse Sturm-Liouville Problem Using Chebyshev Polynomials. Mediterr. J. Math. 2019, 16, 74. [Google Scholar] [CrossRef]
  36. Barcilon, V. On the solution of inverse eigenvalue problems of high orders. Geophys. J. Int. 1974, 39, 143–154. [Google Scholar] [CrossRef] [Green Version]
  37. McLaughlin, J.R. An Inverse Eigenvalue Problem of Order Four. SIAM J. Math. Anal. 1976, 7, 646–661. [Google Scholar] [CrossRef]
  38. McLaughlin, J.R. An Inverse Eigenvalue Problem of Order Four—An Infinite Case. SIAM J. Math. Anal. 1978, 9, 395–413. [Google Scholar] [CrossRef]
  39. Perera, U.; Böckmann, C. Solutions of Direct and Inverse Even-Order Sturm-Liouville Problems Using Magnus Expansion. Mathematics 2019, 7, 544. [Google Scholar] [CrossRef] [Green Version]
  40. Barcilon, V. On the uniqueness of inverse eigenvalue problems. Geophys. J. R. Astron. Soc. 1974, 38, 287–298. [Google Scholar] [CrossRef] [Green Version]
  41. Mathematica; Wolfram Research, Inc.: Champaign, IL, USA, 2016.
  42. Matlab; The MathWorks, Inc.: Natick, MA, USA, 2014.
  43. Taşeli, H. Accurate numerical bounds for the spectral points of singular Sturm–Liouville problems over − < x < . J. Comput. Appl. Math. 2000, 115, 535–546. [Google Scholar] [CrossRef]
  44. Marletta, M.; Pryce, J.D. Automatic solution of Sturm-Liouville problems using the Pruess method. J. Comput. Appl. Math. 1992, 39, 57–78. [Google Scholar] [CrossRef] [Green Version]
Figure 1. (a) Iteratively reconstructed potential (green dashed line) and exact potential (blue solid line) p ( x ) = cos ( x ) (b) Log absolute error in the reconstructed p.
Figure 1. (a) Iteratively reconstructed potential (green dashed line) and exact potential (blue solid line) p ( x ) = cos ( x ) (b) Log absolute error in the reconstructed p.
Mathematics 08 02074 g001
Figure 2. Reconstructed potential (blue) and exact potential (red) p ( x ) = cos ( x ) (a) δ = 0.1 (b) δ = 0.05 .
Figure 2. Reconstructed potential (blue) and exact potential (red) p ( x ) = cos ( x ) (a) δ = 0.1 (b) δ = 0.05 .
Mathematics 08 02074 g002
Figure 3. (a) Reconstructed potential (blue) and exact potential (red) p ( x ) = | 1 x | 0.5 (b) Log absolute error in the reconstructed p.
Figure 3. (a) Reconstructed potential (blue) and exact potential (red) p ( x ) = | 1 x | 0.5 (b) Log absolute error in the reconstructed p.
Mathematics 08 02074 g003
Figure 4. (a) Reconstructed potential (blue) and exact potential (red) p 2 ( x ) = π / 2 x (b) Log absolute error in the reconstructed p 2 .
Figure 4. (a) Reconstructed potential (blue) and exact potential (red) p 2 ( x ) = π / 2 x (b) Log absolute error in the reconstructed p 2 .
Mathematics 08 02074 g004
Figure 5. (a) Reconstructed potential (blue) and exact potential (red) p 1 ( x ) = | x π / 2 | π / 4 (b) Log absolute error in the reconstructed p 1 .
Figure 5. (a) Reconstructed potential (blue) and exact potential (red) p 1 ( x ) = | x π / 2 | π / 4 (b) Log absolute error in the reconstructed p 1 .
Mathematics 08 02074 g005
Table 1. Absolute errors for the first three eigenvalues of Equation (27) from Table 1 of Taşeli [43] and using Magnus method.
Table 1. Absolute errors for the first three eigenvalues of Equation (27) from Table 1 of Taşeli [43] and using Magnus method.
λ lDirichlet ErrorNeumann Error
[43]Magnus Method[43]Magnus Method
14.5 1.70 × 10 8 1.58 × 10 8 1.50 × 10 8 1.68 × 10 8
6.0 2.60 × 10 11 5.80 × 10 10 2.00 × 10 12 5.80 × 10 10
7.5 2.00 × 10 14 2.20 × 10 9 0.00 × 10 0 2.20 × 10 9
54.5 1.88 × 10 6 1.10 × 10 5 1.10 × 10 5 1.19 × 10 5
6.0 3.30 × 10 8 8.02 × 10 9 2.73 × 10 9 8.00 × 10 9
7.5 5.50 × 10 11 3.07 × 10 8 6.00 × 10 12 3.07 × 10 8
94.5 9.99 × 10 1 9.85 × 10 4 9.85 × 10 4 1.12 × 10 3
6.0 6.01 × 10 6 3.06 × 10 8 6.11 × 10 7 2.54 × 10 8
7.5 1.70 × 10 8 1.08 × 10 7 2.00 × 10 9 1.08 × 10 7
Table 2. Absolute errors for the first three eigenvalues of Equation (28) from Table 2 of Taşeli [43] and using Magnus method.
Table 2. Absolute errors for the first three eigenvalues of Equation (28) from Table 2 of Taşeli [43] and using Magnus method.
λ lDirichlet ErrorNeumann Error
[43]Magnus Method[43]Magnus Method
0.019940 1.10 × 10 8 5.60 × 10 12 1.90 × 10 9 1.12 × 10 11
50 6.07 × 10 10 3.17 × 10 11 6.00 × 10 11 3.17 × 10 11
60 1.15 × 10 13 9.39 × 10 11 2.48 × 10 14 9.39 × 10 11
0.059140 3.56 × 10 7 3.78 × 10 10 3.90 × 10 8 3.34 × 10 10
50 2.20 × 10 8 1.12 × 10 10 5.97 × 10 9 1.12 × 10 10
60 5.00 × 10 12 3.37 × 10 10 6.75 × 10 13 3.37 × 10 10
0.097540 7.69 × 10 7 1.85 × 10 8 1.13 × 10 7 1.93 × 10 8
50 6.65 × 10 7 3.83 × 10 10 1.29 × 10 7 3.82 × 10 10
60 2.50 × 10 11 1.15 × 10 9 3.18 × 10 12 1.15 × 10 9
Table 3. Eigenvalues and absolute errors for the first five eigenvalues of Equation (29) using Magnus method.
Table 3. Eigenvalues and absolute errors for the first five eigenvalues of Equation (29) using Magnus method.
kMagnus ValueMagnus Error
1 6.2323550470257 × 10 2 1.76 × 10 4
2 2.7717584125111 × 10 2 6.02 × 10 5
3 1.5581393052633 × 10 2 4.36 × 10 5
4 9.0660619438460 × 10 3 9.34 × 10 4
5 1.6285043066920 × 10 3 5.32 × 10 3
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Perera, U.; Böckmann, C. Solutions of Sturm-Liouville Problems. Mathematics 2020, 8, 2074. https://doi.org/10.3390/math8112074

AMA Style

Perera U, Böckmann C. Solutions of Sturm-Liouville Problems. Mathematics. 2020; 8(11):2074. https://doi.org/10.3390/math8112074

Chicago/Turabian Style

Perera, Upeksha, and Christine Böckmann. 2020. "Solutions of Sturm-Liouville Problems" Mathematics 8, no. 11: 2074. https://doi.org/10.3390/math8112074

APA Style

Perera, U., & Böckmann, C. (2020). Solutions of Sturm-Liouville Problems. Mathematics, 8(11), 2074. https://doi.org/10.3390/math8112074

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop