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  - , etc. Notable among latest works on interpolating with trigonometric function includes , 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 .
 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.
 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 .
The concept of numerical instability and its proof by  , 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  :
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  as a solution to the numerical instability that exist in the use of finite difference schemes.
Since the discovery of this method, researchers like   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
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:
From (3) and (5) we have
From (2) and (4) we have
Also from (3) and (5) we have
For the discrete representation
From (7), (8) and (11)
Let , , ,
, and (16)
, , (19)
Substitute (19) into (17) and (18), we have the integrator in the form (20) and (21) respectively
3. Properties of the Integration Method
3.1. Definition 
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 
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 
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
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 
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
Consider Equation (19), we can also write
Let be defined as a point in the interior of the interval whose points are y and , applying mean value theorem, we have
Taking the absolute value of both sides
If we let
then our Equation (23) turns to
which is the condition for convergence.
3.6. Consistence of the Integration Method
Consider an initial value problem of the form
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
where , then
If , , , and , and then (28) and (29) reduced to
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 .
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
The condition that
is the necessary and sufficient condition that our new method (19) be stable and convergent.
From (29) we have
Applying the mean value theorem as before, we have
Taking the absolute value of both sides
If we let
then our Equation (34) turns to
which is the condition for convergence.
and , , given , then
, 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
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
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
The new standard scheme will be named (New SCH STD) with
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
The hybrid scheme (NEW SCH SIN) by using , step size
The hybrid scheme (NEW SCH EXP) by using , step size
The choice of this denominator function as step size is informed by the works of  .
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 . The choice of appropriate values for variables λ and r can be determined using the conditions set by  . 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  . We can conclude that the discrete model is suitable for the solution of first order ordinary differential equation as proposed.
 Liu, S.J., Chen, Z.L. and Zhu, Y.P. (2015) Rational Quadratic Trigonometric Interpolation Spline for Data Visualization. Mathematical Problems in Engineering, 2015, Article ID: 983120.
 Rana, S.S., Dube, M. and Tiwari, P. (2013) Rational Trigonometric Interpolation and Constrained Control of the Interpolant Curves. International Journal of Computer Applications, 67, 40-44.
 Anguelov, R. and Lubuma, J.M.S. (2003) Nonstandard Finite Difference Method by Nonlocal Approximation. Mathematics and Computers in Simulation, 61, 465-475.