Next Article in Journal
Multifractal Characteristics of Uranium Grade Distribution and Spatial Regularities in a Sandstone-Type Uranium Deposit in Xinjiang, China
Previous Article in Journal
A Fractional (q,q) Non-Extensive Information Dimension for Complex Networks
Previous Article in Special Issue
Solving General Fractional Lane-Emden-Fowler Differential Equations Using Haar Wavelet Collocation Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New Hybrid Block Method for Solving First-Order Differential System Models in Applied Sciences and Engineering

by
Mufutau Ajani Rufai
1,*,
Bruno Carpentieri
1 and
Higinio Ramos
2,3
1
Faculty of Engineering, Free University of Bozen-Bolzano, 39100 Bolzano, Italy
2
Scientific Computing Group, Universidad de Salamanca, Plaza de la Merced, 37008 Salamanca, Spain
3
Department of Mathematics, Escuela Politecnica Superior de Zamora, Campus Viriato, 49022 Zamora, Spain
*
Author to whom correspondence should be addressed.
Fractal Fract. 2023, 7(10), 703; https://doi.org/10.3390/fractalfract7100703
Submission received: 30 August 2023 / Revised: 22 September 2023 / Accepted: 23 September 2023 / Published: 24 September 2023

Abstract

:
This paper presents a new hybrid block method formulated in variable stepsize mode to solve some first-order initial value problems of ODEs and time-dependent partial differential equations in applied sciences and engineering. The proposed method is implemented considering an adaptive stepsize strategy to maintain the estimated error in each step within a specified tolerance. In order to evaluate the performance and usefulness of the proposed technique in real-world applications, several differential problems from applied sciences and engineering, such as the SIR model, Jacobi elliptic function problem, and chemical reactions problems, are solved numerically. The results of numerical simulations in this work demonstrate that the proposed method is more efficient than other existing numerical methods used for comparisons.

1. Introduction

