Next Article in Journal
Pricing Contingent Claims in a Two-Interest-Rate Multi-Dimensional Jump-Diffusion Model via Market Completion
Previous Article in Journal
Comparative Study of Monte Carlo Simulation and the Deterministic Model to Analyze Thermal Insulation Costs
Previous Article in Special Issue
Convergence and Stability Improvement of Quasi-Newton Methods by Full-Rank Update of the Jacobian Approximates
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Block Hybrid Method with Equally Spaced Grid Points for Third-Order Initial Value Problems

by
Salma A. A. Ahmedai Abd Allah
,
Precious Sibanda
,
Sicelo P. Goqo
*,
Uthman O. Rufai
,
Hloniphile Sithole Mthethwa
and
Osman A. I. Noreldin
School of Mathematics, Statistics and Computer Science, University of KwaZulu-Natal, Pietermaritzburg 3209, South Africa
*
Author to whom correspondence should be addressed.
AppliedMath 2024, 4(1), 320-347; https://doi.org/10.3390/appliedmath4010017
Submission received: 23 November 2023 / Revised: 5 February 2024 / Accepted: 19 February 2024 / Published: 1 March 2024
(This article belongs to the Special Issue Contemporary Iterative Methods with Applications in Applied Sciences)

Abstract

:
In this paper, we extend the block hybrid method with equally spaced intra-step points to solve linear and nonlinear third-order initial value problems. The proposed block hybrid method uses a simple iteration scheme to linearize the equations. Numerical experimentation demonstrates that equally spaced grid points for the block hybrid method enhance its speed of convergence and accuracy compared to other conventional block hybrid methods in the literature. This improvement is attributed to the linearization process, which avoids the use of derivatives. Further, the block hybrid method is consistent, stable, and gives rapid convergence to the solutions. We show that the simple iteration method, when combined with the block hybrid method, exhibits impressive convergence characteristics while preserving computational efficiency. In this study, we also implement the proposed method to solve the nonlinear Jerk equation, producing comparable results with other methods used in the literature.

1. Introduction

Differential equations with an appropriate set of conditions are key tools in modeling a plethora of real-world phenomena including engineering problems. It is crucial to have a clear understanding of the behavior and stability of solutions to these differential equations. Many numerical methods have been introduced to find an approximate solution for those cases. In this study, we propose a block hybrid method for the numerical solution of third-order initial value problems (IVPs) of the form
u = f z , u , u , u , z [ a , b ] ,
u ( a ) = u 0 , u ( a ) = u 0 , u ( a ) = u 0 ,
where the prime denotes the derivative with respect to the independent variable z, and f is a continuous linear or nonlinear function in the interval [ a , b ] .
Block hybrid algorithms are a class of numerical methods that blend linear multi-step approaches with interpolation using power series. These techniques were first introduced by Gragg and Stetter [1] and Shampine and Watts [2] and involve the inclusion of an additional point within each step of the formula. This allows for more precise approximations of solutions to differential equations, as well as gives better convergence rates. Since the pioneering work of Gear [3], block methods have been extensively studied and used in the literature to solve initial value problems and boundary value problems. Motsa [4] proposed an overlapping grid hybrid block method with equally spaced and optimal grid points for linear and nonlinear first-order IVPs. He found that the overlapping grid approach gave better performance than the standard non-overlapping grid in terms of reducing the local truncation error. Shateyi [5] applied a block hybrid method with equally spaced grid points to solve linear and nonlinear first-order IVPs. He reported that equally spaced grid points provided high rates of convergence outperforming the fourth-order Runge–Kutta method [6]. Orakwelu [7] developed an implicit block hybrid method for solving first-, second-, and third-order IVPs. He investigated the convergence rates, accuracy, and robustness of implicit block hybrid algorithms. He further investigated the performance of these algorithms when various countable off-points were imposed between grid points in the derivation process. El-Hawary and Mahmoud [8] presented the spline functions with four collocation points to solve second-order IVPs. They showed that the spline collocation scheme is convergent with order seven under certain conditions on the collocation point parameters. They found that the stability properties of the method are analyzed rigorously, and regions of absolute stability are determined based on the parameter values.
Various iterative techniques have been discussed for solving initial value and boundary value problems in ordinary and partial differential equations [9,10]. These iterative methods yield the solution or an approximation of it through a sequence of successive iterations. In the case of initial value problems, these iterative approaches can be formulated either in an integral or differential manner [11]. These approaches linearize the nonlinear equations governing the problem around the previous iteration, leading to linear differential equations at each step [12]. However, the coefficients of these equations may vary with the independent variable, necessitating numerical methods for obtaining the approximation solution.
In this study, we are motivated by the work in [4] and propose a block hybrid method for third-order IVPs with equally spaced grid points. Several strategies have been developed in the literature to solve third-order IVPs. For example, Osa et al. [13] used a multi-step block implicit hybrid method. They reported that the block method performed better than predictor–corrector methods [14] in terms of being more time-efficient, cost-effective, and accurate. Rufai and Ramos [15] modified the block method using variable step sizes instead of equally spaced points. They reported that the technique was efficient in terms of computational time and minimized truncation errors. Orakwelu et al. [16] proposed a single-step hybrid block method for solving third-order ODEs without first converting to an analogous first-order system. They reported that the scheme can be implemented without the use of starting values or predictors, avoiding the necessity for complex subroutines method. Several other authors have utilized the equal step size procedure for solving the third-order IVPs [17,18,19].
Block hybrid methods, used to solve differential equations, are primarily applicable to linear equations. However, when dealing with nonlinear equations, the method necessitates a preliminary step of linearization before its implementation. This prerequisite ensures the effectiveness and accuracy of methods when applied to nonlinear equations. In this study, we developed an intra-step equally spaced grid point block hybrid method that is used in conjunction with the simple iteration method (SIM) to solve third-order IVPs. They demonstrated that the SIM linearization technique generates accurate solutions for nonlinear differential equations [20]. The SIM approach is based on transforming a nonlinear ordinary differential equation into an iterative scheme made up of linear equations, which are then solved using a block hybrid method numerical approach. Linearization methods based on truncated Taylor series approximations are employed to simplify the nonlinear terms of nonlinear differential equations. Relaxation methods based on the assumption that nonlinear terms are known from previous iterations can also be used to convert a nonlinear problem to a linear discretizable problem. Newton-based linearization techniques such as the quasi-linearization method [21], local linearization method [22], Keller-box method [23], and relaxation method [24], have extensively been used to linearize differential equations. All these methods are based on one-term Taylor series expansion and are thereby susceptible to series truncation errors even before errors associated with the numerical method used to solve the linearised problem. In this paper, we propose a method that seeks to circumvent the problem of errors associated with linearization. The proposed method uses ideas akin to those of fixed point iteration to develop iterative schemes, called Simple Iteration Methods (SIMs), for solving nonlinear differential equations. Moreover, SIMs are easy to implement. We show the effectiveness of the proposed method (referred to as the SIM-BHM) through numerical experiments, demonstrating that it gives fast convergence and accurate solutions.

2. Development of the SIM-BHM Iterative Method

This section describes the development of the SIM-BHM with equally spaced points for the numerical solution to third-order IVPs (1) and (2). Before starting the development of the hybrid block method, we need to linearize the function f using the SIM scheme as below:

2.1. The SIM Linearization Scheme

In this section, we present the development of the SIM. If the nonlinear equation f = 0 is rearranged in the form z = g ( z ) , then an iterative method can be written as
z s + 1 = g s , s = 0 , 1 , 2 , 3 ,
where g s = g ( z s ) , and s is the number of iterative steps. This method is called fixed point iteration, the successive substitution method, or the simple iteration method. The simple iteration method is touted as a differential equation equivalent to the fixed point iteration method for solving nonlinear equations. Thus, we begin by reviewing some basic concepts of the fixed point iteration method. The following theorem gives the conditions of the fixed point iteration method.
Theorem 1. (Fixed Point Theorem).
  • If g is continuous on a closed interval [ a , b ] , for all x [ a , b ] , then g has at least one fixed point in [ a , b ] .
  • If, in addition, g ( z ) exists on ( a , b ) and a positive constant K exists with
    | g ( z ) | K < 1 , f o r a l l z ( a , b ) ,
    then there is exactly one fixed point in [ a , b ] .
Assume that c in ( a , b ) is a fixed point for g. Then, if there is an initial guess in ( a , b ) , the sequence
z s + 1 = g s , s = 0 , 1 , 2 ,
converges to the unique fixed point c. For proof of Theorem 1, see [25].

Simple Iteration Method for the ODE

