Next Article in Journal
A Note on a Class of Generalized Parabolic Marcinkiewicz Integrals along Surfaces of Revolution
Next Article in Special Issue
Variational Problems and Applications
Previous Article in Journal
A Complementary Dual of Single-Valued Neutrosophic Entropy with Application to MAGDM
Previous Article in Special Issue
Markov Moment Problem and Sandwich Conditions on Bounded Linear Operators in Terms of Quadratic Forms
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Variational Estimation Methods for Sturm–Liouville Problems

by
Elena Corina Cipu
1,2,* and
Cosmin Dănuţ Barbu
1
1
Department of Applied Mathematics, University Politehnica of Bucharest, Spl. Independentei, No. 313, RO-060042 Bucharest, Romania
2
CiTi, Faculty of Applied Sciences, University Politehnica of Bucharest, Spl. Independentei, No. 313, RO-060042 Bucharest, Romania
*
Author to whom correspondence should be addressed.
Mathematics 2022, 10(20), 3728; https://doi.org/10.3390/math10203728
Submission received: 31 July 2022 / Revised: 28 September 2022 / Accepted: 7 October 2022 / Published: 11 October 2022
(This article belongs to the Special Issue Variational Problems and Applications)

Abstract

:
In this paper, we are concerned with approach solutions for Sturm–Liouville problems (SLP) using variational problem (VP) formulation of regular SLP. The minimization problem (MP) is also set forth, and the connection between the solution of each formulation is then proved. Variational estimations (the variational equation associated through the Euler–Lagrange variational principle and Nehari’s method, shooting method and bisection method) and iterative variational methods (He’s method and HPM) for regular RSL are unitary presented in final part of the paper, which ends with applications.

1. Introduction

Nonlinearities are different from linear type by a function, an operator or a system that is nonlinear or is the case in which only some characteristics of it are known. The existence of the solution and the dependence of conditions for solving some classes of differential equations described by an operator is specified by the general framework of the Sturm–Liouville problem, with parametric conditions at the limit. The general framework of the Sturm–Liouville problem with parametric conditions at the limit is specified in the first part of the paper. The existence of the solution and the dependence of conditions is specified through the connection between the differential operator and Green’s function. Based on the properties of Green’s function, the operator used to analyze the behavior of the solution of the parameters given by the boundary conditions is specified. Variational problems derived from the initial RSLP are outlined with different type conditions in order to estimate the solution.
Let be the operator L = d d x p ( x ) d d x + ρ ( x ) as part of the regular Sturm–Liouville problem (RSL). The Sturm–Liouville (SL) problem expressed by the differential equation and the boundary conditions
a ( x ) d 2 u d x 2 + b ( x ) d u d x + c ( x ) u λ d ( x ) u = 0 ,
B 1 : a 1 u ( a ) + a 2 u ( a ) = 0 , | a 1 | + | a 2 | 0 , a 1 , a 2 R , B 2 : b 1 u ( b ) + b 2 u ( b ) = 0 , | b 1 | + | b 2 | 0 , b 1 , b 2 R
could be written as
L u + λ s ( x ) u = 0 , x ( a , b ) = I , λ R
with p ( x ) = a ( x ) , ρ ( x ) = c ( x ) , s ( x ) = d ( x ) in case b ( x ) = a ( x ) and with integrant factor μ = k e a x b ( t ) a ( t ) d t , p ( x ) = μ a ( x ) , ρ ( x ) = μ c ( x ) , s ( x ) = μ d ( x ) in case b ( x ) a ( x ) (see [1,2]).
The Sturm–Liouville equation is regular in the interval [ a , b ] if the functions verify the condition p ( x ) > 0 and s ( x ) > 0 , x I or s 0 and the operator L : H = L 2 ( I ) C 2 ( I ) L 2 ( I ) is self-adjoint with real eigenvalues and orthogonal eigenfunctions in space L s 2 ( I ) according to the inner product
f , g = a b f g d x f , g L 2 ( I ) ; f , g s = a b s f g d x in L s 2 ( I ) .
and for a given λ , there exist two linearly independent solutions of a R S L equation in the I interval, L 2 ( I ) = f : I R , a b | f ( x ) | 2 d x < .
We denote D ( L ) as the domain of L that is defined by
D ( L ) = y C ( [ a , b ] ) , y L 2 ( I ) , y satisfies B 1 , B 2 for general case D ( L ) = y C ( [ a , b ] ) , ( p y ) C ( [ a , b ] ) , y satisfies B 1 , B 2 for regular case .
The adjoint operator, L , associated to the operator L verifies L f , g = f , L g , f , g H and L is self-adjoint if L = L . Additionally, the operator L is symmetric if L f , g = f , L g , f , g D ( L ) . For the operator defined for SL problems, one obtains
g , L f f , L g = p f g f g | a b , f , g D ( L )
and the condition L f , g = f , L g holds if p f g f g a b = 0 verified in D ( L ) and the Lagrange identity is expressed by
g L f f L g = p f g f g , f , g D ( L ) .
Remark 1. 
L for SL problems is the self-adjoint operator if
p u v u v a b = p ( b ) u ( b ) v ( b ) p ( a ) u ( a ) v ( a ) = 0 .
For example, for p ( b ) = p ( a ) and periodic conditions
u ( a ) = u ( b ) = A , u ( a ) = u ( b ) = B
or antiperiodic conditions
u ( a ) = u ( b ) = A , u ( a ) = u ( b ) = B
the operator L is self-adjoint.
The RSL eigenvalue problem is to find v D ( L ) such that for L v + λ v = 0 with λ , the eigenvalue associated with v is the eigenfunction. For RSL problems, all the eigenvalues are real and positive [3,4,5], and there exists an infinite number of eigenvalues. The sequence of the eigenvalues ( λ n ) n is considered such that λ 0 < λ 1 < < λ n < with lim x λ n = . For each eigenvalue λ n , the corresponding eigenfunction v n is unique up to a constant factor and has exactly n 1 zeros in interval ( a , b ) . The set V = v n n , v n D ( L ) is complete in the D ( L ) space, and the solution of RSL is represented by a generalized Fourier series of the eigenfunction
u = n = 1 c n v n ( x ) .

2. General Framework of SLPs