Ordinary Differential Equations (ODEs) are essential for modeling real-life model problems in applied sciences and engineering. For instance, ODEs have found extensive application in modeling real-world problems in chemical, electrical, and mechanical engineering, physics, chemistry, and mathematical biology. Motivated by the diverse applications of ODEs in modeling physical phenomena, this paper introduces a new numerical method for integrating the following initial value problem (IVP):
v ( t ) = f ( t , v ) , v ( t 0 ) = v 0 , v R d , f [ t 0 , t N ] × R d .
The above equation is defined over an integration interval [ t 0 , t N ] and has initial conditions specified at t = t 0 . It is worth noting that the function f ( t , v ) in Equation (1) satisfies Lipchitz’s condition, which guarantees the existence and uniqueness of the solution of the first-order IVP.
The first-order ODE given in Equation (1) has numerous applications in different engineering and applied sciences. One of the applications is chemical kinetics, where the rate of change of concentration of a chemical species in a reaction can be modeled using a first-order ODE. In this case, the function f ( t , v ) in Equation (1) represents the reaction rate, while the solution v ( t ) represents the concentration of the species at time t. Another application of first-order ODEs can be used to model the behavior of electrical circuits and systems in electrical engineering. In this case, the function f ( t , v ) in Equation (1) represents the current or voltage in the circuit, while the solution v ( t ) represents the state of the circuit or system at time t. In summary, a function f ( t , v ) in Equation (1) represents the rate of change, and the solution v ( t ) denotes the state of the system at time t.
A numerical approximation is important for understanding the behavior of differential systems when analytical solutions are known or unknown; numerical methods have been extensively applied to solving real-life model problems in many fields, including engineering and applied sciences. Numerical approximations play a critical role in developing and validating computational models by comparing approximate results with analytical solutions. Developing numerical techniques for solving differential systems is also a significant research activity since numerical methods provide valuable insights that enable scientists and engineers to make informed decisions and improve their understanding of real-world phenomena. For these reasons, much research is devoted to developing numerical techniques to solve Equation (1).
One of the numerical techniques for solving IVPs in (1) is the Runge–Kutta methods (RKMs). Explicit Runge–Kutta Methods (ERKMs) and Implicit Runge–Kutta Methods (IRKMs) are two main subclasses of RKMs. These methods form a wide range of numerical techniques that offer acceptable accuracy and efficiency in terms of both absolute and relative errors. ERKMs and IRKMs are one-step methods that use a Taylor series expansion to approximate the solution and are easy to implement and computationally efficient. However, explicit Runge–Kutta methods are frequently inappropriate for addressing stiff differential systems. Stiff problems are those where the solution varies rapidly over a short time. ERKMs can be unstable when applied to stiff problems, leading to inaccurate results. Therefore, researchers in numerical analysis often use IRKMs to integrate stiff differential problems. However, unlike ERKMs, IRKMs require solving nonlinear equations at each time step. This system of equations is usually solved using numerical methods such as Newton’s method. IRKMs are more stable than ERKMs and can handle stiff problems more effectively. However, IRKMs are computationally more expensive than ERKMs and require more memory. In [1,2,3], the theory of RK-type methods, including the development, implementation, order condition, stability, convergence, and theoretical analysis of RK methods, are well discussed.
Multistep methods (MMs) are another numerical technique used to integrate (1). Unlike one-step methods, MMs utilize approximations of the solution at multiple mesh points to compute the subsequent approximate solutions, enhancing accuracy but demanding more function evaluations and computational time. Many experts in numerical solutions to differential equations have utilized various MMs to solve (1). The most common MMs for solving (1) are the linear multistep and the Adams–Bashforth methods. MMs can be used to solve stiff and nonstiff differential problems but require more memory and computational resources. Detailed discussions of MMs are reported in [4,5,6,7], including their development, implementation, stability, convergence, and theoretical analysis.
The literature contains many other approximation techniques for solving IVPs and related problems. These techniques are intended to offer satisfactory solutions to IVPs. Some of these approximation techniques comprise embedded Runge–Kutta type methods, general linear methods, Numerov-type methods, spline methods, finite difference, Nyström methods, space-time numerical methods, block methods, and various collocation methods. These methods are extensively discussed in various research papers and books, such as [8,9,10,11,12,13,14,15,16,17,18,19,20,21,22].
The main novelty of this paper is the introduction of a new hybrid block method (NHBM), which is capable of solving efficiently the IVP in (1). The NHBM is developed to overcome the limitations of the high number of function evaluations and stepsizes required by existing methods. The significance of the results presented in this paper is that the NHBM is an accurate and efficient approach for solving (1). The numerical experiments conducted in this study demonstrate that the NHBM outperforms other existing methods in terms of accuracy and efficiency. The NHBM also requires fewer total number of function evaluations and stepsizes, making it a more practical and computationally efficient method for solving (1). Overall, the improved accuracy and efficiency of the NHBM make it a valuable numerical method for researchers and practitioners in various fields, including engineering, physics, and applied mathematics.

2. Derivation of the Proposed Method

To approximate the theoretical solution of the problem defined in Equation (1) on any subinterval [ t n , t n + 1 ] , we can use a polynomial of degree five denoted by P ( t ) and can be expressed as follows:
v ( t ) P ( t ) = j = 0 5 k j t j ,
where coefficients k j are unknown real numbers that satisfy collocation conditions at specific nodes. The v ( t ) derivative is obtained as follows:
v ( t ) P ( t ) = j = 1 5 k j j t j 1 .
To develop the NHBM, we select three hybrid points b 1 , b 2 , and b 3 within the interval [ t n , t n + 1 ] , where 0 < b 1 < b 2 < b 3 < 1 . We then substitute t n into Equation (2) and t n , t n + b 1 , t n + b 2 , t n + b 3 , and t n + 1 into Equation (3) to obtain the following system of equations:
P ( t n ) = v n , P ( t n ) = f n , P ( t n + b 1 ) = f n + b 1 , P ( t n + b 2 ) = f n + b 2 , P ( t n + b 3 ) = f n + b 3 , P ( t n + 1 ) = f n + 1 ,
where the coefficients in Equation (2) are represented by k n , where n ranges from 0 to 5. Additionally, v n + j and f n + j represent the approximations of v ( t n + j ) and v ( t n + j ) = f ( t n + j , y ( t n + j ) ) , respectively.
Once we have determined the values of k n for n = 0 , 1 , 2 , 3 , 4 and 5, we substitute them into Equation (2) and replace the variable t with t n + x h . This substitution allows us to obtain a continuous approximation using the polynomial in Equation (2), which is given by the following equation:
P ( t n + x h ) = v n + h β 0 ( x ) f n + β b 1 ( x ) f n + b 1 + β b 2 ( x ) f n + b 2 + β b 3 ( x ) f n + b 3 + h β 1 ( x ) f n + 1 ,
where β i ( x ) i = 0 , b 1 , b 2 , b 3 , 1 are coefficients that depend on b 1 , b 2 , b 3 .
To obtain the proposed NHBM, we substitute x = 1 , x = 1 5 , 1 2 , and x = 4 5 into Equation (4). After inserting b 1 = 1 5 , b 2 = 1 2 , b 3 = 4 5 and simplifying, we obtain the following formulas:
v n + b 1 = v n + h 40,525 f n + b 1 10,112 f n + b 2 + 4525 f n + b 3 + 20,331 f n 1269 f n + 1 270,000 , v n + b 2 = v n + h 2125 f n + b 1 + 1024 f n + b 2 125 f n + b 3 + 405 f n + 27 f n + 1 6912 , v n + b 3 = v n + 2 h 2300 f n + b 1 + 2816 f n + b 2 + 1175 f n + b 3 + 567 f n 108 f n + 1 16,875 , v n + 1 = v n + 1 432 h 125 f n + b 1 + 128 f n + b 2 + 125 f n + b 3 + 27 f n + 27 f n + 1 .
The above formulas provide the approximations for v ( t n + 1 5 ) , v ( t n + 1 2 ) , v ( t n + 4 5 ) , and v ( t n + 1 ) , respectively.