The approach proposed below is a systematic approach based on a rule used to simplify products of unknown functions and their derivatives. We assume that
f = f ( z , u , δ , γ ) ,
is a nonlinear function in u, δ , and γ . To develop the iteration scheme, Equation (4) is expressed as a sum of its linear as nonlinear components as
f = L ( z , u , δ , γ ) + N ( z , u , δ , γ ) ,
The SIM linearization scheme is developed so that the linear component L can be written as the sum of linear terms (unknown functions) at iteration s + 1
L = L 2 ( z s ) γ s + 1 + L 1 ( z s ) δ s + 1 + L 0 ( z s ) u s + 1
The nonlinear terms are a combination of known (at iteration s) and unknown =functions. In the nonlinear term, the function with a higher derivative is taken as unknown. Provided that the nonlinear component N can be written as the sum of nonlinear terms
N = N 2 ( z s , u s , δ s , γ s ) γ s + 1 + N 1 ( z s , u s , δ s ) δ s + 1 + N 0 ( z s , u s ) u s + 1 + N ( z s )
Substituting Equations (6) and (7) into Equation (5), we obtain
f = L 2 ( z s , u s , δ s , γ s ) γ s + 1 + L 1 ( z s , u s , δ s ) δ s + 1 + L 0 ( z s , u s ) u s + 1 + N ( z s ) ,
where L 2 , L 1 , L 0 , and N are known functions from a previous iteration given by
L 2 = L 2 ( z s ) + N 2 ( z s , u s , δ s , γ s ) , L 1 = L 1 ( z s ) + N 1 ( z s , u s , δ s ) , L 0 = L 0 ( z s ) + N 0 ( z s , u s ) .
Equation (8) is the linearized form of f.
The objective of this paper is to develop an alternative and novel block hybrid method approach that is based on an equally spaced grid (see Figure 1).
We define the block by Υ r and the fixed step length or size is denoted by h, where
h = z r + 1 z r , r = 0 , 1 , 0 , , N 1 ,
where N is the number of blocks in [ a , b ] . For each block Υ r = [ z r , z r + 1 ] , we define m + 1 intra-step points p j as
p j = j m , j = 0 , 1 , 2 , 3 , , m ,
z r + p j = z r + p j h .
In using the BHM, we approximate the solution u ( z ) to Equation (1) using a polynomial of the form
u ( z ) U ( z ) = k = 0 m + 3 a r , k ( z z r ) k .
In addition, we approximate the first and second derivatives as
δ ( z ) U ( z ) = k = 1 m + 3 a r , k k ( z z r ) k 1 ,
γ ( z ) U ( z ) = k = 2 m + 3 a r , k k ( k 1 ) ( z z r ) k 2 ,
respectively, where a r , k are m + 4 unknown coefficients. We apply collocation at the unknown nodes z r + p j and in Equations (10) and (11). We assume that the approximated solutions (12)–(14) satisfy Equations (1) and (2). The coefficients a r , k are obtained from a system of m + 4 equations with m + 4 unknowns generated from
U ( z r + p j ) = f r + p j , j = 0 , 1 , 2 , , m .
U ( z r ) = a r , 0 = u r , r = 0 , 1 , , N 1 .
U ( z r ) = a r , 1 = δ r , r = 0 , 1 , , N 1 .
U ( z r ) = 2 a r , 2 = γ r , r = 0 , 1 , , N 1 .
where f r + p j = f ( z r + p j , u r + p j , δ r + p j , γ r + p j ) , f r = f ( z r , u r , δ r , γ r ) , u r = u ( z r ) , δ r = δ ( z r ) , and γ r = γ ( z r ) . The above system of Equations (15)–(18) can be written in matrix notation as
1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 2 0 0 0 0 0 0 0 0 6 0 0 0 0 0 0 0 6 24 h p 1 60 h 2 p 1 2 120 h 3 p 1 3 ( m + 3 ) ( m + 2 ) ( m + 1 ) h m p 1 m 0 0 0 6 24 h p 2 60 h 2 p 2 2 120 h 3 p 2 3 ( m + 3 ) ( m + 2 ) ( m + 1 ) h m p 2 m 0 0 0 6 24 h p 3 60 h 2 p 3 2 120 h 3 p 3 3 ( m + 3 ) ( m + 2 ) ( m + 1 ) h m p 3 m 0 0 0 6 24 h 60 h 2 120 h 3 ( m + 3 ) ( m + 2 ) ( m + 1 ) h m a r , 0 a r , 1 a r , 2 a r , 3 a r , 4 a r , 5 a r , 6 a r , m + 3
= u r δ r γ r f r f r + p 1 f r + p 2 f r + p 3 f r + 1
We solve the above system of Equation (19) to obtain a r , k . We use Mathematica code to obtain a r , k . The code is as follows
p = T a b l e p j = j m , { j , 0 , m } ; m = L e n g t h [ p ] 1 ; T a b l e [ z r + p o i n t s [ [ j ] ] = ( r + p [ [ j ] ] ) h , { i , m + 1 } ] ; U = s u m [ a k ( z z r ) k , { k , 0 , m + 3 } ] ; d U = s u m [ a k k ( z z r ) k 1 , { k , 1 , m + 3 } ] ; d d U = s u m [ a k k ( k 1 ) ( z z r ) k 2 , { k , 2 , m + 3 } ] ; d d d U = s u m [ a k k ( k 1 ) ( k 2 ) ( z z r ) k 3 , { k , 3 , m + 3 } ] ; i n i t i a l 1 = ( U / . z z r ) = = u r ; i n i t i a l 2 = ( d U / . z z r ) = = δ r ; i n i t i a l 3 = ( d d U / . z z r ) = = γ r ; N e q u a t i o n s = T a b l e [ ( d d d U / . z z r + p [ [ j ] ] ) = = f r + p [ [ j ] ] , { j , 1 , m + 1 } ] ; N a l l e q u a t i o n s = J o i n [ { i n i t i a l 1 } , { i n i t i a l 2 } , { i n i t i a l 3 } , N e q u a t i o n s ] / / S i m p l i f y .
Substituting a r , k into Equations (12)–(14) and evaluating the result at collocation points (10) and (11) yield
u r + p i = u r + h p i δ r + 1 2 h 2 p i 2 γ r + h 3 j = 0 m α i , j f r + p j ,
δ r + p i = δ r + h p i γ r + h 2 j = 0 m μ i , j f r + p j ,
γ r + p i = γ r + h j = 0 m κ i , j f r + p j ,
where α i , j , μ i , l , and κ i , l are known constants. By substituting Equation (8) into the BHM scheme system of Equations (20)–(22), we obtain a system of equations in matrix form
D 11 u r + p s + 1 + D 12 δ r + p s + 1 + D 13 γ r + p s + 1 = R 1 ,
D 21 u r + p s + 1 + D 22 δ r + p s + 1 + D 23 γ r + p s + 1 = R 2 ,
D 31 u r + p s + 1 + D 32 δ r + p s + 1 + D 33 γ r + p s + 1 = R 3 ,
where
D 11 = I m h 3 A L 0 , r + p , D 12 = h 3 A L 1 , r + p , D 13 = h 3 A L 2 , r + p , D 21 = h 2 B L 0 , r + p , D 22 = I m h 2 B L 1 , r + p , D 23 = h 2 B L 2 , r + p , D 31 = h C L 0 , r + p , D 32 = h C L 1 , r + p , D 33 = I m h C L 2 , r + p , R 1 = u r + h P δ r + 1 2 h 2 P 2 γ r + h 3 A 0 F r + h 3 A N r , R 2 = δ r + h P γ r + h 2 B 0 F r + h 2 B N r , R 3 = γ r + h C 0 F r + h C N r .
In Equations (23)–(25), I m is an identity matrix of size ( m × m ), and all other matrices of size ( m × m ) are defined as follows
A = α 1 , 1 α 1 , 2 α 1 , m α 2 , 1 α 2 , 2 α 2 , m α m , 1 α m , 2 α m , m , B = μ 1 , 1 μ 1 , 2 μ 1 , m μ 2 , 1 μ 2 , 2 μ 2 , m μ m , 1 μ m , 2 μ m , m
C = κ 1 , 1 κ 1 , 2 κ 1 , m κ 2 , 1 κ 2 , 2 κ 2 , m κ m , 1 κ m , 2 κ m , m , A 0 = α 1 , 0 0 0 0 α 2 , 0 0 0 0 α m , 0 ,
B 0 = μ 1 , 0 0 0 0 μ 2 , 0 0 0 0 μ m , 0 C 0 = κ 1 , 0 0 0 0 κ 2 , 0 0 0 0 κ m , 0 ,
L 2 , r + p = L 2 ( z r + p 1 s , u r + p 1 s , δ r + p 1 s , γ r + p 1 s ) 0 0 L 2 ( z r + p m s , u r + p m s , δ r + p m s , γ r + p m s ) ,
L 1 , r + p = L 1 ( z r + p 1 s , u r + p 1 s , δ r + p 1 s ) 0 0 0 L 1 ( z r + p 2 s , u r + p 2 s , δ r + p 2 s ) 0 0 0 L 1 ( z r + p m s , u r + p m s , δ r + p m s ) ,
L 0 , r + p = L 0 ( z r + p 1 s , u r + p 1 s ) 0 0 0 L 0 ( s z r + p 2 , u r + p 2 s ) 0 0 0 L 0 ( z r + p m s , u r + p m s ) ,
P = p 1 0 0 0 p 2 0 0 0 p m
whereas vectors of size ( m × 1 ) are defined as follows
u r + p s + 1 = u r + p 1 s + 1 u r + p 2 s + 1 u r + p m s + 1 , δ r + p s + 1 = δ r + p 1 s + 1 δ r + p 2 s + 1 δ r + p m s + 1 , γ r + p s + 1 = γ r + p 1 s + 1 γ r + p 2 s + 1 γ r + p m s + 1 , u r = u r s u r s u r s , δ r = δ r s δ r s δ r s ,
γ r = γ r s γ r s γ r s , F r = L 2 ( z r s , u r s , δ r s , γ r s ) γ r s + L 1 ( z r s , u r s , δ r s ) δ r s + L 0 ( z r s , u r s ) u r s + N ( z r s ) L 2 ( z r s , u r s , δ r s , γ r s ) γ r s + L 1 ( z r s , u r s , δ r s ) δ r s + L 0 ( z r s , u r s ) u r s + N ( z r s ) L 2 ( z r s , u r s , δ r s , γ r s ) γ r s + L 1 ( z r s , u r s , δ r s ) δ r s + L 0 ( z r s , u r s ) u r s + N ( z r s ) ,
N r = N ( z r s ) N ( z r s ) N ( z r s ) ,
We may rewrite System (25) in compact form as follows
D 11 D 12 D 13 D 21 D 22 D 23 D 31 D 32 D 33 u r + p s + 1 δ r + p s + 1 γ r + p s + 1 = R 1 R 2 R 3 .
System (26) may now be written in the form
D U = R ,
where the matrix D is of size ( 3 m × 3 m ). The approximate solutions are found by solving Equation (27) as
U = inv ( D ) R .
where inv is the inverse of the matrix.

