Mathematical Problems in Engineering

Mathematical Problems in Engineering / 2020 / Article

Research Article | Open Access

Volume 2020 |Article ID 6351279 |

Reem Allogmany, Fudziah Ismail, Zanariah Abdul Majid, Zarina Bibi Ibrahim, "Implicit Two-Point Block Method for Solving Fourth-Order Initial Value Problem Directly with Application", Mathematical Problems in Engineering, vol. 2020, Article ID 6351279, 13 pages, 2020.

Implicit Two-Point Block Method for Solving Fourth-Order Initial Value Problem Directly with Application

Academic Editor: Stefan Balint
Received18 Dec 2019
Revised14 Jul 2020
Accepted10 Aug 2020
Published15 Sep 2020


This paper proposes an implicit block method with two-point to directly solve the fourth-order Initial Value Problems (IVPs). The implicit block method is derived by adopting Hermite interpolating polynomial as the basis function, incorporating the first derivative of to enhance the solution’s accuracy. A block formulation is presented to acquire the numerical approximation at two points simultaneously. The introduced method’s basic properties, including order, zero stability, and convergence, are presented. Numerical experiments are carried out to verify the accuracy and efficiency of the proposed method compared with those of the several existing methods. Application in ship dynamics is also presented which yield impressive results for the proposed two-point block method.

1. Introduction

Higher order differential equations (ODEs) can be used to model problems arising from the field of applied sciences and engineering [1]. The problem of static deflection of a uniform beam which can be modelled as a fourth-order initial value problem (IVP) is a good example of a real problem in engineering [2, 3]. The classical way of solving them is by reducing the equation into the system of first-order ODEs, but this process is too rigorous compared with the direct methods [4, 5]. In addition, many researchers have presented direct methods to avoid the reduction effort [615]. To enhance the efficacy of numerical methods, the block method is introduced with the idea of producing simultaneously -point of the approximate solutions at one time step. Kuboye and Omar [7] and Jacob [8] introduced efficient zero-stable numerical block methodologies for fourth-order initial value problems (IVPs). For block multistep method, approximations are used to generate the approximate solutions at the new block. Omar and Kuboye [9] proposed a block method for directly solving general fourth-order IVPs by increasing the step number; Adeyeye and Omar [14] also presented a direct six-step block method. Recently, researchers have been exploring numerical methods with more functional evaluations with the aim of obtaining numerical solutions with very high precision. This is how hybrid methods in nature arises. Jator [16] and Yap and Ismail [10] developed four-step hybrid block methods, and Abdelrahim and Omar [11] generalized three off-step points. Kayode et al. [12] modified the implicit hybrid block method to directly solve the IVPs associated with fourth-order ODEs. Generally, collocation and interpolation techniques are utilized as direct methods. The points need to be collocated and interpolated after which a system of linear equations must be resolved in order to obtain the method’s coefficients. Therefore, we developed a two-point block the implicit method by using a strategy in our proposed method which can be implemented in a straightforward manner. The method replaces the function by interpolating and integrating the polynomial. Our derivation only requires us to do interpolation and integration to obtain the coefficients. The method also involves the fifth derivative of the solution to equation (1), which aims at acquiring better accuracy. In general, the accuracy of a method increases with increase in the order of the method. But, with the idea of incorporation of higher derivative of the solution in the process, higher and better accuracy is achieved without corresponding increase in the order of the method. We are concerned with the development of the numerical solution of IVPs for fourth-order ODEs of the form:where the fifth derivative of the solution to equation (1) is

We assume that in equation (1) is differentiable to a desired order in region and satisfies Lipchitz condition in its second, third, fourth, and fifth terms as follows:for all points , , , and , in the region R. Then, the IVPs in equation (1) have a unique solution in R ([17, 18]).

The paper is structured in the following manner: in Section 2, a discussion is undertaken about the derivation of the block method of order six (I2PBDO6) along with the basic idea on how the block works. The implementation of this method to solve general fourth-order ODEs is proposed in Section 3. Numerical experiments and the application of the ship dynamics problem are given in Section 4 to show the accuracy and efficiency of the proposed method. The study’s conclusion is finally provided in Section 5.