3. Analysis of the Proposed Method

To evaluate the accuracy of the proposed method, we need to determine its local truncation errors and order of accuracy. This is achieved by expressing the NHBM method given in (5) in a different form by rewriting it as:
A ¯ V n = h B ¯ V n ,
where, A ¯ and B ¯ are constant matrices that are obtained from the formulas in Equation (5), and
V n = v n , v n + b 1 , v n + b 2 , v n + b 3 , v n + 1 , V n = v n , v n + b 1 , v n + b 2 , v n + b 3 , v n + 1 .
This new form of the NHBM method is a recurrence relation that relates the values of V n and V n to the matrices A ¯ and B ¯ and the stepsize h. The matrix A ¯ is derived from the coefficients of the polynomial approximation in Equation (5), while the matrix B ¯ is obtained by multiplying the coefficients of the polynomial by the corresponding values of b 1 , b 2 , b 3 , and 1. This matrix form of the NHBM method is useful for numerical implementation and analysis, as it allows us to compute the values of V n and V n recursively using matrix operations.
We then proceed to define the operator L associated with the NHBM method given in Equation (5) as follows:
L [ v ( t n + i ) ; h ] = v ( t n + i h ) F i h , v ( t n ) , v ( t n ) , v ( t n + b 1 h ) , v ( t n + b 2 h ) , v ( t n + b 3 h ) , v ( t n + h ) ,
where the operator definition in (7), involves computing a summation for all i belonging to set i = { 0 , b 1 , b 2 , b 3 , 1 } . We assume that v ( t n ) is sufficiently differentiable. To obtain the local truncation errors of the NHBM method, we expand the term in Equation (5) as a Taylor series about the point t n in powers of h. The following local truncation errors (lte) are obtained by comparing the Taylor series expansion of the exact solution with the Taylor series expansion of the numerical approximation:
l t e = 2 h 6 v ( 6 ) t n 703,125 , h 6 v ( 6 ) t n 1,152,000 , 2 h 6 v ( 6 ) t n 703,125 , h 7 v ( 7 ) t n 5,040,000 T .
Due to the fact that our method is a one-step method, it exhibits zero stability; furthermore, the NHBM has an accuracy greater than one, which means it is consistent with the differential systems under consideration. As a result, the NHBM is convergent since it is both zero-stable and consistent.

4. Error Estimation of the Proposed Method