In this section, we will mention certain conditions that the functions defining the operator L fulfill for different SL or RSL problems.
A second method to find the solution of a RSL problem, different from the generalized Fourier series development, is described using the Green function and two linear independent solutions. The section ends with the analysis of the Fourier equation with different types of boundary conditions.
Remark 2. 
For Sturm–Liouville problems, we consider two types of assumptions that are usually used (see [4])
  • General assumptions:
    (1) 
    p ( x ) C ( [ a , b ] ) , differentiable in x = a , p ( x ) 0 x ( a , b ] with p ( a ) = 0 , p ( a ) 0 .
    For p ( x ) = ( x a ) φ ( x ) , we suppose φ C ( [ a , b ] ) , φ ( x ) 0 , x [ a , b ]
    (2) 
    ρ , f C ( [ a , b ] ) .
  • RSL assumptions
    (1) 
    p ( x ) > 0 and s ( x ) > 0 or s 0 on [ a , b ] ;
    (2) 
    p , ρ , s C ( [ a , b ] ) ;
    (3) 
    a 1 , a 2 , b 1 , b 2 R ;
    (4) 
    p , ρ , s continuously differentiable on [ a , b ] .
Known equations, such as Fourier, Graetz–Nusslet, Collatz and Airy equations, for which RSL assumptions are verified, are given in Table 1, and the first eigenvalues and eigenfunctions are depicted in Figure 1 and Figure 2.
Some other known equations, such as the Legendre differential equation, Chebysev’s differential equation or Bessel equation, must be transformed into the Sturm–Liouville form that we considered in (3). These forms are specified in Table 2. Other SL equations for which general assumptions are fulfilled are exemplified in Table 3.
For an example of a singular SL, a discontinuity on the middle of the interval is considered [ a , b ] , x 0 = ( b a ) / 2 , with ρ ( x ) 0 and p ( x ) = 1 , x 0 , x 0 c 2 , x x 0 , 1 , c 0 , c 1 .
The problem to solve is L u + λ u = 0 , u ( 0 ) = u ( 1 ) = 0 and with u x 0 u x 0 + = 0 , u x 0 c 2 u x 0 + = 0 transmission conditions.
The asymptotic behavior of the eigenvalues leads to λ n n 2 2 π c 1 + c 2 , as n .

2.1. Resolvent Operator and Green Function

This RSL problem is solved using a Green function solution for the resolvent operator R ( λ ) = ( L + λ I ) 1 of the form
R ( λ ) f = φ λ ( x ) ω ( λ ) a x f ( t ) ψ λ ( t ) d t + ψ λ ( x ) ω ( λ ) x b f ( t ) φ λ ( t ) d t
where ψ λ , φ λ are non-trivial classical solutions of ( L + λ I ) f = 0 which satisfy
a 1 φ λ ( a ) + a 2 φ λ ( a ) = 0 b 1 ψ λ ( b ) + b 2 ψ λ ( b ) = 0
A simple normalization that eliminates some complexity can be specified by requiring
φ λ ( a ) = a 2 , φ λ ( a ) = a 1 ψ λ ( b ) = b 2 , ψ λ ( b ) = b 1
Then, the Wronskian ω ( λ ) = v 1 ( x ) v 2 ( x ) v 1 ( x ) v 2 ( x ) of these solutions is a function that depends only on λ :
ω ( λ ) = ( p φ λ ) ψ λ φ λ ( p ψ λ )
Therefore, ω ( λ ) 0 , x [ a , b ] or ω ( λ ) 0 . The Wronskian vanishes if φ λ , ψ λ is a dependent set of functions of x, which is precisely when both functions satisfy ( L + λ I ) h = 0 as well as the specified conditions at x = a and x = b , meaning λ is an eigenvalue of L (see [6]).
For the RSL case, the equation L u = 0 has two linear independent solutions, v 1 and v 2 such that a 1 v 1 ( a ) + a 2 v 1 ( a ) = 0 , b 1 v 2 ( b ) + b 2 v 2 ( b ) = 0 , the Green function G : [ a , b ] × [ a , b ] R ,
G ( x , y ) = v 1 ( y ) v 2 ( x ) / m , a y x b v 1 ( x ) v 2 ( y ) / m , a x y b ,
m = p ( x ) ω ( λ ) has the properties
(i)
G C 1 [ a , b ] 2 , G ( x , y ) = G ( y , x ) and satisfies the boundary conditions according to each variable;
(ii)
G C 2 [ a , b ] 2 M , with L ρ G p ( x ) G x x ( x , y ) p ( x ) G x ( x , y ) + ρ ( x ) G ( x , y ) = 0 over [ a , b ] 2 M and M = { ( x , y ) x = y } ;
(iii)
G x y + , y G x y , y = lim ε 0 ε > 0 G x ( y + ε , y ) G x ( y ε , y ) = 1 p ( x ) , G x is discontinuous on M.
Now let T be the operator T u ( x ) = a b G ( x , y ) u ( y ) d y defined on C [ a , b ] . Using the properties of the Green function G and the continuity of u, one obtains that Tu C 2 [ a , b ] and is the solution of the equation L u = f . The function Tu satisfies the same boundary conditions as u C 2 [ a , b ] , then T ( L u ) ( x ) = u ( x ) , and T is the inverse operator of L. The problem of eigenvalues and eigenfunctions L u + λ u = 0 , B 1 u ( a ) = 0 ; B 2 u ( b ) = 0 becomes T u = μ u , with μ = 1 / λ . For results for the construction of the operator T for fractional SLPs, see [7,8].
Remark 3 
(Rayleigh quotient). The eigenvalues of the operator L are lower bounded by a real constant. The smallest eigenvalue of the SL eigenvalue problem satisfies
λ 0 = min u 0 u D ( L ) L u , u u , u s = min u 0 u D ( L ) p u u a b + a b p u 2 + ρ u 2 d x a b u 2 s d x
and the minimum u 0 is achieved if u 0 is the eigenfunction corresponding to λ 0 .

2.2. Sturm–Liouville Fourier Problems

