Preprint

Article

On the Asymptotic of Solutions of Odd Order Two-Term Differential Equations

Altmetrics

Downloads

Views

Comments

A peer-reviewed article of this preprint also exists.

Yaudat T. Sultanaev,Nur F. Valeev,

Elvira A. Nazirova^*

Yaudat T. Sultanaev,Nur F. Valeev,

Elvira A. Nazirova^*

This version is not peer-reviewed

Submitted:

30 November 2023

Posted:

01 December 2023

You are already at the latest version

Alerts

Abstract

The work is devoted to the development of methods for constructing asymptotic formulas as x→∞ of a fundamental system of solutions of linear differential equations generated by a symmetric two-term differential expression of odd order. The coefficients of the differential expression belong to classes of functions that allow oscillation (for example, those that do not satisfy the classical Titchmarsh-Levitan regularity conditions). As a model equation, the 5th order equation i2p(x)y′′′′′+p(x)y′′′′′+q(x)y=λy, for which various cases of behavior of the coefficients p(x),q(x), is investigated. New asymptotic formulas are obtained for the case when the function h(x)=−1+p−1/2(x)∉L1[1,∞) significantly influences the asymptotics of solutions to the equation. The case when the equation contains a nontrivial bifurcation parameter is studied.

Keywords:

Subject: Computer Science and Mathematics - Mathematics

1. Introduction

The asymptotic behavior as

x \to \infty

of fundamental system of solutions of arbitrary order singular differential equations, in addition to being of independent interest, is an effective method for studying of qualitative spectral characteristics for corresponding differential operators ([1,2,3], etc.). As a rule, in these books differential equations with regular coefficients with regular growth at infinity are investigated. Therefore, the study of the asymptotic behavior of solutions to equations with coefficients from other classes of functions is of particular interest. Such classes of functions were described by us in the papers [4]. Let us also note the works [5,6,7], where differential operators with distribution coefficients were studied.

For example, in the work [7] asymptotic formulas were obtained for the fundamental system of solutions of a to-term equation of even order

{(- 1)}^{n} {(p (x) y^{(n)})}^{(n)} + q (x) y = λ y, x \in [1, \infty),

where the locally summable function p can be represented

p (x) = {(1 + r (x))}^{- 1},

r \in L_{1} [1, \infty),

and q is a generalized function representable for some fixed k,

0 \leq k \leq n

in the form

q = σ^{(k)}

, where

σ \in L_{1} [1, \infty)

k < n

| σ | (1 + | r |) (1 + | σ |) \in L_{1} [1, \infty),

k = n

Since 2014 we have published a series of articles devoted to the study of the asymptotic behavior of solutions to singular ODEs with regularly oscillating coefficients ([4,8,9,10,11,12,13,14]). In this case, a new approach was used for the study, based on a sequence of matrix transformations and the use of Campbell’s identity [16].

The use of this approach made it possible to obtain new asymptotic formulas in different cases. For example, in [4,12,14] new asymptotic formulas were obtained for solutions of the Sturm-Liouville equation

y^{″} + (μ^{2} + \frac{s i n (x^{β})}{x^{α}}) y = 0, 0 < α \leq 1, β > \frac{α}{2} + 1

under some relations between

α, β, μ

. Note that

μ

has the meaning of a bifurcation parameter. By the way this equation is one of equations for testing of new methods for constructing asymptotic formulas (see, for example, [15], p. 160).

Equations of odd order for irregular classes of coefficients (in the Titchmarsh-Levitan sense) have been less studied. In the works [7,11,13] the asymptotics of solutions of odd order equations was studied in the case when the coefficient of the highest derivative is equal or equivalent to unity.

Here we develop an approach that was proposed in [4,8,9,10,11,12,13,14] and can be implemented to study the asymptotic behavior as

x \to \infty

of a fundamental system of solutions of two-term equation of arbitrary odd order of the form

l y = \frac{i}{2} [{(p (x) y^{(n)})}^{(n + 1)} + {(p (x) y^{(n + 1)})}^{(n)}] + q (x) y = λ y, x \geq 1

(1)

for various cases of behavior of the coefficients

p (x), q (x)

2. Transition to the ODE system using quasi-derivatives

