Next Article in Journal
Bidimensional Deterministic Model for Diffusion and Settling of Particles
Next Article in Special Issue
An h-Adaptive Poly-Sinc-Based Local Discontinuous Galerkin Method for Elliptic Partial Differential Equations
Previous Article in Journal
A Class of Janowski-Type (p,q)-Convex Harmonic Functions Involving a Generalized q-Mittag–Leffler Function
Previous Article in Special Issue
Lifting Theorems for Continuous Order-Preserving Functions and Continuous Multi-Utility
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Application of the Block Hybrid Methods to Solve Linear and Non-Linear First Order Differential Equations

Department of Mathematical and Computational Sciences, University of Venda, Thohoyandou 0950, South Africa
Axioms 2023, 12(2), 189; https://doi.org/10.3390/axioms12020189
Submission received: 5 December 2022 / Revised: 3 February 2023 / Accepted: 7 February 2023 / Published: 11 February 2023
(This article belongs to the Special Issue Mathematical Modelling and Applications)

Abstract

:
Block hybrid methods with intra-step points are considered in this study. These methods are implemented to solve linear and nonlinear single and systems of first order differential equations. The stability, convergence, and accuracy of the proposed methods are qualitatively investigated through the absolute and residual error analysis in some selected cases. A number of different numerical examples are tested to demonstrate the efficiency and applicability of the proposed methods. In this study we also implement the proposed methods to solve chaotic systems such as the Glukhvsky–Dolzhansky system, producing very comparable results to those already in the literature.

1. Introduction

Many equations that model physical or dynamical systems in science and engineering are initial value problems. As the name suggests, these are differential equations with prescribed initial conditions that specify the values of the unknown functions at the given specific points in the domain. These initial value problems (IVPs) specify how the systems evolve with time and variations of different parameters (if any), given the initial conditions. There are many initial value problems that we are interested in, that do not have any exact solutions or at least, analytical solutions, that we have figured out yet. In this paper, most of the examples used have exact solutions, since the main objective is the application of the proposed methods.
A plethora of studies have been conducted on block methods since the pioneering work of Shampine and Watts [1] on block implicit one-step methods. Block methods have been developed in order to obtain numerical solutions at more than one point at a time [1]. Brugnano and Trigianite [2] showed that the block methods contain the main and additional methods. Ramos et al. [3] developed a two-step continuous block method with intra-step points through interpolation and collocation. The study observed that the main advantages of the block methods are: (i) overcoming the overlapping pieces of solutions and (ii) that they are self-starting, thus avoiding the use of other methods to obtain starting solutions. Yap et al. [4] proposed the block hybrid collocation method with off-step points. Yap and Ismail [5] proposed the block-hybrid collocation method with three off-step points. Awari [6] studied a class of generalized two-step block hybrid numerical methods. Alharbi and Almatrafi [7] established exact and numerical solitary wave structures to the variant Boussinesq system. Xia [8] studied a generalized Riemann-type hydrodynamical equation and established the existence of a weak solution to the equation in a lower order Sobolev space. Alharbi and Almatrafi [9] presented exact solitary wave and numerical solutions for geophysical KdV equation.
Recently, Ononobgo et al. [10] developed a numerical algorithm for one and two-step hybrid block methods for a numerical solution of first order initial value problems using a method of collocation and Taylor’s series. This then gave a system of nonlinear equations, which was then solved to give a hybrid block method. More detailed work on hybrid block methods can be seen in Gear [11] and Motsa [12]. Yakubu et al. [13] derived the second derivation block hybrid method for the continuous integration of differential systems with the interval of integration.
Ramos and Rufai [14] proposed an implicit two-step block method that incorporates fourth derivatives for solving linear and nonlinear problems. Most recently, Motsa [15] presented a variation of the block methods for integrating systems of initial value problems.
In this study, we apply the newly developed block hybrid linear multi-step methods with off-step points to solve systems of linear and non-linear differential equations. It has been proved that the additional off-step points significantly improve the accuracy of these methods as well as ensuring consistency, zero-stability, and convergence [12]. This happens even if the grid-sizes are kept constant.
Block hybrid methods are used to find the approximate solutions of first order initial value problems of the form
d y d t = f ( t , y ( t ) ) , y ( t 0 ) = y 0 ,
where y 0 is a given initial condition. The one-step method for solving initial value problems ( I V P s ) and parabolic differential equations are considered in the implementations of the block hybrid algorithms. These I V P s are solved over an interval 0 t T , which can be partitioned as 0 = t 0 < t 1 < t 2 < < t N 1 < t N = T . The step size is calculated as h = t n 1 t n , for n = 0 , 1 , 2 , , N 1 . Now, the differential Equation (1) is solved in the N non-overlapping blocks, [ t n , t n + 1 ] , using the known initial condition y ( t n ) for n = 0 , 1 , , N 1 . Then the M 1 points, called the off-steps or intra-step points, are introduced to improve the accuracy of the methods. We remark that M is the number of the collocation points.
The points used in the solution process in each block t n , t n + 1 are t n , t n + p 1 , t n + p 2 , , t t n + p M 1 , t n + p M , where t n + p M = t n + 1 , with p M = 1 . In the block hybrid method, equally spaced intra-step points defined by
t n + p i = t n + i M h = n + i M ,
i = 1 , 2 , , m , are considered [12]. For each interval, say, I n , the first order Equation (1) is solved to obtain a solution, y n + 1 ( t ) . The solution in the first interval for example is denoted by y ( 1 ) ( t 0 ) = y 0 and on the subsequent intervals I n (for n = 1 , 2 , , N 1 ) , as y ( n + 1 ) ( t n ) = y ( n ) ( t n ) .
Motsa [12] showed that the block hybrid method with the intra-step points p 1 , p 2 , , p M 1 has the form
y n + p i = y n + h β i f n + h j = 1 m α i , j f n + p j ,
i = 1 , 2 , , M where α i , j and β i are known constants that depend on the nature of the intra-step points.