3. Analysis of the Developed SIM-BHM

This section provides the fundamental properties of the proposed SIM-BHM.

3.1. Order and Error Constant

Minimizing the truncation error is a fundamental goal in numerical analysis utilized as a measure to improve the accuracy of the method.
Proposition 1.
In the integration block [ z r , z r + 1 ] , for the block hybrid method defined as
u r + p i = u r + h p i δ r + 1 2 h 2 p i 2 γ r + h 3 j = 0 m α i , j f r + p j ,
the local truncation error is given by
L i [ u ( z r ) ; h ] = h m + 4 ( m + 4 ) ! p i m + 4 ( m + 4 ) ( m + 3 ) ( m + 2 ) j = 1 m α i j p i m + 1 u ( m + 4 ) ( z r ) + O ( h m + 5 ) .
Proof. 
Given that u is a sufficiently differentiable function, the local truncation error can be analyzed using the linear operator L defined as
L i [ u ( z r ) ; h ] = C i , 0 u ( z r ) + C i , 1 h u ( z r ) + C i , 2 h 2 u ( z r ) + C i , 3 h 3 u ( z r ) + + C i , q h q u ( q ) ( z r ) ,
where C i , 0 , C i , 1 , , C i , q are constant. The method is said to be order q if
C ^ 0 = C ^ 1 = = C ^ q + 2 = 0 , and C ^ q + 3 0 ,
where C ^ q = [ C 1 , q , C 2 , q , , C m , q ] T . The vector C ^ q + 3 = [ C 1 , q + 3 , C 2 , q + 3 , , C m , q + 3 ] T is the error constant of the method. Using Equation (29), the truncation error of the BHM for third-order can be written in terms of the linear operator as
L i [ u ( z r ) ; h ] = u ( z r + h p i ) u ( z r ) h p i u ( z r ) 1 2 h 2 p i 2 u ( z r ) h 3 j = 0 m α i j u ( z r + p j h ) ,
Expanding the terms u ( z r + h p i ) and u ( z r + h p i ) using Taylor series about z r and then substituting them into Equation (33), we obtain
L i [ u ( z r ) ; h ] = k = 1 K h k + 2 p i k + 2 ( k + 2 ) ! u ( k + 2 ) ( z r ) k = 1 K k h k + 2 k ! j = 0 m α i j p j k 1 u ( k + 2 ) ( z r ) + O ( h K + 1 ) ,
Expanding Equation (34), we obtain
L i [ u ( z r ) ; h ] = k = 1 m + 1 h k + 2 p i k + 2 ( k + 2 ) ! u ( k + 2 ) ( z r ) k = 1 m + 1 k h k + 2 k ! j = 0 m α i j p j k 1 u ( k + 2 ) ( z r ) + k = m + 2 K h k + 2 p i k + 2 ( k + 2 ) ! u ( k + 2 ) ( z r ) k = m + 2 K k h k + 2 k ! j = 0 m α i j p j k 1 u ( k + 2 ) ( z r ) + O ( h K + 1 ) ,
where m is a positive integer. This may be written as
L i [ u ( z r ) ; h ] = k = 3 m + 1 h k k ! p i k k ( k 1 ) ( k 2 ) j = 0 m α i j p j k 3 u ( k ) ( z r ) + k = m + 4 K h k k ! p i k k ( k 1 ) ( k 2 ) j = 0 m α i j p j k 3 u ( k ) ( z r ) + O ( h K + 1 ) ,
Following [4,26], we note that
j = 0 m α i j p j k 3 = p i k k ( k 1 ) ( k 2 ) , for k = 4 , , m + 1 .
Substituting Equation (37) into Equation (36), we obtain
L i [ u ( z r ) ; h ] = h m + 4 ( m + 4 ) ! p i m + 4 ( m + 4 ) ( m + 3 ) ( m + 2 ) j = 1 m α i j p i m + 1 u ( m + 4 ) ( z r ) + O ( h m + 5 ) .
Using Equation (38), the error constant vector is
C ^ m + 4 = [ C 1 , m + 4 , C 2 , m + 4 , , C m , m + 4 ] T ,
where
C i , m + 4 = p i m + 4 ( m + 4 ) ! 1 ( m + 1 ) ! j = 1 m α i j p i m + 1 .
Table 1 shows the truncation error and order for m = 2 , 3 , 4 , 5 , 6 . In the table, we note that increasing the number of intra-step points improves the order of the method.

3.2. Zero Stability

The zero stability of the SIM-BHM pertains to the performance of Equation (20) as the step size h approaches zero [27].
Definition 1.
Zero stability is confirmed if the roots of the characteristic polynomial ρ ( λ ) are such that | λ j | 1 and all roots with | λ j | = 1 have a multiplicity that does not exceed 2.
Using Equations (23)–(25), as h 0 , the matrix system reduces to
I m u r + p A 1 u r 1 + p = 1 0 0 0 1 0 0 0 1 u r + p 1 u r + p 2 u r + p m 0 0 1 0 0 1 0 0 1 u r 1 + p 1 u r 1 + p 2 u r 1 + p m ,
The characteristic polynomial ρ ( λ ) is given by
ρ ( λ ) = det ( λ I m A 1 ) = λ m ( λ 1 ) .
In this case, the block hybrid method derived from m equally spaced intra-step points is zero-stable and consistent for any selection of intra-step points.

3.3. Absolute Stability

Definition 2.
A region is said to have absolute stability or be A-stable if it contains the entire left half-plane.
A region of absolute stability for the method can be defined as
( y ) = { y C : H ( y ) < 1 } .
Applying the test equation [28]
u = λ ^ 3 u , λ < 0 ,
to the new method gives
u r + p i = H ( y ) u r , y = λ ^ 3 h 3 ,
where the H ( y ) is given by
H ( y ) = A ˜ 0 + y p 0 + y 2 p 0 2 + y 3 β 0 A ˜ 1 y 3 α 1 .
The stability function ( y ) can be determined by finding the dominant eigenvalues of the matrix H ( y ) . The region of absolute stability for this method is shown in Figure 2. The stability function ( y ) is to be bounded if ( y ) 1 . The bounded region is illustrated in Figure 3. The BHM is absolute-stable as the stability region contains the entire left half complex plane for m = 2 , 3 , 4 , 5 .

4. Implementation and Computational Procedure of the Proposed SIM-BHM

In this section, we present solutions to third-order IVPs (1)–(2) using the novel SIM-BHM iterative method. In order to evaluate the convergence of the method, we have evaluated the absolute error (AE) and absolute error estimate (AEE) between two consecutive iterations. Suppose U ( z ) , δ ( z ) and γ ( z ) are the approximate solutions. The absolute error is defined as
A E u = | u r + p i U r + p i | , A E u = | u r + p i δ r + p i | , A E u = | u r + p i γ r + p i | .
We compute the absolute error estimate between two consecutive iterations using the following terms
A E E u = | U r + p i s + 1 U r + p i s | , A E E δ = | δ r + p i s + 1 δ r + p i s | , A E E γ = | γ r + p i s + 1 γ r + p i s | .
where s represents the solution at the corresponding iteration. We enforce the following conditional stopping procedure in the SIM-BHM iterative method:
  • If A E E < T o l , the current iteration solutions are deemed satisfactory, and the SIM-BHM proceeds to the next block.
  • If A E E > T o l , the iteration count is incremented, and the SIM-BHM continues within the same block.
