Research

# On the spectral investigation of the scattering problem for some version of one-dimensional Schrödinger equation with turning point

Zaki FA El-Raheem1* and AH Nasser2

Author Affiliations

1 Department of Mathematics, Faculty of Education, Alexandria University, Alexandria, Egypt

2 Department of Mathematics, Faculty of Industrial Education, Helwan University, Cairo, Egypt

For all author emails, please log on.

Boundary Value Problems 2014, 2014:97  doi:10.1186/1687-2770-2014-97

 Received: 30 January 2014 Accepted: 1 April 2014 Published: 6 May 2014

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited.

### Abstract

In this paper we introduce and investigate the eigenvalues and the normalizing numbers as well as the scattering function for some version of the one-dimensional Schrödinger equation with turning point on the half line.

MSC: 58C40, 34L25.

##### Keywords:
initial value problem; the eigenvalues; normalizing numbers; scattering function; asymptotic formula

### 1 Introduction

The solution of many problems of mathematical physics are reduced to the spectral investigation of a differential operator. The differential operator is called regular if its domain is finite and its coefficients are continuous, otherwise it is called a singular differential operator. The Sturm-Liouville theory occupies a central position in the spectral theory of regular operator. During the development of quantum mechanics there was an increase in the interest of spectral theory of singular operators, on which we will restrict our attention. The first basic role in the development of the spectral theory of singular operators dates back to Titchmarsh [1]. He gave a new approach in the spectral theory of singular differential operator of the second order by using contour integration. Also Levitan [2] gave a new method, he obtained the eigenfunction expansion in an infinite interval by taking the limit of a regular case. In the last 35 or so years, due to the needs of mathematical physics, in particular, quantum mechanics, the question of solving various spectral problems with explosive factor has appeared in the study of geophysics and electromagnetic fields; see [3,4]. The spectral theory of differential operators with explosive factor is studied by Tikhonov [5], Gasymov [6]. For earlier results on various aspects of solvability theory of boundary value problems and spectral theory in the half line case, the situation closely related to the principal topic of this paper, we refer, for instance, to [7-10]. Notice that the paper [11] presented an approximate construction of the Jost function for some Sturm-Liouville boundary value problem in the case ρ ( x ) = 1 by means of the collocation method. In the present paper we introduce and investigate the eigenvalues and the normalizing numbers as well as the scattering function for some version of the one-dimensional Schrödinger equation with turning point on the half line as in (1.1), (1.2). In [12,13], and [14] the weight functions introduced are considered as applications of the discontinuous wave speed problem on a non-homogeneous medium as in our case, while the introduction of the weight function ρ ( x ) which is given by (1.3) as ± signs causes an excess of analytical difficulties. In [15] the author studied the spectral property in a finite interval, while in the present work we consider the half line which gives rise both to a continuous and a discrete spectrum; the latter is treated by the scattering function. In [16] the author considered the weight function of the form

