Prerpints.org logo
Preprint
Article

On Solutions of the Third-Order Ordinary Differential Equations of Emden-Fowler Type

Altmetrics

Downloads

81

Views

19

Comments

0

A peer-reviewed article of this preprint also exists.

This version is not peer-reviewed

Submitted:

26 July 2023

Posted:

27 July 2023

You are already at the latest version

Alerts
Abstract
For a linear ordinary differential equation (ODE in short) of the third order, results are presented that supplement the theory of conjugate points and extremal solutions by W. Leighton, Z. Nehari, M. Hanan. It is especially noted the sensitivity of solutions to the initial data, which makes their numerical study difficult. Similar results were obtained for the third-order nonlinear equations of the Emden-Fowler type.
Keywords: 
Subject: Computer Science and Mathematics  -   Applied Mathematics

1. Introduction

Ordinary differential equations are widely used in science, mathematical modeling, physics, biology, and other fields. For students, this is another stage of mathematical learning, after elementary algebra, geometry and different types of equations. Differential equations, considered as models of natural and industrial processes, can catch dynamics, and show the development, and evolution of modeled processes. The theory of ordinary differential equations, when taught to students, starts with linear ones. There is well developed theory, which provides almost exhaustive answers to questions and problems concerning linear ODE with constant coefficients. For equations with coefficients, depending on an independent variable, the theory provides knowledge of the structure of a set of solutions and definitions of the main properties of solutions. When passing to nonlinear ODE, many principal concepts are lost, for instance, the superposition principle and extendability of solutions to infinity (for regular linear equations, without singularities). For practical purposes, numerical analysis often can be performed, and it is sufficient to find solutions. However, sensitive dependence of solutions to the initial data occurs often in practical problems and in theoretical studies of nonlinear ODE.
It can make difficult or even impossible the numerical study. In this article we attract attention to the occurrence of this and similar problems even in the study of relatively simple linear and nonlinear ODE. For this, we have chosen the theory of conjugate points, as developed in the works of W. Leighton, Z. Nehari [8], M. Hanan [4], T. Sherman [10], V. Kondrat’ev [6,7], etc.
Few words about the conjugate points. The first order linear ODE of the type x ′ = a ( t ) x + b ( t ) have no surprises, if the coefficients are continuous functions. The second order linear ODE of the form x ′ ′ + p ( t ) x ′ + g ( t ) x = f ( t ) are more intriguing. They can exhibit oscillatory behavior of solutions, and a lot of problems arise. It is firm, however, that the oscillatory behavior of solutions is intuitively connected to zeros, and the rate of oscillation can be measured by the number of zeros.
When passing to linear equations of order three, and higher, the notion of a zero of a solution becomes more complicated. There are multiple zeros, that is, zeros for a solution, and some of its derivatives. Nevertheless, the rate of oscillation of some classes of linear ODE of higher orders can be measured, introducing the notion of conjugate points. Definitions for the third order linear equations are provided below, following the work by M. Hanan, who, in turn, used the theory of conjugate points for classes of linear equations of order four, developed earlier by W. Leighton and Z. Nehari. The definition of conjugate points for the third order linear equations seems to be tricky since it uses the extremality property of zeros and leads to the concept of an extremal solution. Further analysis of these definitions simplifies the problem and results in the efficient criteria for finding conjugate points, as in Theorem uid7. So conjugate points for some classes of equations of order higher than two, play the same role as ordinary zeros play for the second order equations. So there are conjugate points, and there are the so called extremal solutions, attributed to them. In this article we provide some additional information on the theory of conjugate points for linear ODE of order three. The structure of extremal solutions is revealed, and their remarkable properties are disputed. Namely, the mutual location of extremal solutions is described and, among others, the sensitive dependence of extremal solutions on the initial conditions is discussed.
The second part of the article concerns the nonlinear equations of the Emden-Fowler type, which in some sense, behave similarly to linear equations, discussed in the `linear’ part of the paper. Both the analytical approach and the numerical study are used. Some results in the literature, concerning nonlinear equations, are reminded and their consequences to extremal solutions for the Emden-Fowler type equations are formulated.
Consider the linear equation
x ′ ′ ′ = p ( t ) x ,
where p ( t ) is a positive valued continuous function. This equation was studied by many authors; we will mention the works [4,6,7], books [1,3,5,11].
The theorem by M. Hanan [4] states that there exist special solutions x i ( t ) with
x ( 0 ) = 0 , x ′ ( 0 ) > 0 , x ′ ′ ( 0 ) < 0 ,
which have a double zero at some point η i > 0 . These points form ascending sequence and are called by conjugate points to t = 0 for equations of the form (1). These points can be visualized considering the equation
x ′ ′ ′ = p x , p = constant > 0 .
We recall the structure of a set of solutions to equations of the type (1) and discuss what happens under passage to nonlinear equations
x ′ ′ ′ = p | x | γ x , γ > 0 .

2. Linear equation

W. Leighton and Z. Nehari [8] have investigated oscillatory properties of solutions to linear equations of the form
x ( 4 ) = p ( t ) x ,
where p ( t ) is a continuous positive (or negative) valued function. To measure the rate of oscillation of equation (5) they introduced the notion of conjugate points. Similar task was accomplished by M. Hanan [4] who studied the third order linear differential equations
x ′ ′ ′ = p ( t ) x + q ( t ) x ′ + r ( t ) x ′ ′
with continuous coefficients. The conjugate points for the equation (6) were introduced in the following way. Suppose that equation (6) has a solution x ( t ) which vanishes at t = a and has at least n + 2 zeros a 1 , … , a n + 2 in ( a , + ∞ ) . Then η n is called ([4], p. 920) the n-th conjugate point to t = a with respect to equation (6) if it is a smallest possible value of a n + 2 as x ( t ) ranges over all possible solutions of (6) for which x ( a ) = 0 . Due to this extremality property of a conjugate point the solution which produces the n-th conjugate point is called the n-th extremal solution.
To get description of conjugate points the following characterization of linear equations was introduced.
Definition 1 [4]. Equation (6) is said to be of Class I if any its solution x ( t ) for which x ( a ) = x ′ ( a ) = 0 and x ′ ′ ( a ) > 0 is positive for t < a .
Definition 2 [4]. Equation (6) is said to be of Class II if any its solution x ( t ) for which x ( a ) = x ′ ( a ) = 0 and x ′ ′ ( a ) > 0 is positive for t > a .
There are equations that belong to both classes (for instance, x ′ ′ ′ = 0 ) and there exist equations (for example, x ′ ′ ′ = − x ′ ) that are neither of Class I nor of Class II.
Several criteria for equation (6) to be of Class I or Class II are given in [3,4,11] of which we mention only the simplest one. Namely, equation (1) is of Class I if p ( t ) < 0 and of Class II if p ( t ) > 0 .
The characterization of conjugate points for equations of Class I and Class II was given by M. Hanan ([4] theorems 2.6, 2.7 and 4.4).
Theorem 2.1
Conjugate points η i ( a ) of an equation of Class I are the zeros (which are simple) of the so calledprincipal solution, that is the solution x ( t ) which satisfies the initial conditions
x ( a ) = x ′ ( a ) = 0 , x ′ ′ ( a ) = 1 .
Conjugate points η i ( a ) (if any) of an equation of Class II form ascending sequence
a < η 1 < … < η n < …
The respective extremal solutions x i ( t ) ( i = 1 , … ) have a simple zero at t = a , double zero at η i (that is, x i ( η i ) = x ′ ( η i ) = 0 ) and exactly i − 1 simple zeros in ( a , η i ( a ) .
From now consider equations (1) with positive p ( t ) . They belong to Class II. It is an easy matter to show that nontrivial solutions of (1) with the initial conditions x ( a ) = 0 , x ′ ( a ) ≥ 0 , x ′ ′ ( a ) ≥ 0 or x ( a ) = 0 , x ′ ( a ) ≤ 0 , x ′ ′ ( a ) ≤ 0 do not vanish for t > a . So any extremal solution x i ( t ) satisfies (up to multiplication by − 1 ) the conditions x i ( a ) = 0 ,   x i ′ ( a ) > 0 , x i ′ ′ ( a ) < 0 . It is convenient to attribute the angle ϕ i = arctan x i ′ ′ ( a ) x i ′ ( a ) to the extremal solution.
Studying the properties of these angles have led us to the following result.
Theorem 2.2
For a linear equation of Class II:
1) the angles ϕ i ( a ) corresponding to the extremal solutions x i ( t ) are arranged as
− π / 2 < ϕ 2 < ϕ 4 < … < ϕ 2 i < … < ϕ 2 i + 1 < … < ϕ 3 < ϕ 1 < 0 ;
2) solutions x ( t ) , defined by the initial conditions
x ( a ) = 0 , ϕ ∈ ( ϕ 2 i , ϕ 2 i + 2 ) , ϕ = arctan x ′ ′ ( a ) x ′ ( a ) i = 0 , 1 , 2 , …
have for t > a exactly 2 i + 2 simple zeros, and solutions defined by the initial conditions
x ( a ) = 0 , ϕ ∈ ( ϕ 2 i + 1 , ϕ 2 i − 1 ) , i = 0 , 1 , 2 , …
have for t > a exactly 2 i + 1 simple zeros ( ϕ 0 means − π 2 and ϕ − 1 is set to zero).
Proof. 
We consider extremal solutions x k ( t ) , for which the first order derivative at t = a is positive, and the second order derivative is negative. Taking into account that all zeros of an extremal function x k ( t ) in the interval ( a , η k ) are simple, and x k ( t ) changes sign at each of them, we conclude that x k ( t ) is positive for t > η k , if k is odd, and it is negative for the same t if k is even.
Let us compare two extremal solutions with odd numbers. We wish to show that ϕ 2 k + 1 > ϕ 2 l + l if k < l . Suppose that x 2 k + 1 ′ ( a ) = x 2 l + 1 ′ ( a ) and comparison is made only for the second order derivatives at t = a . This is possible always, since extremal solutions can be multiplied by a constant, and this does not change the corresponding angle ϕ ( a ) . By Theorem uid7, η 2 k + 1 < η 2 l + 1 . For t > η 2 l + 1 both extremal functions are positive. Let us show that the assumption x 2 l + 1 ′ ′ ( a ) > x 2 k + 1 ′ ′ ( a ) or, which is the same, ϕ 2 k + 1 < ϕ 2 l + l , leads to contradiction. Note, that the case x 2 l + 1 ′ ′ ( a ) = x 2 k + 1 ′ ′ ( a ) is excluded, since then both extremal solutions must coincide, by the unique solvability of the Cauchy problems. Consider the difference y = x 2 l + 1 − x 2 k + 1 . The function y ( t ) is a solution of the same equation with the initial conditions y ( a ) = 0 , y ′ ( a ) = 0 , y ′ ′ ( a ) > 0 . Therefore y ( t ) > 0 for t > a , since the equation is of Class II, and any solution with positive y ′ ′ ( a ) is positive to the right of a double zero. However, this contradicts the fact that y ( η 2 l + 1 ) = x 2 l + 1 ( η 2 l + 1 ) − x 2 k + 1 ( η 2 l + 1 ) < 0 .
It can be proved similarly, that ϕ 2 l > ϕ 2 k if l > k .
Let us show now, that ϕ 2 k − 1 > ϕ 2 l for any positive integers k and l . We assume that x 2 k − 1 ′ ( a ) = x 2 l ′ ( a ) and x 2 k − 1 ′ ′ ( a ) < x 2 l ′ ′ ( a ) . Consider the function y = x 2 l − x 2 k − 1 . Since y ( t ) has a double zero at t = a and y ′ ′ ( a ) > 0 , by assumption, y ( t ) is positive for t > a . Two cases are possible. The first one, 2 l > 2 k − 1 . Then η 2 l > η 2 k − 1 . One has that y ( η 2 l ) = x 2 l ( η 2 l ) − x 2 k − 1 ( η 2 l ) < 0 . This is in contradiction with positivity of y ( t ) for t > a . The second possible case, 2 l < 2 k − 1 . Then η 2 l < η 2 k − 1 . One has that y ( η 2 k − 1 ) = x 2 l ( η 2 k − 1 ) − x 2 k − 1 ( η 2 k − 1 ) < 0 . The contradiction with positivity of y ( t ) for t > a is obtained again.
The first statement of the theorem is proved.
Let us pass to the second statement. Consider a solution x ( t ) which vanishes at t = a and is defined by the angle ϕ ∈ ( ϕ 2 k − 1 , ϕ 2 k + 1 ) , k = 1 , 2 , … Let us compare this solution with the extremal solutions x 2 k − 1 and x 2 k + 1 provided that the first derivatives of all three solutions at t = a are equal. The functions u = x 2 k − 1 − x and v = x − x 2 k + 1 are positive for t > a . It follows that x 2 k − 1 ( t ) > x ( t ) > x 2 k + 1 ( t ) for t > a . Denote simple zeros of the function x 2 k − 1 in the interval ( a , η 2 k − 1 ) by τ i and simple zeros of the function x 2 k + 1 by t i . Consider the relative locations of zeros τ i and t j . For even i the function x 2 k − 1 is positive in ( τ i , τ i + 1 ) , and the function x 2 k + 1 is positive in ( t i , t i + 1 ) (we accept that τ 0 = t 0 = a ). Show that τ i < t i < t i + 1 < τ i + 1 for i = 0 , 2 , 4 , … If this is not the case, then there exists an interval ( τ i , τ i + 1 ) with even i such that x 2 k + 1 ( t ) < 0 for t ∈ [ τ i , τ i + 1 ] . Then, by Lemma 2 in [8], there exists a number μ such that the function z ( t ) = x 2 k − 1 ( t ) + μ x 2 k + 1 ( t ) has a double zero in the interval ( τ i , t i + 1 ) . Since the function z ( t ) is a solution of Class II, it does not vanish for t > t i + 1 . On the other hand, the number μ is positive, since otherwise z = x 2 k − 1 + μ x 2 k + 1 > 0 for t ∈ ( τ i , t i + 1 ) , what contradicts the existence of a zero in this interval. Then z ( η 2 k − 1 ) = x 2 k − 1 ( η 2 k − 1 ) + μ x 2 k + 1 ( η 2 k − 1 ) < 0 because the first addend is zero, and x 2 k + 1 ( η 2 k − 1 ) < x 2 k − 1 ( η 2 k − 1 ) < 0 . But z ( η 2 k + 1 ) = x 2 k − 1 ( η 2 k + 1 ) + μ x 2 k + 1 ( η 2 k + 1 ) > 0 since the second addend is zero, and x 2 k − 1 ( t ) is positive for t > η 2 k − 1 . Therefore, the function z ( t ) changes sign in the interval ( η 2 k − 1 , η 2 k + 1 ) , which lies to the right of t = t i + 1 . The obtained contradiction means that any interval ( τ i , t i + 1 ) , where the function x 2 k − 1 is positive, contains a subinterval ( t i , t i + 1 ) , where the function x 2 k + 1 ( t ) is positive. We have proved that zeros of the functions x 2 k − 1 and x 2 k + 1 are arranged as
a < t 1 < τ 1 < τ 2 < t 2 < t 3 < … < τ 2 k − 2 < t 2 k − 2 < t 2 k − 1 < t 2 k < t 2 k + 1 < η 2 k + 1 .
It follows from the inequalities x 2 k − 1 ( t ) < x ( t ) < x 2 k + 1 ( t ) , which are valid for t > a , that in any interval ( t 1 , τ 1 ) , … , ( τ 2 k − 2 , t 2 k − 2 ) , ( t 2 k − 1 , η 2 k − 1 ) there exists a zero of x ( t ) . Thus the function x ( t ) has in the interval ( a , η 2 k − 1 ) totally 1 + ( 2 k − 1 ) = 2 k zeros, and x ( η 2 k − 1 ) < 0 . It follows that there exists an extra zero in the interval ( η 2 k − 1 , t 2 k ) , and the number of zeros of x ( t ) in [ a , t 2 k ) is not less than 2 k + 1 . At the same time, the function x ( t ) cannot have more zeros. If it had more, then it had more by two zeros at least. Then the minimal number of zeros in the interval [ a , t 2 k ] would be 2 k + 3 . Since t 2 k < η 2 k + 1 , this would contradict the choice of η 2 k + 1 as the minimal 2 k + 3 -th zero over all solutions, vanishing at t = a (recall the “extremal” definition of a conjugate point above). The conclusion is that x ( t ) has exactly 2 k + 1 zeros for t > a .
The proof for solutions, defined by the inequalities ϕ ∈ ( ϕ 2 k , ϕ 2 k + 2 ) , can be conducted in a similar way. □
Consider, for instance, the equation
x ′ ′ ′ = x ,
which has a general solution
x ( t ) = C 1 e t + C 2 e − t / 2 cos 3 2 t + C 3 e − t / 2 sin 3 2 t .
Solutions that vanish at t = 0 are
x ( t ) = − e − t / 2 e 3 t / 2 C 2 − C 2 cos 3 2 t − C 3 sin 3 2 t .
In order to find solutions that are zero at t = 0 and have double zeros to the right of t = 0 consider the system (with respect to the unknowns C 2 and C 3 )
x ( η ) = C 2 ( − e η + e − η / 2 cos 3 2 η ) + C 3 e − η / 2 sin 3 2 η = 0 , x ′ ( η ) = C 2 ( − e η − 1 2 e − η / 2 cos 3 2 η − 3 2 e − η / 2 sin 3 2 η ) + C 3 ( − 1 2 e − η / 2 sin 3 2 η + 3 2 e − η / 2 cos 3 2 η ) = 0 .
This system can have nontrivial solutions C 2 and C 3 only if the coefficient determinant is zero. Therefore we get the equation to find the conjugate (to t = 0 ) points η k
3 2 e − η + 3 2 e η / 2 sin 3 2 η − 3 2 e η / 2 cos 3 2 η = 0
The graph (appropriately scaled to fit the area) of the function on the left side is depicted in Figure 1.
Conjugate points to t = 0 of equation x ′ ′ ′ = x coincide with conjugate points to t = 0 of the adjoint equation y ′ ′ ′ = − y . The conjugate points ξ k of the equation y ′ ′ ′ = − y are exactly the zeros (all of them are simple zeros) of a solution y ( t ) that satisfies the initial conditions y ( 0 ) = y ′ ( 0 ) = 0 .
First four extremal solutions (with respect to t = 0 ) for equation (9) are depicted in Figure 2.
The angles ϕ k = arctan x k ′ ′ ( 0 ) x k ′ ( 0 ) that define the initial values for extremal solutions  x k ( t ) are arranged
− π / 2 < ϕ 2 < ϕ 4 < … < ϕ 2 i < … < ϕ 2 i + 1 < … < ϕ 3 < ϕ 1 < 0 .
The respective extremal solutions  x k ( t ) are defined then as
x k ( t ) = tan ϕ k + 1 3 e t + e − 1 2 t 1 − tan ϕ k 3 sin 3 2 t − 1 + tan ϕ k 3 cos 3 2 t
The extremal solutions expressed in terms of conjugate points η
x ( t ) = e − 1 2 η sin 3 2 η e t − e − 1 2 t cos 3 2 t + e − 1 2 t cos 3 2 η − e η e − 1 2 t sin 3 2 t
Below the tan ϕ k are given (rounded to 12 digits) for the first six extremal solutions
tan ϕ 1 = − 0.9969785430602899 ; tan ϕ 2 = − 1.0000131330764497 ; tan ϕ 3 = − 0.9999999430896078 ; tan ϕ 4 = − 1.0000000002466165 ; tan ϕ 5 = − 0.9999999999989313 ; tan ϕ 6 = − 1.0000000000000047 .
The sequences { ϕ 2 i } and { ϕ 2 i + 1 } tend from below and from above respectively to the angle ϕ * = − π 4 . The particular solution defined by the angle ϕ * is
x * ( t ) = 2 3 e − 1 2 t sin 3 2 t .
This is oscillatory solution with equidistant zeros.
Preprints 80623 g003
Figure 3. Angles ϕ k tending to ϕ * (schematically, since real convergence is hardly visualizable because the straight lines are very close to the limit).
Figure 3. Angles ϕ k tending to ϕ * (schematically, since real convergence is hardly visualizable because the straight lines are very close to the limit).
Preprints 80623 g003
The convergence ϕ i → − π 4 is rapid and therefore the extremal solutions are hardly detectable numerically.

3. Phase plane for linear equations

In order to visualize extremal solutions consider the variable change [1]
u 1 = x ′ x − 1 , u 2 = x ′ ′ x − 1 ,
which acts on intervals where x ( t ) ≠ 0 . After standard calculations one arrives at the system
u 1 ′ = u 2 − u 1 2 , u 2 ′ = 1 − u 1 u 2 .
The phase portrait for this system is depicted in Figure 4.

4. Emden Fowler type equations

Consider equations of the type
x ′ ′ ′ = p | x | γ x ,
for γ > 0 . These equations exhibit generally the same behavior as the linear one with the exception that all solutions that eventually tend to infinity x ( t ) → ∞ do so in finite time.
It can be proved that the plane of the initial data ( x ′ ( a ) , x ′ ′ ( a ) ) looks generally as that for the linear case (Figure 3).
Instead of straight lines defined by angles ϕ k in Figure 3 there are curves Γ k that are still arranged in two sequences converging to a specific branch Γ * . This branch contains the initial data (namely, x ′ ( a ) and x ′ ′ ( a ) ) for extendable and oscillatory solutions that have not double zeros.
Generally, the following is true for solutions of (17) that vanish at t = a . There exists a one-dimensional branch of the initial conditions Γ 1 that extends to the 4-th quadrant of the ( x ′ ( a ) , x ′ ′ ( a ) ) -plane ( x ′ ( a ) > 0 , x ′ ′ ( a ) < 0 ) and possesses the property: any solution of
x ′ ′ ′ = p | x | γ x , x ( a ) = 0 , ( x ′ ( a ) , x ′ ′ ( a ) ) ∈ Γ 1
is positive in the interval ( a , η 1 ) , has a double zero at η 1 and is positive for t > η 1 . In the case of a linear equation ( γ = 0 in (17)) this branch is the straight line marked as ϕ 1 in Figure 3. Similar suggestions can be stated for other branches.
Consider equation (17) with p = c o n s t > 0 and γ > 0 .
Theorem 4.1For any t = a there exist branches Γ i ( i = 1 , 2 , … ) of the initial values ( x ′ ( a ) , x ′ ′ ( a ) ) possessing the properties:
1) Γ i locate in the quadrant ( x ′ > 0 , x ′ ′ < 0 ) of the ( x ′ ( a ) , x ′ ′ ( a ) ) -plane;
2) Γ i emanate from the origin ( 0 , 0 ) and extend to infinity;
3) Γ i are ordered as
Γ 2 ≺ Γ 4 ≺ … ≺ Γ 2 i ≺ … ≺ Γ 2 j + 1 ≺ … ≺ Γ 3 ≺ Γ 1
in the meaning that any vertical line x ′ = c o n s t > 0 crosses branches Γ i in indicated order;
4) solutions x i ( t ) of equation (17) subject to the initial conditions x ( a ) = 0 , ( x ′ ( a ) , x ′ ′ ( a ) ) ∈ Γ 2 i + 1 , ( i = 0 , 1 , … ) have exactly i − 1 simple zeros in the interval ( a , η ) and a double zero at t = η (different η for different x);
5) solutions x i ( t ) of equation (17) subject to the initial conditions x ( a ) = 0 , ( x ′ ( a ) , x ′ ′ ( a ) ) ∈ Γ 2 i , ( i = 1 , … ) have exactly i simple zeros in the interval ( a , η ) and a double zero at t = η (different η for different x);
6) solutions x ( t ) , defined by the initial conditions
x ( a ) = 0 , Γ 2 i ≺ ( x ′ ( a ) , x ′ ′ ( a ) ) ≺ Γ 2 i + 2 ,
have for t > a exactly 2 i + 2 simple zeros, and solutions defined by the initial conditions
x ( a ) = 0 , Γ 2 i + 1 ≺ ( x ′ ( a ) , x ′ ′ ( a ) ) ≺ Γ 2 i − 1 , i = 0 , 1 , 2 , …
have for t > a exactly 2 i + 1 simple zeros ( Γ 0 means semi-axis ( x ′ ( a ) = 0 , x ′ ′ ( a ) < 0 ) and Γ − 1 means ( x ′ ( a ) > 0 , x ′ ′ ( a ) = 0 ) ).
Proof. Fix R > 0 and consider a quarter of a circle of radius R in the ( x ′ ( a ) , x ′ ′ ( a ) ) -plane lying in the 4-th quadrant ( x ′ > 0 ,   x ′ ′ < 0 ). Parameterize the arc by the angle ϕ = arctan x ′ ′ ( a ) x ′ ( a ) that takes values in the interval [ 0 , − π 2 ] . Due to Lemma uid28 a solution with the initial conditions x ( a ) = 0 , x ′ ( a ) = R , x ′ ′ ( a ) = 0 or, equivalently φ = 0 , is positive for t > a . A solution with the initial conditions x ( a ) = 0 , x ′ ( a ) = 0 , x ′ ′ ( a ) = − R or, equivalently φ = − π 2 , is negative for t > a .
Denote solutions of equation (17) satisfying the initial conditions x ( a ) = 0 , x ′ 2 ( a ) + x ′ ′ 2 ( a ) = R 2 , x ′ ( a ) > 0 , x ′ ′ ( a ) < 0 by x ( t ; φ ) . Solutions x ( t ; φ ) continuously change together with φ .
Let φ decrease from zero. For small φ solutions x ( t ; φ ) are positive for t > a . Let φ 1 = inf { φ ∈ ( − π / 2 , 0 ) : x ( t ) > 0 f o r t > a i n t h e i n t e r v a l o f d e f i n i t i o n } . Such a value exists since solutions x ( t ; φ ) for φ → − π / 2 are negative. A solution x ( t ; φ 1 ) has a zero, otherwise it is positive in the interval of definition.
This zero (denote it η 1 ( R ) ) must be double zero ( x = 0 ,   x ′ = 0 ). If it is simple ( x ′ ≠ 0 ) then there exists some φ > φ 1 with the same property and this contradicts the definition of φ 1 .
Proceeding in this way we can obtain η 2 and so on for solutions that are eventually positive.
Considering solutions x ( t ; φ ) starting from φ = − π / 2 we obtain solutions with a double zero at η 2 , η 4 and so on. These solutions are negative eventually.
By construction,
η 1 < η 3 < … < η o d d < …
and, similarly,
η 2 < η 4 < … < η e v e n < …
It can be shown that odd and even η alternate. Look at Figures –. If solutions x ( t ; φ 1 ) ,   x ( t ; φ 2 ) and x ( t ; φ 3 ) behave like in Figure then “between” “blue” and “magenta” solutions there exists a solution with the first double zero ξ < η 1 . This contradicts choice of η 1 as the first double zero. Similarly for positions of solutions as in Figure . The only possible location of η -s is as in Figure 6.
Preprints 80623 g005
Figure 5.  
Figure 5.  
Preprints 80623 g005
Preprints 80623 g006
Figure 6.  
Figure 6.  
Preprints 80623 g006
Preprints 80623 g007
Figure 7.  
Figure 7.  
Preprints 80623 g007

5. More on nonlinear equations

Solutions of equations of the form
x ′ ′ ′ = p | x | γ x ,
where continuous coefficient p = p ( t ) > 0 and γ > 0 possess the following properties.
Lemma 5.1Any solution of an equation (19) with the initial conditions
x ( a ) = 0 , x ′ ( a ) ≥ 0 , x ′ ′ ( a ) ≥ 0 , x ′ 2 ( a ) + x ′ ′ 2 ( a ) > 0
remains positive for t > a .
Similarly, any solution with the conditions
x ( a ) = 0 , x ′ ( a ) ≤ 0 , x ′ ′ ( a ) ≤ 0 , x ′ 2 ( a ) + x ′ ′ 2 ( a ) > 0
remains negative for t > a .
Proof. Let (20) hold. Then one of the values x ′ ( a ) or x ′ ′ ( a ) is positive and x ( t ) is positive also in some right neighborhood ( a , b ) of the point t = a . Suppose x ( b ) = 0 . Then x ′ ( b 1 ) = 0 for some b 1 ∈ ( a , b ) and x ′ ′ ( b 2 ) < 0 for some b 2 ∈ ( a , b 1 ) . This is impossible, however, since x ′ ′ ( a ) ≥ 0 and x ′ ′ ′ ( t ) > 0 whenever x ( t ) is positive, namely, x ′ ′ ′ ( t ) > 0 for t ∈ ( a , b ) .
Similarly x ( t ) < 0 for t > a if (21) holds.□
Since we are looking for solutions of (19) that have double zero consider the case
x ( a ) = 0 , x ′ ( a ) > 0 , x ′ ′ ( a ) < 0 .
The case
x ( a ) = 0 , x ′ ( a ) < 0 , x ′ ′ ( a ) > 0
is symmetric.
Due to Theorem uid25 there are solutions with the initial values x ( a ) = 0 , ( x ′ ( a ) , x ′ ′ ( a ) ) ∈ Γ 1 , which have a double zero at some point η > a , and x ′ ′ ( η ) > 0 . Then such a solution is positive for t > η . Solutions with x ( a ) = 0 and ( x ′ ( a ) , x ′ ′ ( a ) ) ∈ Γ i also have a double zero at some point η (individual for any solution) and do not change the sign for t > η . For odd i these solutions are positive and for even i they are negative for t > η . Therefore the following assertion is true.
Lemma 5.2Let 0 < p m i n ≤ p ( t ) ≤ p m a x . Then any solution with x ( a ) = 0 , x ′ ( a ) > 0 which has a double zero at some point η > a has the vertical asymptote.
Proof. Follows from Lemma 4.2 and Lemma 4.3 in [1]. □
Suppose x ( t ) is a solution of the equation (19) with the initial conditions x ( 0 ) = 0 , x ′ ( 0 ) > 0 , x ′ ′ ( 0 ) < 0 and with a double zero at t = η > 0 . For brevity, let us call such a solution extremal solution. In contrast to the linear case, multiples of extremal solutions need not be an extremal one.
Let p = c o n s t > 0 . Any extremal solution may oscillate, having several zeros in the interval ( 0 , η ) and a double zero at t = η . Then it is either eventually positive, or negative, like in Figure 6 and Figure 6. Due to Lemma uid33, any such solution must have the vertical asymptote. The next statement, which is a consequence of Theorem 4.4 in [1], describes the asymptotics of extremal solutions.
Lemma 5.3Any extremal solution has asymptotics C ( t * − t ) − γ 3 ( 1 + o ( 1 ) ) , where
C = 3 γ ( 3 γ + 1 ) ( 3 γ + 2 ) p 3 γ
and t = t * is the vertical asymptote for a solution.
There is another description of the asymptotic behavior of extremal solutions. Let us, following ([1] §4.3), introduce new functions
u 1 ( t ) = x ′ ( t ) x β 1 , u 2 ( t ) = x ′ ′ ( t ) x β 2 , w h e r e β 1 = 1 + γ 3 , β 2 = 1 + 2 γ 3 .
By standard differentiation of u i ( t ) with respect to t and taking into account (17) we can obtain a two-dimensional non-autonomous system of differential equations. Continue with introduction the independent variable z = ∫ η t x γ 3 ( s ) d s for any extremal solution x ( t ) . In view of Lemmas uid33 and uid34, z → + ∞ as t → t * , where t = t * is an asymptote for x ( t ) . Functions u i , considered as functions of the variable z , satisfy the autonomous system
u 1 ′ = u 2 − 3 + γ 3 u 1 2 , u 2 ′ = p − 3 + 2 γ 3 u 1 u 2
Any extremal solution x ( t ) generates a phase trajectory ( u 1 ( z ) , u 2 ( z ) ) , which tends to a single critical point of (25).
Preprints 80623 g008
Figure 8. Phase portrait for system (25), p = 1 , γ = 2 .
Figure 8. Phase portrait for system (25), p = 1 , γ = 2 .
Preprints 80623 g008

6. Conclusion

The structure of a set of extremal solutions, associated with the conjugate points accordingly to the linear theory by M. Hanan [4], is revealed. There are two groups of extremal solutions. Solutions of the first group, after several oscillations, go to + ∞ . Solutions of the second group go to − ∞ . Their initial conditions become extremely close, so they are very difficult to discover numerically. Both groups are separated by a solution with an infinite number of zeros. This solution for p = c o n s t in the equation (1) is unique.
For nonlinear equations of the Emden-Fowler type, the structure of a set of solutions, satisfying the initial conditions x ( 0 ) = 0 ,   x ′ ( 0 ) > 0 , x ′ ′ ( 0 ) < 0 , and having a double zero to the right of t = 0 , is similar to that of linear equations. There are no conjugate points, however. To find numerically the initial conditions for a solution, which has a double zero at some prescribed point, can be a difficult task. Solutions of nonlinear equations also form two groups. Solutions of both groups end respectively at + ∞ , or at − ∞ . Both groups are separated by an oscillatory solution (-ns), which is hard to find. In the two below pictures the graphs of two solutions to the equation x ′ ′ ′ = | x | 2 x are depicted.
Preprints 80623 g009
Figure 9. The solution of (17) with p = 1 , γ = 2 , x ( 0 ) = 0 , x ′ ( 0 ) = 1 , x ′ ′ ( 0 ) = 0.71193167642507565462395 .
Figure 9. The solution of (17) with p = 1 , γ = 2 , x ( 0 ) = 0 , x ′ ( 0 ) = 1 , x ′ ′ ( 0 ) = 0.71193167642507565462395 .
Preprints 80623 g009
Preprints 80623 g010
Figure 10. The solution of (17) with p = 1 , γ = 2 , x ( 0 ) = 0 , x ′ ( 0 ) = 1 , x ′ ′ ( 0 ) = 0.71193167642507565462396 .
Figure 10. The solution of (17) with p = 1 , γ = 2 , x ( 0 ) = 0 , x ′ ( 0 ) = 1 , x ′ ′ ( 0 ) = 0.71193167642507565462396 .
Preprints 80623 g010
They belong to different groups, their behavior is essentially different, but the initial conditions for the second derivative differ only after 22 decimal places.
When looking for solutions of a certain structure (for example, solutions with a given number of zeros), one must be careful, since even for formally simple ODE, they (solutions) can be overlooked when conducting numerical experiments. This applies to ODEs already of the third order and higher. Data on the mutual arrangement of solutions can be used in the study of multiple solutions of boundary value problems [9].

References

  1. I.V. Astashova. Qualitative properties of solutions to quasi-linear ordinary differential equations. Moscow St. University of Economics, Statistics and Informatics, Moscow, 2010 (Russian). [CrossRef]
  2. M. Gera and F. Sadyrbaev. Multiple solutions of a third order boundary value problem. Math. Slovaca, 42 (1992), N 2, 173-180.
  3. M. Greguš. Third Order Linear Differential Equations. Series: Mathematics and its Applications, Vol. 22. Springer, 1987.
  4. M. Hanan. Oscillation criteria for third order linear differential equations. Pacific J. Math., Vol. 11 (1961), 919-944.
  5. I. Kiguradze, T.A. I. Kiguradze, T.A. Chanturia. Asymptotic Properties of Solutions of Nonautonomous Ordinary Differential Equations, Springer Science & Business Media, 2012.
  6. V.A. Kondrat’ev, V. A. Oscillation of solutions of linear equations of third and fourth order. (Russian) Trudy Moskov. Mat. Obsc. 8 1959 259-281. MR104877.
  7. V. A. Kondrat’ev, Oscillatory properties of solutions of the equation y(n)+p(x)y=0. (Russian) Trudy Moskov. Mat. Obsc. 10 1961 419-436. MR141842 34.
  8. W. Leighton and Z. Nehari. On the oscillation of solutions of self-adjoint linear differential equations of the fourth order. Trans. Amer. Math. Soc., 89:325-377, 1961.
  9. F. Sadyrbaev. Multiplicity results for third order two-point boundary value problems. Acta Universitatis Latviensis, ser. Mathematics. Differential Equations(ISSN 1407-2157), vol. 616, 1999, 5-16. https://papers.ssrn.com/sol3/papers.cfm?abstract_id=3148561.
  10. T. Sherman. Conjugate points and simple zeros for ordinary linear differential equations. Trans. Amer. Math. Soc., 146:397-411, 1969. [CrossRef]
  11. C. Swanson. Comparison and oscillation theory of linear differential equations. Academic Press, New York, 196.
Preprints 80623 g001
Figure 1. First five conjugate points η 1 = 4.233207192 , η 2 = 7.859792868 , η 3 = 11.48739599 , η 4 = 15.11499470 , η 5 = 18.74259343 .
Figure 1. First five conjugate points η 1 = 4.233207192 , η 2 = 7.859792868 , η 3 = 11.48739599 , η 4 = 15.11499470 , η 5 = 18.74259343 .
Preprints 80623 g001
Preprints 80623 g002
Figure 2. First four extremal solutions (scaled appropriately to fit the plot) corresponding to the conjugate points η 1 = 4.233207192 , η 2 = 7.859792868 , η 3 = 11.48739599 , η 4 = 15.11499470 .
Figure 2. First four extremal solutions (scaled appropriately to fit the plot) corresponding to the conjugate points η 1 = 4.233207192 , η 2 = 7.859792868 , η 3 = 11.48739599 , η 4 = 15.11499470 .
Preprints 80623 g002
Preprints 80623 g004
Figure 4. Phase portrait for system (16)
Figure 4. Phase portrait for system (16)
Preprints 80623 g004
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.
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.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

facebook logotwitter logolinkedin logo
MDPI Initiatives
Important Links
Subscribe

Disclaimer

Terms of Use

Privacy Policy

© 2024 MDPI (Basel, Switzerland) unless otherwise stated