Here, T o l represents the user-defined tolerance. Within each block, the method iterates until the error A E E falls below T o l . Once the A E E converges to within an acceptable criteria, the procedure advances to the next block or concludes the computation. By applying this conditional structure, we ensure that the accuracy of our solutions is systematically improved within each block, thereby enhancing control and precision in our numerical computations.

Algorithm

To illustrate how to implement this SIM-BHM, we demonstrate an algorithm with the steps provided below:
1.
Define function
  • Input: Initial value problem function f ( z , u , u , u ) , interval [ a , b ] , number of blocks N, tolerance T o l , number of intra-step points m.
  • Output: Approximate solution for u ( z ) .
2.
Linearization scheme
  • Linearize f by using Equation (8).
3.
Discretization
  • Divide interval [ a , b ] into N blocks and determine the step size h.
4.
Collocation points
  • Generate a set of m + 1 fractions in Equation (10).
  • Calculate collocation points within each block using Equation (11).
5.
Approximate solutions
  • Approximate the solutions within each block using a polynomial of degree m + 3 in Equations (12)–(14), which satisfies Equations (1) and (2).
6.
Collocation equations
  • Set up collocation equations and the initial conditions within each block based on Equations (15)–(18).
  • Solve the system obtained from collocation Equations (15)–(18) to obtain coefficients a r , k .
7.
Initialize iteration
  • Initialize iteration s = 1 .
8.
Solve equations
  • Solve the system of equations obtained from collocation Equations (20)–(22) to obtain approximation solutions.
9.
Iterate
  • Repeat Steps 7–10 until convergence or maximum iterations reached:
    (a)
    Use Equation (28) to compute approximation solutions u r + p j s + 1 , δ r + p j s + 1 , and γ r + p j s + 1 at collocation points.
    (b)
    Compute absolute error estimate (AEE) between consecutive iterations using Equation (45).
    (c)
    If A E E u , A E E δ , and A E E γ T o l , accept solutions and move to next block (Step 10).
    (d)
    If A E E u , A E E δ , and A E E γ > T o l , reject solutions and increase iteration s by 1.
10.
Output
  • Output the approximate solutions u r + p j , δ r + p j , and γ r + p j at the block r.

5. Numerical Experimentation

In the next section, we test the method by implementing the SIM-BHM with some specific third-order IVPs from the literature.
Example 1.
Consider the linear third-order IVP [13,17,18,19]
u 3 sin z = 0 , z [ 0 , b ] , u ( 0 ) = 2 , u ( 0 ) = 0 , u ( 0 ) = 1 ,
This IVP has the exact solution: u ( z ) = 3 cos z + 1 2 z 2 2 .
Here, we selected m = 6 , and to implement the SIM-BHM6 scheme, we defined
f = 3 sin z , L 2 = 0 , L 1 = 0 , L 0 = 0 , N = 3 sin z .
In Table 2, we compare the maximum absolute errors of the SIM-BHM6 method with different variants of the hybrid block method [13,17,18,19]. Adesanya et al. [17] used the block method with six collocation points, Areo et al. [18] used the one-twelfth multi-step block method, Skwame et al. [19] used the equally spaced block method with five collocation points, and Osa et al. [13] used the multi-step block method with fifth–fourth collocation points. It is evident that the SIM-BHM6 consistently outperforms the existing block methods in terms of reducing maximum absolute errors. Figure 4 showcases the number of iterations required in each block. The result provides strong evidence of the impressive convergence properties of SIM-BHM6, achieving good accuracy in just two iterations when T o l = 10 18 across all blocks, as shown in Figure 4a,b. It is particularly noteworthy that the method attains a high level of accuracy while maintaining computational efficiency.
Example 2.
Consider the linear third-order IVP [15,29,30]
u 2 u 3 u + 10 u + 34 e 2 z + 16 e 2 z + 10 z 2 6 z 34 = 0 , z [ 0 , b ] , u ( 0 ) = 0 , u ( 0 ) = 0 , u ( 0 ) = 0 ,
The exact solution: u ( z ) = z 2 e 2 z z 2 + 3 .
In this example, we used SIM-BHM6 with
f = 2 u + 3 u 10 u 34 e 2 z 16 exp ( 2 z ) 10 z 2 + 6 z + 34 , L 2 = 2 , L 1 = 3 , L 0 = 10 , N = 34 e 2 z 16 e 2 z 10 z 2 + 6 z + 34 .
Table 3 gives a comparison of the maximum absolute errors using the SIM-BHM6 method with different variants of the block method [15,29,30]. As shown in Table 3, Rufai and Ramos [15] used a variable step-size fourth-derivative block method, Awoyemi1 et al. [29] used a linear multi-step block method with five collocation points, and Allogmany and Ismail [30] used a fourth and fifth derivative block method with three implicit collocation points. We note that SIM-BHM6 outperformed the other block methods in terms of reducing maximum absolute errors. Figure 5 shows that the method achieves good accuracy in two iterations for T o l = 10 10 across all blocks.
Example 3.
Consider the nonlinear third-order IVP [31]
u z u + z 2 u 2 z sin z + cos z z 2 sin 2 z = 0 , z [ 0 , b ] , u ( 0 ) = 0 , u ( 0 ) = 1 , u ( 0 ) = 0 ,
The exact solution of Example 3 is u ( z ) = sin z .
In this example, we selected m = 4 , m = 5 , and m = 6 . We define
f = z u z 2 u 2 + z sin z cos z + z 2 sin 2 z , L 2 = z , L 1 = 0 , L 0 = z u , N = z sin z cos z + z 2 sin 2 z .
In Table 4, we present a comparative analysis of the maximum absolute errors achieved using the SIM-BHM4, SIM-BHM5, and SIM-BHM6 methods with the work in [31]. Adeyeye and Omar [31] employed a block method with equally spaced collocation points. Notably, our findings show enhancement in accuracy when adopting SIM-BHM, particularly as we increase the number of intra-step points and reduce the step sizes. Figure 6a–c provide a visual representation of the number of iterations required across various blocks. In the case of shorter intervals, as in Figure 6a with 10 blocks, we observe that the maximum number of iterations required is five. As shown in Figure 6b, when 100 blocks are utilized, the maximum number of iterations reduces to four. These observations underscore the role of smaller step sizes in reducing the maximum iteration count within the blocks. For larger intervals, as illustrated in Figure 6a with 300 blocks, the maximum number of iterations remains at five. This suggests that, even in scenarios with a larger solution domain, the SIM-BHM maintains its efficient convergence characteristics. Figure 6d provides further evidence of SIM-BHM’s effectiveness, indicating that the method returns a maximum absolute error less than O ( 10 13 ) . This level of precision underscores the robustness and accuracy of the SIM-BHM approach. SIM-BHM6 generally outperforms both SIM-BHM4 and SIM-BHM5.
Example 4.
Consider a nonlinear system of third-order IVPs [15]
u = 1 2 w v e 4 z , u ( 0 ) = 1 , u ( 0 ) = 1 , u ( 0 ) = 1 , z [ 0 , b ] , v = 8 3 u w e 2 z , v ( 0 ) = 1 , v ( 0 ) = 2 , v ( 0 ) = 4 , w = 27 v u , w ( 0 ) = 1 , w ( 0 ) = 3 , w ( 0 ) = 9 ,
with the exact solutions given as u ( z ) = exp ( z ) , v ( z ) = exp ( 2 z ) , and w ( z ) = exp ( 3 z ) .
In the case of a system of equations, the SIM-BHM6 parameters are given by
f 1 = 1 2 w v e 4 z , L 21 = 0 , L 11 = 0 , L 01 = 0 , N 1 = 1 2 w v e 4 z , f 2 = 8 3 u w e 4 z , L 22 = 0 , L 12 = 0 , L 02 = 0 , N 2 = 8 3 u w e 4 z , f 3 = 27 w v , L 23 = 0 , L 13 = 0 , L 03 = 0 , N 3 = 27 w v .
Table 5 shows a comparison of the maximum absolute errors achieved using SIM-BHM6 with the work of Rufai and Ramos [15]. Table 5 illustrates the superior accuracy and computational efficiency of SIM-BHM6 compared to the methodology proposed in [15]. SIM-BHM6 consistently yields more precise approximate solutions across a range of scenarios while maintaining efficient computational performance. Figure 7 provides a representation of the maximum iterations within the computational blocks. Notably, the reduction in user-defined tolerance led to an increase in the number of iterations between the blocks. We observed a maximum of six iterations. Figure 8 gives a comparison of the exact vs. numerical of the displacements in Figure 8a and phase portrait in Figure 8b.
Example 5.
Consider a nonlinear IVP [13] given as
u ( 2 z u + u ) u = 0 , z [ 0 , b ] , u ( 0 ) = 1 , u ( 0 ) = 1 2 , u ( 0 ) = 0 ,
with the exact solution:
u ( z ) = 1 2 + log ( 2 + z 2 z ) .
In this example, we use SIM-BHM4 with
f = ( 2 z u + u ) u , L 2 = 2 z u , L 1 = u , L 0 = 0 , N ( z ) = 0 .
Table 6 provides a comprehensive comparison between the maximum absolute errors achieved using the SIM-BHM4 method and the results obtained by Osa and Olaoluwa [13]. Figure 9 illustrates the number of iterations within the computational blocks for different values of b for a tolerance T o l = 10 10 . The maximum number of iterations required was s = 4 for short and large intervals. The corresponding absolute error profiles for u ( z ) , u ( z ) and u ( z ) are plotted in Figure 10a. Figure 10b shows the exact vs. numerical solution for u ( z ) , u ( z ) , and u ( z ) .