2. Derivation of Method

This section depicts the proposed method’s derivation based on Hermite Interpolating Polynomial , which interpolates at two points. The polynomial has the formwhereand ; is the generalized Lagrange polynomial, , and is a positive integer.

The method computes the approximate solutions and concurrently at and , where becomes the starting point and is the last point in the block with step size . The evaluation solution of at the point will be used as the next iteration’s initial value:where is the fifth derivative of the solution and

The approximate solution at the point can be retrieved by integrating equation (1) once, twice, thrice, and four times, in relation to over the interval . The following formula can be reached:

Let which gives

Then, in equations (9)–(12) will be replaced by Hermite interpolating polynomial .

Taking leads to and and then making changes in the integration limit from to , we havewhere

Evaluating the integrals in equations (13)–(16) gives

Now, we have four formulas, one for the approximate solution and the others for the approximation of the first, second, and third derivatives of the proposed solution at the block’s end point . Therefore, we need to consider the evaluation of at the point over to have a two-point implicit block method. By applying the same technique as for the formula for , we have the following formula at :

2.1. Order of the Method

The equations (15)–(22) can be written in a matrix difference equation to obtain the order of the new method as follows:where , and are defined as

It is possible to define the linear operator that is linked to equation (26) as

Expanding equation (28) in the Taylor series yields

The linear operator and the new method have order if , , which is the error constant. Hence, in our method, and . Thus, it can be concluded that the order of the proposed method is 6.

2.2. Zero Stability

The method (18)–(25) is deemed zero-stable in case the roots of the first characteristic polynomial is found to satisfy .

Now, in our proposed method, we will use the following technique to find the matrix of the first characteristic polynomial.

By replacing equation (18) into equation (22), we have

Also, substituting equations (18) and (19) into equation (23), we get

Substituting equations (18)–(21) into equation (24), we have

By substituting equations (18)–(22) into equation (25), we have

By using the equations (18)–(22) and (30)–(33), the general form of the matrix can be written as , identity matrix,

Then, implies that . Hence,

Therefore, it can be inferred that the proposed method has zero stability.

3. Implementation

In this section, we briefly mention the developed method’s implementation. The approximate solutions in equations (15)–(22) will be estimated using the predictor-corrector schemes. The predictor equations using Taylor method are

Define equation (35) as the initial approximation and apply the derived method in equations (18)–(25) as the corrector. In order to obtain the corrector iteration, the following equations will be used

Then, we evaluate the functions and which we use to compute the following approximate solutions at the second point:

Next, we evaluate the functions and that we will use in the next corrector iteration. Hence, we repeat the procedure given in equations (36) and (37) to perform the next corrector iterations until the end of the interval.

4. Results and Discussion

Within this section, we tested the I2PBDO6 method to solve single as well as a system of initial value challenges in the form of equation (1). We will see how the method of two-point block (order 6) produces comparable results to the existing direct methods of similar characteristics. Most of the existing methods are of higher order in comparison to the proposed method. The following notations will be used in the tables:(i)I2PBDO6: the direct implicit two-point block method proposed in this paper, order 6(ii)MCM: direct maximal order multiderivative collocation method [6](iii)ZSB7: the direct zero-stable block method of order seven [7](iv)ZSDM6: the direct zero-stable method of order six [8](v)IHB8: the direct implicit hybrid block method of order eight [12](vi)FDM6: the direct linear multistep method of order six [16](vii)BHC8: the direct block hybrid collocation method of order eight [10](viii)FSHB8: the direct four-step hybrid block method of order eight [11](ix)SSB8: the direct seven-step block method of order eight [9](x)LIBO7: the direct implicit six-step block method of order seven [14](xi)RKTF5: the direct RKTF method of order five with four-stage [13](xii)RKF5: the RK method with six-stage given in Lambert [19](xiii)h: step size(xiv)FCN: number of function evaluations(xv)NS: number of steps(xvi)MAXE: maximum absolute error

4.1. Problems

Problem 1. Exact solution:

Problem 2. Exact solution:

Problem 3. Exact solution:

Problem 4. Exact solution:

Problem 5. Exact solution:

Problem 6. Exact solution:

Problem 7.