Asymptotic Behaviour of Solutions of Certain Third Order Nonlinear Differential Equations via Phase Portrait Analysis

The global phase portrait describes the qualitative behaviour of the solution set for all time. In general, this is as close as we can get to solving nonlinear systems. The question of particular interest is: For what parameter values does the global phase portrait of a dynamical system change its qualitative structure? In this paper, we attempt to answer the above question specifically for the case of certain third order nonlinear differential equations of the form ( ) 0 x ax g x cx + + + =    . The linear case where ( ) g x bx =   is also considered. Our phase portrait analysis shows that under certain conditions on the coefficients as well as the function g , we have asymptotic stability of solutions.


Introduction
A line connecting the plotted points in their chronological order shows temporal evolution more clearly on the graph.The complete line on the graph (i.e. the sequence of measured values or list of successive iterates plotted on a phase space graph) describes a time path or trajectory [1].A trajectory that comes back upon itself to form a closed loop in phase space is called an orbit [2].
An orbit for a system usually indicates that the dynamical system under consideration is conservative.We also note that each plotted point along any trajectory has evolved directly from the preceding point.As we plot each successive point in phase space, the plotted points migrate around.Orbits and trajectories therefore reflect the movement or evolution of the dynamical system.Thus, an orbit or trajectory moves around in the phase space with time.The trajectory is a neat, concise geometric picture that describes part of the system's history.When drawn on a graph, a trajectory must not always be smooth; instead, it can zigzag all over the phase space, mostly for discrete data [3] [4] [5].
The phase space plot is a world that shows the trajectory and its development.Depending on various factors, different trajectories can evolve for the same system.The phase space plot and such a family of trajectories together are a phase space portrait, phase portrait, or phase diagram.
A phase space with plotted trajectories ideally shows the complete set of all possible states that a dynamical system can ever be in.

The Flow Defined by a Differential Equation
We next describe the notion of the flow of a system of differential equations.We begin with the linear system ( ) 0 , 0 , .
The solution to the initial value problem associated with ( 1) is given by ( ) 0 e .

At x t x =
The set of mappings e :  may be regarded as describing the motion of points 0 n x ∈  along trajectories of (1).This set of mappings is called the flow of the linear system (1).

Remark
The mapping e At t φ = satisfies the following basic properties for all For the nonlinear system ( ) we define the flow t φ and show that it satisfies the above basic properties.Subse- quently we introduce the notion of maximal interval of existence ( ) by ( ) 0 I x since the end points a and b of the maximal interval generally depends on 0 x .

Definition
Let E be an open subset of n  and let ( ) is called a flow of the differential Equation (2).t φ is also referred to as the flow of the vector field ( ) f x .3) It can be shown that the basic properties (i)-(iii) of linear flows are also satisfied by nonlinear flows [6].

Remark
4) The following theorem, provides a method of computing derivatives in coordinates.

Theorem
 is differentiable at 0 x , then the partial derivatives , , 1, , , Thus, if f is a differentiable function, the derivative Df is given by the n n × Jaco- bian matrix.

Definition
An equilibrium x * of the system ( ) is called hyperbolicif all eigenvalues of the Jacobian

( )
Df x * have non-zero real part.

The Hartman-Groβman Theorem
The Hartman-Groβman Theorem [7] is another very important result in the local qualitative theory of ordinary differential equations.The theorem shows that near a hyperbolic equilibrium point 0 x , the nonlinear system ( ) has the same qualitative structure as the linear system with ( ) Throughout this section we shall assume that the equilibrium point 0 x has been translated to the origin.

Definition
Two autonomous systems of differential equations such as (7) and ( 8) are said to be topologically equivalent in a neighborhood of the origin or to have the same qualitative structure near the origin if there is a homeomorphism Φ mapping an open set U containing the origin onto an open set V containing the origin which maps trajectories of (7) in U onto trajectories of (8) in V and preserves their orientation by time in the sense that if a trajectory is directed from 1 x to 2 x in U, then its image is directed from ( ) If the homeomorphism Φ preserves the parameteri- zation by time, then the systems ( 7) and ( 8) are said to be topologically conjugate in a neighborhood of the origin.

Example
Consider the linear systems x Ax =  and y By x Rx Φ =, where Then one can easily check that 1 B RAR − = , and letting ( ) is a solution of the first system through 0 x , then At Bt The phase plane portraits of the two systems are shown in Figure 2. It clearly shows that the mapping ( ) x Rx Φ = is simply a rotation through 45˚ and thus it is a ho- meomorphism.

Theorem (Hartman-Groβman)
Let E be an open subset of n  containing the origin, suppose ( ) x I ∈ and i.e.Φ maps trajectories of the nonlinear system ( ) near the origin and preserves the parametrization by time.