6. On the Application of the SIM-BHM to Solve Nonlinear Jerk Equations

The concept of jerk is particularly relevant in the study of motion and control systems, especially in engineering and physics [32,33]. In physics, the jerk equation models particle motion under varying forces, providing valuable insights into dynamic systems [34]. The equation is used to model systems where the rate of change in acceleration is important, such as in vibration analysis, control systems, and motion planning. Additionally, jerk analysis finds application in signal processing and control theory, facilitating the filtering of noisy signals and the detection of anomalies, contributing to technological advancements. Higher derivatives of motion, such as snap, are also important in motion control and can be experienced in everyday life, such as on trampolines and roller coasters [35].
We consider a class of nonlinear jerk equation IVP containing a third derivative of position with respect to time that describes the rate of change in acceleration, and it is given by
u + β 1 u + β 2 ( u ) 3 + β 3 u 2 u + β 4 u u u + β 5 u ( u ) 2 = 0 , z [ 0 , b ] ,
u ( 0 ) = 0 , u ( 0 ) = β 0 , u ( 0 ) = 0 ,
where the parameters β 0 , β 1 , β 2 , β 3 , β 4 , and β 5 are constant. Since there is no analytic solution to Equation (46), we compared the results of the solution profiles to approximated solutions that have been previously reported in the literature [36,37,38,39].
  • Case 1:  β 1 = β 4 = β 5 = 0 and β 1 = β 2 = 1 . Then, the nonlinear jerk Equation (46) takes the form
    u + u 2 u + ( u ) 3 = 0 , z [ 0 , b ] ,
    u ( 0 ) = 0 , u ( 0 ) = β 0 , u ( 0 ) = 0 ,
    Gottlieb [38] employed the harmonic balance method (HBM) approach approximation solution to solve the nonlinear Equation (48). He found that HBM yields a good approximated solution of the period and displacement amplitude of oscillations for a range of values of initial velocity, given as
    u ( z ) = β 0 Ω sin Ω z , Ω = 1 2 2 3 β 0 2 + 9 β 0 4 + 16 β 0 2
    where Ω is the angular frequency. To solve Equation (48) numerically, we use SIM-BHM7 with
    f = u 2 u ( u ) 3 , L 2 = 0 , L 1 = u 2 ( u ) 2 , L 0 = 0 , N ( z ) = 0 .
    Case 1 is solved using SIM-BHM7 with h = 0.01 and T o l = 10 14 in the domain of integration [ 0 , 20 ] . The displacement and phase trajectory are given in Figure 11a and Figure 11b, respectively. A comparison of the approximated analytical vs. numerical solution is depicted by plotting the two results on the same graph. The corresponding iterations per block and the absolute error profiles for the displacement are plotted in Figure 12a and Figure 12b, respectively, for different values of the initial velocity β 0 .
  • Case 2:  β 2 = β 3 = β 5 = 0 , β 1 = 1 and β 4 = 1 . For the Case 2 values of the parameters, the nonlinear jerk Equation (46) is given as
    u + u u u u = 0 , z [ 0 , b ] ,
    u ( 0 ) = 0 , u ( 0 ) = β 0 , u ( 0 ) = 0 .
    Mirzabeigy and Yildirim [39] employed the modified differential transform method (MDTM) to obtain approximate periodic solutions to Equation (50), given as
    u ( z ) = β 0 Ω sin Ω z + β 0 96 Ω 3 ( 48 Ω 9 β 0 2 48 ) sin Ω z β 0 2 sin 3 Ω z , Ω = 1 2 β 0 2 + 4 .
    In Case 2, we employ SIM-BHM3 with
    f = u + u u u , L 2 = u u , L 1 = 1 , L 0 = 0 , N ( z ) = 0 .
    To determine the accuracy of SIM-BHM3, the numerical results are compared with Kashkari and Alqarni [40]. In their work, they implemented a two-step hybrid block method (TSHBM) with a polynomial of degree 6. Table 7 gives a comparison of the SIM-BHM for β 0 = 0.2 and β 0 = 0.3 . It can be observed that SIM-BHM3 gives similar results to MDTM and TSHBM within a fast CPU time.
Figure 13a depicts the number of iterations required per block for different values of the initial velocity β 0 = 0.2 , 0.4 , 0.6 . It can be observed that SIM-BHM3 converges within five iterations. Figure 13b illustrates the phase portrait for different values of the initial velocity.
  • Case 3: β 2 = β 3 = β 4 = 0 , β 1 = 1 , and β 5 = 1 .
    u + u + u ( u ) 2 = 0 , z [ 0 , b ] ,
    u ( 0 ) = 0 , u ( 0 ) = β 0 , u ( 0 ) = 0 .
    The periodic solutions to Equation (52) are given as [39]
    u ( z ) = β 0 Ω sin Ω z + β 0 96 Ω 3 [ ( 48 Ω 9 β 0 2 48 ) sin Ω z β 0 2 sin 3 Ω z + ( 12 Ω 3 β 0 2 + 48 Ω 48 Ω 3 ) z cos Ω z ] , Ω = 2 1 4 β 0 2 .
    Here, we set
    f = u u ( u ) 2 , L 2 = u u , L 1 = 1 , L 0 = 0 , N ( z ) = 0 .
    The displacement, velocity, acceleration, and phase portrait trajectories are given in Figure 14. A comparison of the exact (MDTM) vs. numerical (SIM-BHM3) solutions is depicted by plotting the two results on the same graph. The corresponding phase portrait and iterations per block are plotted in Figure 15a and Figure 15b, respectively, for different values of β 0 .
  • Case 4: β 1 = β 2 = β 3 = β 4 = β 5 = 1 . For the Case 4 values of the parameters, the nonlinear jerk Equation (46) is given as
    u + u + ( u ) 3 + u 2 u + u u u + u ( u ) 2 = 0 , z [ 0 , b ] ,
    u ( 0 ) = 0 , u ( 0 ) = β 0 , u ( 0 ) = 0 .
    With linearization, we obtain
    f = u ( u ) 3 u 2 u u u u u ( u ) 2 , L 2 = u u u u , L 1 = 1 ( u ) 2 u 2 , L 0 = 0 , N ( z ) = 0 .
    Equation (54) is solved using SIM-BHM7 with h = 0.01 in the domain of integration [ 0 , 40 ] for different values of the initial velocity β 0 . The displacement and iterations per block utilized to solve Equation (54) are plotted in Figure 16 for different values of the initial velocity.
    Figure 17 depicts phase portraits for Case 4 when varying β 0 .
  • Case 5: β 1 = β 2 = β 3 = β 4 = β 5 = 4 . For the Case 5 values of the parameters, the nonlinear jerk Equation (46) is given as
    u + 4 u + ( u ) 3 + u 2 u + u u u + u ( u ) 2 = 0 , z [ 0 , b ] ,
    u ( 0 ) = 0 , u ( 0 ) = β 0 , u ( 0 ) = 0 .
    With linearization, we obtain
    f = 4 u 4 ( u ) 3 4 u 2 u 4 u u u 4 u ( u ) 2 , L 2 = 4 u u 4 u u , L 1 = 4 4 ( u ) 2 4 u 2 , L 0 = 0 , N ( z ) = 0 .
    Equation (56) is solved using SIM-BHM7 with h = 0.01 and b = 40 for different values of the initial velocity. Figure 18 illustrates the solution u ( z ) and the iteration required. Figure 19 shows the phase portrait.