Consider the operator L = d d x p ( x ) d d x + ρ ( x ) and nonhomogeneous equation L u + λ s ( x ) u = f ( x ) with functions p ( x ) smooth, ρ ( x ) positive, and also (i) p ( a ) = 0 or p ( b ) = 0 or both or (ii) the interval I is infinite. In this section, we study the Fourier problem
u + α u = f ( x ) , x ( a , b ) , BVP conditions   ( 2 ) and   ( 3 )
with different type of boundary conditions. For example, in the case of Dirichlet conditions u ( a ) = u ( b ) = 0 , the operator L = d 2 d x 2 + α in space C 0 ( I ) is self-adjoint.
In Example 1, we study the case for α = 0 . Additionally, for the case α = n 2 > 0 , general solutions of the homogeneous equation are v n ( x ) = A exp ( n x ) + B exp ( n x ) , and for boundary conditions u ( a ) = u ( b ) = 0 , the RSL solution is v = 0 .
In Example 2, for α < 0 , we consider different cases, shown in Table 4 and Table 5.
Example 1. 
Let us consider the RSL equation u ( x ) = f x , with general solution v ( x ) = m x + n for the homogeneous equation and associated Green’s function
G ( x , y ) = x ( b y ) / b , 0 x y y ( b x ) / b , y x < b .
Using the superposition principle, the solution of the problem defined for u ( 0 ) = A , u ( b ) = B is u ( x ) = v 1 ( x ) + v 2 ( x ) , v 1 ( x ) = ( b x ) A + x B and
v 2 ( x ) = 0 b G ( x , y ) f ( y ) d y = 1 b ( b x ) 0 x y f ( y ) d y x 1 y ( b y ) f ( y ) d y .
Changing the boundary conditions in the previous problem, we now consider
u ( x ) = f x , B 1 : u ( 0 ) u ( 0 ) = 0 , B 2 : u ( b ) + u ( b ) = 0 .
Solving the initial value problem u ( x ) = f x , B 1 : u ( 0 ) = A , u ( 0 ) = A , one finds the solution u ( x ) = A ( 1 + x ) 0 b ( x y ) f ( y ) d y and boundary condition B 2 leads to u ( x ) = 0 b G ( x , y ) f ( y ) d y with G ( x , y ) = ( 1 + x ) ( b + 1 y ) / ( b + 2 ) , x < y ( 1 + y ) ( b + 1 x ) / ( b + 2 ) , y < x .
Example 2. 
For α = n 2 and a = 0 , b = π , general solutions of the equation are v n ( x ) = A cos ( n x ) + B sin ( n x ) with λ n = n 2 the eigenvalues and for u ( 0 ) = u ( π ) = 0 the eigenfunctions are v n = sin ( n x ) . The general solution is a Fourier series:
u ( x ) = n = 1 B n sin ( n x ) , B n = f ( x ) , sin ( n x ) sin ( n x ) , sin ( n x ) = 2 π 0 π f ( x ) sin ( n x ) d x
  • Case 1.1
According to Table 4, for Cases 1.1, we consider a 2 = 0 and B 1 is u ( 0 ) = 0 and the eigenfunctions v n = sin ( λ n x ) . The eigenvalues corresponding to Cases 1.1a and 1.1b are λ n = n π b 2 and λ n = ( 2 n + 1 ) π 2 b 2 , respectively. For problems P1c
P 1 c : u ( x ) = λ u + f ( x ) in ( 0 , b ) ; u ( 0 ) = 0 , b 1 u ( b ) + b 2 u ( b ) = 0
the general solution is u ( x ) = n = 1 c n v n ( x ) with the eigenvalues determined by the equation tan ( λ n b ) = b 2 b 1 λ n . The determination of the first eigenvalues is graphically presented in Figure 3.
  • Case 1.2
The eigenfunctions corresponding to Case 1.2 are v n = cos ( λ n x ) with the eigenvalues determined by the equation tan ( λ n b ) = b 1 b 2 λ n .
  • Case 2
The eigenfunctions corresponding to Case 2.1 are v n = λ n cos ( λ n x ) + sin ( λ n x ) with the eigenvalues determined by the equation tan ( λ n b ) = ( b 1 + b 2 ) λ n b 1 b 2 λ n . If b has the form ( 2 n + 1 ) π 2 · b 2 b 1 , then λ n = b 1 b 2 is the eigenvalue for the problem. In Figure 4a, the determination of the first eigenvalues is graphically presented as the roots of the function tan ( x ) + ( b 1 + b 2 ) x b 1 b b 2 x with notation x = λ n b and in Figure 4b, the corresponding eigenfunctions are plotted.
For Case 2.2, the eigenfunctions are v n = λ n a 1 cos ( λ n x ) + sin ( λ n ) , and the eigenvalues are the solutions of the nonlinear equation tan ( λ n b ) = ( a 1 b 1 ) λ n a 1 + λ n .
Example 3. 
The conditions can be considerably weakened with respect to continuity and differentiability. In some cases, changes of variables, dependent and independent, may transform a problem from singular to regular; see [4].
For construction of the solutions, the Dirac function is used. Green’s function verifies d d x p ( x ) d G ( x , y ) d x + ρ ( x ) G ( x , y ) = δ ( x y ) , and expresses the response under homogeneous boundary conditions to a forcing function consisting of a concentrated unit of inhomogeneity at x = y .
For the problem u ( x ) = λ u ( x ) + f ( x ) in ( 0 , b ) ; u ( 0 ) = u ( b ) = 0 , λ = n 2 , the solution is u ( x ) = 0 b G ( x , y ) f ( y ) d y using the Green function
G ( x , y ) = sin ( n x ) sin ( n ( b y ) ) n sin ( n b ) , 0 x < y sin ( n ( b x ) ) sin ( n y ) n sin ( n b ) , y < x b .
That leads to the representation
u ( x ) = s i n ( n x ) n sin ( n b ) 0 x sin ( n ( b y ) ) f ( y ) d y sin ( n ( b x ) ) n sin ( n b ) x b sin ( n y ) f ( y ) d y .
Remark 4. 
Using the definition of the norm convergence, namely: “A sequence ( φ n ) n in L s 2 ( I ) converges to φ L s 2 ( I ) if lim n φ n φ s = 0 , i.e., φ n φ in L s 2 norm”, some sequences δ n ( x ) could be used instead of δ ( x ) in order to obtain the Green function.
Starting from the definition δ ( x y ) = 0 , x y , x = y and use some properties (see [6,9]):
δ is symmetric with δ ( a x ) = 1 | a | δ ( x ) , δ ( x ) = lim n δ n ( x ) , δ n ( x ) = n π e n 2 x 2 , δ n ( x ) = sin 2 ( n x ) n π x 2 or δ n ( x ) = n π 1 + n 2 x 2 , also
δ x 2 a 2 = 1 | 2 a | [ δ ( x + a ) + δ ( x a ) ] a n d + f ( y ) δ ( x y ) d x = f ( x ) .

3. Variational RSL Problems