2. The Solution Method

Following Motsa [12], the block hybrid method for first order differential equations is based on approximating the exact solution y ( t ) of the linear or non-linear differential Equation (1) by
y ( t ) Y ( t ) = k = 0 M + 1 c n , k ( t t n ) k ,
where c n , i are the unknown coefficients in the [ t n , t n + 1 ] block. With the coefficients c n , i being obtained from a system of M + 2 equations with M + 2 unknowns generated from
Y ˙ n + p i = f ( t n + p i , y n + p i ) , i = 0 , 1 , 2 , , M
Y ( t n ) = c n , 0 = y n , n = 0 , 1 , , N 1 ,
where the dot denotes differentiation with respect to time t . As clearly explained in Motsa [12], the unknown constants c n , k are generated through Mathematica code. The code is as follows:
M = 2 ; p o i n t s = T a b l e p i = i / M , { i , θ , M } ; T a b l e [ t n + p i = ( n + p i ) × h , { i , θ , M } ] ; Y = S u m [ c n , k × ( t t n ) k , { k , θ , M } ] ; e q u a t i o n s = T a b l e [ ( D [ Y , t ] / t t n + p i ) f n + p i , { i , θ , M } ] / / s i m p l y ; u n k n o w n s = T a b l e [ c n , i , { i , θ , M + 1 } ] ; i n i t i a l = c n , θ y n ; A l l e q u a t i o n s = J o i n [ e q u a t i o n s , { i n i t i a l } ] ; s o l u t i o n s = S o l v e [ A l l e q u a t i o n s , u n k n o w n s ]
Upon solving the system of equations that arise from Equation (5), we can have c n , 0 = y n , and c n , 1 = f n for all values M. These c n , k coefficients for k 2 are given for example, when M = 2 ,
c n , 2 = 3 f n + f n + 1 4 f n + 1 2 2 h , c n , 3 = 2 f n + f n + 1 2 f n + 1 2 3 h 2 .
Details of M 3 can be found in the study by Motsa [12]. We then obtain the block hybrid method equation by substituting the expressions for c n , k in the approximation Y ( t ) and by evaluating the result at the collocation points t n + p i for i = 1 , 2 , , M . When the nodes are equally spaced, say, { 0 , 1 2 , 1 } with M = 2 , we get
y n + 1 2 = y n + 1 24 h 5 f n f n + 1 + 8 f n + 1 2 , y n + 1 = y n + 1 6 h f n + f n + 1 + 4 f n + 1 2 .
The above equation can be written in matrix form as follows:
y n + 1 2 y n + 1 = y n y n + h 5 f n 24 f n 6 + h 1 3 1 24 2 3 1 6 f n + 1 2 f n + 1 .
The block hybrid method with the intra-step points p 1 , p 2 , p 3 , , p m 1 , which in general is given by (3).
Equation (3), can be written in matrix as
Y n + m = Y n + h A m F n + m + h B m F n ,
where the respective coefficient matrices are defined as
A m = α 1 , 1 α 1 , 2 , , α 1 , m α 2 , 1 α 2 , 2 α 2 , m α m , 1 α m , 2 α m , m , B m = β 1 0 0 0 β 2 0 0 0 β m ,
and the column vectors are:
Y n + m = y n + p 1 , y n + p 2 , , y n + p m T , F n + m = f n + p 1 , f n + p 2 , , f n + p m T ,
Y n = y n , y n , , y n T , F n = f n , f n , , f n T .
The coefficient matrices A m and B m , when M = 1 , 2 , 3 , 4 , 5 are
A 1 = 1 2 , B 1 = 1 2 ,
A 2 = 1 3 1 24 2 3 1 6 , B 2 = 1 24 1 6 1 6 1 3 ,
A 3 = 19 72 5 72 1 72 4 9 1 9 0 3 8 3 8 1 8 , B 3 = 1 72 13 216 17 216 4 27 1 9 4 27 7 24 5 24 5 24 ,
A 4 = 323 1440 11 120 53 1440 19 2880 31 90 1 15 1 90 1 360 51 160 9 40 21 160 3 320 16 40 2 15 16 45 7 90 , B 4 = 71 2880 17 640 19 480 263 5760 11 90 2 45 13 180 31 360 9 40 51 640 33 320 81 640 59 180 7 60 7 60 31 180 ,
A 5 = 1427 7200 133 1200 241 3600 173 7200 3 800 43 150 7 225 7 225 1 75 1 450 219 800 57 400 57 400 21 800 3 800 64 225 8 75 64 225 14 225 0 25 96 25 144 25 144 25 96 19 288 , B 5 = 187 7200 113 9000 809 36,000 61 2250 1073 36,000 229 2250 4 2250 89 2250 113 2250 127 2250 723 4000 63 2000 9 160 147 2000 333 4000 292 1125 52 1125 76 1125 22 225 124 1125 487 1440 11 180 113 1440 47 360 193 1440 .
To validate the block hybrid method results obtained in this study, we consider the equation:
y = y t 2 + 1 , 0 t 2 , y ( 0 ) = 0.5 ,
with the exact solution
y ( t ) = ( t + 1 ) 2 0.5 e t .
Table 1 clearly shows the high accuracy of the block hybrid method even with very high step sizes. The results from the fourth order Ruge–Kutta method were generated with the step size h = 0.02 , whilst those of the BHM were generated with the step size h = 0.2 .
We note that the block hybrid method gives a consistently superior result, which can be clearly observed in Table 2.