We considered five cases of nonlinear jerk Equation (46) with different parameter combinations for Equations (48)–(56). This allowed us to test the SIM-BHM approach on various forms of the jerk equation. For Cases 1–3, known analytical solutions from previous works were available to validate our numerical results. The SIM-BHM-produced displacement profiles, velocities, accelerations, and phase portraits are in agreement with these established solutions. In Cases 4 and 5, there were no analytical solutions. The SIM-BHM generated physically realistic oscillatory behaviors as the initial velocity parameter β 0 was varied. This indicates that the method can reliably solve these types of jerk equations numerically. Across all test cases, the SIM-BHM converged rapidly, typically requiring only 2–8 iterations per block. Even for large solution domains and higher β 0 values producing larger oscillations, the method maintained its fast convergence. The displacement profiles, phase portraits, and plots of the iterations clearly illustrate the accuracy and consistency of the SIM-BHM approach. Comparisons with previous works validated the high precision of the numerical solutions obtained. In conclusion, the SIM-BHM proves to be an effective technique for solving an important class of nonlinear jerk equations. It generates solutions in close agreement with the known results. The method exhibits robust computational performance and accuracy across a wide range of problem scenarios.
The numerical experiments consistently showed the SIM-BHM achieving higher accuracy when a smaller step size h was used. For example, reducing h from 10 1 to 10 2 led to errors 1–2 orders of magnitude smaller across test problems. This implies that the method exhibits the expected property of increased accuracy for smaller discretization. Larger solution domains were also handled effectively. Cases 1–3 of the jerk equations covered intervals up to b = 20 with good precision. Cases 4 and 5 were solved over the wider range of b = 40 while maintaining accuracy. The SIM-BHM modeled problems requiring integration over large physical spaces. Additional support comes from Example 5, where accuracy improved for the longer interval of b = 1.5 compared to b = 1 . These examples provide strong evidence that SIM-BHM’s accuracy increases predictably with smaller h and can be systematically extended to large b values through domain discretization.
Figure 11. Displacement and phase trajectory for using SIM-BHM7 when T o l = 10 14 , b = 20 , h = 0.01 , and N = 2000 .
Figure 11. Displacement and phase trajectory for using SIM-BHM7 when T o l = 10 14 , b = 20 , h = 0.01 , and N = 2000 .
Appliedmath 04 00017 g011
Figure 12. Iterations per block and absolute error for using the SIM-BHM7 when β 0 = 5 , T o l = 10 14 , b = 20 , h = 0.01 , and N = 2000 .
Figure 12. Iterations per block and absolute error for using the SIM-BHM7 when β 0 = 5 , T o l = 10 14 , b = 20 , h = 0.01 , and N = 2000 .
Appliedmath 04 00017 g012
Figure 13. Displacement and phase trajectory for using the SIM-BHM3 when T o l = 10 18 , b = 20 , h = 0.01 , and N = 2000 .
Figure 13. Displacement and phase trajectory for using the SIM-BHM3 when T o l = 10 18 , b = 20 , h = 0.01 , and N = 2000 .
Appliedmath 04 00017 g013
Figure 14. Displacement, velocity, acceleration, and phase portrait for using the SIM-BHM3 when β 0 = 0.2 , T o l = 10 18 , b = 20 , h = 0.01 , and N = 200 .
Figure 14. Displacement, velocity, acceleration, and phase portrait for using the SIM-BHM3 when β 0 = 0.2 , T o l = 10 18 , b = 20 , h = 0.01 , and N = 200 .
Appliedmath 04 00017 g014
Figure 15. Displacement and phase trajectory for using SIM-BHM3 when T o l = 10 18 .
Figure 15. Displacement and phase trajectory for using SIM-BHM3 when T o l = 10 18 .
Appliedmath 04 00017 g015
Figure 16. Displacement and iterations per block for using SIM-BHM7 when T o l = 10 18 , b = 40 , and N = 4000 .
Figure 16. Displacement and iterations per block for using SIM-BHM7 when T o l = 10 18 , b = 40 , and N = 4000 .
Appliedmath 04 00017 g016
Figure 17. Phase trajectory for using the SIM-BHM7 when T o l = 10 18 , b = 40 , and N = 4000 .
Figure 17. Phase trajectory for using the SIM-BHM7 when T o l = 10 18 , b = 40 , and N = 4000 .
Appliedmath 04 00017 g017
Figure 18. Displacement and iterations per block for using SIM-BHM7 when T o l = 10 18 , b = 40 , and N = 4000 .
Figure 18. Displacement and iterations per block for using SIM-BHM7 when T o l = 10 18 , b = 40 , and N = 4000 .
Appliedmath 04 00017 g018
Figure 19. Phase portraits for using SIM-BHM7 when T o l = 10 18 , b = 40 , and N = 4000 .
Figure 19. Phase portraits for using SIM-BHM7 when T o l = 10 18 , b = 40 , and N = 4000 .
Appliedmath 04 00017 g019

7. Conclusions

We introduced a novel SIM-BHM iterative method for solving third-order IVPs with a fixed step size. The primary objective of the SIM-BHM is to obtain accurate solutions to the IVPs. To evaluate the convergence of the method, we evaluated the absolute error and the absolute error estimate between consecutive iterations. We applied the SIM-BHM to a variety of IVPs, both linear and nonlinear. We compared its performance with existing block methods. The results consistently demonstrated the superior performance of the SIM-BHM in terms of returning the lowest maximum absolute errors, highlighting its accuracy and efficiency. The SIM-BHM offers a versatile approach by employing a selective intra-step method and ensuring convergence. The method gives highly accurate solutions with good computational performance. The main findings with regard to this numerical method are as follows:
  • Smaller intra-step sizes generally lead to smaller truncation errors.
  • Increasing the number of collection points improves the accuracy of SIM-BHM.
  • For large intervals, the SIM-BHM gave more accurate approximations than some existing BHMs, for example, those in [13,15,17,18,19,29,30,31].
  • The SIM-BHM gives robust computational performance using only a few iterations.
The SIM-BHM method was applied to solve nonlinear jerk equations. Jerk equations are an important class of third-order problems with applications in fields like engineering and physics. Five cases of jerk equations with varying parameter combinations were tested. Comparisons to known approximate solutions from previous works showed excellent agreement for Cases 1–3. The method generated physically realistic oscillatory behaviors for Cases 4 and 5. Across all jerk equation problems, the SIM-BHM maintained its fast convergence, typically requiring only a few iterations per block even for large solution domains and higher oscillation amplitudes. The accuracy of the numerical solutions was validated through comparisons. This demonstrates the SIM-BHM provides an effective technique for reliably solving nonlinear jerk equations.
In future work, the following would allow for improving the performance of the SIM-BHM:
  • Implementing adaptivity with variable step size would allow the method to better capture solutions with varying behaviors over the domain. This could enhance accuracy.
  • Investigating the use of Legendre–Gauss–Lobatto and Gauss–Lobatto collocation points, which have been shown to improve accuracy.
  • Applying optimally spaced intra-step points within each block, as studies have observed non-uniform point distributions lead to higher precision than uniform spacing [7].
  • Applying the overlapping domain approach, where blocks intersect, as this has been demonstrated to reduce errors versus standard non-overlapping grids [4,41].
  • Extending the SIM-BHM framework to other problem types, such as higher-order initial value problems, boundary value problems, singular differential equations, and fractional differential equations, to evaluate performance on broader classes of equations.

Author Contributions

Conceptualization, S.A.A.A.A.A.; Investigation, S.A.A.A.A.A. and U.O.R.; Writing—original draft, S.A.A.A.A.A.; Writing—review & editing, P.S., S.P.G., H.S.M. and O.A.I.N.; Supervision, P.S., S.P.G. and O.A.I.N. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the University of KwaZulu-Natal.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

Abbreviations

The following abbreviations are used in this manuscript:
AEAbsolute error
AEEAbsolute error estimate
BHMBlock hybrid method
HBMHarmonic balance method
MDTMModified differential transform method
ODEOrdinary differential equation method
TSHBMTwo-step hybrid block method
Max AELargest absolute error
SIMSimple iteration method

