Abstract

For any 𝑛, the contour integral 𝑦=cosh𝑛+1𝑥𝐶(cosh(𝑧𝑠)/(sinh𝑧sinh𝑥)𝑛+1𝑑𝑧,𝑠2=𝜆, is associated with differential equation 𝑑2𝑦(𝑥)/𝑑𝑥2+(𝜆+𝑛(𝑛+1)/cosh2𝑥)𝑦(𝑥)=0. Explicit solutions for 𝑛=1 are obtained. For 𝑛=1, eigenvalues, eigenfunctions, spectral function, and eigenfunction expansions are explored. This differential equation which does have solution in terms of the trigonometric functions does not seem to have been explored and it is also one of the purposes of this paper to put it on record.

1. Introduction

When one considers eigenfunction expansions associated with second-order ordinary differential equations, as Titchmarsh does in his book [1], one is concerned with solutions of the equation 𝑑2𝑦(𝑥)𝑑𝑥2+𝑞(𝑥)𝑦(𝑥)=𝜆𝑦(𝑥)(1.1)along with certain boundary conditions, and one tends to say that the only case in which one can solve this equation explicitly in elementary terms for all 𝜆 is the case 𝑞(𝑥)=0, when the solutions are of course trigonometric functions.

Now in fact this is not true, and there is in particular one problem which does not seem to have been explored, and it is the purpose of this paper to put it on record. Here is the problem: 𝑑2𝑦(𝑥)𝑑𝑥2+𝜆+𝑛(𝑛+1)sech2𝑥𝑦(𝑥)=0,(1.2) which can be solved explicitly in elementary terms when 𝑛 is integral. The explicit solution was known to Kamke [2], but Kamke does not anyway explore the consequences for eigenfunction expansions nor does Titchmarsh discuss this problem, although he does discuss problems close to it, for example, 𝑑2𝑦(𝑥)𝑑𝑥2+𝜈𝜆+214sec2𝑥𝑦(𝑥)=0,(1.3) on (𝜋/2,𝜋/2), which leads, when 𝜈=𝑛, to an expansion in series involving associated Legendre functions.

It is perhaps worth remarking how our interest in this problem arises. In [3] there is the question of travelling waves and steady solutions for a discrete reaction-diffusion equation of the type 𝑢𝑛=𝑢𝑛+12𝑢𝑛+𝑢𝑛1+𝑓(𝑢𝑛),(1.4) where the function 𝑓 is “bistable”. That is, there exist three numbers 𝑈1, 𝑈2, 𝑈3, 𝑈1<𝑈2<𝑈3, such that 𝑓(𝑈1)=𝑓(𝑈2)=𝑓(𝑈3)=0(1.5) with 𝑓<0 in (𝑈1,𝑈2) and 𝑓>0 in (𝑈2,𝑈3). A prototypical 𝑓 would be 𝑓(𝑢)=𝐴sin(𝑢)+𝐹,(1.6) where 𝐴 and 𝐹 are positive constants and 𝐹<𝐴, so that 𝑈1=sin1𝐹𝐴,𝑈2=𝜋+sin1𝐹𝐴,𝑈3=2𝜋+sin1𝐹𝐴.(1.7) Such equations arise in a number of different applications, for example, in dislocation theory where 𝑢𝑛 is the displacement of the 𝑛th atom in some material, or in neurobiology where 𝑢𝑛 is typically the electric potential of the 𝑛th nerve cell, and in both these applications the interest is in monotonic solutions 𝑢𝑛 with 𝑢𝑛𝑈1 as 𝑛, 𝑢𝑛𝑈3 as 𝑛 [36].

The basic question is whether there exist such solutions with the form of a travelling wave, 𝑢(𝑛𝑐𝑡), 𝑐0, or of a steady solution or standing wave, where 𝑐=0, and there is an important distinction between these two cases. For a travelling wave, 𝑐0, 𝑢𝑛 is clearly a function of the continuous variable 𝑡, and indeed because of (1.4), a differentiable function of 𝑡. This leads to the difference-differential equation 𝑢𝑡(𝑥,𝑡)=𝑢(𝑥+1,𝑡)2𝑢(𝑥,𝑡)+𝑢(𝑥1,𝑡)+𝑓𝑢(𝑥,𝑡).(1.8) If, however, 𝑐=0, then, as in [3], we have to study the purely difference equation 𝑢(𝑥+1,𝑡)2𝑢(𝑥,𝑡)+𝑢(𝑥1,𝑡)+𝑓𝑢(𝑥,𝑡)=0,(1.9) and the solutions may be discontinuous since there is nothing that now connects values of 𝑢(𝑥) with values of 𝑢(𝑥+𝛿) for |𝛿|<1. It is best therefore to think of the solution of (1.9) as a number of (monotonic) sequences 𝑢𝛼𝑛 indexed by 𝛼, each satisfying 𝑢𝑛+12𝑢𝑛+𝑢𝑛1+𝑓(𝑢𝑛)=0.(1.10) The simplest case would be that there is just one such sequence (modulo the translation 𝑛𝑛+𝑘, 𝑘 integral), but it is possible that there may be a finite number, or even a partial or total continuum.

In view of applications, where the distance between atoms or nerve cells is small, it is more natural to think of (1.8) in the form 𝑢𝑡(𝑥,𝑡)=𝑢(𝑥+𝜖,𝑡)2𝑢(𝑥,𝑡)+𝑢(𝑥𝜖,𝑡)+𝑓𝑢(𝑥,𝑡),(1.11) where 𝜖 is small and represents the distance between atoms or nerve cells. A tempting approximation is then 𝑢𝑡(𝑥,𝑡)𝜖2𝑢𝑥𝑥,(𝑥,𝑡)+𝑓𝑢(𝑥,𝑡)(1.12) and in order to make sense of the scaling, in [3] the authors introduced a factor 𝜖2 in front of 𝑓. This therefore leads to a comparison between the solutions of 𝑢𝑛=𝑢𝑛+12𝑢𝑛+𝑢𝑛1+𝜖2𝑓(𝑢𝑛𝑢)=0,(1.13)𝑡(𝑥,𝑡)=𝑢𝑥𝑥(𝑥,𝑡)+𝑓(𝑢).(1.14) For the continuous diffusion problem, the answer is both simple and well known [4, 5].

Given a function 𝑓 that is bistable, there is just one possible wave-speed 𝑐, and this value of 𝑐 is 0, that is, there is a steady solution if and only if 𝑈3𝑈1𝑓(𝑢)𝑑𝑢=0.(1.15) (The proof is a simple phase plane argument, and 𝑐=0 implies (1.15) follows by multiplying (1.14) by 𝑢 and integrating.)

The solution in the discrete case is however different, as discussed in [6]. There may continue to be steady solutions where (1.15) no longer holds. Consider specifically the case (1.6), so that 𝑢𝑛=𝑢𝑛+12𝑢𝑛+𝑢𝑛1𝜖2sin𝑢𝑛+𝐹.(1.16) The case corresponding to (1.15) is 𝐹=0, but the authors, in [3], have shown that for 𝐹 sufficiently small, say |𝐹|<𝐹crit, there exist precisely two steady solutions of (1.16), and 𝐹crit, which of course depends on 𝜖, can be evaluated for small 𝜖. Specifically, 𝐹crit𝐵𝑒𝜋2/𝜖,(1.17) where the constant 𝐵=64𝜋𝜋0(sin2(𝑠)/𝑠)𝑑𝑠 is given. For |𝐹|>𝐹crit, the solutions move and the equation has travelling wave solutions instead of steady solutions.

In order to prove results such as (1.17), one has to regard (1.16) as a singular perturbation of the steady continuous-diffusion equation 𝑢𝑥𝑥sin𝑢=0,(1.18) for which the solution (satisfying 𝑢()=0, 𝑢()=2𝜋) is 𝑈=4tan1𝑒𝑥. If we linearize (1.18) about 𝑈, we obtain 𝜙cos(𝑈)𝜙=0.(1.19) But multiplying (1.18) by 𝑈 and integrating lead to (1/2)(𝑈)2=1cos𝑈, so that since 𝑈=2sech(𝑥), we have cos𝑈=12sech2(𝑥). The linearization (1.19) thus becomes 𝜙+(1+2sech2𝑥)𝜙=0, which is of course (1.2) with 𝜆=1 and 𝑛=1. Thus the selfadjoint operator 𝑇 given (in 𝐿2(,)) by 𝑇𝜙=𝜙2sech2(𝑥)𝜙 has an eigenvalue at 1, with eigenfunction 𝑈 (differentiation of (1.18) shows that 𝑈 satisfies (1.19)). This fact, together with the additional fact that the spectrum of 𝑇 is continuous above 0 (since sech2(𝑥)𝐿(0,)), is highly relevant to the work in [3] and led to our interest more generally in the spectral problem (1.2).

The explicit solution for any 𝑛 using contour integrals different from what Kamke did is known to [7]. For more information on this problem one can see [7, 8].

2. Preliminaries

We want to know expansion of an arbitrary function 𝑓(𝑥) in terms of eigenfunctions. So one needs to know the following. Let 𝜃(𝑥,𝜆) and 𝜙(𝑥,𝜆) be the solutions of (1.1) such that 𝜙(0,𝜆)=sin𝛼,𝜙(0,𝜆)=cos𝛼,𝜃(0,𝜆)=cos𝛼,𝜃(0,𝜆)=sin𝛼,(2.1) where 𝛼 is real. 𝑊𝑥(𝜙,𝜃)=𝑊0(𝜙,𝜃)=1. The general solution of (1.1) is of the form 𝜓(𝑥,𝜆)=𝜃(𝑥,𝜆)+𝑚(𝜆)𝜙(𝑥,𝜆)𝐿2(0,).(2.2) The spectrum is defined by means of the function 𝑘(𝜆)=lim𝛿0𝜆0{𝑚(𝑢+𝑖𝛿)}𝑑𝑢,(2.3) which exists for all real 𝜆 and 𝑘(𝜆) is a nondecreasing function. The expansion of a function 𝑓(𝑥) in terms of the spectral function depends on the following lemmas taken from [1].

Lemma 2.1. Without detailing, let the interval be (0,): 1𝑓(𝑥)=𝜋0𝜙(𝑥,𝜆)𝑑𝑘(𝜆)0𝜙(𝑡,𝜆)𝑓(𝑡)𝑑𝑡.(2.4) If 𝑚(𝜆) has poles, then 𝑓(𝑥)=𝑁=0𝜙𝑁(𝑡)0𝜙𝑁(𝑡,𝜆)𝑓(𝑡)𝑑𝑡.(2.5)

Lemma 2.2. Without detailing, let the interval be (,). If 𝑞(𝑥) is an even function, then 𝑚1(𝜆)=𝑚2(𝜆). So the expansion formula is 1𝑓(𝑥)=𝜋𝜃(𝑥,𝜆)𝑑𝜉1𝜃(𝑦,𝜆)𝑓(𝑦)𝑑𝑦+𝜋𝜙(𝑥,𝜆)𝑑𝜁𝜙(𝑦,𝜆)𝑓(𝑦)𝑑𝑦,(2.6) where 𝜉(𝜆)=lim𝛿0𝜆01𝑚1(𝑢+𝑖𝛿)𝑚2(𝑢+𝑖𝛿)𝑑𝑢,𝜉1(𝜆)=2𝑚2,(𝜆)𝜁(𝜆)=lim𝛿0𝜆0𝑚1(𝑢+𝑖𝛿)𝑚2(𝑢+𝑖𝛿)𝑚1(𝑢+𝑖𝛿)𝑚2(𝑢+𝑖𝛿)𝑑𝑢,𝜁1(𝜆)=2𝑚2.(𝜆)(2.7)

3. Main Results

We are now dealing with (1.2) in the case where 𝑛 is integral. Without loss of generality, we may suppose 𝑛0, but since 𝑛=0 reduces (1.2) to the simple trigonometric case, we are in fact interested only in 𝑛>0. We first prove that a solution is given by 𝑦=cosh𝑛+1𝑥𝐶cosh(𝑧𝑠)(sinh𝑧sinh𝑥)𝑛+1𝑑𝑧,𝑠2=𝜆,(3.1) where the contour 𝐶 is taken round the point 𝑧=𝑥 and no other zero of sinh𝑧sinh𝑥. This is slight variant of a form which Titchmarsh uses in his discussion of (1.3). The proof below will show (3.1), being continuous at least formally, to be a solution of (1.2) where 𝑛 is not an integer, but the difficulty then is to choose a suitable contour, since the integrand has a branch point at 𝑧=𝑥.

Remark 3.1. We also remark that it is obvious that we can express the solution (3.1) equivalently ignoring some multiplicative constants as 𝑦(𝑥)=cosh𝑛+1𝑥𝑑𝑛(cosh𝑥𝑑𝑥)𝑛𝐶cosh(𝑧𝑠)sinh𝑧sinh𝑥𝑑𝑧.(3.2)

Theorem 3.2. The contour integral (3.1) satisfies the differential equation (1.2).

Proof. We see that 𝑦=(𝑛+1)tanh(𝑥)𝑦+(𝑛+1)cosh𝑛+2𝑥𝐶cosh(𝑧𝑠)(sinh𝑧sinh𝑥)𝑛+2𝑦𝑑𝑧,(3.3)+𝑛(𝑛+1)sech2(𝑥)𝑦=(𝑛+1)cosh𝑛+1𝑥𝐶cosh(𝑧𝑠)(𝑛+1)sinh2𝑧+sinh𝑧sinh𝑥+𝑛+2(sinh𝑧sinh𝑥)𝑛+3𝑑𝑧.(3.4) Integrating (3.1) by parts, 𝑦(𝑥)=𝑛+1𝑠2𝐶cosh(𝑧𝑠)(𝑛+1)sinh2𝑧+sinh𝑧sinh𝑥+𝑛+2(sinh𝑧sinh𝑥)𝑛+3𝑑𝑧,(3.5) so that 𝜆𝑦=(𝑛+1)cosh𝑛+1𝑥𝐶cosh(𝑧𝑠)(𝑛+1)sinh2𝑧+sinh𝑧sinh𝑥+𝑛+2(sinh𝑧sinh𝑥)𝑛+3𝑑𝑧.(3.6) Comparing (3.4) and (3.6), we see that 𝑦+𝑛(𝑛+1)sech2(𝑥)𝑦=𝜆𝑦, so that 𝑦(𝑥) satisfies (1.2), as required.

Remark 3.3. We now point out that the factor cosh(𝑧𝑠) played little part in the argument. Certainly, the argument would have washed equally well if we had replaced cosh(𝑧𝑠) by sinh(𝑧𝑠): 𝑦2(𝑥)=cosh𝑛+1𝑥𝐶sinh(𝑧𝑠)(sinh𝑧sinh𝑥)𝑛+1𝑑𝑧.(3.7)

Theorem 3.4. The contour integral (3.7) satisfies the differential equation (1.2).

Proof. Proof is the same as the above theorem. So we omit it.

Remark 3.5. Furthermore, once the integrands have poles at 𝑧=𝑥, the solution can be evaluated by calculating the relevant residues. For example, in the trivial case 𝑛=0, when we should recover the trigonometric functions, the residues of cosh(𝑧𝑠)sinh𝑧sinh𝑥(3.8) are cosh(𝑥𝑠),cosh𝑥(3.9) so that the solution (3.7) becomes multiples of cos(𝑥𝜆) (similarly sin(𝑥𝜆)), as we expect.
We can generalize Theorems (3.2) and (3.4) by defining the following operator:𝑇𝑓(𝑥)=cosh𝑛+1𝑥𝐶𝑓(𝑧)(sinh𝑧sinh𝑥)𝑛+1𝑑𝑧,(3.10)where 𝑓 is a differentiable function as long as one can pick up residue.

Corollary 3.6. If 𝑓(𝑧)=cosh(𝑧𝑠)(sinh(𝑧𝑠)), then we obtain Theorems (3.2) and (3.4). The operator 𝑇 is also linear.

4. The Explicit Solution Given by Residues for 𝑛=1

We now require the residues ofcosh(𝑧𝑠)(sinh𝑧sinh𝑥)2,sinh(𝑧𝑠)(sinh𝑧sinh𝑥)2.(4.1) Since cosh(𝑧𝑠)(sinh𝑧sinh𝑥)2={cosh(𝑥𝑠)+(𝑧𝑥)𝑠sinh(𝑥𝑠)+}{1(𝑧𝑥)tanh(𝑥)+}(𝑧𝑥)2cosh2,(𝑥)(4.2) we see that the residue at 𝑧=𝑥 is𝑠sinh(𝑥𝑠)tanh(𝑥)cosh(𝑥𝑠)cosh2,(𝑥)(4.3) so that one solution is𝑦1=𝜆sin(𝑥𝜆)+tanh(𝑥)cos(𝑥𝜆).(4.4) By examining the residue of the second equation of (4.1), we see that a second solution is𝑦2=𝜆cos(𝑥𝜆)tanh(𝑥)sin(𝑥𝜆).(4.5)

Remark 4.1. The solution can also be obtained from (3.2). For we have already seen, from our brief discussion of the case 𝑛=0, that the integral in (3.2) is just a multiple of cos(𝑥𝜆)/cosh𝑥 (or of sin(𝑥𝜆)/cosh𝑥 if we replace cosh(𝑧𝑠) by sinh(𝑧𝑠)) hence in the first case, (3.2) gives a multiple of 𝑑cosh(𝑥)𝑑𝑥cos(𝑥𝜆)cosh(𝑥)=𝜆sin(𝑥𝜆)tanh(𝑥)cos(𝑥𝜆),(4.6)in accordance with (4.4).

Remark 4.2. Wronskian 𝑊(𝑦1(𝑥),𝑦2(𝑥))=𝜆(𝜆+1). We now have two linearly independent solutions.

Remark 4.3. The general solution is 𝑦(𝑥)=𝑐1𝑦1(𝑥)+𝑐2𝑦2(𝑥).

Lemma 4.4. 𝑦1(𝑥) is an odd function but 𝑦2(𝑥) is an even function.

Proof. 𝑦1(𝑥)=𝑦1(𝑥),𝑦2(𝑥)=𝑦2(𝑥).(4.7)

5. Eigenvalues and Eigenfunctions for 𝑛=1 When 𝑦(0)=𝑦(𝑏)=0

Theorem 5.1. Eigenvalues with associated boundary conditions 𝑦(0)=𝑦(𝑏)=0 are the zeros of 𝜆tan(𝑏𝜆)+tanh(𝑏)=0;(5.1) furthermore, one and only one eigenvalue lies in the interval 1𝑘2𝜋<𝑏1𝜆<𝑘+2𝜋(5.2) for every integral 𝑘0.

Proof. One can see from Remark 4.3 that if 𝑦(0)=0, then 𝑐2=0. So that 𝑦(𝑏)=0 implies 𝜆sin(𝑏𝜆)+tanh(𝑏)cos(𝑏𝜆)=0,𝜆tan(𝑏𝜆)=tanh(𝑏).(5.3) One can see immediately that the eigenvalues belong to the interval (5.2).
To prove the second part we use the following strategy. Multiply (5.1) by 𝑏 and set 𝑥=𝑏𝜆, then denote (𝑥)=𝑥tan(𝑥)+𝑏tanh(𝑏). So (𝑥)=sec2(𝑥){(1/2)sin(2𝑥)+𝑥}. Notice that (𝑥) is an even and does not intersect 𝑥-axis, where 1/2<𝑥<1/2. If 𝑥<1/2, then (𝑥)<0 and if 𝑥>1/2, then (𝑥)>0. So the monotonicity of (𝑥) implies that (𝑥) has only one zero belonging to the interval (5.2) for every integral 𝑘0.

Remark 5.2. So 𝑦𝑘 is the eigenfunction in the form of (4.4): 𝑦𝑘=𝜆𝑘sin(𝑥𝜆𝑘)+tanh(𝑥)cos(𝑥𝜆𝑘).(5.4)

Corollary 5.3. One can orthonormalize the eigenfunctions.

Proof. 𝑏0𝑦2𝑘𝑑𝑥=4𝜆𝑘tanh(𝑏)cos2(𝑏𝜆𝑘)+2𝜆𝑘(𝜆𝑘+1)𝑏(𝜆𝑘1)sin(2𝑏𝜆𝑘)4𝜆𝑘.(5.5) The orthonormalized eigenfunctions denoted by Φ(𝑥,𝜆𝑘), Φ(𝑥,𝜆𝑘)=(𝜆𝑘sin(𝑥𝜆𝑘)+tanh(𝑥)cos(𝑥𝜆𝑘))/𝑏0𝑦2𝑘𝑑𝑥.

Remark 5.4. An arbitrary function 𝑓(𝑥) in terms of eigenfunctions follows: 𝑓(𝑥)=𝑘=0𝐶𝑘Φ(𝑥,𝜆𝑘),where𝐶𝑘=𝑏0𝑓(𝑥)Φ(𝑥,𝜆𝑘)𝑑𝑥.(5.6)

6. Eigenvalues and Eigenfunctions for 𝑛=1 When 𝑦(0)=𝑦(𝑏)=0

Theorem 6.1. The eigenvalues with 𝑦(0)=𝑦(𝑏)=0 are the zeros of (6.1); furthermore, there exists one and only one eigenvalue lying in the interval (5.2) for every integral 𝑘0.

Proof. One can see from Remark 4.3 that if 𝑦(0)=0, then either 𝑐1=0 or 𝜆=1. If 𝜆=1, then the associated eigenfunction is zero. So this is useless. Hence, 𝑐1=0. 𝑦(𝑏)=0 implies 𝜆tan(𝑏𝜆)+sech2(𝑏)tan(𝑏𝜆)+𝜆tanh(𝑏)=0.(6.1)So it is obvious that the zeros (eigenvalues) belong to the interval (5.2).
Set 𝑥=𝑏𝜆. (𝑏2𝜆=𝑥2𝜆=𝑥2/𝑏2). Equation (6.1) is denoted by (𝑥): 𝑥(𝑥)=2𝑏2+sech2𝑥(𝑏)tan(𝑥)+𝑏tanh(𝑏).(6.2) It is enough to show that (𝑥) is monotonic: 𝑥(𝑥)=2+𝑏2sech2(𝑏)+𝑥sin(2𝑥)𝑏2cos2+1(𝑥)𝑏tanh(𝑏).(6.3) We see that (𝑥)>0 everywhere. We therefore conclude that (𝑥) is monotonic.

Remark 6.2. The associated eigenfunctions are 𝑦𝑘=𝜆𝑘cos(𝑥𝜆𝑘)tanh(𝑥)sin(𝑥𝜆𝑘).(6.4)

Corollary 6.3. One can orthonormalize the eigenfunctions.

Proof. 𝑏0𝑦2𝑘𝑑𝑥=4𝜆𝑘tanh(𝑏)sin2(𝑏𝜆𝑘)+2𝜆𝑘(𝜆𝑘+1)𝑏+(𝜆𝑘1)sin(2𝑏𝜆𝑘)4𝜆𝑘.(6.5) The orthonormalized eigenfunctions denoted by Ψ(𝑥,𝜆𝑘), Ψ(𝑥,𝜆𝑘)=(𝜆𝑘cos(𝑥𝜆𝑘)tanh(𝑥)sin(𝑥𝜆𝑘))/𝑏0𝑦2𝑘𝑑𝑥.

Remark 6.4. Therefore, an arbitrary 𝑓(𝑥) in terms of orthonormalized eigenfunctions is 𝑓(𝑥)=𝑘=0𝐶𝑘Ψ(𝑥,𝜆𝑘),𝐶𝑘=𝑏0𝑓(𝑥)Ψ(𝑥,𝜆𝑘)𝑑𝑥.(6.6)

7. Spectral Function 𝑚(𝜆) over (0,) and Expansion

Now let 𝜃(𝑥,𝜆) and 𝜙(𝑥,𝜆) be the solutions of 𝑑2𝑦𝑑𝑥2+𝜆+2sech2(𝑥)𝑦=0,(7.1) which satisfy (2.1); so that 𝜙(𝑥,𝜆)=cos𝛼{𝜆sin(𝑥𝜆)+tanh𝑥cos(𝑥𝜆)}+𝜆+1sin𝛼{𝜆cos(𝑥𝜆)tanh𝑥sin(𝑥𝜆)}𝜆,(7.2)𝜃(𝑥,𝜆)=sin𝛼{𝜆sin(𝑥𝜆)+tanh𝑥cos(𝑥𝜆)}+𝜆+1cos𝛼{𝜆cos(𝑥𝜆)tanh𝑥sin(𝑥𝜆)}𝜆.(7.3) Now we need to find 𝑚(𝜆). This suggests that 𝜓(𝑥,𝜆)=𝜃(𝑥,𝜆)+𝑚(𝜆)𝜙(𝑥,𝜆)𝐿2(0,).(7.4) To get demanded result, one needs to find the asymptotics of 𝜃(𝑥,𝜆) and 𝜙(𝑥,𝜆) as 𝑥 and 𝜆>0: 𝜙(𝑥,𝜆){𝜆cos𝛼+𝜆(𝜆+1)sin𝛼+𝑖(𝜆cos𝛼(𝜆+1)sin𝛼)}𝑒𝑖𝑥𝜆2𝜆(𝜆+1)=𝑀1(𝜆)𝑒𝑖𝑥𝜆,{𝜃(𝑥,𝜆)𝜆sin𝛼+𝜆(𝜆+1)cos𝛼+𝑖(𝜆sin𝛼(𝜆+1)cos𝛼)}𝑒𝑖𝑥𝜆2𝜆(𝜆+1)=𝑀(𝜆)𝑒𝑖𝑥𝜆.(7.5) Finally, we must arrange the linear combination so that the terms 𝑒𝑖𝑥𝜆 cancel. That is, 𝑀(𝜆)+𝑚(𝜆)𝑀1(𝜆)=0. Hence, 𝑚(𝜆)=𝑖𝜆(𝜆+1)(𝜆2+𝜆+1)sin𝛼cos𝛼𝜆cos2𝛼+(𝜆+1)2sin2𝛼,𝑚(𝜆)=𝜆(𝜆+1)𝜆cos2𝛼+(𝜆+1)2sin2𝛼when𝜆>0,0when𝜆<0.(7.6) We see that the spectrum is continuous for 𝜆>0. But we have a point spectrum for 𝜆<0. So that an arbitrary 𝑓(𝑥) is a linear combination of integrand and series. The spectral function is calculated from (2.3). Hence, 𝑑𝑘(𝜆)=𝜆(𝜆+1)𝜆cos2𝛼+(𝜆+1)2sin2𝛼when𝜆>0,0when𝜆<0.(7.7)(1) In particular, if 𝛼=0, then𝑑𝑘(𝜆)=𝜆+1𝜆when𝜆>0,0when𝜆<0.(7.8) So that in this case there is no eigenvalue for 𝜆<0. From Lemma 2.1, one can see the expansion of1𝑓(𝑥)=𝜋0𝜆sin(𝑥𝜆)+tanh𝑥cos(𝑥𝜆)𝜆𝑑𝜆0𝜆sin(𝑦𝜆)+tanh𝑦cos(𝑦𝜆)𝜆+1𝑓(𝑦)𝑑𝑦.(7.9) Similarly, if 𝛼=𝜋/2, then 𝑚(𝜆)=𝜆𝜆+1when𝜆>0,0when𝜆<0..(7.10)(2) Hence, there exits only one eigenvalue at 𝜆=1 and the corresponding eigenfunction 𝜙(𝑥,1) where 𝜙(𝑥,𝜆) is (7.2). So the spectrum calculated by (2.3) is𝑑𝑘(𝜆)=𝜆𝜆+1when𝜆>0,0when𝜆<0..(7.11) Therefore, from Lemma 2.1, an expansion of function 𝑓(𝑥) in terms of eigenfunctions and spectral function follows:1𝑓(𝑥)=𝜋0𝜆cos(𝑥𝜆)tanh𝑥sin(𝑥𝜆)×𝜆+1𝑑𝜆0𝜆cos(𝑦𝜆)tanh𝑦sin(𝑦𝜆)𝜆𝑓(𝑦)𝑑𝑦+𝑐1sech(𝑥),(7.12) where 𝑐1 is a constant. So we have proved the following theorem regarding the nature of 𝑚(𝜆).

Theorem 7.1. If 𝛼=0, there is no eigenvalue where 𝜆<0. If 𝛼=𝜋/2, then there exists only one eigenvalue at 𝜆=1 and the corresponding eigenfunction is 𝜙(𝑥,1)=𝑠𝑒𝑐(𝑥).

Finally, one can ask what is the range of 𝛼 in the case of 𝜆<0 and 𝑛=1.

Theorem 7.2. If 0<𝛼<𝜋/2, then there are precisely two eigenvalues except at 𝛼=𝜋/4. Hence there are two eigenfunctions, namely, 𝜙(𝑥,𝜆1) and 𝜙(𝑥,𝜆2).

Proof. We check the zeros of both the numerator and denominator of (7.6). After working out the algebra, we see that the zeros of the numerator and the denominator are 𝜇1=(2+tan2𝛼)+tan𝛼4+tan2𝛼+12,𝜇2=(2+tan2𝛼)tan𝛼4+tan2𝛼+12,𝜆1=(2tan2𝛼+1)+4tan2𝛼+12tan2𝛼,𝜆2=(2tan2𝛼+1)4tan2𝛼+12tan2𝛼.(7.13) If 𝑚(𝜆) has poles, then they are the eigenvalues. If so, the eigenfunctions are 𝜙(𝑥,𝜆1) and 𝜙(𝑥,𝜆2). Finally, there is one thing to be proved in this case 𝛼=𝜋/4. To do this, we expand the zeros around 𝛼=𝜋/4 and at the end 𝛼𝜋/4: 𝜆1=3+52+20125𝜋10𝛼4+400+3525𝜋100𝛼42𝜋+𝑂𝛼43,𝜆2=352+20+125𝜋10𝛼4+4003525𝜋100𝛼42𝜋+𝑂𝛼43.(7.14) Similarly, 𝜇1=3+62+76126𝜋𝛼4+288+1186𝜋36𝛼42𝜋+𝑂𝛼43,𝜇2=362+76126𝜋𝛼4+2881186𝜋36𝛼42𝜋+𝑂𝛼43.(7.15) That completes the proof.

8. Expansion over (,)

Now consider the interval (,) instead of (0,). From (2.7) the following can be calculated: 𝜉(𝜆)=𝜆(𝜆+1)2𝜆sin2𝛼+2(𝜆+1)2cos2𝛼𝜁when𝜆>0,0when𝜆<0,(𝜆)=(𝜆+1)𝜆2𝜆cos2𝛼+2(𝜆+1)2sin2𝛼when𝜆>0,0when𝜆<0.(8.1) Now one can use (7.2), (7.3), and (8.1) gets the following expansion from Lemma (2.2): 1𝑓(𝑥)=𝜋0𝜆cos(𝑥𝜆)tanh𝑥sin(𝑥𝜆)2(𝜆+1)𝑑𝜆𝜆cos(𝑦𝜆)tanh𝑦sin(𝑦𝜆)𝜆+𝑓(𝑦)𝑑𝑦0𝜆sin(𝑥𝜆)+tanh𝑥cos(𝑥𝜆)2𝜆𝑑𝜆𝜆sin(𝑦𝜆)+tanh𝑦cos(𝑦𝜆)𝜆+1𝑓(𝑦)𝑑𝑦+𝑐1sech(𝑥),(8.2) where 𝑐1 is a constant.