For optimal performance, variable stepsize methods are preferred. Here, we utilize an approach to estimate local errors at step endpoints. This involves employing a lower-order formula to approximate v ( t n + 1 ) from the NHBM method, then comparing this to an approximation with a smaller stepsize. The difference helps estimate NHBM’s error in each step, aiding stepsize adjustment for accuracy within tolerance. This balances larger steps for smooth solutions and smaller ones for rapid changes, enhancing accuracy while saving computational effort. In order to determine a reasonable estimate of the local error using the proposed NHBM, we use the following fourth-order method:
v n + 1 * = v n + h 54 25 f n + b 1 + 4 f n + b 2 + 25 v n + b 3 ,
with ( l t e ) = 17 h 5 v ( 5 ) t n 4800 + O h 6 . The following E R E S T formula estimates local error:
E R E S T =   v n + 1 v n + 1 * .
In order to continue with our calculations, it is necessary to ensure that the absolute value of E R E S T is equal to or less than the specified absolute tolerances ( A T ). If this condition is met, we accept the results and proceed by increasing the stepsize. However, if E R E S T exceeds A T , we must reject the results and repeat the computations using a different strategy to improve accuracy. This computation involves adjusting the stepsize using the following formula:
h v a r = Λ × h o l d A T E R E S T p + 1 ,
where A T denotes the absolute tolerance, and 0 < Λ < 1 is an adjustment factor that prevents stepsize failures during computation. This formula allows us to adjust the stepsize based on the magnitude of the error to achieve the desired accuracy level while minimizing computational effort. An additional p is also included in (9), which denotes the accuracy of the lower-order method presented in (8).

5. Implementation

Outlined steps for the proposed NHBM to integrate problem (1) using the strategy in Section 4 are as follows:
  • Formulation: Rewrite Equation (5) as F ( v ) = 0 , with unknowns represented by V ˜ = ( v b 1 , v b 2 , v b 3 , , v 1 ) .
  • Newton’s Method: Since NHBM is implicit, use Newton’s Method (NM) to solve nonlinear equations iteratively: V ˜ i + 1 = V ˜ i J i 1 F i , where J is the Jacobian matrix of F .
  • Initial Values: Start NM with v n + j = v n + ( j h ) f n for each iteration, with n = 0 , 1 , , N 1 and j = b 1 , b 2 , b 3 , 1 .
  • Apply the NHBM to solve first-order IVPs as a system of d equations: v = f ( t , v T ) , v ( t 0 ) = v 0 , where v = ( v 1 , , v d ) T .
  • Solution and Stopping: Solve resulting system of 4 d nonlinear equations for d-dimensional IVPs. Utilize NM similarly for scalar IVPs. Implement a stopping criterion ( | V ˜ i + 1 V ˜ i | 2 × AT ) and set maximum iteration limit to 120 for NM execution.

6. Numerical Simulations

In this section, we demonstrate the numerical performance of the NHBM on several real-life models of the form (1). The formula used to measure the accuracy of the proposed method is given by:
M A E = max j = 0 , 1 , , N v ( t j ) v j .
where v ( t j ) represents the exact solution, while v j denotes the solution obtained by the NHBM at every time point represented by the discrete grid t j . This formula is used to compare the accuracy of the NHBM with existing numerical methods. Throughout this paper, different methods and performance measures are represented using the following abbreviations:
  • NHBM denotes the new hybrid block method defined in (5).
  • OOBM denotes one-step optimized block method in [23].
  • ode15s is the numerical solver available in MATLAB version R2021b.
  • RADAUIIA denotes an implicit Radau method of the fifth order in [1].
  • IGRKM represents the implicit Gauss–Legendre method of the sixth order in [4].
  • THBM is the two-step hybrid block method in [20].
  • FE denotes the total number of function evaluations.

6.1. Example 1

We first consider the SIR model problem from Kashkari and Syam [23], which arises in mathematical biology. This SIR model describes infectious disease dynamics over time. It involves the following three connected differential equations representing susceptible S ( t ) , infected I ( t ) , and recovered R ( t ) individuals:
d S d t = Δ ( 1 S ) α I S ,
d I d t = Δ I θ I + α I S ,
d R d t = Δ R + θ I ,
where Δ , θ , and α are positive parameters. Additionally, we introduce the function v ( t ) , which is defined as the sum of the susceptible, infected, and recovered individuals at time t. This function is given by:
v ( t ) = S ( t ) + I ( t ) + R ( t ) .
By adding Equation (10), Equation (11) and (12), we obtain the following expression for the derivative of v ( t ) :
v ( t ) = Δ ( 1 v ( t ) ) .
For a particular closed population, assuming an initial condition of v ( 0 ) = 0.5 and Δ = 0.5 , we can rewrite the above equation as:
v ( t ) = 0.5 ( 1 v ( t ) ) , v ( 0 ) = 0.5 t [ 0 , 1 ] , v ( t ) = 1 0.5 exp ( 0.5 t ) .
The results presented in Table 1 provide evidence of the effectiveness and better performance of the proposed method. To further demonstrate the accuracy of the method, we include graphical representations in Figure 1 and Figure 2 of the absolute errors of the NHBM, as well as a comparison of the exact and discrete solutions of the NHBM using A T = 10 8 , h 0 = 10 2 and t [ 0 , 60 ] . The plots clearly indicate that the NHBM produces highly accurate numerical solutions that closely match the exact solutions. These results highlight the effectiveness of the proposed method in accurately approximating solutions to the SIR model in mathematical biology.