We define problem P 1 as follows:
p u + ρ u + λ s u = f , on ( a , b )
B 1 u ( a ) : a 1 u ( a ) + a 2 u ( a ) = 0 B 2 u ( b ) : b 1 u ( b ) + b 2 u ( b ) = 0
and the set V = v C 1 ( [ a , b ] ) , v piecewise continuous on [ a , b ] , B 1 v ( a ) , B 2 v ( b ) verified . For v D ( L ) , we have a b p u v d x + a b ( ρ + λ s ) u v d x = a b f v d x , v V (see [10,11]).
Variational problem (VP 1 ) associated to the problem P 1 is as follows: find u V such that a ( u , v ) = l v , v V with
a ( u , v ) = a b p u v + ( ρ + λ s ) u v d x , u , v V ; l v = a b f v d x , v V .
The functional F : V R , F v = 1 2 a ( u , v ) l v , v V expresses the difference between the internal elastic energy and the load potential.
Lemma 1. 
(i) 
l : V R is linear;
(ii) 
Let a : V × V R with λ positive eigenvalue of RSL ( P 1 ), then a ( u , v ) is bilinear functional, positive and symmetric.
Proof. 
(i) Let v 1 , v 2 V and α , β R , then l ( α v 1 + β v 2 ) = α l v 1 + β l v 2 is the result that is obtained from the properties of the scalar product.
(ii) Let u 1 , u 2 , v V and α , β R then
a α u 1 + β u 2 , v = a b p α u 1 + β u 2 v + ( ρ + λ s ) α u 1 + β u 2 v d x = α a b p u 1 v + ( ρ + λ s ) u 1 v d x + β a b p u 2 v + ( ρ + λ s ) u 2 v d x = α a u 1 , v + β a u 2 , v
Let u , v 1 , v 2 V and α , β R , then
a u , α v 1 + β v 2 = a b p u α v 1 + β v 2 + ( ρ λ s ) u α v 1 + β v 2 d x = α a b p u v 1 + ( ρ + λ s ) u v 1 d x + β a b p u v 2 + ( ρ + λ s ) u v 2 d x = α a u , v 1 + β a u , v 2
Let u V , then
a ( u , u ) = a b p u 2 + ( ρ + λ s ) u 2 d x = p a b u 2 d x + a b ( ρ + λ s ) u 2 d x
The weight function p ( x ) in [ a , b ] is positive and in RSL( P 1 ) conditions, ρ + λ s is a positive function in [ a , b ] accordingly a ( u , u ) 0 , u V , and hence a ( · , · ) is positive. Let u , v V , then
a ( u , v ) = a b p u v + ( ρ + λ s ) u v d x = a ( v , u )
Consequently, a ( · , · ) is symmetric. □
Minimization problem ( M P 1 ) associated to ( V P 1 ) is as follows.
Find u V such that F u = min v V F v with
F v = 1 2 a b p ( v ) 2 + ( ρ + λ s ) v 2 d x a b f v d x
Theorem 1. 
(1) 
u V is the solution of ( V P 1 ) if u is solution of ( M P 1 ) ;
(2) 
u H = C 1 ( [ a , b ] ) C 2 ( ( a , b ) ) , u solution of ( V P 1 ) , then u solution of ( P 1 ) .
Proof. 
(1) (i) Let u V be the solution of ( V P 1 ) , then a ( u , v ) = l v , v V . For any w V , denoting v = w u V , we have
F w F u = 1 2 a ( u + v , u + v ) l ( u + v ) F u = F u + F v + a ( u , v ) F u = a ( u , v ) + F ( v ) = a ( u , v ) l v + l u + F ( v ) = 1 2 a ( v , v ) 0
meaning that
F u = min w V F w i . e . u solution of ( M P 1 )
(1) (ii) Let u V be the solution of ( M P 1 ) , then F u = min w V F w ; therefore, F w F u 0 , w V .
Using w = u + t v , t R , one finds F ( u + t v ) F u 0 , t R meaning
1 2 a ( u + t v , u + t v ) l ( u + t v ) 1 2 a ( u , u ) + l u 0 , t R , u , v V
1 2 a ( u , u ) + t a ( u , v ) + t 2 2 a ( v , v ) l u t l v 1 2 a ( u , u ) + l u 0 , t R , u , v V
1 2 a ( v , v ) t 2 + a ( u , v ) l v t 0 t R
Using the positivity of the term a ( v , v ) , one finds
[ a ( u , v ) l v ] 2 0 a ( u , v ) = l v v V
meaning the u solution of ( V P 1 ) .
(2) Let u V H solution of ( V P 1 ) , then
a b p u v + ( ρ + λ s ) u v d x = a b f v d x a b p u v d x + a b ( ρ + λ s ) u v d x = a b f v d x p u v a b a b p u v d x + a b ( ρ + λ s ) u v d x = a b f v d x p u v a b + a b p u + ( ρ + λ s ) u f v d x = 0 , v V
In case of a self-adjoint operator for L, such as in the case of periodic or antiperiodic boundary conditions, we have
p u v a b = p ( b ) u ( b ) v ( b ) p ( a ) u ( a ) v ( a ) = 0 , u , v V
and one obtains a b φ v d x = 0 , v V for φ ( x ) = ( p u ) + ( ρ + λ s ) u C ( a , b ) that is φ 0 over interval ( a , b ) .
This means that the u solution of P 1 also verifies the boundary conditions. □
The theorem proved above transfers the search space of the solution u of the problem P 1 to the search space for the solution of the problem M P 1 , where the existence is ensured through the Lax–Milgram theorem for a coercive quadratic form, even more general from the Lions–Stampacchia theorem, where a ( · , · ) is of a symmetric positive bilinear form.

4. Variational Approaches for VP of RSL

4.1. Nehari Variational Method

Let u C 2 ( [ a , b ] ) and F : [ a , b ] × C 2 ( [ a , b ] ) × C ( [ a , b ] ) [ a , b ] be a function that has continuous second-order derivatives with respect to all of its arguments. According to the Euler–Lagrange variational principle, a necessary condition for the functional
J ( u ) = a b F x , u , u d x ,
to be stationary at u is that u is a solution of the Euler–Lagrange equation (see [3])
F u d d x F u = 0 , a x b
with the Dirichlet conditions u ( a ) = A and u ( b ) = B .
For a nonlinear RSL problem,
u = f ( x , u 2 ) , f H , x ( a , b ) = I ; u ( a ) = A and u ( b ) = B ,
for H = f C ( [ a , b ] × [ 0 , ] ) | f verifies ( ip 1 ) , ( ip 2 ) and boundary conditions .
( ip 1 ) : f ( x , y ) > 0 for y > 0 ( ip 2 ) : ν , y ν f ( x , y ) is a non - decreasing function of y [ 0 , )
For the problem (22), looking for the extremum value of the functional
J ( u ) = a b u 2 0 u 2 ( x ) f ( x , y ) d y d x ,
for the set V a , b = u | u C ( [ a , b ] ) , u piecewise continuous in [ a , b ] , u ( a ) = u ( b ) = 0 the functional J ( u ) is not bounded. Using the Nehari method, a new condition on the function f and u is imposed:
a b u 2 d x = a b u 2 f ( x , u 2 ) d x ,
which is satisfied by the solutions of (22).
Let us consider the set V a , b = u | u V a , b verifies ( 24 ) . For I given, let μ ( a , b ) = inf u V J ( u ) . Then u V a , b such that J ( u ) = μ ( a , b ) and also (see [12]), for u V a , b with J ( u ) = μ ( a , b ) , w = | u | V a , b is a positive solution of (22).
The function μ ( a , b ) is continuous with respect to both arguments and
μ ( a , b ) = inf c < d [ a , b ] μ ( c , d ) with lim b a μ ( a , b ) = .
Remark 5. 
For a partition Δ n : a = x 0 < x 1 < < x n 1 < x n = b of the interval I, over each subinterval [ x i ; x i + 1 ] , consider u i V x i , x i + 1 normalized with the Nehari condition and
F o r x [ x i ; x i + 1 ] : u ( x ) = ( 1 ) i | u i ( x ) | , J ( u ) = μ n 1 ( x 1 , x 2 , , x n 1 )
F o r Δ n given : μ n 1 ( x 1 , x 2 , , x n 1 ) = i = 1 n μ ( x i 1 , x i ) ,
then the solution u ( x ) is in V a , b and is vanishing n 1 times over interval I. Additionally, if | u k ( x k ) | | u k + 1 ( x k ) | then μ n 1 ( x 1 , x 2 , , x n 1 ) is not a minimum of i = 1 n μ ( x i 1 , x i ) .