Local Truncation Error

Given a function y ( t ) , which is sufficiently differentiable, we can express the block hybrid methods following Motsa [17], in terms of the linear i , with
i [ y ( t n ) ; h ] = y ( t n + h p i ) y ( t n ) h β i y ( t n ) h j = 1 M α i , j y ( t n + h p j ) , i = 1 , 2 , , M .
We can then expand the terms of y ( t n + h p i ) and y ( t n + h p i ) using Taylor series about t n to obtain
i [ y ( t n ) ; h ] = C i , 0 y ( t n ) + C i , 1 h y ( t n ) + C i , 2 h 2 y ( t n ) + + C i , q h q y ( q ) ( t n ) ,
where C i , 0 , C i , 1 , C i , q are constants. The method is said to be an order of q if
C 0 ^ = C 1 ^ = = C q + 1 ^ = 0 and C q + 2 ^ 0 .
Upon expanding Equation (16) using Taylor series we get:
i [ y ( t n ) ; h ] = k = 1 m p i k k ! h k y ( k ) ( t n ) h β i y ( t n ) k = 1 m p i k k ! j = 1 M α i , j p j k 1 y ( k ) ( t n ) + O ( h m + 1 ) ,
where m is a positive integer. The above Equation (18) can be expanded to obtain:
i [ y ( t n ) ; h ] = k = 1 M + 1 p i k k ! h k y ( k ) ( t n ) h β i y ( t n ) k = 1 M + 1 p i k k ! j = 1 M α i , j p j k 1 y ( k ) ( t n ) + h M + 2 ( M + 2 ) ! p i M + 2 ( M + 2 ) j = 1 M α i , j p j M + 1 y ( M + 2 ) ( t n ) + O ( h M + 3 )
Motsa [17] showed that from numerical evaluation the sum of the first three terms of Equation (19) is zero. Therefore the truncation error for the block hybrid with M + 1 nodes in the interval [ t n , t n + 1 ] is
i [ y ( t n ) ; h ] = h M + 2 ( M + 2 ) ! p i M + 2 ( M + 2 ) j = 1 M α i , j p j M + 1 y ( M + 2 ) ( t n ) + O ( h M + 3 ) .
Thus, the order of the block hybrid method for say 3 n o d e s will be 5, which is the least order of applying this method. We also remark that as M increases so does the order of convergence. This then increases the accuracy of the method, as will be observed later in this study.

3. Applications

Randomly selected numerical experiments are studied to show the performance and viability of the proposed methods.

3.1. Numerical Linear Examples

We consider a general linear first order equation of the form:
y ˙ = f ( t , y ) = ψ ( t ) + ϕ ( t ) y , t 0 < t < t T ,
where ψ and ϕ are known functions of t . Following Motsa [12], the block hybrid scheme for solving linear equations is given by
y n + p i = y n + h β i ψ n + ϕ n y n + h j = 1 n α i , j ψ n + p j + ϕ n + p j y n + p j , i = 1 , 2 , 3 , , M .
The above Equation (22) can be expressed in matrix form as:
A n Y n = B n ,
with
A n = 1 0 0 0 1 0 0 0 1 h A m ϕ n + p 1 0 0 0 ϕ n + p 2 0 0 0 ϕ n + p m .
B n = y n y n y n + h A m ψ n + p 1 ψ n + p 2 ψ n + p m + h B m ψ n + ϕ n y n ψ n + ϕ y n ψ n + ϕ n y n ,
and
Y n = y n + p 1 , y n + p 2 , , y n + p M T .

3.1.1. Linear Example 1

In this section we are going to randomly selected examples to demonstrate the efficiency of the block-hybrid method. First, we consider the following example;
y = t e 3 t 2 y , 0 t 1 , y ( 0 ) = 0 ,
with an exact solution given as
y ( t ) = 1 5 t e 3 t 1 25 e 3 t + 1 25 e 2 t .
In this case, we have ψ ( t ) = t e 3 t and ϕ ( t ) = 2 . As can be clearly observed in Figure 1, the accuracy of the method greatly increases as M increases. Figure 2 shows that there is an excellent agreement between the exact and the approximate solutions. This can be seen also on the errors of order 10 15 in Figure 1.

3.1.2. Linear Example 2

The example below was randomly selected to demonstrate the robustness of the proposed method.
y = 2 2 t y t 2 + 1 , 0 t 1 , y ( 0 ) = 0 ,
with exact solution given as
y ( t ) = 2 t + 1 t 2 + 1 .
In example 2, ψ ( t ) = 2 t 2 + 1 and ϕ ( t ) = 2 t t 2 + 1 . In Figure 3, we observe that the accuracy is up to the 10 13 error levels, which is highly commendable. This is observed further in Figure 4, which shows an excellent agreement between the exact and approximate solutions.

3.1.3. Linear Example 3

The example below was chosen to illustrate the efficiency of the BHM.
y = 4 t y t 2 + 1 + 2 ( 1 + t ) 2 , 0 t 6 , y ( 0 ) = 1 ,
with exact solution given as
y ( t ) = 2 t + 1 ( t 2 + 1 ) 2 .
In this example, ψ ( t ) = 2 ( 1 + t ) 2 and ϕ ( t ) = 4 t 1 + t 2 . The exact and approximate solutions are depicted in Figure 5. We observe that there is an excellent agreement between the two sets of solutions. This can also be perceived in Figure 6 which depicts a very high level of accuracy.

3.1.4. Linear Example 4