6.2. Example 2

Furthermore, we solve Robertson’s stiff problem, which is a well-known differential system that arises in chemical sciences. This system describes the kinetics of an autocatalytic reaction and is known for its stiffness, making it challenging to solve numerically. This problem involves a set of three nonlinear differential equations, which are given by the following equations:
v 1 ( t ) = 10 4 v 2 ( t ) v 3 ( t ) v 1 ( t ) 25 , v 1 ( 0 ) = 1 , v 2 ( t ) = v 1 ( t ) 25 3 × 10 7 v 2 ( t ) 2 10 4 v 2 ( t ) v 3 ( t ) , v 2 ( 0 ) = 0 , v 3 ( t ) = 3 × 10 7 v 2 ( t ) 2 , v 3 ( 0 ) = 0 , 0 t 40 ,
where, v 1 ( t ) , v 2 ( t ) , and v 3 ( t ) represent the concentrations of three chemical species at time t. The system is subject to initial conditions v 1 ( 0 ) = 1 , v 2 ( 0 ) = 0 , and v 3 ( 0 ) = 0 , and is defined over the time interval 0 t 40 . To evaluate the accuracy of our method, we compare the numerical results with a reference solution obtained from [24]. Specifically, we use the maximum absolute errors to compare the obtained numerical results with the reference solution at the endpoint t N = 40 . The reference solution is given by:
v 1 ( t N ) = 0.7158270687194135 , v 2 ( t N ) = 9.185534764558135 × 10 6 ,
v 3 ( t N ) = 0.28416374574582 .
Using the above reference solution, we determine the accuracy of our method and gain insights into the behavior of Robertson’s stiff problem.
In Example 2, we compare the accuracy of the proposed NHBM method with two other commonly used methods, namely ode15s and RADAUIIA. The results are presented in Table 2, which clearly shows that the NHBM method outperforms the other two methods in terms of accuracy. By providing a more accurate and efficient approach for stiff differential systems, our study contributes to the development of new numerical methods for solving challenging problems in chemistry and chemical engineering.

6.3. Example 3

Let us consider a Jacobi elliptic function problem that arises in physics. The problem is defined by the following system of differential equations:
v 1 ( t ) = v 2 ( t ) v 3 ( t ) , v 1 ( 0 ) = 0 , v 2 ( t ) = v 1 ( t ) v 3 ( t ) , v 2 ( 0 ) = 1 , v 3 ( t ) = 1 2 v 1 ( t ) v 2 ( t ) , v 3 ( 0 ) = 1 , 0 t 50 .
The exact solution to this problem is given in terms of Jacobi elliptic functions S N , C N , D N , which can be expressed by series expansions. Specifically, the exact solution is given by:
v 1 ( t ) = J a c o b i S N [ t , 1 / 2 ] , v 2 ( t ) = J a c o b i C N [ t , 1 / 2 ] , v 3 ( t ) = J a c o b i D N [ t , 1 / 2 ] ,
where S N , C N , D N are Jacobi elliptic functions.
In Example 3, we compared the proposed NHBM method with two existing numerical methods, namely ode15s and RADAUIIA. Our comparison results are presented in Table 3, which clearly demonstrates that the NHBM method outperforms the other two methods in terms of FE and MAE. Our study contributes to developing a new efficient numerical method for solving challenging problems in physics by providing a more accurate and efficient approach for Jacobi elliptic differential systems.

6.4. Example 4