4.2. Variational Estimations for RSL

In the following, two variational estimation methods are presented, the shooting method and bisection method, consisting in solving the variational equations associated to the problem given.
  • Shooting method:
( P 2 ) p u + ( q + λ s ) u = 0 , x [ a , b ] u ( a ) = 0 , u ( b ) = 0
For λ eigenvalue and u λ ( x ) , the corresponding eigenfunction u λ ( a ) 0 and y = u λ u λ ( a ) is the normalized eigenfunction with y ( a ) = 1 , which is the solution for the variational equation associated to (28) with the initial value conditions:
( V I P ) p y + ( q + λ s ) y = 0 , x [ a , b ] y ( a ) = 0 , y ( a ) = 1
with y ( b ) = 0 (see [4]).
Algorithm of the shooting method:
Step 1
Determine an interval of an eigenvalue and make a guess;
Step 2
Solve VIP ( P 2 ) and find the eigenfunction u = u λ ( x ) ;
Step 3
If u λ ( b ) = 0 or u λ ( b ) < ε given, then Stop.
Else, find λ the root of u λ ( b ) = 0 in a given interval and update λ .
GO TO Step 1.
  • Bisection method
For SL eigenvalue problem ( P 3 ) with functions and constants satisfying RSL assumptions (1)–(4):
( P 3 ) p u + ( q + λ s ) u = 0 , x [ a , b ] B 1 u ( a ) : a 1 u ( a ) + a 2 u ( a ) = 0 B 2 u ( b ) : b 1 u ( b ) + b 2 u ( b ) = 0
the related variational initial value problem V I P 3 is
( V I P 3 ) p u + ( q + λ s ) u = 0 , , x [ a , b ] u ( a ) = a 2 a 1 2 + a 2 2 ; u ( a ) = a 1 a 1 2 + a 2 2
For the eigenvalue λ denoting u λ , the corresponding eigenfunction has u λ ( 0 ) 0 , and y = u λ is the normalized u λ eigenfunction such that a 1 y ( a ) + a 2 y ( a ) = 0 . In this case, λ is the eigenvalue for P 3 if F ( λ ) : = B 2 y ( b ) = b 1 y ( b ) + b 2 y ( b ) = 0 .
The function w λ ( x ) = u λ λ satisfies the variational initial value problem ( V I V P 3 ) .
( V I V P 3 ) p w + ( q + λ s ) w + s w = 0 w ( a ) = 0 , w ( a ) = 0
and F ( λ ) = b 1 u λ ( b ) + b 2 u λ ( b ) is a continuously differentiable function on λ with F ( λ ) = b 1 w λ ( b ) + b 2 w λ ( b ) 0 .
Remark 6. 
Under RSL assumptions (1)–(4), if λ is the eigenvalue of ( P 3 ) and y = u λ is the corresponding normalized eigenfunction, then there exists ( λ inf , λ sup ) containing λ such that F ( λ inf ) F ( λ sup ) < 0 and the approximate sequence ( λ n ) n is convergent, λ n λ and y n = u λ n are the corresponding eigenfunctions obtained by solving ( V I V P 3 ) such that y n y and y n y .
For instance, the solution of the problem
( E P 3 ) p ( x ) u + ρ u ( x ) = f ( x ) , x ( 0 , 1 ) u ( 0 ) = 0 , u ( 1 ) = 0
with p , ρ , f C ( [ 0 , 1 ] ) verifying (RSL) conditions (1)–(4) is obtained solving the associated
( V I E P 3 ) p y + ρ y = f ( x ) , x ( 0 , 1 ) y ( 0 ) = 0 , y ( 1 ) = s
The solution of ( E P 3 ) is determined such that u s ( x ) = u p ( x ) + s v ( x ) with u s ( 1 ) = 0 and u p is a particular solution of L y = f and v satisfies
L v = 0 , v ( 0 ) = 0 ; v ( 0 ) = 1 .

4.3. Iterative Variational Methods for RSL

Among analytical estimation methods, the variational iteration method (VIM or He’s methods, see [13]) and homotopy perturbation method (HPM) (see [14,15]) are considered to find approximations for the nonlinear equation
L u + λ s ( x ) u = f ( x , u , u ) , x ( a , b ) = I , λ R ,
under different boundary conditions (Dirichlet, Neumann or general case B 1 u ( a ) , B 2 u ( b ) ) .

4.3.1. He’s Variational Method (VIM)

For the nonlinear Equation (35), we define N as the nonlinear operator such that (35) becomes
L u + N u = g ( x ) , x ( a , b ) = I ,
and the correction functional for the general Lagrange multiplier method is
u n + 1 ( x ) = u n ( x ) + 0 x μ t , x , λ L u n ( t ) + N u ˜ n ( t ) g ( t ) d t ,
with u ˜ n considered as restricted variation, δ u ˜ n = 0 , and μ t , x , λ a Lagrange multiplier determined through the calculus of variations from (37); see [10,16].
δ u n + 1 ( x ) = δ u n ( x ) + δ 0 x μ t , x , λ L u n ( t ) + N u ˜ n ( t ) g ( t ) d t .

4.3.2. Homotopy Perturbation Method (HPM)