References

  1. Gragg, W.B.; Stetter, H.J. Generalized multistep predictor-corrector methods. J. Assoc. Comput. Mach. 1964, 11, 188–209. [Google Scholar] [CrossRef]
  2. Shampine, L.F.; Watts, H. Block implicit one-step methods. Math. Comput. 1969, 23, 731–740. [Google Scholar] [CrossRef]
  3. Gear, C.W. Hybrid methods for initial value problems in ordinary differential equations. J. Soc. Ind. Appl. Math. Ser. B Numer. Anal. 1965, 2, 69–86. [Google Scholar] [CrossRef]
  4. 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]
  5. Shateyi, S. On the Application of the Block Hybrid Methods to Solve Linear and Non-Linear First Order Differential Equations. Axioms 2023, 12, 189. [Google Scholar] [CrossRef]
  6. Burden, R.; Faires, J. Numerical Analysis, 9th ed.; Cengage Learning: Boston, MA, USA, 2011. [Google Scholar]
  7. Orakwelu, M.G. Generalised Implicit Block Hybrid Algorithms for Initial Value Problems. Ph.D. Thesis, University of KwaZulu-Natal, Peitermaritzburg, South Africa, 2019. [Google Scholar]
  8. El-Hawary, H.; Mahmoud, S. On some 4-point spline collocation methods for solving second-order initial value problems. Appl. Numer. Math. 2001, 38, 223–236. [Google Scholar] [CrossRef]
  9. Bellman, R.E.; Kalaba, R.E. Quasilinearization and Nonlinear Boundary-Value Problems; Elsevier: New York, NY, USA, 1965. [Google Scholar]
  10. Mandelzweig, V.; Tabakin, F. Quasilinearization approach to nonlinear problems in physics with application to nonlinear ODEs. Comput. Phys. Commun. 2001, 141, 268–281. [Google Scholar] [CrossRef]
  11. Kelley, C.T. Iterative Methods for Linear and Nonlinear Equations; SIAM: Philadelphia, PA, USA, 1995. [Google Scholar]
  12. Ramos, J.I. On the variational iteration method and other iterative techniques for nonlinear differential equations. Appl. Math. Comput. 2008, 199, 39–69. [Google Scholar] [CrossRef]
  13. Osa, A.L.; Olaoluwa, O.E. A fifth-fourth continuous block implicit hybrid method for the solution of third order initial value problems in ordinary differential equations. Appl. Comput. Math 2019, 8, 50–57. [Google Scholar]
  14. Butcher, J.C. Numerical Methods for Ordinary Differential Equations; John Wiley & Sons: Hoboken, NJ, USA, 2016. [Google Scholar]
  15. Rufai, M.A.; Ramos, H. A variable step-size fourth-derivative hybrid block strategy for integrating third-order IVPs, with applications. Int. J. Comput. Math. 2022, 99, 292–308. [Google Scholar] [CrossRef]
  16. Orakwelu, M.G.; Otegbeye, O.; Mambili-Mamboundou, H. A class of single-step hybrid block methods with equally spaced points for general third-order ordinary differential equations. J. Niger. Soc. Phys. Sci. 2023, 5, 1484. [Google Scholar] [CrossRef]
  17. Adesanya, A.O.; Udoh, D.M.; Ajileye, A. A new hybrid block method for the solution of general third order initial value problems of ordinary differential equations. Int. J. Pure Appl. Math. 2013, 86, 365–375. [Google Scholar] [CrossRef]
  18. Areo, E.; Omojola, M. One-Twelveth Step Continuous Block Method. Int. J. Pure Appl. Math. 2017, 114, 165–178. [Google Scholar] [CrossRef]
  19. Skwame, J.; Dalatu, P.; Sabo, J.; Mathew, M. Numerical application of third derivative hybrid block methods on third Order Initial Value Problem of ordinary differential equations. Int. J. Stat. Appl. Math. 2019, 4, 90–100. [Google Scholar]
  20. Otegbeye, O.; Ansari, M.S. A finite difference based simple iteration method for solving boundary layer flow problems. AIP Conf. Proc. 2022, 2435, 020055. [Google Scholar]
  21. Parsaeitabar, Z.; Nazemi, A. A Third-degree B-spline Collocation Scheme for Solving a Class of the Nonlinear Lane—Emden Type Equations. Iran. J. Math. Sci. Inform. 2017, 12, 15–34. [Google Scholar]
  22. Motsa, S.S. A new spectral local linearization method for nonlinear boundary layer flow problems. J. Appl. Math. 2013, 2013, 423628. [Google Scholar] [CrossRef]
  23. Keller, H.B.; Cebeci, T. Accurate numerical methods for boundary-layer flows. II: Two dimensional turbulent flows. AIAA J. 1972, 10, 1193–1199. [Google Scholar] [CrossRef]
  24. Zubik-Kowal, B.; Vandewalle, S. Waveform relaxation for functional-differential equations. SIAM J. Sci. Comput. 1999, 21, 207–226. [Google Scholar] [CrossRef]
  25. Burden, A.M.; Burden, R.L.; Faires, J.D. Numerical Analysis, 10th ed.; Cengage Learning: Boston, MA, USA, 2016. [Google Scholar] [CrossRef]
  26. Butcher, J.C. Implicit Runge–Kutta processes. Math. Comput. 1964, 18, 50–64. [Google Scholar] [CrossRef]
  27. Lambert, J.D. Computational Methods in Ordinary Differential Equations; Wiley: London, UK, 1973. [Google Scholar]
  28. Dahlquist, G. Convergence and stability in the numerical integration of ordinary differential equations. Math. Scand. 1956, 4, 33–53. [Google Scholar] [CrossRef]
  29. Awoyemi, D.; Kayode, S.; Adoghe, L. A five-step P-stable method for the numerical integration of third order ordinary differential equations. Am. J. Comput. Math. 2014, 4, 119–126. [Google Scholar] [CrossRef]
  30. Allogmany, R.; Ismail, F. Implicit three-point block numerical algorithm for solving third order initial value problem directly with applications. Mathematics 2020, 8, 1771. [Google Scholar] [CrossRef]
  31. Adeyeye, O.; Omar, Z. Solving third order ordinary differential equations using one-step block method with four equidistant generalized hybrid points. Int. J. Appl. Math. 2019, 49, 1–9. [Google Scholar]
  32. Kyriakopoulos, K.J.; Saridis, G.N. Minimum jerk path generation. In Proceedings of the 1988 IEEE International Conference on Robotics and Automation, Philadelphia, PA, USA, 24–29 April 1988; pp. 364–369. [Google Scholar]
  33. Hayati, H.; Eager, D.; Pendrill, A.M.; Alberg, H. Jerk within the context of science and engineering—A systematic review. Vibration 2020, 3, 371–409. [Google Scholar] [CrossRef]
  34. Eichhorn, R.; Linz, S.J.; Hänggi, P. Transformations of nonlinear dynamical systems to jerky motion and its application to minimal chaotic flows. Phys. Rev. E 1998, 58, 7151. [Google Scholar] [CrossRef]
  35. Eager, D.; Pendrill, A.M.; Reistad, N. Beyond velocity and acceleration: Jerk, snap and higher derivatives. Eur. J. Phys. 2016, 37, 065008. [Google Scholar] [CrossRef]
  36. Gottlieb, H. Simple nonlinear jerk functions with periodic solutions. Am. J. Phys. 1998, 66, 903–906. [Google Scholar] [CrossRef]
  37. Liu, C.S.; Chang, J.R. The periods and periodic solutions of nonlinear jerk equations solved by an iterative algorithm based on a shape function method. Appl. Math. Lett. 2020, 102, 106151. [Google Scholar] [CrossRef]
  38. Gottlieb, H. Harmonic balance approach to periodic solutions of non-linear jerk equations. J. Sound Vib. 2004, 271, 671–683. [Google Scholar] [CrossRef]
  39. Mirzabeigy, A.; Yildirim, A. Approximate periodic solution for nonlinear jerk equation as a third-order nonlinear equation via modified differential transform method. Eng. Comput. 2014, 31, 622–633. [Google Scholar] [CrossRef]
  40. Kashkari, B.S.; Alqarni, S. Two-step hybrid block method for solving nonlinear jerk equations. J. Appl. Math. Phys. 2019, 7, 1893. [Google Scholar] [CrossRef]
  41. Mkhatshwa, M.P. Overlapping Grid Spectral Collocation Methods for Nonlinear Differential Equations Modelling Fluid Flow Problems. Ph.D. Thesis, University of KwaZulu-Natal, Peitermaritzburg, South Africa, 2020. [Google Scholar]