For our last example, we consider a singularly perturbed parabolic time-dependent partial differential equation that arises in chemical engineering. This problem is a reaction–diffusion process in a catalytic pellet, where the concentration of the reactant is modeled by v ( x , t ) , x [ 0 , 1 ] and t [ 0 , T ] . The reaction rate is assumed to be proportional to the concentration of the reactant, and the rate of diffusion is assumed to be small. The problem is given as follows:
ϵ v x x v t = k v 2 , ( x , t ) D = [ 0 , 1 ] × [ 0 , T ] , v ( x , 0 ) = ( 1 x ) 2 , x [ 0 , 1 ] , v ( 0 , t ) = 0 , t [ 0 , T ] , v ( 1 , t ) = 0 , t [ 0 , T ] ,
where k is a positive constant (reaction rate) and ϵ is a small positive constant (diffusion rate). Boundary conditions are at x = 0 and x = 1 , with an initial condition at t = 0 . However, the exact solution is unknown.
To solve Example 4 using the NHBM, the second-order spatial derivative is discretized using the method of lines. The following centered difference approximation is used to discretize v x x :
v x x ( x i , t ) v i + 1 2 v i + v i 1 Δ x 2 ,
where v i = v ( x i , t ) , Δ x = ( x N + 1 x 0 ) N + 1 , and N represents the internal number of spatial nodes.
To evaluate the accuracy of our proposed method, we compare NHBM with the IGRKM method. We used NDSolve, a numerical method for solving differential equations implemented in Mathematica, as the reference solution to determine the MAE. Our comparison results are presented in Table 4, which clearly demonstrates that our proposed NHBM method outperforms the IGRKM method in terms of accuracy. To further demonstrate the effectiveness of our approach, we present surface plots of the proposed NHBM and NDSolve solutions for Example 4 in Figure 3. The results in Table 4 and Figure 3 distinctly showcase the proposed method’s efficacy in accurately approximating solutions to the reaction–diffusion process that arises in the chemical engineering domain.

7. Concluding Remarks

This paper proposes a new numerical method that utilizes a variable stepsize formulation to solve differential systems in (1) efficiently. We also discuss a technique to control the stepsize of the method to enhance solution accuracy. The efficacy of the proposed approach is validated through its application to four distinct model differential problems within applied sciences and engineering, and the results are compared to existing solutions. The variable stepsize implementation used in this paper produces better performance, as evidenced by the computed results. We conclude that the NHBM can be a valuable numerical method for solving first-order differential model problems in various fields.

Author Contributions

The proposed method was developed and analyzed by M.A.R.; The experiment section and figure plotting were carried out by M.A.R. and H.R.; Conceptualization, M.A.R., B.C. and H.R.; methodology, M.A.R., H.R. and B.C.; software, M.A.R. and H.R.; validation, M.A.R., H.R. and B.C.; formal analysis, M.A.R., H.R. and B.C.; investigation, M.A.R., H.R. and B.C.; writing—original draft preparation, M.A.R. and B.C.; writing—review and editing, H.R., M.A.R. and B.C.; supervision, H.R., M.A.R. and B.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Acknowledgments

The first two authors sincerely thank Provincia Autonoma di Bolzano/Alto Adige—Ripartizione Innovazione, Ricerca, Universita e Musei for their invaluable support in providing an exceptional research environment and good research facilities. We are also very grateful to the anonymous reviewers for their valuable comments that helped improve the article’s quality.

Conflicts of Interest

The authors declare that they have no conflict of interest.