For the nonlinear Equation (35) we define the operators L and N for q [ 0 , 1 ]
L [ Φ ( x , q ) ] = d d x p ( x ) d d x Φ ( x , q ) ,
N [ Φ ( x , q ) ] = d d x p ( x ) d d x Φ ( x , q ) + ( ρ + λ s ) ( x ) Φ ( x , q ) f ( x , Φ ( x , q ) , Φ x ( x , q ) ) ,
given by the maximum order of derivation from the equation and by the form of the equation; see [14,17]. We write the zero-order equation associated with the initial equation:
( 1 q ) L Φ ( x , q ) u 0 ( x ) = h q N [ Φ ( x , q ) ]
with h as a nonzero parameter, u 0 as a first analytical approximation of the function u with conditions
Φ ( x , 0 ) = u 0 ( x ) ; Φ ( x , 1 ) = u ( x ) , x [ a , b ] .
where u 0 is a initial function that verifies the boundary conditions B 1 u ( a ) : a 1 u ( a ) + a 2 u ( a ) = 0 , B 2 u ( b ) : b 1 u ( b ) + b 2 u ( b ) = 0 could be obtained from polynomial approximation developing the function f.
We develop Φ ( x , q ) by a Taylor series in the vicinity of the origin in relation to the second variable
Φ ( x , q ) = u 0 ( x ) + 1 u m ( x ) q m ; u m ( x ) = 1 m ! m Φ x m ( x , q ) q = 0
A good choice for h (in relation to the error obtained compared to the initial equation) leads to u ( x ) = u 0 ( x ) + m = 1 u m ( x ) .
The equation of order m:
  • case m = 1 L u m ( x ) = h N [ u m 1 ]
  • case m 2 L f m ( x ) f m 1 ( x ) = h N [ u m 1 ]
with boundary conditions B 1 u m ( a ) , B 2 u m ( b ) .
For the approximation of order 1, we have
p ( x ) u 1 ( x ) = h p ( x ) u 0 ( x ) + ( ρ ( x ) + λ s ( x ) ) u 0 ( x ) f ( x , u 0 ( x ) , u 0 ( x ) ) ε 0 ( x )
with ε 1 ( x , h ) = N [ u 1 ( x ) ] .
The parameter h at step 1 is chosen such that the value of max x I ε 1 ( x , h ) is the smallest possible and becomes the next value of ε 1 , but also could be taken as h = 1 . Iteratively, for m 2
L u m + λ s ( x ) u m = f ( x , u m 1 , u m 1 ) , x I ; B 1 u m ( a ) , B 2 u m ( b )
with start condition u 0 being known and stop condition max x u m u m 1 < ϵ .

4.3.3. Applications

Let us consider a nonlinear RSL, such as the following problem:
u ( x ) + λ u ( x ) = f x , u , u , x [ 0 , b ] ; B 1 u ( 0 ) , B 2 u ( b ) .
The corresponding correction functional (37) to the problem (46) for the variational iteration method leads to the general Lagrange multiplier
for λ > 0 , μ t , x , λ = 1 2 α e α ( t x ) e α ( x t ) = 1 α sinh α ( t x ) , α = λ
for λ < 0 , μ t , x , λ = 1 α sin α ( t x ) , α = λ
and for f x , u , u = g ( x ) , one finds
u n + 1 ( x ) = u n ( x ) + 1 α 0 x sinh α ( t x ) u n ( t ) + α 2 u n ( t ) g ( t ) d t , λ > 0 ,
u n + 1 ( x ) = u n ( x ) 1 α 0 x sin α ( t x ) u n ( t ) + α 2 u n ( t ) g ( t ) d t , λ < 0
with u 0 ( x ) = A + B x . Particularly, for g ( x ) = x , the first step leads to
u 1 ( x ) = u 0 ( x ) + 1 α 0 x sinh α ( t x ) λ A + ( λ B 1 ) t d t , λ = α 2 ,
u 1 ( x ) = u 0 ( x ) 1 α 0 x sin α ( t x ) λ A + ( λ B 1 ) t d t , λ = α 2
from where
u 1 ( x ) = 2 A + 2 B x + α 2 B 1 α 3 sinh ( α x ) A cosh ( α x ) x α 2 , λ = α 2 ,
u 1 ( x ) = 2 A + α 2 B 1 α 3 sin ( α x ) A cos ( α x ) + x α 2 , λ = α 2 ,
For n > 1 , we have
u n + 1 ( x ) = u n ( x ) + 1 α 0 x sinh α ( t x ) u n ( t ) + α 2 u n ( t ) t d t , λ > 0 ,
u n + 1 ( x ) = u n ( x ) 1 α 0 x sin α ( t x ) u n ( t ) + α 2 u n ( t ) t d t , λ < 0
Constants A and B are determined from the boundary condition imposed to the last function u n computed, and for u ( x ) = u n ( x ) , boundary conditions are imposed, resulting in a system for the constants A , B . Thus, the solution of the problem is obtained.
In the case of using HPM, Equation (44), with λ = α 2 , for the first step becomes
u 1 ( x ) = h u 0 ( x ) + λ u 0 ( x ) x = h ε 0 ( x )
and for case m 2
u m u m 1 ( x ) = h u m 1 ( x ) + λ u m 1 ( x ) x = ε m ( x , h ) .
One obtains
u 1 ( x ) = h λ A x 2 2 + ( λ B 1 ) x 3 3 ,
u 2 = h ( 1 h ) ( λ A + ( λ B 1 ) x ) + h x h 2 λ 2 A x 2 2 + λ ( λ B 1 ) x 3 3 = ε 1 ( x , h ) .
from where
u 2 = h ( 1 h ) ( λ A x 2 2 + ( λ B 1 ) x 3 6 ) + h x 3 6 h 2 λ 2 A x 4 24 + λ ( λ B 1 ) x 5 603 .
The solution will be u ( x ) = u 0 ( x ) + u 1 ( x ) + u 2 ( x ) + .
Additionally, if we consider the nonlinearity in (46) through function f x , u , u = g ( x ) u , that is, for g ( x ) = x 2 , x ( l , l ) we are in the harmonic oscillator case ( λ < 0 ), then the following correction functional appears
u n + 1 ( x ) = u n ( x ) + 1 α 0 x sinh α ( t x ) u n ( t ) + ( α 2 g ( t ) ) u n ( t ) d t , λ > 0 ,
u n + 1 ( x ) = u n ( x ) 1 α 0 x sin α ( t x ) u n ( t ) + ( α 2 g ( t ) ) u n ( t ) d t , λ < 0
as an eigenfunction of the equation Hermite polynomials appears.
For B 1 u ( l ) : u ( l ) = 0 ; B 2 u ( l ) : u ( l ) = 0 and ( λ < 0 )
u 1 ( x ) = A + B x 1 α 0 x sin α ( t x ) t 2 ( A + B t ) + α 2 ( A + B t ) d t ,
u n + 1 ( x ) = u n ( x ) 1 α 0 x sin α ( t x ) u n ( t ) + ( α 2 + t 2 ) u n ( t ) d t ,
u 1 ( x ) = h u 0 ( x ) + ( x 2 + α 2 ) u 0 ( x ) = h ε 0 ( x )
and for case m 2
u m u m 1 ( x ) = h u m 1 ( x ) + ( x 2 + α 2 ) u m 1 ( x ) = ε m ( x , h ) .
The two methods are fast convergent methods.
Variational iteration methods, such as VIM and HPM, could be used also for nonlinear propagation problems in which the temporal variable is considered, for example, for the coupled pseudo-parabolic equation, or the one-dimensional coupled Burgers equation numerically studied in [18]. The nonlinear coupled Burgers equations are also studied in [19] as an application of EOHAM (extension optimal homotopy asymptotic method) in which homotopy is combined with perturbation techniques. The Newell–Whitehead–Segel equation (NWSE) was also studied using the VIM technique and He’s polynomials [20].