This example was randomly selected to illustrate the robustness of the proposed method.
y = 2 y + t 3 e 2 t 0 t 1 , y ( 0 ) = 1 ,
with exact solution given as
y ( t ) = e 2 t 4 ( t 4 + 4 ) .
In this example, ψ ( t ) = t 3 e 2 t and ϕ ( t ) = 2. Figure 7 and Figure 8 depict, respectively, the errors and the comparisons of the exact and approximate solutions. Moreover, we observe the robustness of the numerical method in solving these types of first order initial value problems.

3.2. Numerical Nonlinear First Order Differential Equations

In this subsection, we examine some few selected non-linear first order differential equations to demonstrate the strength of the methods under discussion. The quasi-linearization method is used to linearize the equations first. The non-linear first order differential equations are first linearized to enable us to apply the BHMs. We consider a non-linear first order differential equation of the form
y = ( t ) y + ( t , y ) ,
where, ( t ) is a known function of t, and ( t , y ) is a nonlinear function of y . In this study, we will use the quasi-linearization (QLM) iteration method, developed by Bellman and Kalaba [18]. The QLM approach is based on Taylor series expansion of the nonlinear term ( t , y ) . In this method, we assume that the difference between the current and previous iteration ( y r + 1 y r ) is small. We have, ( t , y ) ( t , y r ) + y ( y r + 1 y r ) . This gives the linearized approximation Equation (35)
y r + 1 = ( t ) y r + 1 + ( t , y r ) + y y r + 1 y r .
We then apply the block hybrid method scheme with
ϕ ( t ) = ( t ) + y , ψ ( t ) = ( t , y r ) y r y .
We discuss some few randomly selected examples to demonstrate the strength of the block-hybrid method.

3.2.1. Nonlinear Example 1 (Riccati Equation)

This example was conveniently selected to demonstrate the robustness of the proposed method. The riccati equation though nonlinear its exact solution can easily be found.
y = 1 + y t + y 2 t 2 , 1 t 3 , y ( 1 ) = 0 ,
with the exact solution given as
y ( t ) = t tan ( ln t ) .
Following Equation (36), we deduce that, ( t ) = 1 t and ( t , y ) = 1 + y 2 t 2 with y = 2 y t 2 . Therefore,
y ˙ r + 1 = 1 t + 2 y r t 2 y r + 1 + 1 y r 2 t 2 ,
with ϕ ( t ) = 1 t + 2 y r t 2 , ψ ( t ) = 1 y r 2 t 2 . We then apply the block-hybrid method to solve the resultant linear equations and the results as depicted in Figure 9 and Figure 10. Again, we see that the current method gives an accuracy as high as 10 14 , which is highly commendable.

3.2.2. Nonlinear Example 2

This example was randomly selected to demonstrate the robustness of the proposed BHM.
y = y 3 + t y 4 6 , 0 t 1 , y ( 0 ) = 2 ,
with the exact solution given as
y ( t ) = 2 ( 4 t 4 + 5 e t ) 1 3 .
We have ( t ) = 1 3 and ( t , y r ) = t y 4 6 , with y r = 2 t y r 3 3 .
Hence,
y ˙ r + 1 = 1 3 + 2 t y r 3 3 y r + 1 t y r 4 2 .
We have, ϕ ( t ) = 1 3 + 2 t y r 3 3 and ψ ( t ) = t y r 4 2 . Upon application of the BHMs on Equation (43) we obtain the results, which are depicted in Figure 11 and Figure 12.

3.3. Nonlinear Systems of First Order Equations

Lastly, we apply the block-hybrid methods to non-linear systems of N equations of the form:
y ˙ 1 = f 1 ( t , y 1 , y 2 , , y N ) = 1 ( t , y 2 , y 3 , , y N ) y 1 + 1 ( t , y 1 , y 2 , , y N ) , y ˙ 2 = f 2 ( t , y 1 , y 2 , , y N ) = 2 ( t , y 1 , y 3 , , y N ) y 2 + 2 ( t , y 1 , y 2 , , y N ) , y k ˙ = f k ( t , y 1 , y 2 , , y N ) = k ( t , y 1 , y 3 , y k 1 , y k + 1 , , y N ) y k + 2 ( t , y 1 , y 2 , , y N ) , y N ˙ = f N ( t , y 1 , y 2 , , y N ) = N ( t , y 1 , y 3 , , y N 1 ) y N + N ( t , y 1 , y 2 , , y N ) ,
where k is the component of the non-linear function that is a coefficient to y k in the k - th equation and k ( t ) is the remaining component which may or may not be a non-linear function for k = 1 , 2 , , N . At each iteration, denoted by r + 1 , the decoupling iterative scheme takes the form:
y ˙ k , r + 1 = k ( t , y 1 , r + 1 , y 2 , r + 1 , , y k 1 , r + 1 , y k + 1 , r , y N , r ) y k , r + 1 + k ( t , y 1 , r + 1 , y 2 , r + 1 , , y k 1 , r + 1 , y k + 1 , r , y N , r ) ,
for each k = 1 , 2 , , N . The block-hybrid method for a linear first order Equation (37) can now be applied to the above Equation (45) with ϕ k = k ( t , y 1 , r + 1 , y 2 , r + 1 , , y k 1 , r + 1 , y k + 1 , r , y N , r ) and ψ k = k ( t , y 1 , r + 1 , y 2 , r + 1 , , y k 1 , r + 1 , y k + 1 , r , y N , r ) , for each k = 1 , 2 , , N . The scheme is now called the Relaxation Block Hybrid Method (Motsa [12]), for solving these systems of equations. We then develop a quasilinearization scheme by applying the linearization sequentially in y k only to obtain
y ˙ k , r + 1 = k ( t , y 1 , r + 1 , y 2 , r + 1 , , y k 1 , r + 1 , y k + 1 , r , y N , r ) y k , r + 1 + k ( t , y 1 , r + 1 , y 2 , r + 1 , , y k 1 , r + 1 , y k + 1 , r , y N , r ) + k y k ( y r + 1 y k , r )
The BHMs are then applied with ϕ k = k + k y k , ψ k = k y k , r k y k . The BHM implemented with the quasilinearization (46) can be referred to as the Local Quazilinearization Block Hybrid Method (LQBHM) [12].