References

  1. Butcher, J.C. Numerical Methods for Ordinary Differential Equations; John Wiley & Sons: New York, NY, USA, 2008. [Google Scholar]
  2. Lambert, J.D. Numerical Methods for Ordinary Differential Systems: The Initial Value Problem; John Wiley & Sons: New York, NY, USA, 1991. [Google Scholar]
  3. Ascher, U.M.; Petzold, L.R. Computer Methods for Ordinary Differential Equations and Differential-Algebraic Equations; Society for Industrial and Applied Mathematics (SIAM): Philadelphia, PA, USA, 1998. [Google Scholar]
  4. Gustafsson, S. Numerical Methods for Differential Equations: An Introduction to Scientific Computing; Springer: Berlin/Heidelberg, Germany, 2017. [Google Scholar]
  5. Brugnano, L.; Trigiante, D. Solving Differential Problems by Multistep Initial and Boundary Value Methods; Gordon and Breach Science Publishers: Langhorne, PA, USA, 1998. [Google Scholar]
  6. Atkinson, K.E. Elementary Numerical Analysis, 2nd ed.; John Wiley & Sons: New York, NY, USA, 1993. [Google Scholar]
  7. Atkinson, K.E.; Han, W.; Stewart, D. Numerical Solution of Ordinary Differential Equations; John Wiley and Sons: Hoboken, NJ, USA, 2009. [Google Scholar]
  8. Hairer, E.; Norsett, S.P.; Wanner, G. Solving Ordinary Differential Equations I, Nonstiff Problems; Springer: Berlin, Germany; New York, NY, USA, 1993. [Google Scholar]
  9. Mehdizadeh, K.M.; Shokri, A.; Wang, Y.; Bazm, S.; Navidifar, G.; Khakzad, P. Qualitatively Stable Schemes for the Black–Scholes Equation. Fractal Fract. 2023, 7, 154. [Google Scholar] [CrossRef]
  10. Piao, X.; Bu, S.; Kim, D.; Kim, P. An embedded formula of the Chebyshev collocation method for stiff problems. J. Comput. Phys. 2017, 351, 376–391. [Google Scholar] [CrossRef]
  11. Kovalnogov, V.N.; Fedorov, R.V.; Simos, T.E. New hybrid symmetric two step scheme with optimized characteristics for second order problems. J. Math. Chem. 2018, 20, 1–29. [Google Scholar] [CrossRef]
  12. Rufai, M.A.; Ramos, H. Numerical solution of third-order boundary value problems by using a two-step hybrid block method with a fourth derivative. Comput. Math. Methods 2021, 3, e1166. [Google Scholar] [CrossRef]
  13. Rufai, M.A.; Ramos, H. Numerical solution of Bratu’s and related problems using a third derivative hybrid block method. Comput. Appl. Math. 2020, 39, 322. [Google Scholar] [CrossRef]
  14. Hoffman, J.D. Numerical Methods for Engineers and Scientists, 3rd ed.; Marcel Dekker, Inc.: New York, NY, USA, 2001; Chapter 8. [Google Scholar]
  15. Anastassi, Z.A.; Kosti, A.A.; Rufai, M.A. A Parametric Method Optimised for the Solution of the (2+1)-Dimensional Nonlinear Schrödinger Equation. Mathematics 2023, 11, 609. [Google Scholar] [CrossRef]
  16. Rufai, M.A.; Ramos, H. Numerical Solution for Singular Boundary Value Problems Using a Pair of Hybrid Nyström Techniques. Axioms 2021, 10, 202. [Google Scholar] [CrossRef]
  17. Medvedeva, M.A.; Simos, T.E.; Tsitouras, C. Variable step-size implementation of sixth-order Numerov-type methods. Math. Methods Appl. Sci. 2020, 43, 1204–1215. [Google Scholar] [CrossRef]
  18. Rufai, M.A. An efficient third derivative hybrid block technique for the solution of second-order BVPs. Mathematics 2022, 10, 3692. [Google Scholar] [CrossRef]
  19. Rufai, M.A.; Ramos, H. An adaptive one-point second-derivative Lobatto-type method for solving efficiently differential systems. Int. J. Comput. Math. 2022, 99, 1687–1705. [Google Scholar]
  20. Ajileye, G.; Amoo, S.A.; Ogwumu, O.D. Hybrid Block Method Algorithms for Solution of First Order Initial Value Problems in Ordinary Differential Equations. J. Appl. Comput. Math. 2019, 7, 2. [Google Scholar]
  21. Uddin, M.; Ali, H. The Space–Time Kernel-Based Numerical Method for Burgers’ Equations. Mathematics 2018, 6, 212. [Google Scholar] [CrossRef]
  22. Hairer, E.; Wanner, G. Solving Ordinary Differential Equations II; Volume 14 of Springer Series in Computational Mathematics; Springer: Berlin, Germany, 2010. [Google Scholar]
  23. Kashkari, B.S.; Syam, M.I. Optimization of one step block method with three hybrid points for solving first-order ordinary differential equations. Results Phys. 2019, 12, 592–596. [Google Scholar] [CrossRef]
  24. Mazzia, F.; Iavernaro, F. Test Set for Initial Value Problems Solvers; University of Bari: Bari, Italy, 2003. [Google Scholar]