5. Conclusions

In the first part of the paper, definitions and results are presented connected to regular and singular Sturm–Liouville problems. Some types of direct singular SLPs were solved in [5,21] and a study of the inverse SLP algorithm was made. We defined in a different manner the SLP, and different boundary conditions were considered. All the figures were made using Matlab codes, the academic versions.
In the core of the paper, the variational formulation (VP) through a bilinear functional positive and symmetric is made. The minimization problem (MP) is also outlined through the functional of energy, and the equivalence of the formulations under some conditions imposed for RSL problems is proved.
Variational estimations are in the final part of the paper through the construction of the solution trough variational equations associated to the problem, such as the shooting method and bisection method, or using a sequential analytical approximate solution that is constructed according to the accuracy established. Here, we present He’s variational method and the homotopy method. In the closing part, a is taken into account and the sequentiality of the transition from one step to another is specified for both methods. Al-Khaled et al. (see [22]) solve numerically a SLP using the general Sinc–Galerkin and Newton method but for different types of boundary conditions. In the paper, He’s method, the Adomian method and Lagrange multiplier for special ODEs were given in detail, numerical results being obtained for Duffing and Titchmarch equations. We considered our applications the interval ( 0 , b ) and general conditions B 1 u ( 0 ) , B 2 u ( b ) for a linear and a nonlinear case of f x , u , u .
In [23] spectral problems of the nonlocal SLP with an integral B 2 u ( b ) were studied. Kernel of the operator, properties of the first eigenvalue and oscillation properties of eigenfunctions to the nonlocal problem were expressed. Additionally, the solution of the Cauchy problem for the SL equation on a star graph was constructed in [24].
For fractional differential equations, VIM could also be a very powerful instrument, in which Equations (36) and (37) are written using the Caputo fractional derivative, see [25]. This is our intention for the new study.
Nonlinear RSL problems could appear in the case of non-Newtonian fluid flows. Variational estimation methods are efficient techniques for finding analytical approximate solutions for a class of problems and also for optimal problems when looking for a minimum, using the functional of energy.

Author Contributions

Conceptualization, E.C.C.; methodology, E.C.C. and C.D.B.; software, E.C.C. and C.D.B.; validation, E.C.C. and C.D.B.; formal analysis, E.C.C.; investigation, E.C.C. and C.D.B.; writing—original draft preparation, E.C.C. and C.D.B.; writing—review and editing, E.C.C. and C.D.B.; supervision, E.C.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors want to thank to the referees who allowed us to improve ourselves and our article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dautray, R.; Lions, J.L. Mathematical Analysis and Numerical Methods for Science and Technology; Springer: Berlin/Heidelberg, Germany, 2000. [Google Scholar]
  2. Fowler, A.C. Mathematical Models in the Applied Sciences; Cambridge University Press: Cambridge, UK, 1997. [Google Scholar]
  3. Tyn Myint-U, T.; Debnath, L. Linear Partial Differential Equations for Scientists and Engineers, 4th ed.; Springer: Berlin, Germany, 2007; ISBN 0-8176-4393-1. [Google Scholar]
  4. Guenther, R.B.; Lee, J.W. Sturm-Liouville Problems. Theory and Numerical Implementation; CRC Press: Boca Raton, FL, USA, 2019; ISBN 9781138345430. [Google Scholar]
  5. Perera, U.; Böckmann, C. Solutions of Sturm-Liouville Problems. Mathematics 2020, 8, 2074. [Google Scholar] [CrossRef]
  6. Hassana, A.A. Green’s Function Solution of Non-Homogenous Regular Sturm-Liouville Problem. J. Appl. Comput. Math. 2017, 6, 362. [Google Scholar] [CrossRef] [Green Version]
  7. Klimek, M. Spectrum of Fractional and Fractional Prabhakar Sturm–Liouville Problems with Homogeneous Dirichlet Boundary Conditions. Symmetry 2021, 13, 2265. [Google Scholar] [CrossRef]
  8. Razdan, A.K.; Ravichandran, V. Fundamentals of Partial Differential Equations; Springer: Singapore, 2022. [Google Scholar] [CrossRef]
  9. Guo, H.; Qi, J. Sturm–Liouville Problems Involving Distribution Weights and an Application to Optimal Problems. J. Optim. Theory Appl. 2020, 184, 842–857. [Google Scholar] [CrossRef]
  10. Altintan, D.; Ugur, O. Variational iteration method for Sturm-Liouville differential equations. Comput. Math. Appl. 2008, 58, 322–328. [Google Scholar] [CrossRef] [Green Version]
  11. Johnson, C. Numerical Solution of Partial Differential Equations by the Finite Element Method; Cambridge University Press: Cambridge, UK, 1987; ISBN 0521345146. [Google Scholar]
  12. Chen, C.N. A survey of nonlinear Sturm-Liouville equations. In Sturm-Liouville Theory: Past and Present; Amrein, W.O., Hinz, A.M., Eds.; Pearson D.P.: Basel, Switzerland, 2005; pp. 201–216. [Google Scholar]
  13. He, J.H.; Wu, X.H. Variational iteration method: New development and applications. Comput. Math. Appl. 2007, 54, 881–894. [Google Scholar] [CrossRef] [Green Version]
  14. Hayat, T.; Khan, M.; Ayub, M. On the non-linear flows with slip boundary condition. Z. Angew. Math. Phys. 2005, 56, 1012–1029. [Google Scholar] [CrossRef]
  15. Yusufoglu, E.; Bekir, A. On the extended tanh method applications of nonlinear equations. Int. J. Nonlinear Sci. 2007, 4, 10–16. [Google Scholar]
  16. Neamaty, A.; Darzi, R. Comparison between the Variational Iteration Method and the Homotopy Perturbation Method for the Sturm-Liouville Differential Equation. Bound. Value Probl. 2010, 2010, 317369. [Google Scholar] [CrossRef] [Green Version]
  17. Zhang, T.-T.; Jia, L.; Wang, Z.-C. Analytic Solution for Steady Slip Flow between Parallel Plates with Micro-Scale Spacing. Chin. Phys. Lett. 2008, 25, 180. [Google Scholar] [CrossRef]
  18. Nadeem, M.; Yao, S. Solving system of partial differential equations using variational iteration method with He’s polynomials. Int. J. Math. Comput. Sci. 2019, 19, 203–211. [Google Scholar] [CrossRef]
  19. Fiza, M.; Chohan, F.; Ullah, H.; Islam, S.; Bushnaq, S. An extension of the optimal homotopy asymptotic method with applications to nonlinear coupled partial differential equations. J. Math. Comput. Sci. 2019, 19, 218–229. [Google Scholar] [CrossRef] [Green Version]
  20. Nadeem, M.; Yao, S.; Parveen, N. Solution of Newell-Whitehead-Segel equation by variational iteration method with He’s polynomials. J. Math. Comput. Sci. 2020, 20, 21–29. [Google Scholar] [CrossRef] [Green Version]
  21. Herron, I.H. Solving singular boundary value problems for ordinary differential equations. Caribb. J. Math. Comput. Sci. 2013, 15, 1–30. [Google Scholar]
  22. Al-Khaled, K.; Hazaimeh, A. Comparison Methods for Solving Non-Linear Sturm–Liouville Eigenvalues Problems. Symmetry 2020, 12, 1179. [Google Scholar] [CrossRef]
  23. Liu, Z.; Qi, J. The Properties of Eigenvalues and Eigenfunctions for Nonlocal Sturm–Liouville Problems. Symmetry 2021, 13, 820. [Google Scholar] [CrossRef]
  24. Kanguzhin, B.; Aimal Rasa, G.H.; Kaiyrbek, Z. Identification of the Domain of the Sturm–Liouville Operator on a Star Graph. Symmetry 2021, 13, 1210. [Google Scholar] [CrossRef]
  25. Nagdy, A.; Hashem, K.H. Numerical solutions of nonlinear fractional differential equations by variational iteration method. J. Nonlinear Sci. Appl. 2021, 14, 54–62. [Google Scholar] [CrossRef]