3.3.1. Nonlinear Systems, Example 1

We consider the following Lorenz system given by
y 1 ˙ = σ ( y 2 y 1 ) , y 2 ˙ = y 1 y 3 + ρ y 1 y 2 , y 3 ˙ = y 1 y 2 c y 3 ,
with initial conditions y 1 ( 0 ) = 1 , y 2 ( 0 ) = 5 , y 3 ( 0 ) = 10 . The BHM parameters:
ϕ 1 = σ , ψ 1 = σ y 2 , ϕ 2 = 1 , ψ 2 = y 1 y 3 + ρ y 1 , ϕ 3 = c , ψ 3 = y 1 y 2 .
We select some few graphical examples to demonstrate how powerful the proposed method is. We can clearly observe in Figure 13, Figure 14 and Figure 15 that the accuracy of the method increases as the number of iterations increases. An accuracy up to 10 14 is achieved by this method. In Figure 16, we observe the chaotic nature of the solutions. We observe that solutions exhibit irregular oscillations that persist as t , but never repeat exactly. The motions are further aperiodic.
Figure 17, Figure 18, Figure 19, Figure 20, Figure 21 and Figure 22 depict phase portraits of the system when varying the system parameters. In this study we are not going to go into the physics and explanations of the portraits. Our focus is to generate accurate solutions, since the Lorenz equations have been the subject of hundreds of research papers, see for example Fang and Hao [19], and Hao et al. [20].
In Figure 23, Figure 24 and Figure 25, we display the 3 D evolutions of the three trajectories when varying the three parameters. Fundamentally, the trajectories seem to exhibit similar shapes when varying the three parameters.
In Figure 26, Figure 27 and Figure 28, we display the phase portraits of y 2 and y 3 when varying the three parameters.

3.3.2. Nonlinear Systems, Example 2

Further, we consider the Glukhovsky–Dolzhanksy system (see Garashchuk et al. [21]) given by
y 1 ˙ = σ ( y 1 y 2 ) a y 2 y 3 , y 2 ˙ = r y 1 y 2 y 1 y 3 , y 3 ˙ = b y 3 + y 1 y 2 ,
where σ , r , b are the physical parameters. This system has an additional nonlinear term compared to the Lorenz system. This extra nonlinear term leads to essential differences in the analytical structure and dynamics of the system. This system describes the following physical processes: convective fluid motion in a ellipsoidal rotating cavity, a rigid body rotation in a resisting medium, the forced motion of a gyrostat, and a convective motion in harmonically oscillating horizontal fluid. Selected phase portraits are depicted in Figure 29, Figure 30, Figure 31 and Figure 32. Figure 29 shows that very chaotic patterns are initially exhibited and thereafter the flows exhibit somehow steady states as t .

3.3.3. Nonlinear Systems Example 3

Consequently, the study considers a dissipative chaotic system with no equilibrium given by
y 1 ˙ = y 2 , y 2 ˙ = y 1 + y 2 y 3 , y 3 ˙ = y 1 2 4 y 2 2 + 1 ,
This dissipative structure/system is characterized by the spontaneous appearance of symmetry breaking (anisotropy) and the formation of complex and chaotic structure Brogliato et al. [22]. In these structures, interacting particles exhibit long range correlations. The current system (49), is similar to a thermodynamically open system that is operating out of the thermodynamic equilibrium in an area with which it exchanges energy and matter. Figure 33, Figure 34, Figure 35, Figure 36 and Figure 37 display the chaotic behaviors of these dissipative chaotic systems.

3.3.4. Nonlinear Systems Example 4

Consequently, we demonstrate the applicability of the proposed method using the five-dimensional non-dissipative Lorenz model given as
y 1 ˙ = σ ( y 2 y 1 ) , y 2 ˙ = y 1 y 3 + ρ y 1 y 2 , y 3 ˙ = y 1 y 2 y 1 y 4 b y 3 , y 4 ˙ = y 1 y 3 2 y 1 y 5 c y 4 , y 5 ˙ = 2 y 1 y 4 4 b y 5 .
Figure 38, Figure 39, Figure 40, Figure 41 and Figure 42 depict the errors, and we can clearly see that the level of accuracy for this proposed method is very high. The last few selected Figure 43, Figure 44, Figure 45, Figure 46, Figure 47, Figure 48, Figure 49, Figure 50, Figure 51, Figure 52, Figure 53 and Figure 54 display the selected phase portraits when varying the parameters.

4. Conclusions

This study proposed the application of newly developed block hybrid linear multi-step methods with off-step points for solving linear and nonlinear single and systems of differential equations. Numerical results for the current methods are excellent and compare very well with the exact solutions. The high levels of convergence indicate that the HBMs are very good candidates to solve high order systems of nonlinear systems of equations. We consequentially remark that the study observed that the numerical approximations converged quickly after very few iterations, even with very many collocation points and small step sizes. We also observed that the block hybrid methods are far superior to some classical numerical methods such as the Runge–Kutta methods. This great accuracies of the block hybrid methods and their user friendliness will go a long way in their applications in more complex models. Unfortunately, at the moment BHMs can only be applied to linear equations. Nonlinear equations have to be linearized first.