Main Results
In [8] Okereke demonstrated very clearly the veracity of the Hartman-Groβman theorem by considering the simulation of the nonlinear and linearized system of ordinary differential equations in terms of their phase portrait analysis.Consider the nonlinear system; The equilibria of the above system is obtained by setting Solving the above equations we obtain the equilibria as (0, 0) and ( ) 3, .
To obtain the linearization at the origin, we begin by computing the Jacobian: ( ) Evaluating the Jacobian at the first equilibrium gives ( ) and therefore the linearization of our system at (0, 0) is =− < , we immediately see that the origin is a sad- dle for the linearized system.Evaluating the Jacobian at the second equilibrium gives ( ) and therefore the linearization of our system at ( ) = > , thus the equilibrium point is a centre for the linearized system.
In the simulation which follows we will consider only the nontrivial equilibrium point ( ) 3, .

MathCAD Simulation
a) The given nonlinear system MathCAD [9] format as follows.Solution matrix is given in Figure 3.
0 : 5 t = − Initial value of independent variable.b) The phase portrait of the linearized system near the origin ( )  can be recast in MathCAD format as follows.Solution matrix is given in Figure 5.
0 : 5 t = − Initial value of independent variable.Figure 5. Solution matrix for the system The solution profiles are depicted in Figures 6(a)-(c).

Observation
The phase portraits of the nonlinear system near ( ) 3, and linearized system about the origin, show stability but not asymptotic stability.This is because the graph is a centre, and as a result, we conclude that the system is conservative.In each case we see that the phase portraits for the nonlinear and linearized system are topologically the same near the equilibrium point ( )

Phase Portrait Analysis for Stability of Third Order ODE
In this section we consider a third order linear equation which is equivalent to the system where a, b, c are all positive constants.
We study the asymptotic properties of the above system with the help of MathCAD simulation.The constants a, b, c are chosen such that ( ) 0 x ab c − > .
0 : 4 t = Initial value of independent variable.
1 : 15 t = Terminal value of independent variable.
Vector of initial values.

The General Nonlinear Third Order ODE
We now consider the more general nonlinear third order ODE given by ( ) ( ) ( ) ( ) We have the following theorem.
x ax bx cx p t + + + =    (14) We have the following result following immediately from the above theorem.

Corollary
Given that Then every solution ( )

Conclusion
In this study, we investigated the stability analysis of certain third order linear and nonlinear ordinary differential equations.We employed the method of phase portrait analysis.We showed, using simulation that the Hartman-Groβman Theorem is verified, for a second order linearized system as an example, approximates the nonlinear system preserving the topological features.In the case of the third order nonlinear system ( ) 0 x ax g x cx + + + =    , we stated appropriate theorems guaranteeing asymptotic stability of solutions.For the linear case where ( ) , our phase portrait analysis shows that under certain conditions on the coefficients as well as the function g , we have asymptotic stability of solutions.
be the solution of the initial value problem (3) defined on its maximal interval of existence ( ) 0 I x .Then for ( ) 0 t I x ∈ , the set of mappings t φ defined by

1 )
We can think of the initial point as being fixed and let curve or trajectory of the system (2) through the point 0x E ∈ .Naturally the mapping is identified simply with its graph in I E × and a trajectory is visualized as a motion along a curve Γ through the point 0x E ∈ of the phase space n  (Figure1(a)).On the other hand, if we think of the point 0 x as varying throughout K E ⊂ , then the flow of the differential Equation (2), as the motion of all points in the set K (Figure1(b)).

Figure 2 .
Figure 2. (a) Phase portrait of x Ax =  ; (b) Phase portrait of y By =  .
the solution of the second system passing through 0 Rx .In other words Φ maps trajectories of the first system onto trajectories of the second system and preserving the parametrization, since e e .
with zero real part.Then there exists a homeomorphism Φ of an open set U containing the origin onto an open set V containing the origin such that for each 0 x U ∈ , there is an open interval 0 I ⊂  containing zero such that for all 0 0 Terminal value of independent variable.
values.The solution profiles are depicted in Figures4(a)-(c).
Figure 4. (a) Trajectory of ( ) 1 x t ; (b) Trajectory of ( ) 2 x t ; (c) Phase portrait of system near Terminal value of independent variable.
, Y t t N D Solution matrix.0 : t S = Independent variable values.
First solution function values.
Second solution function values.
Figure 6.(a) Trajectory of ( ) 1 x t of linearized system; (b) Trajectory of ( ) 2 x t of linearized values.The solution matrix for the above system is given in Figure7, while the solution profiles are depicted in Figures8(a)-(d).

4. 5 . Remark 1 )
We note that (1) and (2) are the well known Routh-Hurwitz conditions[6] for the asymptotic stability of the following third-order homogeneous linear differential equaclearly satisfied and from the simulation (Figures8(a)-(c)) we can see the truth in the limit conditions (15).

Figure 8 (
d) depicts a spiral sink in the simulation, and this further stresses the asymptotic nature of the solutions.