On Trigonometric Numerical Integrator for Solving First Order Ordinary Differential Equation ()
1. Introduction
1.1. Formulation of the Interpolating Function
Finite difference schemes have been in the forefront of the methods of using discrete models to approximate the solution of ordinary differential equations. Among the techniques used in building finite difference, scheme is the use of interpolation which requires the design of a basis function that is adequately differentiable in the domain of the numerical integration. Such basis function is then used to create a discrete version of the differential equation involved. This method has been used in the works of [1] - [6], etc. Notable among latest works on interpolating with trigonometric function includes [7], who constructed a new piecewise rational quadratic trigonometric spline with four local positive shape parameters in each subinterval is to visualize some given planar data. The order of approximation of the developed interpolating function was found to
.
[8] developed a new method for smooth rational cubic trigonometric interpolation based on values of function which is being interpolated. This rational cubic trigonometric spline is used to constrain the shape of the interpolant in such a way that the uniqueness of the interpolating function for the given data would be replaced by uniqueness of the interpolating curve for the given data.
[9] constructed a new quadratic trigonometric B-spline with control parameters to address the problems related to two dimensional digital image interpolation. The newly constructed spline is then used to design an image interpolation scheme. Most of these works concentrated on application of trigonometric splines to digital imaging. We are interested in a general use of trigonometric interpolation functions for creating discrete models for the solution of ordinary differential Equation.
1.2. Nonstandard Modeling Techniques
The need for the nonstandard method came up due to some shortcomings of the standard method, in which the qualitative properties of the exact solutions are not usually transferred to the numerical solution. These shortcomings may create a lot of problems, which may affect the stability properties of the standard approach [10].
The concept of numerical instability and its proof by [10] [11], has led to the establishment of five major modeling rules proposed for the construction of difference schemes that will exhibit numerical stability. He also used these techniques to derive numerical models that are exact schemes for some classes of ordinary differential equations.
A finite difference scheme is called nonstandard finite difference method, if at least one of the following conditions is met [12] :
a) In the discrete derivative, the traditional denominator is replaced by a non-negative function
such that,
as
;
b) Non-linear terms that occur in the differential equation are approximated in a non-local way, i.e. by a suitable function of several points of the mesh. The concept of nonstandard finite difference schemes was proposed by [10] as a solution to the numerical instability that exist in the use of finite difference schemes.
Since the discovery of this method, researchers like [5] [12] among others have suggested ways of building numerically reliable schemes using the nonstandard modeling rules. In this work we will create a scheme using the both of these two techniques above.
2. Derivation of the Scheme
Let assume that a solution of a differential equation can be represented by a function
(1)
where
and Q are simulation parameters and A and B are arbitrary constants.
Let a first order ordinary differential equation possess a real valued solution and be differentiable in its domain several times, then from (1) we can write:
(2)
(3)
(4)
(5)
From (3) and (5) we have
(6)
From (2) and (4) we have
(7)
Also from (3) and (5) we have
(8)
For the discrete representation
(9)
(10)
and
(11)
From (7), (8) and (11)
(12)
Or
(13)
Let
,
,
,
Then
(14)
Or
(15)
,
,
,
and
(16)
,
,
,
(17)
,
(18)
Let
,
,
,
,
(19)
Substitute (19) into (17) and (18), we have the integrator in the form (20) and (21) respectively
(20)
(21)
3. Properties of the Integration Method
3.1. Definition [13]
Any algorithm for solving a differential equation in which the approximation
to the solution at
can be calculated iff
and h are known is called a one step method. It is a common practice to write the functional dependence
on the quantities
and h in the form
where
is the incremental function.
3.2. Definition [1]
A numerical scheme with an incremental
is said to be consistent with the initial value problem
if the incremental function is identically zero at
when
.
3.3. Theorem [13]
Let the incremental function of the scheme defined in the one step scheme above be continuous and jointly as a function of its arguments in the region defined by
and
,
where
and let there exists a constant L such that
for all
and
in the region just defined then the relation
is a necessary condition for the convergence of the new scheme.
3.4. Theorem [1]
Let
and
denote two different numerical solution of the differential equation with the initial condition specified a
and
respectively such that
,
.
If the two numerical estimates are generated by the integration scheme, we have
The condition that
is the necessary and sufficient condition for the stability and convergence of the schemes.
3.5. Proof of Convergence of the Integration Method
The increment function
can be written in the form
(22)
Consider Equation (19), we can also write
(23)
Let
be defined as a point in the interior of the interval whose points are y and
, applying mean value theorem, we have
We define
,
,
Therefore
(24)
Taking the absolute value of both sides
(25)
If we let
then our Equation (23) turns to
(26)
which is the condition for convergence.
3.6. Consistence of the Integration Method
Consider an initial value problem of the form
(27)
Having an integrator of the form
which can be obtained using (17) and (18) above and applying the rule in Section 1.2 a).
The renormalized nonstandard form of the schemes will be
,
(28)
,
(29)
where
,
then
If
,
,
,
and
,
and
then (28) and (29) reduced to
(30)
It is a known fact that a consistent method has order of at least one. Therefore, the new numerical integrator is consistent since Equations (28) and (29) can be reduced to (30) when
.
3.7. Stability Analysis of the Integration Method
We shall establish the stability analysis of the integrator by considering the theorem established by [14].
Let
and
denote two different numerical solutions of initial value problem of ordinary differential Equation (25) with the initial conditions specified as
and
respectively, such that
,
. If the two numerical estimates are generated by the integrator (19). From the increment function (26), we have
(31)
(32)
The condition that
(33)
is the necessary and sufficient condition that our new method (19) be stable and convergent.
Proof
From (29) we have
(34)
Then let
(35)
and
(36)
Therefore,
(37)
Applying the mean value theorem as before, we have
(38)
We define
,
,
Therefore
(39)
Taking the absolute value of both sides
(40)
If we let
then our Equation (34) turns to
(41)
which is the condition for convergence.
and
,
, given
, then
(42)
and
, for every
(43)
Then we conclude that our method (29) is stable and convergent.
Note: Similar arguments can be used to proof the stability, convergence and consistency of the scheme
.
4. The Implementation of the Integration Method
4.1. Application of the Finite Difference Schemes to a Deferential Equation I
We derive a scheme for the first equation thus
(44)
(45)
(46)
(47)
Since we are integrating from zero and
where
,
,
,
.
4.2. Application of the Finite Difference Schemes to Logistic Model
We now derive a scheme for the first model
(48)
(49)
(50)
(51)
Since we are integrating from zero we will use
,
where
,
,
,
.
4.3. The Hybrid Nonstandard Schemes
The new scheme will be obtained by substituting the derivatives and applying it to our scheme
(52)
The new standard scheme will be named (New SCH STD) with
,
(53)
For each of the examples the hybrid schemes will be obtained by applying the renormalization techniques.
The hybrid scheme (NEW SCH h) is obtained by using
, step size
(54)
The hybrid scheme (NEW SCH SIN) by using
, step size
,
(55)
The hybrid scheme (NEW SCH EXP) by using
, step size
,
(56)
The choice of this denominator function as step size is informed by the works of [5] [12].
4.4. Experimentation and Result
The following are the 3D graphs obtained from the schemes when applied to the two models (Figures 1-14). We have used same parameters, step size, denominator functions and simulation parameters ∝ and Q to test the two differential equation models.
4.4.1. The Graph for the Method of Differential Equation Model I
Graph of the Schemes with parameters: h = 0.001, ∝ = (4.25 - 14), Q = 36, λ = 0.26.
Figure 1. Solution curves for the schemes of Model Ι.
Figure 2. Graph of absolute error for the schemes in Figure 1.
Graph of the Schemes with parameters: h = 0.0001, ∝ = (4.25 - 14), Q = 36, λ = 0.26.
Figure 3. Solution curves for the schemes of model ΙΙ.
Figure 4. Graph of absolute Error for schemes in Figure 3.
Figure 5. Graph of absolute Error for schemes in Figure 3.
4.4.2. The Graph for the Schemes of the Logistic Model II
Graph of the Schemes with parameters: h = 0.001, ∝ = (6.39 - 7.3), Q = 36, λ = 0.0026.
Figure 6. Solution curves for the schemes of Model IΙ.
Figure 7. Graph of absolute error for the schemes in Figure 6.
Figure 8. Graph of absolute error for the standard scheme in Figure 6.
Figure 9. Error curve for the standard scheme in Figure 6.
Graph of the Schemes of model II with parameters: h = 0.0001, ∝ = (6.39 - 7.3), Q = 36, λ = 0.0016.
Figure 10. Solution curves for the schemes of model ΙΙ.
Figure 11. Solution curves for the standard schemes of model ΙΙ.
Figure 12. Graph of absolute Error for the hybrid Nonstandard schemes in Figure 10.
Figure 13. Graph of absolute Error for the Standard schemes in Figure 11.
Figure 14. Graph of absolute Error for the Standard schemes in Figure 11.
5. Discussion and Conclusion
The derived simulation models have been tested with the control parameters ∝ and Q. We also applied the Nonstandard method by modifying denominator function φ, which also provides for parameters λ and r that can be chosen to obtain iteratively assigned step size as denominator. The discrete model worked for the tested differential equation. The solution curves of the hybrid schemes follow the analytical solutions of the respective equations monotonically are shown in Figure 1, Figure 3, Figure 5, Figure 10 and Figure 12. The numerical properties of the schemes like linear stability, convergence and consistency have been proved analytically. During the course of simulating the equations, we varied these control parameters to obtain family of curves that are very close to the analytic solution. The scheme NEW SCH h and NEW SCH STD are the standard schemes because we maintain the
as step size in both cases. The result showed that the two schemes do not possess smooth monotonicity properties (see Figure 6 and Figure 11) and diverge if
(see Figure 14). These two schemes possess the highest absolute error of deviation from the analytic solution (see Figure 8, Figure 9, Figure 13 and Figure 14). These results confirm the good qualities of Nonstandard modeling techniques [15]. The choice of appropriate values for variables λ and r can be determined using the conditions set by [5] [12]. The graph of Absolute error in Figure 4, Figure 7 and Figure 12 has demonstrated these qualities. The Nonstandard schemes of Model II produced absolute errors that are very close to zero. It is also observed that these two hybrid schemes are exact schemes for the logistic equation when
because the result of the schemes is the same with the analytic solution (see Figure 12). This work can be extended by add higher degree polynomial to the basis function for a possible improved performance [14] [16]. We can conclude that the discrete model is suitable for the solution of first order ordinary differential equation as proposed.