Figure 1. Absolute errors for Example 1.
Figure 1. Absolute errors for Example 1.
Fractalfract 07 00703 g001
Figure 2. Exact solution (blue line) and NHBM solution (red points), for Example 1.
Figure 2. Exact solution (blue line) and NHBM solution (red points), for Example 1.
Fractalfract 07 00703 g002
Figure 3. Surface plots of the proposed NHMB and NDSolve for Example 4 with h 0 = 10 2 , ϵ = 10 2 , N = 29 , k = 0.5 .
Figure 3. Surface plots of the proposed NHMB and NDSolve for Example 4 with h 0 = 10 2 , ϵ = 10 2 , N = 29 , k = 0.5 .
Fractalfract 07 00703 g003
Table 1. Comparison of absolute errors for Example 1 with h 0 = 10 2 .
Table 1. Comparison of absolute errors for Example 1 with h 0 = 10 2 .
tNHBMOOBMTHBM
0.1000 4.4409 × 10 16 1.9984 × 10 15 1.2180 × 10 13
0.2000 8.8180 × 10 16 3.8856 × 10 15 1.3999 × 10 13
0.3000 1.2213 × 10 15 5.4401 × 10 15 1.1849 × 10 12
0.4000 1.6653 × 10 15 6.9944 × 10 15 1.5389 × 10 12
0.5000 1.7764 × 10 15 8.2157 × 10 15 1.1100 × 10 12
0.6000 2.1094 × 10 15 9.5479 × 10 15 5.2702 × 10 12
0.7000 2.3315 × 10 15 1.0547 × 10 14 2.1090 × 10 12
0.8000 2.4425 × 10 15 1.1324 × 10 14 1.2979 × 10 11
0.9000 2.6645 × 10 15 1.2213 × 10 14 3.0823 × 10 11
1.0000 2.7756 × 10 15 1.2879 × 10 14 4.1219 × 10 11
Table 2. Comparison of numerical solutions for Example 2.
Table 2. Comparison of numerical solutions for Example 2.
h 0 ATMethodsFEMAE
10 2 10 9 NHBM275 3.9830 × 10 11
10 2 10 9 ode15s351 2.6289 × 10 9
10 2 10 9 RADAUIIA504 5.2187 × 10 10
10 3 10 10 NHBM420 2.5190 × 10 12
10 3 10 10 ode15s475 1.1331 × 10 10
10 3 10 10 RADAUIIA735 6.0419 × 10 12
Table 3. Comparison of numerical solutions for Example 3.
Table 3. Comparison of numerical solutions for Example 3.
h 0 ATMethodsFEMAE
10 1 10 4 NHBM405 7.5079 × 10 5
10 1 10 4 ode15s646 8.8513 × 10 3
10 1 10 4 RADAUIIA1115 51.0737 × 10 3
10 2 10 5 NHBM630 7.5703 × 10 6
10 2 10 5 ode15s771 1.0373 × 10 3
10 2 10 5 RADAUIIA1323 2.6296 × 10 4
Table 4. Comparison of numerical solutions for Example 4 using h 0 = 10 2 , ϵ = 10 2 , N = 29 , k = 0.5 .
Table 4. Comparison of numerical solutions for Example 4 using h 0 = 10 2 , ϵ = 10 2 , N = 29 , k = 0.5 .
h 0 ATMethodsFEMAE
10 2 10 3 NHBM40 5.8394 × 10 7
10 2 10 3 IGRKM204 4.0275 × 10 3
10 2 10 4 NHBM55 6.5863 × 10 8
10 2 10 4 IGRKM576 1.2977 × 10 4
10 2 10 5 NHBM75 6.9609 × 10 9
10 2 10 5 IGRKM1740 4.1942 × 10 5
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

Rufai, M.A.; Carpentieri, B.; Ramos, H. A New Hybrid Block Method for Solving First-Order Differential System Models in Applied Sciences and Engineering. Fractal Fract. 2023, 7, 703. https://doi.org/10.3390/fractalfract7100703

AMA Style

Rufai MA, Carpentieri B, Ramos H. A New Hybrid Block Method for Solving First-Order Differential System Models in Applied Sciences and Engineering. Fractal and Fractional. 2023; 7(10):703. https://doi.org/10.3390/fractalfract7100703

Chicago/Turabian Style

Rufai, Mufutau Ajani, Bruno Carpentieri, and Higinio Ramos. 2023. "A New Hybrid Block Method for Solving First-Order Differential System Models in Applied Sciences and Engineering" Fractal and Fractional 7, no. 10: 703. https://doi.org/10.3390/fractalfract7100703

APA Style

Rufai, M. A., Carpentieri, B., & Ramos, H. (2023). A New Hybrid Block Method for Solving First-Order Differential System Models in Applied Sciences and Engineering. Fractal and Fractional, 7(10), 703. https://doi.org/10.3390/fractalfract7100703

Article Metrics

Back to TopTop