A Five-Step P-Stable Method for the Numerical Integration of Third Order Ordinary Differential Equations ()
1. Introduction
Linear multistep methods (LMM) for solving first order initial value problems (ivps) is of the form
(1)
where and are uniquely determined and
Conventionally, they are used to solve higher order ordinary differential equations by first reducing them to a system of first order. This approach has been extensively discussed in [1] -[5] . However the method of reducing to a system first order has some serious drawback which includes wastage of human effort and computer time [6] .
The LMM in (1) generates discrete schemes which are used to solve first order odes. Various forms of this LMM have been developed [1] -[4] . Other researchers have introduced the continuous LMM using the continuous collocation and interpolation technique. This has led to the development of continuous LMM of form
(2)
and are expressed as continuous functions of t and are at least differentiable once.
The introduction of continuous collocation methods as against the discrete schemes enhances better global error estimation and ability to approximate solution at all interior points [6] - [10] . In this study, we shall develop continuous multistep collocation method for the solution of third order ordinary differential equations using power series as the basis function.
Power Series Collocation
In [6] [8] [9] , some continuous LMM of Type (2) were developed using power series of form
: (3)
In [10] Chebyshev polynomial function of the form
(4)
where are some Chebyshev function used to develop continuous LMM.
The use power series as basis function for derivation of continuous LMM are based on the property of analytic function that given the Taylor’s polynomial of the form
(5)
The approximate function reduces to as
In this study we proposed the polynomial function of the form in [7] :
(6)
which is of Type (3) to develop a continuous LMM for the solution of initial value problem of the form:
(7)
This paper is organized as follows: Section 1 consists of introduction and background of study; Section 2, we derive a continuous approximation to for exact solution, and specific methods; section 3 consists of the analysis and implementation followed by numerical examples.
2. Derivation of the Method
Consider the third order differential Equation (7), we proposed an approximate solution of the form:
(8)
where.
The derivative of (8) up to the third order yield
(9)
And where are the parameters to be determined. By substituting (8) and (9) into (7) we have
(10)
Collocating (10) at and interpolating (8) at we obtained the system
of equations given below
(11)
(12)
The above equations are solved to obtain the values of which when substituted into Equation (6) yield a method of the form in Equation (13)
(13)
The continuous polynomial obtained when the values of are substituted into (6) and simplified is as follows
(14)
Evaluating (14) at the following discrete method is obtained
(15)
3. Analysis and Implementation of the Method
3.1. Basic properties of the Method
The method (15) is a specific member of the conventional LMM which can expressed as
(16)
Following [1] [2] , we define the local truncation error associated with (16) by the difference operator
(17)
where is assumed to have continuous derivatives of sufficiently high order. Therefore expanding (23) in Taylor series about the point to obtain the expression
(18)
where the are defined as
, , ,
In the sense of [1] , we say that the method (20) is of order p and error constant if
Using the concept above, the method (19) has order and error constant given by
3.2. Zero-stability of the 5-Step Method
Considering the first characteristics polynomial of the method of Equation (15) given as
Putting implying that is a factor. Therefore solving the polynomial it is found that is also a factor of the polynomial and. The other roots which are called spurious roots are and.
3.3. Region of absolute Stability of the 5-Step scheme
Applying the boundary locus method, we have that
In the spirit of Lambert (1973),
By letting and substituting this into the express above to yield
At and for at an interval of we have that. The method is therefore said to be p-stable.
3.4. Implementattion
Single step method can be used to solve higher order ordinary differential equations directly without the need to first reducing it to an equivalent system of first order.
Consider the initial value problem in (7). For our method of order, Taylor series expansion is used to calculate.
and their first, second, third derivatives up to order.
Then the known values of and are substituted into the differential equations. Next the differential equation is differentiated to obtain the expression for higher derivatives using partial differentiation as follows
where
and
where p is the order of the method.
3.5. Numerical Experiments
Our methods of order were used to solve some initial value problems of both general and special nature using Taylor’s series. Our results were compared with the results of other researchers in this area as seen in table 1. In table 2 and table 3, the accuracy of our method is seen in the small error values.
The following initial value problems were used as our test problems:
3.6. Problem 1
Exact solution:.
3.7. Problem 2
Exact solution:.
3.8. Problem 3
,
Exact solution:.
Table 1. showing the result of test problem 1.
Table 2. Showing the result of test problem 2.
Table 3. showing the result of test problem 3.
4. Discussion of Result
We have developed and implemented our methods using Taylor series of the same order as the schemes that we developed. Some special and general third order initial value problems (ivps) were used to test the efficiency of our methods. Our method was found to be zero stable, consistent and convergent. The better accuracy of our method can be shown from the numerical examples.