Figure 1. Illustration of the equally spaced grid points.
Figure 1. Illustration of the equally spaced grid points.
Appliedmath 04 00017 g001
Figure 2. Region of absolute stability for m = 2 , 3 , 4 , 5 .
Figure 2. Region of absolute stability for m = 2 , 3 , 4 , 5 .
Appliedmath 04 00017 g002
Figure 3. Region of absolute stability for m = 2 , 3 , 4 , 5 .
Figure 3. Region of absolute stability for m = 2 , 3 , 4 , 5 .
Appliedmath 04 00017 g003aAppliedmath 04 00017 g003b
Figure 4. Iterations per block using the SIM-BHM6.
Figure 4. Iterations per block using the SIM-BHM6.
Appliedmath 04 00017 g004
Figure 5. Iterations per block using SIM-BHM6.
Figure 5. Iterations per block using SIM-BHM6.
Appliedmath 04 00017 g005
Figure 6. Iterations per block and absolute error for SIM-BHM when m = 4 , m = 5 , and m = 6 .
Figure 6. Iterations per block and absolute error for SIM-BHM when m = 4 , m = 5 , and m = 6 .
Appliedmath 04 00017 g006
Figure 7. Iterations per block using SIM-BHM6 when T o l = 10 8 and T o l = 10 12 .
Figure 7. Iterations per block using SIM-BHM6 when T o l = 10 8 and T o l = 10 12 .
Appliedmath 04 00017 g007
Figure 8. Numerical solution vs. exact solution using SIM-BHM6 when T o l = 10 12 , h = 10 1 , and b = 6 .
Figure 8. Numerical solution vs. exact solution using SIM-BHM6 when T o l = 10 12 , h = 10 1 , and b = 6 .
Appliedmath 04 00017 g008
Figure 9. Iterations per block using the SIM-BHM4 when T o l = 10 10 .
Figure 9. Iterations per block using the SIM-BHM4 when T o l = 10 10 .
Appliedmath 04 00017 g009
Figure 10. Absolute error and numerical solution vs. exact solution using the SIM-BHM4 when T o l = 10 10 , h = 10 1 , and b = 1.5 .
Figure 10. Absolute error and numerical solution vs. exact solution using the SIM-BHM4 when T o l = 10 10 , h = 10 1 , and b = 1.5 .
Appliedmath 04 00017 g010
Table 1. Truncation errors for the SIM-BHM.
Table 1. Truncation errors for the SIM-BHM.
m p i Truncation Errorq
2 1 2 , 1 u ( 6 ) z r h 6 9216 , u ( 6 ) z r h 6 1440 3
3 1 3 , 2 3 , 1 47 u ( 7 ) z r h 7 22044960 , 19 u ( 7 ) z r h 7 1377810 , u ( 7 ) z r h 7 30240 4
4 1 4 , 1 2 , 3 4 , 1 139 u ( 8 ) z r h 8 2642411520 , u ( 8 ) z r h 8 2949120 , 243 u ( 8 ) z r h 8 293601280 , u ( 8 ) z r h 8 645120 5
5 1 5 , 2 5 , 3 5 , 4 5 , 1 9809 u ( 9 ) z r h 9 7087500000000 , 491 u ( 9 ) z r h 9 55371093750 , 1917 u ( 9 ) z r h 9 87500000000 , 1136 u ( 9 ) z r h 9 27685546875 , 149 u ( 9 ) z r h 9 2268000000 6
6 1 6 , 1 3 , 1 2 , 2 3 , 5 6 , 1 4001 u ( 10 ) z r h 10 109709829734400 , 199 u ( 10 ) z r h 10 857108044800 , 29 u ( 10 ) z r h 10 267846264400 , 29 u ( 10 ) z r h 10 26784626400 , 7625 u ( 10 ) z r h 10 4388393189376 , u ( 10 ) z r h 10 391910400 7
Table 2. Maximum absolute error comparison using SIM-BHM6.
Table 2. Maximum absolute error comparison using SIM-BHM6.
bhAdesanya et al. [17]Areo and Omojola [18]Skwame et al. [19]Osa and Olaoluwa [13]SIM-BHM6
Max AE1 10 1 2.2736   ×   10 10 5.8703   ×   10 15 3.5652   ×   10 14 1.3678   ×   10 13 4.4409   ×   10 16
CPU ----0.024886
Max AE5-----1.0658   ×   10 14
CPU ----0.029308
Max AE1 10 2 ----7.7716   ×   10 16
CPU ----0.035961
Max AE5 10 2 ----1.2434   ×   10 14
CPU ----0.096864
Table 3. Maximum absolute error comparison for SIM-BHM6 when T o l = 10 10 .
Table 3. Maximum absolute error comparison for SIM-BHM6 when T o l = 10 10 .
bhRufai and Ramos [15]Awoyemi1 et al. [29]Allogmany and Ismail [30]SIM-BHM6
Max AE1 10 1 2.08069 × 10 11 3.73600 × 10 5 2.84944 × 10 9 9.78773 × 10 13
CPU ---0.025860
Max AE4 10 1 4.50239 × 10 9 -- 3.5040 × 10 10
CPU ---0.031220
Max AE4 10 2 --- 5.2776 × 10 12
CPU ---0.092856
Table 4. Maximum absolute error comparison for SIM-BHM4, SIM-BHM5, and SIM-BHM6.
Table 4. Maximum absolute error comparison for SIM-BHM4, SIM-BHM5, and SIM-BHM6.
bh Tol Adeyeye and Omar [31]SIM-BHM4SIM-BHM5SIM-BHM6
Max AE1 10 1 10 18 4.644063 × 10 13 4.640732 × 10 13 5.173639 × 10 14 5.551115 × 10 16
CPU -0.0131250.0121830.013210
Max AE1 10 2 10 18 -1.4433 × 10 15 1.7764 × 10 15 9.9920 × 10 16
CPU -0.0289990.0317100.031207
Max AE3 10 2 10 18 -4.1911× 10 14 2.1955 × 10 14 1.8152 × 10 14
CPU -0.0834010.0910060.094574
Table 5. Maximum absolute error comparison using the SIM-BHM6.
Table 5. Maximum absolute error comparison using the SIM-BHM6.
Methodbh Tol Max AE u Max AE v Max AE w CPU
Rufai and Ramos [15]3 10 1 10 8 4.36535 × 10 8 3.08415 × 10 8 2.43322 × 10 8 0.172608
SIM-BHM6 4.4973 × 10 11 4.6361 × 10 11 1.6977 × 10 11 0.094197
Rufai and Ramos [15]6 10 1 10 8 8.81911 × 10 4 2.95244 × 10 5 2.71811 × 10 9 0.172608
SIM-BHM6 2.5891 × 10 6 1.1437 × 10 7 7.1468 × 10 10 0.122777
Rufai and Ramos [15]3 10 1 10 12 2.38150 × 10 11 1.50130 × 10 11 2.66831 × 10 12 0.395336
SIM-BHM6 2.1786 × 10 11 3.0987 × 10 11 1.2338 × 10 11 0.090950
Rufai and Ramos [15]6 10 1 10 12 4.53908 × 10 7 2.87265 × 10 8 4.84118 × 10 10 0.56434
SIM-BHM6 1.8995 × 10 6 8.5346 × 10 8 5.7199 × 10 10 0.126770
Table 6. Maximum absolute errors comparison for SIM-BHM4 when T o l = 10 10 .
Table 6. Maximum absolute errors comparison for SIM-BHM4 when T o l = 10 10 .
bhOsa and Olaoluwa [13]SIM-BHM4CPU
Max | u ( z ) U ( z ) | Max | u ( z ) U ( z ) | Max | u ( z ) δ ( z ) | Max | u ( z ) γ ( z ) |
1 10 2 2.6241 × 10 12 1.7764 × 10 15 3.7748 × 10 15 8.5487 × 10 15 0.042517
1.5 10 2 -1.5765 × 10 14 1.7342 × 10 13 3.3595 × 10 13 0.051358
Table 7. Comparison of implementation the SIM-BHM3 when h = 0.125 , b = 2 , and T o l = 10 18 .
Table 7. Comparison of implementation the SIM-BHM3 when h = 0.125 , b = 2 , and T o l = 10 18 .
z β 0 = 0.2 β 0 = 0.3
MDTM [39]TSHBM [40]SIM-BHM3MDTM [39]TSHBM [40]SIM-BHM3
0.000000000
0.1250.0249340.0249340.0249349507235620.0374020.0374020.037402433678462
0.2500.0494800.0494800.0494808958400670.0742200.0742200.074221538790320
0.3750.0732530.0732530.0732552602789680.1098780.1098790.109884305379430
0.5000.0958810.0958810.0958881619185530.1438140.1438140.143837971511998
0.6250.1170080.1170080.1170283721650120.175490.175490.175559287806927
0.7500.13630.1363000.1363488523828760.2043970.2043990.204562874433873
0.8750.1534520.1534530.1535517672555720.230070.2300730.230408498532534
1.0000.1681910.1681910.1683729048351160.2520880.2520930.252707185491022
1.1500.180280.1802810.1805854616806170.2700880.2700950.271126160102852
1.2500.1895250.1895260.1900031769620890.2837690.2837780.285392683210093
1.3750.1957770.1957780.1964828192080380.2929020.2929130.295296898220800
1.5000.1989350.1989370.1999260419197730.2973320.2973440.300693825250087
1.6250.1989480.1989490.2002806290841050.2969840.2969960.301504637583959
1.7500.1958150.1958160.1975411493383350.2918640.2918760.297717328253788
1.8750.1895870.1895890.1917490298290540.2820590.2820070.289386829338821
2.0000.1803650.1803650.1829920500661850.2677340.2677450.276634591040122
CPU (s)--0.020688--0.019756
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

Ahmedai Abd Allah, S.A.A.; Sibanda, P.; Goqo, S.P.; Rufai, U.O.; Sithole Mthethwa, H.; Noreldin, O.A.I. A Block Hybrid Method with Equally Spaced Grid Points for Third-Order Initial Value Problems. AppliedMath 2024, 4, 320-347. https://doi.org/10.3390/appliedmath4010017

AMA Style

Ahmedai Abd Allah SAA, Sibanda P, Goqo SP, Rufai UO, Sithole Mthethwa H, Noreldin OAI. A Block Hybrid Method with Equally Spaced Grid Points for Third-Order Initial Value Problems. AppliedMath. 2024; 4(1):320-347. https://doi.org/10.3390/appliedmath4010017

Chicago/Turabian Style

Ahmedai Abd Allah, Salma A. A., Precious Sibanda, Sicelo P. Goqo, Uthman O. Rufai, Hloniphile Sithole Mthethwa, and Osman A. I. Noreldin. 2024. "A Block Hybrid Method with Equally Spaced Grid Points for Third-Order Initial Value Problems" AppliedMath 4, no. 1: 320-347. https://doi.org/10.3390/appliedmath4010017

APA Style

Ahmedai Abd Allah, S. A. A., Sibanda, P., Goqo, S. P., Rufai, U. O., Sithole Mthethwa, H., & Noreldin, O. A. I. (2024). A Block Hybrid Method with Equally Spaced Grid Points for Third-Order Initial Value Problems. AppliedMath, 4(1), 320-347. https://doi.org/10.3390/appliedmath4010017

Article Metrics

Back to TopTop