Figure 1. (a) Fourier equation; (b) Graetz–Nusselt equation.
Figure 1. (a) Fourier equation; (b) Graetz–Nusselt equation.
Mathematics 10 03728 g001
Figure 2. (a) Collatz problem; (b) Airy problem.
Figure 2. (a) Collatz problem; (b) Airy problem.
Mathematics 10 03728 g002
Figure 3. The eigenvalues determination and the corresponding eigenfunctions example (15).
Figure 3. The eigenvalues determination and the corresponding eigenfunctions example (15).
Mathematics 10 03728 g003
Figure 4. The eigenvalues determination and the corresponding eigenfunctions example 2.1.
Figure 4. The eigenvalues determination and the corresponding eigenfunctions example 2.1.
Mathematics 10 03728 g004
Table 1. Examples of regular Sturm–Liouville problems.
Table 1. Examples of regular Sturm–Liouville problems.
EquationAssociated FunctionsName
u + λ u = f , x ( a , b ) p 1 , s 1 , ρ 0 Fourier
( x u ) = 2 λ x ( 1 x 2 ) u , x ( 0 , 1 ) p ( x ) = x , s ( x ) = 2 x ( 1 x 2 ) Graetz–Nusselt
x 6 u + ( 3 / 4 ) x 4 u = λ u , x ( 1 , 2 ) p ( x ) = x 6 , s ( x ) = 3 4 x 4 Collatz
u + x u = λ x u , x ( 0 , 1 ) p ( x ) = x , s ( x ) = 2 x ( 1 x 2 ) Airy
Table 2. Examples of differential equations and their SL form.
Table 2. Examples of differential equations and their SL form.
TypeEquationSturm–Liouville Form
Legendre u 2 x 1 x 2 + μ 1 x 2 u = 0 , ( ( 1 x 2 ) u ) + μ u = 0
Chebysev ( 1 x 2 ) u x u + n 2 u = 0 ( 1 x 2 u ) = n 2 1 x 2 u
Bessel x 2 u + x u + ( λ 2 n 2 ) y = 0 , ( x u ) + n 2 x u = λ 2 x u
Table 3. Examples of Sturm–Liouville problems.
Table 3. Examples of Sturm–Liouville problems.
( sin ( x ) u ) + cos ( x ) u = λ x 2 u , x ( 0 , b ) p = sin ( x ) , ρ = cos ( x ) s = x 2
( x 2 u ) x sin ( x ) u = λ x u , x ( 0 , b ) p = x , s = 2 x ( 1 x 2 )
Table 4. Examples of Sturm–Liouville problems with a = 0 , b > 0 , x ( 0 , b ) , case 1: a 1 · a 2 = 0 .
Table 4. Examples of Sturm–Liouville problems with a = 0 , b > 0 , x ( 0 , b ) , case 1: a 1 · a 2 = 0 .
case 1.1a a 1 0 , a 2 = 0 b 1 0 , b 2 = 0 u ( 0 ) = 0 , u ( b ) = 0
case 1.1b a 1 0 , a 2 = 0 b 1 = 0 , b 2 0 u ( 0 ) = 0 , u ( b ) = 0
case 1.1c a 1 0 , a 2 = 0 b 1 0 , b 2 0 u ( 0 ) = 0 , b 1 u ( b ) + b 2 u ( b ) = 0
case 1.2 a 1 = 0 , a 2 0 b 1 = 0 , b 2 0 u ( 0 ) = 0 , b 1 u ( b ) + b 2 u ( b ) = 0
Table 5. Examples of Sturm–Liouville problems with a = 0 , b > 0 , x ( 0 , b ) , case 2: a 1 · a 2 0 .
Table 5. Examples of Sturm–Liouville problems with a = 0 , b > 0 , x ( 0 , b ) , case 2: a 1 · a 2 0 .
case 2.1 a 1 0 , a 2 = a 1 b 1 = 0 , b 2 0 u ( 0 ) = 0 , b 1 u ( b ) + b 2 u ( b ) = 0
case 2.2 a 1 > 0 , a 2 = 1 b 1 = 0 , b 2 0 a 1 u ( 0 ) u ( 0 ) = 0 , b 1 u ( b ) + b 2 u ( b ) = 0
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Cipu, E.C.; Barbu, C.D. Variational Estimation Methods for Sturm–Liouville Problems. Mathematics 2022, 10, 3728. https://doi.org/10.3390/math10203728

AMA Style

Cipu EC, Barbu CD. Variational Estimation Methods for Sturm–Liouville Problems. Mathematics. 2022; 10(20):3728. https://doi.org/10.3390/math10203728

Chicago/Turabian Style

Cipu, Elena Corina, and Cosmin Dănuţ Barbu. 2022. "Variational Estimation Methods for Sturm–Liouville Problems" Mathematics 10, no. 20: 3728. https://doi.org/10.3390/math10203728

APA Style

Cipu, E. C., & Barbu, C. D. (2022). Variational Estimation Methods for Sturm–Liouville Problems. Mathematics, 10(20), 3728. https://doi.org/10.3390/math10203728

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