Funding

The research was funded by the University of Venda.

Data Availability Statement

Data is contained within the article.

Acknowledgments

The author would like to acknowledge the reviewers for their respective suggestions, Sandile Motsa for his advice and input and Ndivhuwo Ndou for his input.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Shampine, L.F.; Watts, H.A. Block Implicit One-Step Methods. Math. Comput. 1969, 23, 731–740. [Google Scholar] [CrossRef]
  2. Brugnano, L.; Trigiante, D. Solving Differential Problems by Multistep Initial and Boundary Value Methods; Gordon and Breach Science Publishers: Amsterdam, The Netherlands, 1998. [Google Scholar]
  3. Ramos, H.; Kalogiratou, Z.; Monovasilis, T.H.; Simos, T.E. An optimized two-step hybrid block method for solving general second order initial-value problems. Numer. Algorithms 2016, 72, 1089–1102. [Google Scholar] [CrossRef]
  4. Yap, L.K.; Ismail, F.; Senu, N. An Accurate Block Hybrid Collocation Method for Third Order Ordinary Differential Equations. J. Appl. Math. 2014, 2014, 549597. [Google Scholar] [CrossRef]
  5. Yap, L.K.; Ismail, F. Block Hybrid Collocation Method with Application to Fourth Order Differential Equations. Math. Probl. Eng. 2015, 2015, 561489. [Google Scholar] [CrossRef]
  6. Awari, Y.S. Some generalized two-step block hybrid Numerov method for solving general second order ordinary differential equations without predictors. Sci. World J. 2015, 12, 12–18. [Google Scholar]
  7. Albarbi, A.R.; Almatrafi, M.B. Exact and Numerical SolitaryWave Structures to the Variant Boussinesq System. Symmetry 2020, 12, 1473. [Google Scholar] [CrossRef]
  8. Xia, S. Existence of a weak solution to a generalized Riemanntype hydrodynamical equation. In Applicable Analysis; Taylor & Francis: Abingdon, UK, 2022. [Google Scholar] [CrossRef]
  9. Albarbi, A.R.; Almatrafi, M.B. Exact solitary wave and numerical solutions for geophysical KdV equation. J. King Saud Univ. Sci. 2022, 34, 102087. [Google Scholar] [CrossRef]
  10. Ononogbo, C.B.; Airemen, I.E.; Ezurike, U.J. Numerical Algorithm for One and Two-Step Hybrid Block Methods for the Solution of First Order Initial Value Problems in Ordinary Differential Equations. Appl. Eng. 2022, 6, 13–23. [Google Scholar] [CrossRef]
  11. Gear, C.N. Hybrid methods for initial value problems in Ordinary Differential Equations. SIAM J. Numer. Anal. 1964, 2, 69–86. [Google Scholar] [CrossRef]
  12. Motsa, S.S. Block hybrid methods. In Proceedings of the 13th Annual Workshop on Computational Mathematics and Modelling, University of KwaZulu-Natal, Pietermaritzburg Campus, Durban, South Africa, 5–9 July 1964. [Google Scholar]
  13. Yakubu, D.G.; Shokri, A.; Kumleng, G.M.; Marian, D. Second Derivative Block Hybrid Methods for the Numerical Integration of Differential Systems. Fractal Fract. 2022, 6, 386. [Google Scholar] [CrossRef]
  14. Ramos, H.; Rufai, M.A. A two-step hybrid block method with fourth derivatives for solving third-order boundary value problems. J. Comput. Appl. Math. 2022, 404, 113419. [Google Scholar] [CrossRef]
  15. Motsa, S. Overlapping Grid-Based Optimized Single-Step Hybrid Block Method for Solving First-Order Initial Value Problems. Algorithms 2022, 15, 427. [Google Scholar] [CrossRef]
  16. Burden, R.L.; Faires, J.D. Numerical Analysis, 9th ed.; Brooks/Cole, Cengage Learning: San Francisco, CA, USA, 2011; pp. 302–314. [Google Scholar]
  17. Motsa, S.S. Hybrid block methods for IVPs using Mathematica. In Proceedings of the 14th Annual Workshop on Computational Mathematics and Modelling, University of KwaZulu-Natal, Pietermaritzburg Campus, Durban, South Africa, 4–8 July 2022. [Google Scholar]
  18. Bellman, R.E.; Kalaba, R.E. Quasilinearization and Nonlinear Boundary-Value Problems; Elsevier: New York, NY, USA, 1965. [Google Scholar]
  19. Fang, H.P.; Hao, B.L. Symbolic dynamics of the Lorenz equations. Chaos Solitons Fractals 1966, 7, 217–246. [Google Scholar] [CrossRef]
  20. Hao, B.; Liu, J.; Zheng, W. Symbolic dynamics analysis of the Lorenz equations. Phys. Rev. 1998, 57, 5378. [Google Scholar] [CrossRef]
  21. Garashchuk, I.R.; Kudryashov, N.A.; Sinelshchikov, D.I. On the analytical properties and some exact solutions of the Glukhovsky-Dolzhansky system. J. Phys. Conf. Ser. 2017, 788, 012013. [Google Scholar] [CrossRef]
  22. Brogliato, B.; Lozano, R.; Maschke, B.; Egeland, O. Dissipative Systems Analysis and Control, 2nd ed.; Theory and Applications; Springer: London, UK, 2007. [Google Scholar]
