Dear Label_Name,

Greetings! I hope you are well.

I am writing on behalf of an international publishing house, Science and Education Publishing, USA (http://www.sciepub.com).

Science and Education Publishing (SciEP) is a multidisciplinary open access publisher of Journals covering the fields of Science, Technology and Medicine. All SciEP publications are subject to high-quality peer review, editorial and production processes and are freely available online.

Welcome to the April, 2013 edition of SciEP newsletter, providing monthly summaries of articles on science & education related topics of interest to scientist and professionals. The mission of our journals is to provide an open access forum for publishing new findings on issues pertinent to scientific and academic field. We also invite you to submit original research papers, review articles, short reports, or case studies are considered for publication in our Multidisciplinary Journals.

American Journal of Public Health Research
DOI: 10.12691/ajphr-1-1-3

The Implementation of a Tele-Homecare System with Service Oriented Architecture and HL7 Message Transmission Standard
Ching-Sung Wang

American Journal of Medical and Biological Research
DOI: 10.12691/ajmbr-1-1-2

The Influence of Ethanol on Pyruvate Kinases Activity in Vivo, in Vitro, in Silico
Lelevich Sergey Vladimirovich, Khrustalev Vladislav Victorovich, Barkovsky Eugene Victorovich, Shedogubova Tatyana Aleksandrovna

American Journal of Medical Sciences and Medicine
DOI: 10.12691/ajmsm-1-1-1

Clinical and Histological Features in Wegener Granulomatosis
Li Cavoli Gioacchino, Passantino Rita, Bono Luisa, Tortorici Calogera, Ferrantelli Angelo, Giammarresi Carlo, Rotolo Ugo

American Journal of Infectious Diseases and Microbiology
DOI: 10.12691/ajidm-1-1-2

The Curse of Lymphatic Filariasis: Would the Continual Use of Diethylcarbamazine Eliminate this Scourge in Papua New Guinea?
Anne Lanham, Lillian Mwanri

We have released issues of some journals on February 2013.
American Journal of Public Health Research
American Journal of Medical Sciences and Medicine
American Journal of Medical and Biological Research
American Journal of Pharmacological Sciences
American Journal of Clinical Medicine Research
American Journal of Nursing Research
American Journal of Cancer Prevention
American Journal of Medicine Studies
American Journal of Infectious Diseases and Microbiology
American Journal of Cardiovascular Disease Research
American Journal of Biomedical Research
American Journal of Mechanical Engineering
American Journal of Microbiological Research
Journal of Preventive Medicine
American Journal of Food and Nutrition
American Journal of Hypertension Research
American Journal of Sports Science and Medicine
Journal of Cancer Research and Treatment
American Journal of Zoological Research
Journal of Biomedical Engineering and Technology
American Journal of Microbiological Research
Biomedical Science and Engineering

We are pleased to announce that submissions are now being accepted for peer review for all the journals across all our subject areas in SciEP.

Please submit your manuscript via our online submission system:
http://www.sciepub.com/author/signin

To unsubscribe from our mailing list, please click on the following link: unsubscribe

I am looking forward to hearing from you.

Kind regards,

David Rain
Email:join@sciepub.com
Science and Education Publishing (http://www.sciepub.com)

Solving the Boundary Value Problems for Second-Order Ordinary Differential Equations using a Class of Obrechkoff-Type Block Methods

Taiwo E. Fayode, Bola T. Olabode, Emmanuel A. Areo, Ezekiel O. Omole, Oluwasegun M. Ibrahim*

University of Texas at Austin, USA

Abstract: Employing a combination of collocation and interpolation techniques, this research introduces a novel set of Obrechkoff-type methods designed to address second-order boundary value problems (BVPs) characterized by Neumann and Dirichlet boundary conditions. The methodology involves placing the derivative function equation at all grid points and interpolating the basis function at only two locations, resulting in the development of a series of continuous multistep techniques with variable step numbers. For numerical implementation, our approach employs block mode techniques. We investigate the order, consistency, stability, and convergence of these algorithms to ensure their robustness and reliability. To evaluate the efficacy and precision of the proposed method, comprehensive testing is conducted using Obrechkoff-type problems. Notably, the numerical solutions demonstrate improved performance compared to traditional methods, highlighting the potential of our approach to deliver enhanced accuracy in solving second-order BVPs.

Keywords: Boundary value problems, interpolation, collocation, differential equations, numerical methods.

Received May 13, 2025; Revised June 15, 2025; Accepted June 23, 2025

Citation: Taiwo E. Fayode, Bola T. Olabode, Emmanuel A. Areo, Ezekiel O. Omole, Oluwasegun M. Ibrahim, Solving the Boundary Value Problems for Second-Order Ordinary Differential Equations using a Class of Obrechkoff-Type Block Methods. American Journal of Applied Mathematics and Statistics, 13 (2). 30-37. January 2025. doi:10.12691/ajams-13-2-2

1. Introduction

Mathematics is widely utilized and often considered a prerequisite for progress in various fields, such as biology, economics, chemical kinetics, circuit theory, and others. Physical phenomena in these disciplines are frequently examined through the lens of mathematical models, which often yield equations involving derivatives of unidentified functions in one or more variables commonly known as differential equations.Interestingly, differential equations arising from the modeling of physical events often lack analytical solutions. They form the basis for modeling a diverse range of problems across sciences, technology, engineering, social sciences, and more. Consequently, the development of numerical techniques to approximate solutions becomes imperative. Various numerical methods, including finite difference techniques, finite element techniques, and finite volume techniques, have been devised to address the nature and type of the differential equations at hand.A boundary value problem comprises a differential equation along with additional constraints referred to as boundary conditions. A Dirichlet boundary condition, also known as a first-type boundary condition, specifies the values of the functions themselves. On the other hand, a Neumann boundary condition provides details about the normal derivative of the function at the boundary. A solution to a boundary value problem is a solution to the associated differential equation that also satisfies the specified boundary conditions. The main goal of this study is to numerically solve the boundary value problems of the second-order ordinary differential equations (ODEs) of the type

[picture] | (1)

Where are continuous functions and m is the system dimension. There are numerous engineering and applied science fields that deal with these second-order boundary value problems. Numerous authors have worked to solve equation (1) utilizing numerical approaches for resolving boundary value problems, which can essentially be divided into the shooting method, finite element method, and finite difference method. The technique utilized can be divided into two categories: linear multistep methods (LMMs) and one-step (single-step) techniques; see [1,2,3,4,5,6,7]. Evidently, approaches of the Obreckhoff-type are effective for both start and boundary value problems. The well-known Numerov method is the simplest of the Obreckhoff-type approaches. For second-order initial value problems, authors have created a number of Obreckhoff-type approaches with higher orders in the past. Consider the category of second-order boundary value problem with Dirichlet boundary conditions in equation (1),

[picture] | (2)

and Neumann boundary conditions

[picture] | (3)

Numerous branches of science, engineering, and medicine encounter second-order boundary-value problems. For instance, they can be observed in the bending of Euler-Bernoulli beams Fernandez-Sez et al. [8] and Li et al. [9], curved cantilever beams Ghuku and Saha [10] homogenous-heterogenous flow of fluids over surfaces, Malik and Khan [11], and Hayat et al. [12] among other phenomena. Numerous techniques, including the shooting technique Filipov et al. [13] , Block Nystrom type method Jator and Manathunga [14] , Quartic B-spline method Goh et al. [15]. Boubaker polynomials expansion strategy Koak et al. [16], and the Adomian decomposition method, Aly et al. [17]. The resolution of second-order initial value problems given as follows has been studied by a number of scholars.

(4)

In reality, problems of kind (4) can be found in circuit theory, biology, celestial mechanics, chemical kinetics, and celestial mechanics. The second-order problem has been widely studied; see Ibrahim and Ikhile [1], Ibrahim and Ikhile [2], Ibrahim and Ikhile [3], Ibrahim and Ikhile [4] and other references herein.

Obrechkoff [18] introduced a series of techniques for solving first-order IVPs in 1942 as

(5)

Classical Obreckhoff methods were listed by Lambert [27]. By enhancing and extending Lambert’s work on stability as Boundary Value Methods (BVMs), Ghelardoni [20], expanded its scope, while some BVPs were considered in Ibrahim et al. [21]; Ibrahim and Lawrence [22] . The step length was set at 2, and Usmani’s formula [23] was used as the starting technique.

(6)

Lambert [24] observed that, for large , it is impossible to satisfy the zero-stability of (5), as the error constant decreases more rapidly with an increase in 1 than with . As a result, its stability domain appears to be limited. Investigators have since conducted studies in this area. Shokri [26] separately created higher-order L-stable Obreckhoff methods with varied step sizes that are precise and effective and produce a bigger stability domain. A fascinating thing about Obreckhoff methods is that they can be extended to solve higher-order ordinary differential equations. Considering this, numerous authors have expanded the Obreckhoff techniques (6) to

(7)

where the aforementioned equation successfully and reliably resolves second-order initial value problems. The inclusion of higher derivatives makes this family of approaches stand out as distinct. The -step Obrechkoff method’s general form using the derivative is provided by Lambert [27]. These methods, known as Obreckhoff methods, can be particularly effective with the first few derivatives. Although Obreckhoff’s original work primarily addressed numerical quadrature, it appears that Milne [28] was the first to promote the application of the Obreckhoff formula. However, before implementing the method described in (6), additional one-step numerical methods like the Euler or Trapezoidal method must be introduced. This then can be translated to the block method, Ajayi et al. [25].

2. Derivation of the Block Method

The Block Method is derived with the intention of approximating the exact solution y(x) in the partition for of a power series polynomial of the form,

(8)

where the numbers of collocation points and interpolating points, respectively, are and . The second and fourth derivatives of (8) are thus

(9)

(10)

The interpolation, collocation, and collocation equations in equations (8), (9), and (10) will be used to develop the suggested approaches.

The system of equations arriving from (8), (9) and (10) will be expressed in the form of matrix equation as given below.

(11)

where

(12)

is resolved for the parameters s using the Gaussian elimination approach. Equation (8) requires the parameters to be swapped in order to produce continuous schemes. The approach is then determined by evaluating the continuous schemes at the non-interpolating point.

This results in the equations below

(13)

Combining (13) into (11) in matrix form gives

noting

[picture] | (14)

Solving (14) for the values of parameters aj′ s we obtain

(15)

A linear multistep technique with a continuous coefficient of this type is produced by substituting (15) into (8) to yield

(16)

The coefficients of (16) were obtained as follows, using the transformation

(17)

where

(18)

Evaluating (16) at , gives the discrete scheme below

(19)

The first derivatives of (16) gives

(20)

where

(21)

Evaluating (21) at and 1 gives the following first derivative schemes:

(22)

(23)

(24)

2.1. Block form of One-step Obreckhoff-type Method

Combining (19) and (23) give the block method below

(25)

Solving (25) using matrix inversion gives,

(26)

Writing out (26) explicitly gives

(27)

(28)

Substituting (27) and (28) into (22) and (24) gives

(29)

(30)

3. Analysis of the Method

This section presents the analysis of the properties of the method as follows.

3.1. Local Truncation Error and Order

The developed system order was investigated using the methodology described in Lambert (1973) and Fatunla (1991). Equation (31) displays the error constants as having order and

(31)

While the block method’s order was looked at, it has order , and the error constants are represented in Equation (32) as

(32)

3.2. Consistency

If a linear multistep approach converges with an order greater than 2, it is considered to be consistent. That is, ().

3.3. Zero Stability

For any well-behaved starting value problem, a linear multistep approach is zero-stable if all roots of the first characteristics polynomial lies in the unit disk, and any roots on the unit circle are simple,

(33)

Equation (33) produces when set to zero, proving that the approach is zero stable.

3.4. Convergence

Consistency and zero stability are the required and sufficient conditions for our method to be convergent.

3.5. Region of Absolute Stability

As defined by Jator and Manathunga [14] and Ajayi et al. [25], the region of absolute stability of the established schemes is taken into consideration.

4. Numerical Application on Dirichlet Boundary Conditions

In order to evaluate the effectiveness of the produced methods, we consider the following problems and their sources, Majid et al. [29] and Chew et al. [30].

Problem 1: One-Step Obreckhoff Block Method (1SOBM):

Exact solution:

Problem 2: Direct Adams-Moulton Method (DAMD)

Exact Solution:

Table 1. Showing the comparison of error in the method (1SOMBM) with Chew et al. for Problem 1

Table 2. Showing the comparison of error in the method (1SOMBM) with Chew et al. for Problem 2

Table 3. Showing the exact solutions, the computed results and the error from the 1SOMBM for Problem 1

Table 4. Showing the exact solutions, the computed results and the error from the 1SOMBM for Problem 2

5. Summary and Conclusion

Through the development, rigorous analysis, and effective application of a sophisticated class of Obrechkoff-type block methods have achieved direct solutions for second-order boundary value problems 1 and 2. Our investigation delved into the comparison between Dirichlet and Neumann boundary conditions for second-order boundary value problems within ordinary differential equations. Leveraging continuous coefficients and ensuring exceptional consistency, our collocation approach has given rise to implicit methods of notable accuracy. Notably, these methods demonstrate a remarkable feature: they are implemented seamlessly without the necessity for the development of predictors, eliminating the need for additional methods to generate starting values. The comparison profiles for problems 1 and 2 are visually presented in Figures 1 and 2 above, providing a clear representation of the methods’ performance.

What sets our methods apart is their high order, accuracy, and low error constants. This intrinsic superiority contributes to the methods’ exceptional proficiency, allowing them to achieve a level of perfection that surpasses conventional approaches in Chew et al. [30]. The culmination of these attributes establishes our Obrekhoff-type block methods as an impressive and reliable solution for addressing second-order boundary value problems [31].

Funding

This work did not receive any funding.

Conflict of Interest

The authors declare that there is no conflict of interest regarding the publication of this paper.

Availability of Data and Material

No specific data or unique material was used for this work

[1]   | Ibrahim, O. M, and Ikile, M.N.O. (2020a): A generalized family of symmetric multistep methods with minimal phase-lag for initial value problems in ordinary differential equations. Mediterranean Journal of Mathematics 17, 1-30. https:// link.springer.com/ article/10. 1007/s00009-020-01507-5.

[2]  
Ibrahim, O.M., and Ikile, M.N.O. (2020b): Inverse hybrid linear multistep methods for solving the second order initial value problems in ordinary differential equations. International Journal of Applied and Computational Mathematics, 6(6), 158. https:// link.springer.com/article/10.1007/s40819-020-00910-6.

[3]   | Ibrahim, O.M., and Ikile, M.N.O (2017a): On the construction of high accuracy symmetric super-implicit hybrid formulas with phase-lag properties. Transaction of the Nigerian Association of Mathematical Physics 4, 101-108. https:// www.researchgate.net/ publication/ 377330540.
[4]   | Ibrahim, O.M., and Ikile, M.N.O. (2017a): Highly stable super-implicit hybrid methods for special second-order IVPs. American Journal of Applied Scientific Research, 3(3), 21-27.
[5]   | Omole, E.O., Adeyefa, E.O., Ayodele, V.I., Shokri, A., Wang, Y. (2023) Ninth-order Multistep Collocation Formulas for Solving Models of PDEs Arising in Fluid Dynamics: Design and Implementation Strategies. Axioms 2023, 12, 891.
[6]   | E. Aeao and M. T. Omojola (2015): A new one-twelfth step continuous block method for the solution of modeled problems of ordinary differential equation, American Journal of Mathematics. 5(04), 447-450.
[7]   | Rufai, M.A., Shokri, A., Omole, E.O. A (2023): One-Point Third-Derivative Hybrid Multistep Technique for Solving Second-Order Oscillatory and Periodic Problems. Hindawi J. Math. 2023, 2023, 2343215.
[8]   | Fernandez-Saez, J., Zaera, R., Loya, J. A., & Reddy, J. (2016). Bending of Euler-Bernoulli beams using Eringen’s integral formulation: a paradox resolved. International Journal of Engineering Science, 99, 107-116.
[9]   | Li, S. R., Cao, D. F., & Wan, Z. Q. (2013). Bending solutions of FGM Timoshenko beams from those of the homogenous Euler–Bernoulli beams. Applied Mathematical Modelling, 37(10-11), 7077-7085.
[10]   | Ghuku S, K.N Saha (2016): A theoretical and experimental study on geometric nonlinearity of initially curved cantilever beams. Engineering Science and Technology, an International Journal 19 (1): 135-146.
[11]   | Malik R, M. Khan (2018): Numerical study of homogenous heterogeneous reactions in Sisko fluid flow past a stretching cylinder. Results in Physics 8 (2018) 415 121.
[12]   | Hayat T, K. Muhammed, A. Alsaedi, and S. Asghar (2018): Numerical study for melting heat transfer and homogenous-heterogenous reactions in flow involving carbon nanotubes, Results in Physics 8 (2018) 415 121.
[13]   | Filipov, S. M., Gospodinov, I. D., & Faragó, I. (2017). Shooting-projection method for two-point boundary value problems. Applied Mathematics Letters, 72, 10-15.
[14]   | Jator, S. N., & Manathunga, V. (2018). Block Nyström type integrator for Bratu’s equation. Journal of Computational and Applied Mathematics, 327, 341-349.
[15]   | Goh, J., Majid, A. A., & Ismail, A. I. M. (2012). A quartic B-spline for second-order singular boundary value problems. Computers & Mathematics with Applications, 64(2), 115-120.
[16]   | Koçak, H., Yıldırım, A., Zhang, D. H., & Mohyud-Din, S. T. (2011). The Comparative Boubaker Polynomials Expansion Scheme (BPES) and Homotopy Perturbation Method (HPM) for solving a standard nonlinear second-order boundary value problem. Mathematical and Computer Modelling, 54(1-2), 417-422.
[17]   | Aly E. H., Ebaid A, Rach R. (2012): Advances in the Adomian decomposition method for solving two-point nonlinear Boundary Value Problems with Neumann boundary conditions. Computers and Mathematics with Applications. 2012; 63(6): 1056–65.
[18]   | Obrechkoff N (1942): On Mechanical quadrature (Bulgarian French summary). Spisanie Bulgar Akad. Nauk. 65, 191-289.
[19]   | Lambert J.D. (1973): Computational methods in ordinary differential equations. New York.
[20]   | Ghelardoni. P. Marzulli (1995): Stability of some boundary value methods for IVPs, Ap plied Numerical Mathematics 18 (1) (1995) 144-153.
[21]   | Ibrahim, O.M., and Ikile, M.N.O (2017): Spectral Collocation Method for the Numerical Solution of Ordinary Differential Equations. Transaction of the Nigerian Association of Mathematical Physics 4, 95-100. https://www.researchgate.net/ profile/Oluwasegun-Ibrahim/publication/319763706.
[22]   | Ibrahim, O.M., and Lawrence, P. W. (2019): Spectral rectangular collocation formula: an approach for solving oscillatory initial value problems and/or boundary value problems in ordinary differential equations. Turkish J. Anal. Number Theory, 7(1), 11-17.
[23]   | Usmani R.A (1966): Boundary value techniques for the numerical solution of certain initial value problems in ordinary differential equations J. ACM 13 (1966) 287 295.
[24]   | Lambert J. D. and Mitchel A. R. (1962):On the solution of y’= f(x,y) by a class of high-accuracy differential formulae of low order, Z. Angew. Math. Phys. 13, 223-232.
[25]   | Ajayi, S. A., Muka, K. O., and Ibrahim, O. M. (2019): A family of stiffly stable second derivative block methods for initial value problems in ordinary differential equations. Earthline Journal of Mathematical Sciences, 1(2), 221-239.
[26]   | Shokri A (2013): The new class of implicit 1-stable hybrid Obrechkoff method for the numerical solution of first-order initial value problems. Computer Physics Communications 184 (3) (2013) 529-531.
[27]   | Lambert J.D (1973): Computational methods in ordinary differential equations. John Wiley.
[28]   | Milne, W.E. (1949): A note on the numerical integration of differential equations, J. Res. Nat. Bur. Standards, 43, 537-542 (1949).
[29]   | Majid Z.A., M.M. Hasni, and N. Senu (2013): “Solving second-order linear dirichlet and Neumann boundary value problems by block method,” IAENG International Journal of Applied Mathematics, vol. 43, no. 2, pp. 71-76.
[30]   | Chew K. T., Zanariah A. M., Mohamed S. and Norazak S. (2012): Solving Linear Two-Point Boundary Value Problems by Direct Adams Moulton Method, Applied Mathematical Sciences, Vol. 6, no. 99, 4921-4929.

[31]  
Ibrahim, O. M. (2017). High order symmetric super-implicit hybrid LMMs with minimal phase-lag: Solving differential equations using the numerically stable Super-Implicit Hybrid formula. LAMBERT Academic Publishing. ISBN 978-6202023795. https:// www.amazon.ca/ Symmetric-Super-Implicit-Hybrid-Minimal-Phase-Lag/dp/6202023791.