Let us write equation (1) in the form of a system of ordinary differential equations of the first order. To do this, we will use the apparatus of quasi-derivatives (for more details see [17]). Let us define the functions

q_{n} (x) \in L_{1, l o c} [1, \infty)

so that

q_{n}^{(n)} (x) = q (x)

and introduce into consideration quasi-derivatives defined by the following formulas

\{\begin{matrix} z_{1} = y, z_{n + 2} = \sqrt{p} z_{n + 1}^{'} - i q_{n} z_{1} \\ z_{2} = z_{1}^{'}, z_{n + 3} = z_{n + 1}^{'} + i C_{n}^{1} q_{n} z_{2} \\ . . . . . . . . \\ z_{n} = z_{n - 1}^{'}, z_{2 n} = z_{2 n - 1}^{'} + i {(- 1)}^{n - 1} C_{n}^{n - 2} q_{n} z_{n - 1} \\ z_{n + 1} = \sqrt{p} z_{n}^{'}, z_{2 n + 1} = z_{2 n}^{'} + i {(- 1)}^{n} C_{n}^{n - 1} q_{n} z_{n} . \end{matrix}

(2)

Then equation (1) is equivalent to the relation

z_{2 n + 1}^{'} = λ z_{1} - i {(- 1)}^{n + 1} \frac{q_{n}}{\sqrt{p}} z_{n + 1} .

Let us introduce the column vector

z = c o l u m n (z_{1}, z_{2}, . . ., z_{2 n + 1})

and write equation (1) as a system of ODEs

z^{'} = S z,

where

S (x, λ)

is the Shin-Zettl matrix [17].

S = (\begin{matrix} 0 & 1 & 0 & . . . & 0 & 0 & 0 & 0 & . . . & 0 \\ 0 & 0 & 1 & . . . & 0 & 0 & 0 & 0 & . . . & 0 \\ . . . . . . . . . . . . . . \\ 0 & 0 & 0 & . . . & 1 & 0 & 0 & 0 & . . . & 0 \\ 0 & 0 & 0 & . . . & 0 & \frac{1}{\sqrt{p}} & 0 & 0 & . . . & 0 \\ \frac{i q_{n}}{\sqrt{p}} & 0 & 0 & . . . & 0 & 0 & \frac{1}{\sqrt{p}} & 0 & . . . & 0 \\ 0 & - i n q_{n} & 0 & . . . & 0 & 0 & 0 & 1 & . . . & 0 \\ . . . . . . . . . . . . . . . \\ 0 & 0 & 0 & . . . & 0 & {(- 1)}^{n - 1} i n q_{n} & 0 & 0 & . . . & 1 \\ - i λ & 0 & 0 & . . . & 0 & 0 & \frac{{(- 1)}^{n} i q_{n}}{\sqrt{p}} & 0 & . . . & 0 \end{matrix}),

here the non-zero elements of the matrix

S (x, λ)

are given by the formulas

s_{k j} = 1, j = 1 + k, k = \bar{1, n - 1}, k = \bar{n + 2, 2 n}, s_{n, n + 1} = s_{n + 1, n + 2} = \frac{1}{\sqrt{p}},

s_{n + 1, 1} = \frac{i q_{n}}{\sqrt{p}}, s_{n + k, k} = {(- 1)}^{k - 1} i C_{n}^{k - 1} q_{n}, k = \bar{2, n},

s_{2 n + 1, 2 n - 1} = \frac{{(- 1)}^{n} i q_{n}}{\sqrt{p}}, s_{2 n + 1, 1} = - i λ .

Note that from the relation

q_{n}^{(n)} (x) = q (x)

the function

q_{n} (x)

is determined up to a polynomial of order

n - 1

. However, the fundamental system of solutions of the equation (1) does not depend on the choice of integration constants, which follows directly from the formulas (2). Conditions for choosing the coefficients of the polynomial will be formulated for each case under study.

Further, in order to avoid complicated formulas, we will limit ourselves to considering the 5th order two-term equation

l y = \frac{i}{2} [{(p (x) y^{″})}^{‴} + {(p (x) y^{‴})}^{″}] + q (x) y = λ y, x \geq 1 .

(3)

Using formulas (2), we introduce quasi-derivatives

\{\begin{matrix} z_{1} = y \\ z_{2} = z_{1}^{'} \\ z_{3} = \sqrt{p} z_{2}^{'} \\ z_{4} = \sqrt{p} z_{3}^{'} - i q_{2} z_{1} \\ z_{5} = z_{4}^{'} + 2 i q_{2} z_{2}, \end{matrix}

then equation (3) is equivalent to the relation

z_{5}^{'} = - i λ z_{1} + \frac{i q_{2}}{\sqrt{p}} z_{3}

and can be written as a system of ordinary differential equations

z^{'} = (\begin{matrix} 0 & 1 & 0 & 0 & 0 \\ 0 & 0 & 1 / \sqrt{p} & 0 & 0 \\ i q_{2} / \sqrt{p} & 0 & 0 & 1 / \sqrt{p} & 0 \\ 0 & - 2 i q_{2} & 0 & 0 & 1 \\ - i λ & 0 & i q_{2} / \sqrt{p} & 0 & 0 \end{matrix}) z,

where

z = c o l u m n (z_{1}, z_{2}, z_{3}, z_{4}, z_{5}) .

Let the function

p (x)

admit the representation

\frac{1}{\sqrt{p (x)}} = 1 + h (x), L_{1, l o c} [1, \infty) .

Let us write the last system of equations in the form

z^{'} = (L_{0} + h (x) L_{1} + i q_{2} (x) D_{0} + i h (x) q_{2} (x)) D_{1}) z,

(4)

L_{0} = (\begin{matrix} 0 & 1 & 0 & 0 & 0 \\ 0 & 0 & 1 & 0 & 0 \\ 0 & 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 0 & 1 \\ - i λ & 0 & 0 & 0 & 0 \end{matrix}), L_{1} = (\begin{matrix} 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & 1 & 0 & 0 \\ 0 & 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 \end{matrix}),

D_{0} = (\begin{matrix} 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 \\ 1 & 0 & 0 & 0 & 0 \\ 0 & - 2 & 0 & 0 & 0 \\ 0 & 0 & 1 & 0 & 0 \end{matrix}), D_{1} = (\begin{matrix} 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 \\ 1 & 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & 1 & 0 & 0 \end{matrix}) .

3. Construction of asymptotic formulas

3.1. Case 1

Let the following conditions be satisfied:

h (x), q_{2} (x) \in L_{1} [1, \infty), \tilde{h} (x) \in L_{1, l o c} [1, \infty) .

For example, these conditions are true for

h (x) = \frac{1}{x^{γ}}, γ > 1; q (x) = x^{α} sin x^{β}, α > 0, β > \frac{α + 3}{2} .

Let the constant matrix T reduce the matrix

L_{0}

to diagonal form. Let us make a replacement

z = T u, T^{- 1} L_{0} T = Λ, μ_{k}^{5} = - i λ, k = \bar{1, 5},

Λ = (\begin{matrix} μ_{1} & 0 & 0 & 0 & 0 \\ 0 & μ_{2} & 0 & 0 & 0 \\ 0 & 0 & μ_{3} & 0 & 0 \\ 0 & 0 & 0 & μ_{3} & 0 \\ 0 & 0 & 0 & 0 & μ_{3} \end{matrix}), T = (\begin{matrix} 1 & 1 & 1 & 1 & 1 \\ μ_{1} & μ_{2} & μ_{3} & μ_{4} & μ_{5} \\ μ_{1}^{2} & μ_{2}^{2} & μ_{3}^{2} & μ_{4}^{2} & μ_{5}^{2} \\ μ_{1}^{3} & μ_{2}^{3} & μ_{3}^{3} & μ_{4}^{3} & μ_{5}^{3} \\ μ_{1}^{4} & μ_{2}^{4} & μ_{3}^{4} & μ_{4}^{4} & μ_{5}^{4} \end{matrix}) .

Then system (5) will take the form

u^{'} = (Λ + h (x) T^{- 1} L_{1} T + i q_{2} (x) T^{- 1} D_{0} T + i q_{2} (x) h (x) T^{- 1} D_{1} T) u .

(5)

Obviously, due to the imposed conditions, system (5) satisfies the conditions of Lemma 1 in [3], p. 284 and is L-diagonal, which means we can write out asymptotic formulas as

x \to \infty

for the fundamental system solutions of this system.

z_{k} (x, λ) = T \cdot u_{k} (x, λ) = e^{μ_{k} x} \cdot T \cdot (e_{k} + o (1)), k = \bar{1, 5},

where

e_{k}

are unit basis vectors.

3.2. Case 2

Let us set

\tilde{h} (x) = q_{2} (x) h (x)

. Let the following conditions be satisfied:

q_{2} (x) \notin L_{1} [1, \infty), h (x), q_{3} (x), \tilde{h} (x) \in L_{1} [1, \infty) .

(6)

These conditions are true for

h (x) = \frac{1}{x^{γ}}, γ > 1; q (x) = x^{α} sin x^{β}, α > 0, \frac{α + 3}{2} \geq β > \frac{α + 4}{3} .

Following the approach outlined in the paper [9], we make a replacement in system (4)

z = e^{i q_{3} (x) D_{0}} u .

(7)

We get

u^{'} = e^{- i q_{3} (x) D_{0}} (L_{0} + h (x) L_{1} + i \tilde{h} (x) D_{1}) e^{i q_{3} (x) D_{0}} u .

(8)

Let us apply Campbell’s identity to transform the right-hand side of (8) to

e^{- i q_{3} (x) D_{0}} L_{0} e^{i q_{2} (x) D_{0}} = L_{0} - i q_{3} (x) [D_{0}, L_{0}] + \frac{i^{2} q_{3}^{2} (x)}{2!} [D_{0}, [D_{0}, L_{0}]] -

- \frac{i^{3} q_{3}^{3} (x)}{3!} [D_{0}, [D_{0}, [D_{0}, L_{0}]]] + . . .,

here

[A, B] = A B - B A

is a matrix commutator.

Below we will use the following obvious consideration: if the matrix A is nilpotent, then a nonzero sequence of matrix commutators of the form

[A, [A, . . ., [A, B]] . . .]

is finite.

Note that the matrix

D_{0}

is nilpotent. By sequentially calculating the commutators on the right side of the last relation, we obtain that all terms, starting from the fourth, are equal to zero, and non-zero terms can be calculated

[D_{0}, L_{0}] = (\begin{matrix} 0 & 0 & 0 & 0 & 0 \\ - 1 & 0 & 0 & 0 & 0 \\ 0 & 3 & 0 & 0 & 0 \\ 0 & 0 & - 3 & 0 & 0 \\ 0 & 0 & 0 & 1 & 0 \end{matrix}), [D_{0}, [D_{0}, L_{0}]] = (\begin{matrix} 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 \\ 5 & 0 & 0 & 0 & 0 \\ 0 & 5 & 0 & 0 & 0 \end{matrix}) .

Similar calculations can be carried out for the remaining terms on the right side of (8)

e^{- i q_{3} (x) D_{0}} h (x) L_{1} e^{i q_{3} (x) D_{0}} =

h (x) L_{1} - i q_{3} h (x) [D_{0}, L_{1}] + \frac{h (x) i^{2} q_{3}^{2} (x)}{2!} [D_{0}, [D_{0}, L_{1}]] + . . . .,

[D_{0}, L_{1}] = (\begin{matrix} 0 & 0 & 0 & 0 & 0 \\ - 1 & 0 & 0 & 0 & 0 \\ 0 & 2 & 0 & 0 & 0 \\ 0 & 0 & - 2 & 0 & 0 \\ 0 & 0 & 0 & 1 & 0 \end{matrix}), [D_{0}, [D_{0}, L_{1}]] = (\begin{matrix} 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 \\ 4 & 0 & 0 & 0 & 0 \\ 0 & 4 & 0 & 0 & 0 \end{matrix}) .

Due to

[D_{0}, D_{1}] = 0

the representation is true

e^{- i q_{3} (x) D_{0}} \tilde{h} (x) D_{1} e^{i q_{3} (x) D_{0}} = \tilde{h} (x) D_{1} .

Then the equation (8) can be rewritten as

u^{'} = (L_{0} + h L_{1} + \tilde{h} D_{1} - i q_{3} [D_{0}, L_{0}] + \frac{i^{2} q_{3}^{2}}{2!} [D_{0}, [D_{0}, L_{0}]] - i q_{3} h [D_{0}, L_{1}] + \frac{i^{2} h q_{3}^{2}}{2!} [D_{0}, [D_{0}, L_{1}]]) u .

(9)

Due to the imposed conditions on the functions

h (x), q (x)

the last system can be written as:

u^{'} = (L_{0} + D (x)) u,

where

D (x)

is a matrix whose elements belong to

L_{1} [1, \infty)

. Just as in case 1, let’s make the replacement

u = T v,

then

v^{'} = (Λ + T^{- 1} D (x) T) v .

(10)

The system (10) satisfies the conditions of Lemma 1 in [3] and is L-diagonal, which means, taking into account (7), we can write asymptotic formulas for

x \to \infty

for its fundamental system of solutions

z_{k} (x, λ) = e^{μ_{k} x} \cdot e^{i q_{3} (x) D_{0}} \cdot T \cdot (e_{k} + o (1)), k = \bar{1, 5} .

where

e_{k}

are unit vectors.

Remark 1. Let us note the importance of the resulting equation (9). Imposing various conditions on the coefficients of this equation

h (x), q_{3} (x) h (x), q_{3}^{2} (x) h (x), \tilde{h} (x), q_{3} (x)

and

q_{2}^{3} (x)

, different from the conditions (6) , one can obtain different asymptotics of the fundamental system of solutions with nontrivial properties.

3.3. Case 3

Let us define the function

h_{1} (x)

so that

h_{1}^{'} (x) = h (x)

. Let us now consider the case when

h (x), q_{2} (x) \notin L_{1} [1, \infty), q_{3} (x), h_{1} (x), \tilde{h} (x) \in L_{1} [1, \infty) .

For example, these conditions are true for

h (x) = \frac{1}{x^{γ}}, 0 < γ < 1; q (x) = x^{α} sin x^{β}, 2 > α > 0, \frac{α + 3}{2} \geq β > \frac{α + 3 - γ}{2} .

Just as in case 2, let’s make a replacement in system (4)

z = e^{h_{1} (x) L_{1}} u .

(11)

Then system (4) will take the form

u^{'} = e^{- h_{1} (x) L_{1}} (L_{0} + i q_{2} (x) D_{0} + i \tilde{h} (x) D_{1}) e^{h_{1} (x) L_{1}} u .

(12)

Let us apply Campbell’s identity to transform the right-hand side of (12)

e^{- h_{1} (x) L_{1}} L_{0} e^{h_{1} (x) L_{1}} = L_{0} - h_{1} (x) [L_{1}, L_{0}] + \frac{h_{1}^{2} (x)}{2!} [L_{1}, [L_{1}, L_{0}]] - \frac{h_{1}^{3} (x)}{3!} [L_{1}, [L_{1}, [L_{1}, L_{0}]]] + . . .

Note that the matrix

L_{1}

[L_{1}, L_{0}] = (\begin{matrix} 0 & 0 & - 1 & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 & 1 \\ 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 \end{matrix}), [L_{1}, [L_{1}, L_{0}]] = (\begin{matrix} 0 & 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 0 & 1 \\ 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 \end{matrix}) .

Similar calculations can be carried out for the remaining terms on the right side of (12)

e^{- h_{1} (x) L_{1}} i q_{2} (x) D_{0} e^{h_{1} (x) L_{1}} = i q_{2} (x) D_{0} - i q_{2} (x) h_{1} (x) [L_{1}, D_{0}] + \frac{i q_{2} (x) h_{1}^{2} (x)}{2!} [L_{1}, [L_{1}, D_{0}]] -

- \frac{i q_{2} (x) h_{1}^{3} (x)}{3!} [L_{1}, [L_{1}, [L_{1}, D_{0}]]] + \frac{i q_{2} (x) h_{1}^{4} (x)}{4!} [L_{1}, [L_{1}, [L_{1}, [L_{1}, D_{0}]]]] .

[L_{1}, D_{0}] = (\begin{matrix} 0 & 0 & 0 & 0 & 0 \\ - 1 & 0 & 0 & 0 & 0 \\ 0 & 2 & 0 & 0 & 0 \\ 0 & 0 & - 2 & 0 & 0 \\ 0 & 0 & 0 & 1 & 0 \end{matrix}), [L_{1}, [L_{1}, D_{0}]] = (\begin{matrix} 0 & 0 & 0 & 0 & 0 \\ 0 & - 2 & 0 & 0 & 0 \\ 0 & 0 & 4 & 0 & 0 \\ 0 & 0 & 0 & - 2 & 0 \\ 0 & 0 & 0 & 0 & 0 \end{matrix}),

[L_{1}, [L_{1}, [L_{1}, D_{0}]]] = (\begin{matrix} 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & 6 & 0 & 0 \\ 0 & 0 & 0 & - 6 & 0 \\ 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 \end{matrix}), [L_{1}, [L_{1}, [L_{1}, [L_{1}, D_{0}]]]] = (\begin{matrix} 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & - 12 & 0 \\ 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 \end{matrix}),

e^{- h_{1} (x) L_{1}} i \tilde{h} (x) D_{1} e^{h_{1} (x) L_{1}} = i \tilde{h} (x) D_{1} - i h_{1} (x) \tilde{h} (x) [L_{1}, D_{1}],

[L_{1}, D_{1}] = (\begin{matrix} 0 & 0 & 0 & 0 & 0 \\ 1 & 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & - 1 & 0 \end{matrix}) .

Due to the imposed conditions on the functions

h (x), q (x)

the last system can be written as:

u^{'} = (L_{0} + i q_{2} (x) D_{0} + D (x)) u .

where

D (x)

is a matrix whose elements belong to

L_{1} [1, \infty)

. Unlike case 2, the resulting system is not yet L-diagonal. Let’s make one more transformation

u = e^{i q_{3} (x) D_{0}} v .

(13)

Then the last system will take the form

v^{'} = e^{- i q_{3} (x) D_{0}} (L_{0} + i q_{2} (x) D_{0} + D (x)) e^{i q_{3} (x) D_{0}} v .

(14)

Let us apply Campbell’s identity to transform the right-hand side of (14). Just as in case 2, taking into account the nilpotency of the matrix

D_{0}

and sequentially calculating all the necessary matrix commutators, we obtain the following form of system (14)

v^{'} = (L_{0} + \tilde{D} (x)) v .

Here the matrix

\tilde{D} (x)

is defined by the expression

\tilde{D} (x) = - i q_{3} (x) [D_{0}, L_{0}] + \frac{i^{2} q_{3}^{2}}{2!} [D_{0}, [D_{0}, L_{0}]] + e^{- i q_{3} (x) D_{0}} D (x) e^{i q_{3} (x) D_{0}},

which, obviously, due to the conditions imposed above on the functions

h (x), q (x)

is a matrix with elements summable over

[1, \infty)

Next we make the replacement

v = T s,

then

s^{'} = (Λ + T^{- 1} \tilde{D} (x) T) s .

(15)

System (15) satisfies the conditions of Lemma 1 in [3] and is L-diagonal, which means, taking into account (11),(13), we can write out asymptotic formulas as

x \to \infty

for its fundamental system of solutions

z_{k} (x, λ) = e^{μ_{k} x} \cdot e^{h_{1} (x) L_{1}} \cdot e^{i q_{3} (x) D_{0}} \cdot T \cdot (e_{k} + o (1)), k = \bar{1, 5},

where

e_{k}

are unit vectors.

Summarizing cases 1-3, we find that we have proven the following theorem:

Theorem 1.

Let one of the next conditions be satisfied:

h (x), q_{2} (x) \in L_{1} [1, \infty),

h (x), q_{3} (x), \tilde{h} (x) \in L_{1} [1, \infty)

q_{3} (x), h_{1} (x), \tilde{h} (x) \in L_{1} [1, \infty)

then the asymptotic formulas as

x \to \infty

for the fundamental of equation (3) are valid

y_{i} (x, λ) = e^{μ_{j} x} \cdot (1 + o (1)) .

In point of fact we obtained the asymptotic formulas as

x \to \infty

for vector-function, we may also write down the asymptotic formulas for quasi-derivatives of solutions.

3.4. Counterexample

Let us how that the conditions of Theorem 1 are essential.

Let

h (x), q_{2} (x) \notin L_{1} [1, \infty), q_{3} (x), h_{1} (x) \in L_{1} [1, \infty), \tilde{h} (x) \in L_{1, l o c} [1, \infty) .

In the same way as for case 3, we will make sequential transformations

z = e^{h_{1} (x) L_{1}} u, u = e^{i q_{3} (x) D_{0}} w,

which will bring the system of equations (4) to the form

w^{'} = (L_{0} + i h (x) q_{2} (x)) D_{1}) w + F (x) w,

(16)

where

F (x) \in L_{1} [1, \infty)

The last system of equations allows for a large variety in the asymptotic behavior of the as

x \to + \infty

and can be the subject of a separate study.

We will limit ourselves to considering a model example on which we will demonstrate an unusual property of equations with oscillating coefficients, namely the influence of the algebraic structure of the coefficients of the equation on the asymptotic behavior of the solutions.

Put

h (x) = a sin (e^{x}), q_{2} (x) = sin (k e^{x}),

from which

\tilde{h} (x) = h (x) q_{2} (x) = a sin (e^{x}) sin (k e^{x}) = \frac{1}{2} a [cos ((k - 1) e^{x}) + cos ((k + 1) e^{x})] .

Consider two cases:

k = \pm 1

and

k \neq \pm 1

. Let

k \neq \pm 1

. Define the function

{\tilde{h}}_{1} (x)

so that

{\tilde{h}}_{1}^{'} (x) = \tilde{h} (x)

. Note that in this case

{\tilde{h}}_{1} (x) \in L_{1} [1, \infty)

In the system (16) we put

w = e^{i {\tilde{h}}_{1} (x) D_{1}} v,

then for

v

we obtain the system

v^{'} = (L_{0} + F_{1} (x)) v,

(17)

where

F_{1} (x) \in L_{1} [1, \infty)

and which can easily be reduced to an L-diagonal system. Consequently, the main term of the asymptotics of the fundamental system of solutions (17) , as above, will be determined

z_{j} (x, λ) = e^{μ_{j} x} \cdot T \cdot (e_{j} + o (1)), j = \bar{1, 5} .

Let

k = \pm 1

. Note that now

h (x) q_{2} (x) = a sin (e^{x}) sin (k e^{x}) =

= \frac{1}{2} a [cos ((k - 1) e^{x}) + cos ((k + 1) e^{x})] = \frac{1}{2} a (1 + cos (2 e^{x}))

whence it follows that

h (x) q_{2} (x) \notin L_{1} [1, \infty)

. Let us denote

σ (x) = \frac{i}{2} a cos (2 e^{x})

and represent the system (17) in the form:

w^{'} = (L_{0} + \frac{i}{2} a D_{1} + σ (x) D_{1}) w + F (x) w,

(18)

where

F (x) \in L_{1} [1, \infty),

and the matrix

L_{0} + \frac{i}{2} a D_{1}

is constant. Considering

σ_{1} (x) = - \frac{i}{2} a \int_{x}^{+ \infty} cos (2 e^{ξ}) d ξ \in L_{1} [1, \infty),

we will make a replacement in the system (18)

w = e^{σ_{1} (x) D_{1}} v

and again using the technique described above, we obtain the system

v^{'} = (L_{0} + \frac{i}{2} a D_{1}) w + F_{1} (x) w .

(19)

Here, taking into account that

σ_{1} (x) \in L_{1} [1, \infty)

we have

F_{1} (x) \in L_{1} [1, \infty)

. Let the matrix

\hat{T}

reduce the matrix

L_{0} + \frac{i}{2} a D_{1}

to diagonal form,

{\hat{μ}}_{j} j = \bar{1, 5}

- eigenvalues of the matrix

L_{0} + \frac{i}{2} a D_{1}

. Let’s make a replacement

v = \hat{T} s, {\hat{T}}^{- 1} (L_{0} + \frac{i}{2} a D_{1}) \hat{T} = \hat{Λ} .

Then the system (18) will take the form

s^{'} = (\hat{Λ} + {\hat{T}}^{- 1} F_{1} (x) \hat{T}) s .

The resulting system is equivalent to the L-diagonal system, then, as above, the fundamental system of solutions equation (4) can be represented as

z_{j} (x, λ) = e^{{\hat{μ}}_{j} x} \cdot \hat{T} \cdot (e_{j} + o (1)), j = \bar{1, 5} .

Thus, when the numerical coefficient k passes through the points

k = \pm 1

, the asymptotics of the fundamental system of solutions of equation (3) undergoes a qualitative change. In other words, the points

k = \pm 1

are bifurcation points for the system (4) and corresponding equation (3).

Such bifurcation points for differential equations and (systems of equations) with regularly oscillating coefficients are typical and were recently noted by us in works devoted to the study of the asymptotic behavior of the Sturm-Liouville equation with an oscillating potential ([4,14]).

4. Discussion

The results obtained have important applications in the spectral theory of differential operators generated by the left side of equation (1). In particular, they make it possible to calculate the deficiency indices of the corresponding minimal differential operator.

The authors intend to investigate this issue in the future. In addition, we will be interested in the qualitative nature of the spectrum of such operators.

Funding

The studies of E.A. Nazirova and Ya.T. Sultanaev was funded by the Russian Science Foundation, project no. 23-21-00225.

References

Eastham, M.S.P. The Asymptotic Solution of Linear Differential Systems, Applications of the Levinson Theorem. Clarendon Press, Oxford, 1989.
Fedoryuk, M.V. Asymptotic methods for linear ordinary differential equations, (Russian) Nauka, Moscow, 1983.
Naimark, M. A. Linear differential operators, Harrap, 1967.
Nazirova, E.A.; Sultanaev, Ya.T.; Valeeva, L.N., On a Method for Studying the Asymptotics of Solutions of Sturm–Liouville Differential Equations with Rapidly Oscillating Coefficients. Math. Notes 2022, 112:6, pp.1059–1064.
Konechnaja, N.N.; Mirzoev, K.A.; Shkalikov, A.A., On the Asymptotic Behavior of Solutions to Two-Term Differential Equations with Singular Coefficients. Math. Notes 2023, 104:2, pp.244-–252.
Mirzoev, K.A.; Konechnaja, N.N., Asymptotics of solutions to linear differential equations of odd order. Vestnik Moskov. Univ. Ser.1. Mat. Mekh. 2020 75:1 pp. 22–-26.
Konechnaja, N.N.; Mirzoev, K.A.; Shkalikov, A.A., Asymptotics of Solutions of Two-Term Differential Equations. Math. Notes 2023, 113:2, pp.228-–242.
Nazirova, E.A.; Sultanaev, Ya.T.; Valeev, N.F., Distribution of the eigenvalues of singular differential operators in space of vector-functions. Trans. Moscow Math. Soc. 2014, 75, pp.89-102.
Nazirova, E.A.; Sultanaev, Ya.T.; Valeev, N.F., On a new approach for studying asymptotic behavior of solutions to singular differential equations. Trans. Ufa Math. J. 2015, 3, pp.9–14.
Myakinova, O.V.; Sultanaev, Ya.T.; Valeev, N.F., On the Asymptotics of Solutions of a Singular nth-Order Differential Equation with Nonregular Coefficients. Math. Notes 2018, 104:4, pp.606-611.
Valeev, N.F.; Nazirova, E.A.; Sultanaev, Ya.T., On a Method for Studying the Asymptotics of Solutions of Odd-Order Differential Equations with Oscillating Coefficients, Math. Notes 2021. 109:6, pp.980–-985.
Nazirova, E.A.; Sultanaev, Ya.T.; Valeev, N.F., Construction of asymptotics for solutions of Sturm-Liouville differential equations in classes of oscillating coefficients,(Russian). Vestnik Moskov. Univ. Ser. 1. Mat. Mekh. 2023, 5, pp. 61–-65.
Nazirova, E.A.; Sultanaev, Ya.T., On Deficiency Indices of Singular Differential Operator of Odd Order in Degenerate Case Azerbaijan Journal of Mathematics 2019, 9:1, pp. 125-136.
Nazirova, E.A.; Sultanaev, Ya.T.; Valeev, N.F., The new asymptotics for solutions of the Sturm-Liouville equation. Proceedings of the Institute of Mathematics and Mechanics, 2023, . 49:2, pp. 253-–258.
Bellman, R. Stability theory of differential equations. New York - Toronto - London : Mc- Graw - Hill, 1953.
Rossmann, W. Lie Groups – An Introduction Through Linear Groups. Oxford University Press, 2006.
Everitt, W.N.; Marcus, L. Boundary value problem and symplectic algebra for ordinary differential and quasi-differential operators, AMS. Mathematical Surveys and Monographs 1999, 60, pp. 1–60.

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.