Figure 1. Errors for example 1.
Figure 1. Errors for example 1.
Axioms 12 00189 g001
Figure 2. Graphs of the exact and approximate solutions for example 1.
Figure 2. Graphs of the exact and approximate solutions for example 1.
Axioms 12 00189 g002
Figure 3. Errors for example 2.
Figure 3. Errors for example 2.
Axioms 12 00189 g003
Figure 4. Graphs of the exact and approximate solutions for example 2.
Figure 4. Graphs of the exact and approximate solutions for example 2.
Axioms 12 00189 g004
Figure 5. Graphs of the exact and approximate solutions for example 3.
Figure 5. Graphs of the exact and approximate solutions for example 3.
Axioms 12 00189 g005
Figure 6. Errors for example 3.
Figure 6. Errors for example 3.
Axioms 12 00189 g006
Figure 7. Errors for example 4.
Figure 7. Errors for example 4.
Axioms 12 00189 g007
Figure 8. Graphs of the exact and approximate solutions for example 4.
Figure 8. Graphs of the exact and approximate solutions for example 4.
Axioms 12 00189 g008
Figure 9. Errors for nonlinear example 1 as M is increasing.
Figure 9. Errors for nonlinear example 1 as M is increasing.
Axioms 12 00189 g009
Figure 10. Graphs of the exact and approximate solutions for example 1.
Figure 10. Graphs of the exact and approximate solutions for example 1.
Axioms 12 00189 g010
Figure 11. Errors for non-linear example 2.
Figure 11. Errors for non-linear example 2.
Axioms 12 00189 g011
Figure 12. Graphs of the exact and approximate solutions for example 2.
Figure 12. Graphs of the exact and approximate solutions for example 2.
Axioms 12 00189 g012
Figure 13. Errors for the Lorenz system.
Figure 13. Errors for the Lorenz system.
Axioms 12 00189 g013
Figure 14. Errors for the Lorenz system.
Figure 14. Errors for the Lorenz system.
Axioms 12 00189 g014
Figure 15. Errors for the Lorenz system.
Figure 15. Errors for the Lorenz system.
Axioms 12 00189 g015
Figure 16. Plots of the solutions of y 1 , y 2 , y 3 .
Figure 16. Plots of the solutions of y 1 , y 2 , y 3 .
Axioms 12 00189 g016
Figure 17. Phase portraits of y 1 and y 2 when varying ρ .
Figure 17. Phase portraits of y 1 and y 2 when varying ρ .
Axioms 12 00189 g017
Figure 18. Phase portraits of y 1 and y 2 when varying σ .
Figure 18. Phase portraits of y 1 and y 2 when varying σ .
Axioms 12 00189 g018
Figure 19. Phase portraits of y 1 and y 2 when varying c.
Figure 19. Phase portraits of y 1 and y 2 when varying c.
Axioms 12 00189 g019
Figure 20. Phase portraits of y 1 and y 3 when varying ρ .
Figure 20. Phase portraits of y 1 and y 3 when varying ρ .
Axioms 12 00189 g020
Figure 21. Phase portraits of y 1 and y 3 when varying σ .
Figure 21. Phase portraits of y 1 and y 3 when varying σ .
Axioms 12 00189 g021
Figure 22. Phase portraits of y 1 and y 3 when varying c.
Figure 22. Phase portraits of y 1 and y 3 when varying c.
Axioms 12 00189 g022
Figure 23. Phase portraits of y 1 , y 2 and y 3 when varying ρ .
Figure 23. Phase portraits of y 1 , y 2 and y 3 when varying ρ .
Axioms 12 00189 g023
Figure 24. Phase portraits of y 1 , y 2 and y 3 when varying σ .
Figure 24. Phase portraits of y 1 , y 2 and y 3 when varying σ .
Axioms 12 00189 g024
Figure 25. Phase portraits of y 1 , y 2 and y 3 when varying c.
Figure 25. Phase portraits of y 1 , y 2 and y 3 when varying c.
Axioms 12 00189 g025
Figure 26. Phase portraits of y 2 and y 3 when varying ρ .
Figure 26. Phase portraits of y 2 and y 3 when varying ρ .
Axioms 12 00189 g026
Figure 27. Phase portraits of y 2 and y 3 when varying σ .
Figure 27. Phase portraits of y 2 and y 3 when varying σ .
Axioms 12 00189 g027
Figure 28. Phase portraits of y 2 and y 3 when varying c.
Figure 28. Phase portraits of y 2 and y 3 when varying c.
Axioms 12 00189 g028
Figure 29. Solution profiles of the Glukhovsky–Dolzhanksy system of equations.
Figure 29. Solution profiles of the Glukhovsky–Dolzhanksy system of equations.
Axioms 12 00189 g029
Figure 30. Phase portraits of y 1 and y 2 .
Figure 30. Phase portraits of y 1 and y 2 .
Axioms 12 00189 g030
Figure 31. Phase portraits of y 1 and y 3 .
Figure 31. Phase portraits of y 1 and y 3 .
Axioms 12 00189 g031
Figure 32. Phase portraits of y 2 and y 3 .
Figure 32. Phase portraits of y 2 and y 3 .
Axioms 12 00189 g032
Figure 33. Solution profiles of the chaotic system.
Figure 33. Solution profiles of the chaotic system.
Axioms 12 00189 g033
Figure 34. Phase portraits of y 1 and y 2 .
Figure 34. Phase portraits of y 1 and y 2 .
Axioms 12 00189 g034
Figure 35. Phase portraits of y 1 and y 3 .
Figure 35. Phase portraits of y 1 and y 3 .
Axioms 12 00189 g035
Figure 36. Phase portraits of y 2 and y 3 .
Figure 36. Phase portraits of y 2 and y 3 .
Axioms 12 00189 g036
Figure 37. Phase portraits of y 1 , y 2 and y 3 .
Figure 37. Phase portraits of y 1 , y 2 and y 3 .
Axioms 12 00189 g037
Figure 38. Convergence graph of y 1 .
Figure 38. Convergence graph of y 1 .
Axioms 12 00189 g038
Figure 39. Convergence graph of y 2 .
Figure 39. Convergence graph of y 2 .
Axioms 12 00189 g039
Figure 40. Convergence graph of y 3 .
Figure 40. Convergence graph of y 3 .
Axioms 12 00189 g040
Figure 41. Convergence graph of y 4 .
Figure 41. Convergence graph of y 4 .
Axioms 12 00189 g041
Figure 42. Convergence graph of y 5 .
Figure 42. Convergence graph of y 5 .
Axioms 12 00189 g042
Figure 43. Phase portraits of y 1 and y 2 when varying ρ .
Figure 43. Phase portraits of y 1 and y 2 when varying ρ .
Axioms 12 00189 g043
Figure 44. Phase portraits of y 1 and y 3 when varying ρ .
Figure 44. Phase portraits of y 1 and y 3 when varying ρ .
Axioms 12 00189 g044
Figure 45. Phase portraits of y 1 and y 4 when varying ρ .
Figure 45. Phase portraits of y 1 and y 4 when varying ρ .
Axioms 12 00189 g045
Figure 46. Phase portraits of y 1 and y 5 when varying ρ .
Figure 46. Phase portraits of y 1 and y 5 when varying ρ .
Axioms 12 00189 g046
Figure 47. Phase portraits of y 1 and y 2 when varying σ .
Figure 47. Phase portraits of y 1 and y 2 when varying σ .
Axioms 12 00189 g047
Figure 48. Phase portraits of y 1 and y 3 when varying σ .
Figure 48. Phase portraits of y 1 and y 3 when varying σ .
Axioms 12 00189 g048
Figure 49. Phase portraits of y 2 and y 3 when varying σ .
Figure 49. Phase portraits of y 2 and y 3 when varying σ .
Axioms 12 00189 g049
Figure 50. Phase portraits of y 2 and y 4 when varying σ .
Figure 50. Phase portraits of y 2 and y 4 when varying σ .
Axioms 12 00189 g050
Figure 51. Phase portraits of y 1 and y 5 when varying σ .
Figure 51. Phase portraits of y 1 and y 5 when varying σ .
Axioms 12 00189 g051
Figure 52. Phase portraits of y 1 and y 4 when varying σ .
Figure 52. Phase portraits of y 1 and y 4 when varying σ .
Axioms 12 00189 g052
Figure 53. Phase portraits of y 3 and y 5 when varying σ .
Figure 53. Phase portraits of y 3 and y 5 when varying σ .
Axioms 12 00189 g053
Figure 54. Phase portraits of y 4 and y 5 when varying σ .
Figure 54. Phase portraits of y 4 and y 5 when varying σ .
Axioms 12 00189 g054
Table 1. Maximum Errors for HBM (with h = 0.2 ) and RK4 (with h = 0.02 ).
Table 1. Maximum Errors for HBM (with h = 0.2 ) and RK4 (with h = 0.02 ).
tExact SolutionBHM-ErrorsRK4-Errors
00.500
0.20.829298620919917 2 × 10 15 5.7 × 10 10
0.41.214087651179369 4 × 10 15 1.2 × 10 9
0.61.648940599804753 7 × 10 15 2.0 × 10 9
0.82.127229535753778 1.2 × 10 14 2.9 × 10 9
1.02.640859085770495 1.8 × 10 14 3.9 × 10 9
1.23.179941538631752 2.5 × 10 14 5 × 10 9
1.43.732400016577701 5.2 × 10 14 6.3 × 10 9
1.64.283483787802496 5.2 × 10 14 7.8 × 10 9
1.84.815176267793600 7.5 × 10 14 9.5 × 10 9
2.05.305471950534773 9.9 × 10 14 1.1 × 10 8
Table 2. Comparisons of errors of the results generated by the BHM and those obtained by Burden and Faires [16].
Table 2. Comparisons of errors of the results generated by the BHM and those obtained by Burden and Faires [16].
t 0 ExactBHMErrorAdams 4 th Order Predictor-Corrector MethodError (Burden and Faires [16])
0.00.50000000.500000000.50000000
0.20.82929860.829298600.82929330.0000053
0.41.21408771.214087701.214077620.0000114
0.61.64894061.648940601.64892200.0000186
0.82.12722952.127229502.12720560.0000239
1.02.64085912.640859102.64082860.0000305
1.23.17994153.179941503.17990260.0000389
1.43.73240003.732400003.73235050.0000495
1.64.28348384.283483804.28342080.0000630
1.84.81517634.815176304.81509640.0000799
2.05.30547205.305472005.30537070.0001013
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.

Share and Cite

MDPI and ACS Style

Shateyi, S. On the Application of the Block Hybrid Methods to Solve Linear and Non-Linear First Order Differential Equations. Axioms 2023, 12, 189. https://doi.org/10.3390/axioms12020189

AMA Style

Shateyi S. On the Application of the Block Hybrid Methods to Solve Linear and Non-Linear First Order Differential Equations. Axioms. 2023; 12(2):189. https://doi.org/10.3390/axioms12020189

Chicago/Turabian Style

Shateyi, Stanford. 2023. "On the Application of the Block Hybrid Methods to Solve Linear and Non-Linear First Order Differential Equations" Axioms 12, no. 2: 189. https://doi.org/10.3390/axioms12020189

APA Style

Shateyi, S. (2023). On the Application of the Block Hybrid Methods to Solve Linear and Non-Linear First Order Differential Equations. Axioms, 12(2), 189. https://doi.org/10.3390/axioms12020189

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop