Abstract
In this article, small modification to the Modified Euler Method is proposed. Stability and consistency were tested to determine the end result, and some numerical results were presented, and the CPU time was compared again, and it is recognized that the proposed method is more reliable and compatible with higher efficiency.
1. Introduction
An initial value problem (IVP) is a differential equation with an initial condition that specifies the value of an unknown function at a specific point in the domain. Differential equations are commonly used in physics, chemistry, biology, and economics in science and engineering to solve a variety of physical problems [1, 3].
There are a variety of analytical methods for determining the solution of differential equations. Analytical methods, on the other hand, may be unable to solve such complicated or complex differential equations in some cases. The solution to the difficult differential equations [1–3, 10] is obtained using numerical methods. Numerical methods are extremely useful tools for rapidly solving complex problems when used in conjunction with computer programming.
Many researchers have developed numerical methods for solving ordinary differential equations (ODEs) with initial value problems (IVPs). Many authors have attempted to solve initial value problems (IVP) using a variety of methods, including Euler’s method, Modified Euler Method (MEM), and Improved Modified Euler Method (IMEM), and Improved Euler’s Method. The slope of the function over the interval is estimated by Euler’s procedure using the line tangent to the function at the beginning of the interval. The tangent lines to the solution curve at both ends of the interval are considered in Improved Euler’s (or Heun’s) technique. Some have tried to improve these precision methods, while others have improved them for greater accuracy, stability, and consistency [6–9].
Numerical methods have been improved on occasion to increase performance in accordance with our requirements. By assuming the tangent slope as an average of the arithmetic mean and contra-harmonic mean, [1] proposed to improve the Improved Euler’s (or Heun’s) method. [10] proposed a hybrid numerical method that combines the Modified Euler method, the Improved Euler’s method, and the 2nd-order contra harmonic mean method to solve initial value problems. [2] performed a study on Improving the Improved Modified Euler Method for Better Performance on Autonomous IVP.
The Modified Euler Method and its modification are used in this paper to solve ordinary differential equations in initial value problems. The numerical outcomes are highly promising.
2. The Methodology of the Proposed Method
Consider a first-order ordinary differential equation with initial value problem (IVP)
The simplest and most known numerical method to solve this IVP is Euler’s method which is used to evaluate differential equation involving initial value problem using slope at the start point.
According to [3], the numerical method to solve the IVP using Modified Euler’s Method (MEM) is given by which is used to solve initial value problem using slope at the midpoint.
According to [2], the numerical method to solve the IVP using Improved Modified Euler’s Method (IMEM) is given by which solves ODE involving IVP using slope of the tangent at the midpoint but the value of (4) is updated and hence it better than (3).
Now, the new proposed method assumes: which uses the slope of the tangent at the midpoint but the value of (5) is updated twice as compared to (4).
3. Stability and Consistency Analysis of the Proposed Method
The stability analysis of the proposed method indicated in (5) can be obtained using Dahlquist’s test problem [5]
To check the stability of the proposed method, substituting (6) into (5), we obtain , where
Hence, the stability region of the proposed method is the region indicated in Figure 1 below.
The shaded region in Figure 1 above is the region of stability of the proposed method. Figure 2 refers to the stability regions of the stated methods: Euler method (region enclosed by blue curve) and proposed method (region enclosed by red curve).
In order to check the consistency of the initial value problem (IVP), the numerical formula expressed by [5], Equation (1) can be written in the form of and will be consistent if .
In order to check the consistency of the proposed method of Equation (5), we have and then
which shows that the proposed method is consistent.
4. Examples, Numerical Results, and Discussion
The numerical solutions of for the initial value problem (IVP) using MEM, IMEM, and the proposed method were discussed in this section, and three examples were used to illustrate the new proposed method, and Tables 1–8 show the numerical solutions, absolute errors made, exact solutions, and the average number of function evaluations-CPU time [11], and the results are computed using MATLAB software [4].
Example 1. Solve the initial value problem (IVP) which is reported in [1] using the MEM, IMEM, and the proposed method with step sizes and and compare the errors.
Example 2. Solve the initial value problem (IVP) which is reported in [2].
Example 3. Solve reported in [2].
Solution: For Example 1, Table 1 below shows results using MEM, IMEM, and the new proposed method with step sizes and while Table 2 shows the absolute error made using MEM, IMEM, and the new proposed method to find approximate values of the solution of the initial value problem.
For Example 2, Table 3 below shows results using MEM, IMEM, and the new proposed method with step sizes and while Table 4 shows the absolute error made using MEM, IMEM, and the new proposed method to find approximate values of the solution of the initial value problem.
For Example 3, Table 5 below shows results using MEM, IMEM, and the new proposed method with step sizes and while Table 6 shows the absolute error made using MEM, IMEM, and the new proposed method to find approximate values of the solution of the initial value problem.
Table 7 below shows exact solutions for Examples 1, 2, and 3.
5. Conclusion
A new proposed Modified Euler Method is shown in this paper to solve ordinary differential equations (ODEs) with initial value problems (IVPs). Stability and consistency were evaluated and found to be stable and compatible with the new proposed method. Comparison between MEM, IMEM, and the proposed approach was done. The numerical solutions of the MEM, IMEM, and the new proposed method with step lengths and are shown in Tables 1, 3, and 5. The absolute errors of the numerical solutions obtained in Tables 1, 3, and 5 are shown in Tables 2, 4, and 6, respectively. Table 7 shows the exact solution for each numerical example. Table 8 indicates the CPU time (even though it is machine dependent) and shows that the proposed method takes less time on average to perform the functions evaluation. This shows that the proposed method is the fastest method while IMEM is the slowest method. Hence, the results indicate that the new proposed method better than MEM and IMEM. Future work is necessary to consider the practical applications of this proposed method like solving time fractional diffusion equation and others.
Data Availability
No data were used other than on the table.
Conflicts of Interest
The authors declare that they have no conflicts of interest.