ρ ( x ) = { α 2 ; Im α 0 , 0 x a < π , 1 ; a < x π ,

and the spectra were both continuous and discrete as in our problem. We must notice that the result of this paper is a starting point in calculating the regularized trace formula and solving the inverse scattering problem, which will be investigated later on.

Consider the initial value problem

y + q ( x ) y = μ ρ ( x ) y , 0 x < , (1.1)

y ( 0 ) h y ( 0 ) = 0 , h > 0 , (1.2)

where

ρ ( x ) = { 1 ; 0 x 1 , 1 ; 1 < x < , (1.3)

q ( x ) is a finite real valued function which satisfies

0 ( 1 + x ) | q ( x ) | d x < ,

and μ is a complex spectral parameter. To study the eigenvalues of (1.1)-(1.2), we first consider the case when q ( x ) 0 and h = 0 .

For q ( x ) 0 and h = 0 problem (1.1)-(1.2) takes the form

y = μ ρ ( x ) y , 0 x , (1.4)

y ( 0 ) = 0 , (1.5)

μ = λ 2 . (1.6)

From now on we consider Im λ 0 because according to (1.6) μ covers all the complex plane. Denote by φ o ( x , λ ) the solution of (1.4) with the initial conditions φ o ( 0 , λ ) = 1 , φ o ( 0 , λ ) = 0 . According to (1.3), (1.4) is equivalent to the two equations

y = λ 2 y , 0 x < 1 , y = λ 2 y , 1 < x . (1.7)

It is easy to see that

φ o ( x , λ ) = { cosh λ x ; 0 x < 1 , a o ( λ ) e i λ x + b o ( λ ) e i λ x ; 1 < x < , (1.8)

where a o ( λ ) , b o ( λ ) are calculated from the requirements φ o ( 1 0 , λ ) = φ o ( 1 + 0 , λ ) and φ o ( 1 0 , λ ) = φ o ( 1 + 0 , λ ) , so that (1.8) takes the form

φ o ( x , λ ) = { cosh λ x ; 0 x < 1 , e i λ 2 ( cosh λ i sinh λ ) e i λ x + e i λ 2 ( cosh λ + i sinh λ ) e i λ x ; 1 < x < . (1.9)

For Im λ = 0 , the function φ o ( x , λ ) does not belong to L 2 ( 0 , ) also, for Im λ > 0 , e i λ x 0 as x whereas e i λ x as x , so that it is convenient to consider

e i λ ( cosh λ + i sinh λ ) = 0 (1.10)

as the equation of the eigenvalues μ o = λ o 2 .

From this we have λ o = ( n + 1 4 ) π i , n = 0 , 1 , ± 1 , ± 2 , or

μ n o = ( n + 1 4 ) 2 π 2 , n = 1 , 2 , . (1.11)

Together with the solution φ o ( x , λ ) of (1.4) we introduce the second solution f o ( x , λ ) , which is known as the Jost solution. This solution is defined by the condition

f o ( x , λ ) e i λ x , x . (1.12)

With the aid of (1.7), we have

f o ( x , λ ) = { c o ( λ ) cosh λ x + d o ( λ ) sinh λ x , 0 x < 1 , e i λ x , 1 < x < ,

where the coefficients c o ( λ ) , d o ( λ ) are calculated from the requirements f o ( 1 0 , λ ) = f o ( 1 + 0 , λ ) and f o ( 1 0 , λ ) = f o ( 1 + 0 , λ ) , and the solution becomes

f o ( x , λ ) = { e i λ ( cosh λ i sinh λ ) cosh λ x + e i λ ( i cosh λ sinh λ ) sinh λ x , 0 x < 1 , e i λ x , 1 < x < . (1.13)

It should be noted, here, that the equation of the eigenvalues can be obtained, also, from the condition that the solution f o ( x , λ ) L 2 ( 0 , ; ρ ) ; this condition implies that f o ( 0 , λ ) = 0 , which is the same as (1.10).

Now for q ( x ) 0 , h 0 we denote by f ( x , λ ) the solution of (1.1) which satisfies the condition

f ( x , λ ) e i λ x , x , < λ < .

For x > 1 , (1.1) takes the form y + q ( x ) y = λ 2 y , and in the following, we study its solution and the related spectrum. From [4] this solution has the following representation:

f ( x , λ ) = e i λ x + x K ( x , t ) e i λ t d t , 1 < x < , (1.14)

where Im λ 0 , K ( x , x ) = 1 2 0 x q ( t ) d t , 1 < x < .

For 0 x 1 , the solution f ( x , λ ) has the form

f ( x , λ ) = { a ( λ ) φ ( x , λ ) + b ( λ ) θ ( x , λ ) , 0 x 1 , e i λ x + x K ( x , t ) e i λ t d t , 1 < x < , (1.15)

where φ ( x , λ ) , θ ( x , λ ) is the fundamental system of solutions of (1.1) subject to the initial conditions

φ ( 0 , λ ) = 1 , φ ( 0 , λ ) = h , θ ( 0 , λ ) = 0 , θ ( 0 , λ ) = 1 , (1.16)

where the coefficients a ( λ ) , b ( λ ) are calculated from the requirements f ( 1 0 , λ ) = f ( 1 + 0 , λ ) , f ( 1 0 , λ ) = f ( 1 + 0 , λ ) , from which

a ( λ ) = f ( 1 , λ ) θ ( 1 , λ ) f ( 1 , λ ) θ ( 1 , λ ) , b ( λ ) = f ( 1 , λ ) φ ( 1 , λ ) f ( 1 , λ ) φ ( 1 , λ ) . (1.17)

Further, (1.1), for 0 x 1 , takes the form y + q ( x ) y = λ 2 y , and the fundamental system of solution of this follows from [[4], p.18] by the representation

φ ( x , λ ) = sinh λ x λ + 0 x B ( x , t ) sinh λ t λ d t , A ( x , x ) = 1 2 0 x q ( t ) d t , A ( x , 0 ) = 0 , A ( 0 , 0 ) = h , (1.18)

θ ( x , λ ) = cosh λ x + 0 x A ( x , t ) cosh λ t d t , B ( x , x ) = 1 2 0 x q ( t ) d t , B t | t = 0 = 0 , B ( 0 , 0 ) = 1 . (1.19)

Now we find the characteristic equation of the eigenvalues of (1.1)-(1.2). Since the solution (1.15) belongs to L 2 ( 0 , ) , Im λ > 0 it follows that, for μ = λ 2 to be an eigenvalue, it must satisfy the initial condition (1.2), namely

f ( 0 , λ ) h f ( 0 , λ ) = 0 . (1.20)

From (1.15) and (1.16) we have

f ( 0 , λ ) h f ( 0 , λ ) = b ( λ ) = f ( 1 , λ ) φ ( 1 , λ ) f ( 1 , λ ) φ ( 1 , λ ) . (1.21)

In the following lemmas we study some properties of the eigenvalues of problem (1.1)-(1.2).

Lemma 1.1Under the conditions q ( x ) > 0 ( 0 < x < ), the roots of (1.20), for Im λ > 0 , are simple and lie only on the imaginary axis.

Proof Let λ o , where Im λ o > 0 , be a zero of the function f ( 0 , λ ) h f ( 0 , λ ) , so that

f ( 0 , λ o ) h f ( 0 , λ o ) = 0 . (1.22)

We prove that λ o = i τ o , τ o > 0 . Since f ( x , λ o ) is a solution of (1.1) we have

f ( x , λ o ) + q ( x ) f ( x , λ o ) = λ o 2 ρ ( x ) f ( x , λ o ) , (1.23)

multiplying both sides of this by f ( x , λ o ) ¯ and integrating both sides from 0 to ∞, we have

0 f ( x , λ o ) f ( x , λ o ) ¯ d x + 0 q ( x ) | f ( x , λ o ) | 2 d x = λ o 2 0 ρ ( x ) | f ( x , λ o ) | 2 d x .

Integrating the first integral by parts and using (1.22), (1.15) we obtain

λ o 2 = 0 { | f ( x , λ o ) | 2 + q ( x ) | f ( x , λ o ) | 2 } d x + h | f ( 0 , λ o ) | 2 0 ρ ( x ) | f ( x , λ o ) | 2 d x , (1.24)

where 0 ρ ( x ) | f ( x , λ o ) | 2 d x 0 , from which we deduce that λ o 2 is real and hence λ o is pure imaginary. We turn now to the proof that the roots are simple from (1.22), this is carried out by proving that f ( 0 , λ ) h f ( 0 , λ ) = 0 implies [ f ˙ ( 0 , λ ) h f ˙ ( 0 , λ ) ] 0 , where ‘dot’ denotes differentiation with respect to λ.

Integrating the difference [ f ˙ ( x , λ ) × (1.23) ] [ f ( x , λ ) × d d λ (1.23) ] with respect to x from 0 to ∞ and using (1.20) we get after some calculation that

f ( 0 , λ ) [ f ˙ ( 0 , λ ) h f ˙ ( 0 , λ ) ] = 2 λ 0 ρ ( x ) f 2 ( x , λ ) d x . (1.25)

We prove the reality of f ( x , λ ) .

For x > 1 , λ = i τ the function f ( x , λ ) = e τ x + x K ( x , t ) e τ t d t is real because reality of K ( x , t ) comes from the reality of q ( x ) .

To prove that, for 0 x < 1 , we observe that φ and θ are real. Let λ = i τ ; since φ ( x , λ ) is a solution of (1.1)-(1.2), we have

φ ( x , λ ) + q ( x ) φ ( x , λ ) = τ 2 ρ ( x ) φ ( x , λ ) , φ ( 0 , λ ) = 1 , φ ( 0 , λ ) = h . (1.26)

Taking the conjugate of (1.26) we have

φ ( x , λ ) ¯ + q ( x ) φ ( x , λ ) ¯ = τ 2 ρ ( x ) φ ( x , λ ) ¯ , φ ( 0 , λ ) ¯ = 1 , φ ( 0 , λ ) ¯ = h . (1.27)

It is clear, from (1.26) and (1.27), that φ ( x , λ ) = φ ( x , λ ) ¯ . In a similar way we can prove that θ ( x , λ ) is also real so that the solution f ( x , λ ) for 0 x < 1 is real from which we have f 2 ( x , λ ) = | f ( x , λ ) | 2 and (1.25) takes the form

f ( 0 , λ ) [ f ˙ ( 0 , λ ) h f ˙ ( 0 , λ ) ] = 2 λ 0 ρ ( x ) | f ( x , λ ) | 2 d x . (1.28)

From (1.28) we see that d d λ [ f ( 0 , λ ) h f ( 0 , λ ) ] 0 , which completes the proof. □

Remark 1 For Im λ n > 0 and f ( 0 , λ n ) h f ( 0 , λ n ) = 0 , the function f ( 0 , λ n ) is the eigenfunction of problem (1.1)-(1.2) that corresponds to the negative eigenvalues μ n = λ n 2 = χ n 2 .

Lemma 1.2For all Re λ 0 the function f ( 0 , λ ) h f ( 0 , λ ) does not tend to zero, i.e.

f ( 0 , λ ) h f ( 0 , λ ) 0 , Re λ 0 , < λ < . (1.29)

Proof Since the function f ( x , λ ) is the solution of (1.1), f ( x , λ ) is also a solution, and it can be shown that these two solutions are linearly independent and their Wronskian is

W [ f ( x , λ ) , f ( x , λ ) ] = 2 i λ , (1.30)

so that W [ f ( x , λ ) , f ( x , λ ) ] 0 , for Re λ 0 , so that f ( x , λ ) and f ( x , λ ) is a fundamental system of solutions of (1.1). In particular, putting x = 0 into (1.30) we have

f ( 0 , λ ) f ( 0 , λ ) f ( 0 , λ ) f ( 0 , λ ) = 2 i λ . (1.31)

To prove that f ( 0 , λ ) h f ( 0 , λ ) 0 , Re λ 0 , < λ < , assume to the contrary i.e. f ( 0 , λ ) h f ( 0 , λ ) = 0 , Re λ 0 , < λ < . From (1.31) and (1.20) we reach to contradiction to the assumption, and, consequently, we deduce that f ( 0 , λ ) h f ( 0 , λ ) 0 , Re λ 0 , < λ < . Notice that f ( x , λ ) ¯ = f ( x , λ ) . □

Lemma 1.3For all Re λ 0 the following equality holds:

2 i λ φ ( x , λ ) f ( 0 , λ ) h f ( 0 , λ ) = f ( x , λ ) S ( λ ) f ( x , λ ) , (1.32)

where φ ( x , λ ) is the solution of problem (1.1)-(1.2) and the function

S ( λ ) = f ( 0 , λ ) h f ( 0 , λ ) f ( 0 , λ ) h f ( 0 , λ ) (1.33)

satisfies the properties

S ( λ ) ¯ = S ( λ ) , | S ( λ ) | = 1 , < λ < . (1.34)

It should be noted here that the function S ( λ ) defined by (1.33) is called the scattering function of problem (1.1)-(1.2) and the function f ( 0 , λ ) h f ( 0 , λ ) is called the denominator of S ( λ ) .

Proof As mentioned before (1.30) for all Re λ 0 , f ( x , λ ) and f ( x , λ ) is a fundamental system of solutions of (1.1)-(1.2), so that any linear combination of them is again a solution of (1.1)-(1.2):

φ ( x , λ ) = A ( λ ) f ( x , λ ) + B ( λ ) f ( x , λ ) , (1.35)

where A ( λ ) , B ( λ ) are calculated from the initial conditions φ ( 0 , λ ) = 1 , φ ( 0 , λ ) = h in the form

A ( λ ) = f ( 0 , λ ) h f ( 0 , λ ) 2 i λ , B ( λ ) = f ( 0 , λ ) h f ( 0 , λ ) 2 i λ . (1.36)

Substituting (1.36) into (1.35) we arrive at the required formula (1.32). Further, since f ( x , λ ) ¯ = f ( x , λ ) , it follows from (1.33) that

S ( λ ) = f ( 0 , λ ) h f ( 0 , λ ) f ( 0 , λ ) h f ( 0 , λ ) = f ( 0 , λ ) h f ( 0 , λ ) ¯ f ( 0 , λ ) h f ( 0 , λ ) ,

from which we have

| S ( λ ) | = | f ( 0 , λ ) h f ( 0 , λ ) ¯ | | f ( 0 , λ ) h f ( 0 , λ ) | = 1 ,

and

S ( λ ) ¯ = ( f ( 0 , λ ) h f ( 0 , λ ) f ( 0 , λ ) h f ( 0 , λ ) ) ¯ = f ( 0 , λ ) h f ( 0 , λ ) f ( 0 , λ ) h f ( 0 , λ ) = S ( λ ) .

□

### 2 The asymptotic formulas of eigenvalues and normalizing numbers

The eigenvalues μ = λ 2 of problem (1.1)-(1.2) are the roots of the equation

f ( 0 , λ ) h f ( 0 , λ ) = 0 , Im λ > 0 . (2.1)

In the following we prove that (2.1) has an infinite number of roots and find their asymptotic formula. From (1.15), (1.17), (1.18), and (1.19) we have

f ( 0 , λ ) h f ( 0 , λ ) = b ( λ ) = f ( 1 , λ ) φ ( 1 , λ ) f ( 1 , λ ) φ ( 1 , λ ) . (2.2)

Now, we calculate the asymptotic formula of f ( 1 , λ ) , f ( 1 , λ ) , φ ( 1 , λ ) and φ ( 1 , λ ) . Integrating (1.15) by parts we have, for x 1 , Im λ > 0 ,

f ( x , λ ) = e i λ x K ( x , x ) i λ e i λ x x e i λ x i λ K ( x , x ) t d t , (2.3)

f ( x , λ ) = i λ e i λ x K ( x , x ) e i λ x + x K ( x , x ) x e i λ x d t . (2.4)

Similarly from (1.18) we have

φ ( x , λ ) = cosh λ x + A ( x , x ) sinh λ x λ 1 λ 0 x A ( x , t ) t sinh λ t d t , (2.5)

φ ( x , λ ) = λ sinh λ x + A ( x , x ) cosh λ x + 0 x A ( x , t ) x cosh λ t d t . (2.6)

The following group of inequalities follows from (2.3)-(2.6):

f ( 1 , λ ) = e i λ + O ( e Im λ λ ) , (2.7)

f ( 1 , λ ) = i λ e i λ + O ( e Im λ ) , (2.8)

φ ( 1 , λ ) = cosh λ + O ( e | Re λ | λ ) , (2.9)

φ ( 1 , λ ) = λ sinh λ + O ( e | Re λ | ) . (2.10)

Substituting (2.7)-(2.10) into (2.2), we obtain

f ( 0 , λ ) h f ( 0 , λ ) = i λ e i λ [ cosh λ + i sinh λ ] + O ( e Im λ + | Re λ | ) , λ 0 , Im λ > 0 ; (2.11)

comparing (1.10) and (2.11) we see that f ( 0 , λ ) h f ( 0 , λ ) and f o ( 0 , λ ) = e i λ [ cosh λ + i sinh λ ] have the same number of zeros inside the quadratic contour Γ n where { Γ n : | Re λ | π ( n 1 4 ) , 0 < Im λ π ( n 1 4 ) } , but since f o ( 0 , λ ) has exactly n zeros, namely λ k o = i π ( k 1 4 ) , k = 1 , 2 , , n , f ( 0 , λ ) h f ( 0 , λ ) has an infinite number of zeros, as n , with limiting point at infinity. Denote by λ n the zeros of f ( 0 , λ ) h f ( 0 , λ ) = 0 , so that, by the Rouche theorem, we have

λ n = i ( n + 1 4 ) + ε n . (2.12)

To make (2.12) more accurate, we must refine (2.11). With the aid of Lemma 1.1, λ n lies on the imaginary axis, so that it is sufficient to know the asymptotic of f ( 0 , λ ) h f ( 0 , λ ) for small λ. Let λ = i τ , τ > 0 , we find the asymptotic formula of f ( 0 , i τ ) h f ( 0 , i τ ) for τ . From (2.3), (2.4), (2.5), and (2.6), we have

f ( 1 , i τ ) = e τ + K ( 1 , 1 ) e τ τ + o ( e τ τ ) , 7 f ( 1 , i τ ) = τ e τ K ( 1 , 1 ) e τ + o ( e τ ) , φ ( 1 , i τ ) = cos τ + A ( 1 , 1 ) sin τ τ + o ( 1 τ ) , φ ( 1 , i τ ) = τ sin τ + A ( 1 , 1 ) cos τ + o ( 1 ) , (2.13)

substituting (2.13) into f ( 0 , i τ ) h f ( 0 , i τ ) = 0 , and putting λ n = i τ n we have

cos τ n sin τ n + A ( 1 , 1 ) + K ( 1 , 1 ) τ n ( cos τ n + sin τ n ) + o ( e τ n τ n ) = 0 , (2.14)

and from this and by virtue of the inequality | cos τ n | δ > 0 n, we have

1 tan τ n + α τ n + β τ n tan + o ( 1 τ n ) = 0 , where  α = A ( 1 , 1 ) + K ( 1 , 1 ) , β = A ( 1 , 1 ) K ( 1 , 1 ) . (2.15)

From (2.12), it is easy to see that

τ n = ( n + 1 4 ) + ε n , tan τ n = 1 + 2 ε n + O ( 1 n 3 ) , 1 τ n = 1 n π + O ( 1 n ) . (2.16)

The estimation of ε n follows from (2.15) and (2.16) in the form

ε n = 1 + β α n π + o ( 1 n ) . (2.17)

Therefore

τ n = π ( n + 1 4 ) 1 + c o n π + o ( 1 n ) , c o = 1 π 0 1 q ( t ) d t . (2.18)

Finally

λ n = i [ π ( n + 1 4 ) 1 + c o n π + o ( 1 n ) ] , c o = 1 π 0 1 q ( t ) d t . (2.19)

Definition (The normalizing numbers)

The numbers

a n = def 0 ρ ( x ) | f ( x , λ n ) | 2 d x (2.20)

are called the normalizing numbers of problem (1.1)-(1.2) (notice that f ( x , λ n ) are the eigenfunctions of problem (1.1)-(1.2) corresponding to the eigenvalues λ n ). From (1.28) and the reality of f ( x , λ n ) , we have

a n = def 0 ρ ( x ) | f ( x , λ n ) | 2 d x = [ f ˙ ( 0 , λ n ) h f ˙ ( 0 , λ n ) ] f ( 0 , λ n ) 2 λ n . (2.21)

To evaluate the asymptotic formula of a n we evaluate the asymptotic formula of the right hand side of (2.21). From (1.15), (1.17) we have

[ f ˙ ( 0 , λ n ) h f ˙ ( 0 , λ n ) ] f ( 0 , λ n ) = [ b ˙ ( λ ) + ( 1 + h ) a ˙ ( λ ) ] a ( λ ) , (2.22)

where dots and dashes denote the differentiation with respect to λ and x, respectively, a ( λ ) and b ( λ ) are given by (1.17)

a ( λ ) = f ( 1 , λ ) θ ( 1 , λ ) f ( 1 , λ ) θ ( 1 , λ ) , b ( λ ) = f ( 1 , λ ) φ ( 1 , λ ) f ( 1 , λ ) φ ( 1 , λ ) ,

from which it follows that

a ˙ ( λ ) = f ˙ ( 1 , λ ) θ ( 1 , λ ) + f ( 1 , λ ) θ ˙ ( 1 , λ ) f ˙ ( 1 , λ ) θ ( 1 , λ ) f ( 1 , λ ) θ ˙ ( 1 , λ ) , b ˙ ( λ ) = f ˙ ( 1 , λ ) φ ( 1 , λ ) + f ( 1 , λ ) φ ˙ ( 1 , λ ) f ˙ ( 1 , λ ) φ ( 1 , λ ) f ( 1 , λ ) φ ˙ ( 1 , λ ) . (2.23)

From (1.18), using integration by parts and then putting x = 1 , λ = i τ , we obtain

φ ( 1 , i τ ) = cos τ + A ( 1 , 1 ) sin τ τ + A t ( 1 , 1 ) cos τ τ 2 + o ( 1 τ 2 ) , φ ( 1 , i τ ) = τ sin τ + A ( 1 , 1 ) cos τ + A t ( 1 , 1 ) sin τ τ + o ( 1 τ ) , φ ˙ ( 1 , i τ ) = sin τ + A ( 1 , 1 ) cos τ τ [ A ( 1 , 1 ) + A t ( 1 , 1 ) ] sin τ τ 2 + o ( 1 τ 2 ) , φ ˙ ( 1 , i τ ) = τ cos τ [ 1 + A ( 1 , 1 ) ] sin τ A t ( 1 , 1 ) cos τ τ + A t ( 1 , 1 ) sin τ τ 2 + o ( 1 τ 2 ) . (2.24)

From (1.19), carrying out a similar calculation with respect to θ, we obtain

l θ ( 1 , i τ ) = sin τ τ B ( 1 , 1 ) cos τ τ 2 + o ( 1 τ 2 ) , θ ( 1 , i τ ) = cos τ + B ( 1 , 1 ) sin τ τ + o ( 1 τ ) , θ ˙ ( 1 , i τ ) = cos τ τ [ 1 + B ( 1 , 1 ) ] sin τ τ 2 + 2 B ( 1 , 1 ) cos τ τ 3 + o ( 1 τ 3 ) , θ ˙ ( 1 , i τ ) = sin τ + B ( 1 , 1 ) cos τ τ B ( 1 , 1 ) sin τ τ 2 + o ( 1 τ 2 ) . (2.25)

With the aid of (1.15), similar expressions can be calculated with respect to f ( 1 , i τ ) :

f ( 1 , i τ ) = e τ + K ( 1 , 1 ) e τ τ + o ( e τ τ ) , f ( 1 , i τ ) = τ e τ K ( 1 , 1 ) e τ + o ( e τ ) , f ˙ ( 1 , i τ ) = e τ K ( 1 , 1 ) e τ τ K ( 1 , 1 ) e τ τ 2 + o ( e τ τ 2 ) , f ˙ ( 1 , i τ ) = τ e τ e τ + K ( 1 , 1 ) e τ + o ( e τ ) . (2.26)

From (2.21) and (2.22), the normalizing numbers a n can be written in the form

a n = [ b ˙ ( λ ) + ( 1 + h ) a ˙ ( λ ) ] a ( λ ) 2 λ n . (2.27)

We substitute (2.23), (2.24), (2.25), and (2.26) into (2.27), λ n = i τ n , and we find

a n = e 2 τ n { cos τ n sin τ n + cos 2 τ n + α + β 2 τ n cos τ n sin τ n + β 2 τ n sin 2 τ n + α 2 τ n cos 2 τ n + α 1 τ n cos τ n sin τ n + β 1 τ n cos 2 τ n + o ( 1 τ n ) } , (2.28)

where α = 1 + 2 K ( 1 , 1 ) + 2 A ( 1 , 1 ) , β = 1 + 2 A ( 1 , 1 ) , α 1 = B ( 1 , 1 ) + K ( 1 , 1 ) , and β 1 = K ( 1 , 1 ) B ( 1 , 1 ) . Further, from (2.16) and (2.17) we have

1 τ n = 1 n π [ 1 + o ( 1 τ n ) ] , cos 2 τ n = 3 2 β α n π + o ( 1 τ n ) , sin 2 τ n = 1 2 β α n π + o ( 1 τ n ) , cos τ n sin τ n = 1 2 [ 1 + o ( 1 τ n ) ] . (2.29)

By substituting from (2.29) into (2.28) we obtain the required asymptotic formula for a n :

a n = e 2 τ n { 2 + c 1 n + o ( 1 n ) } , (2.30)

where

c 1 = 8 h + 4 4 π + 3 8 π 0 1 q ( t ) d t + 9 8 π 1 q ( t ) d t .

### Competing interests

The authors declare that they have no competing interests.

### Authors’ contributions

The two authors typed read and approved the final manuscript also they contributed to each part of this work equally.

### Acknowledgements

We are indebted to an anonymous referee for a detailed reading of the manuscript and useful comments and suggestions, which helped us improve this work. This work is supported by the Research Support Unit of Alexandria University.

### References

1. Everitt, WN: Charles Sturm and the development of Sturm-Liouville theory in the years 1900 to 1950. Sturm-Liouville Theory, pp. 271–331. Birkhäuser, Basel (2005)

2. Levitan, BM: Inverse Sturm-Liouville Problems, Nauka, Moscow (1984) English transl., VNU Sci. Press, Utrecht (1987)

3. Alimov, SA: The work of A.N. Tikhonov on inverse problems for the Sturm-Liouville equation. Usp. Mat. Nauk. 31(6(192)), 84–88 English transl., Russ. Math. Surv. 31(6), 87-92 (1976) (1976)

4. Marchenko, VA: Sturm-Liouville Operators and Their Applications, Naukova Dumka, Kiev (1977) English transl., Birkhäuser (1986)

5. Tikhonov, AN: On the uniqueness of the solution of electric investigation problem. Doklady of the USSR. 69(6), 797–800 (in Russian) (1949)

6. Gasymov, MG: Direct and inverse problems of spectral analysis for a class of equations with discontinuous coefficients. Novosibirsk. (1977) (in Russian)

7. Aliev, AR, Mirzoev, SS: On boundary value problem solvability theory for a class of high-order operator-differential equations. Funkc. Anal. Prilozh.. 44(3), 63–65 English transl., Funct. Anal. Appl. 44(3), 209-211 (2010) (2010)

English transl., Funct. Anal. Appl. 44(3), 209-211 (2010)

Publisher Full Text

8. Mirzoev, SS, Aliev, AR, Rustamova, RA: Solvability conditions for boundary-value problems for elliptic operator-differential equations with discontinuous coefficient. Mat. Zametki. 92(5), 789–793 English transl., Math. Notes 92(5), 722-726 (2012) (2012)

English transl., Math. Notes 92(5), 722-726 (2012)

Publisher Full Text

9. Mirzoev, SS, Aliev, AR, Rustamova, RA: On the boundary value problem with the operator in boundary conditions for the operator-differential equation of second order with discontinuous coefficients. J. Mat. Fiz. Anal. Geom.. 9(2), 207–226 (2013)

10. Aliev, AR: On the generalized solution of the boundary-value problem for the operator-differential equations of the second order with variable coefficients. J. Mat. Fiz. Anal. Geom.. 2(1), 87–93 (2006)

11. Aliev, AR, Gasymova, SG, Gasymova, DG, Ahmadzadeh, ND: Approximate construction of the Jost function by the collocation method for Sturm-Liouville boundary value problem. Azerb. J. Math.. 3(2), 45–61 (2013)

12. Freiling, G, Yurko, VA: Inverse Sturm-Liouville Problems and Their Applications, Nova Science Publishers, New York (2001)

13. Efendiev, RF, Orudzhev, HD: Inverse wave spectral problem with discontinuous wave speed. J. Mat. Fiz. Anal. Geom.. 6(3), 255–265 (2010)

14. Efendiev, RF: Spectral analysis for one class of second-order indefinite non-self-adjoint differential operator pencil. Appl. Anal.. 90(12), 1837–1849 (2011). Publisher Full Text

15. El-Raheem, ZFA, Nasser, AH: On the spectral property of a Dirichlet problem with explosive factor. Appl. Math. Comput.. 138(2-3), 355–374 (2003). Publisher Full Text

16. Gasymov, MG, Kakhramanov, AS, Petrosyan, SK: On the spectral theory of linear differential operators with discontinuous coefficients. Akad. Nauk Azerbaĭdzhan. SSR Dokl.. 43(3), 13–16 (in Russian) (1987)