Next Article in Journal
The Forecasting of a Leading Country’s Government Expenditure Using a Recurrent Neural Network with a Gated Recurrent Unit
Previous Article in Journal
Fine-Grained Forward Secure Firmware Update in Smart Home
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Effective Approach Based on Generalized Bernstein Basis Functions for the System of Fourth-Order Initial Value Problems for an Arbitrary Interval

1
Department of Mathematics, University of Sargodha, Sargodha 40100, Pakistan
2
Software Engineering Programme, Faculty of Computing and Informatics, Universiti Malaysia Sabah, Jalan UMS, Kota Kinabalu 88400, Sabah, Malaysia
3
Data Technologies and Applications (DaTA) Research Lab., Faculty of Computing and Informatics, Universiti Malaysia Sabah, Jalan UMS, Kota Kinabalu 88400, Sabah, Malaysia
4
Creative Advanced Machine Intelligence (CAMI) Research Centre, Universiti Malaysia Sabah, Jalan UMS, Kota Kinabalu 88400, Sabah, Malaysia
5
Department of Mathematics, The Islamia University of Bahawalpur, Bahawalpur 63100, Pakistan
*
Authors to whom correspondence should be addressed.
Mathematics 2023, 11(14), 3076; https://doi.org/10.3390/math11143076
Submission received: 8 June 2023 / Revised: 8 July 2023 / Accepted: 10 July 2023 / Published: 12 July 2023

Abstract

:
The system of ordinary differential equations has many uses in contemporary mathematics and engineering. Finding the numerical solution to a system of ordinary differential equations for any arbitrary interval is very appealing to researchers. The numerical solution of a system of fourth-order ordinary differential equations on any finite interval [ a , b ] is found in this work using a symmetric Bernstein approximation. This technique is based on the operational matrices of Bernstein polynomials for solving the system of fourth-order ODEs. First, using Chebyshev collocation nodes, a generalised approximation of the system of ordinary differential equations is discretized into a system of linear algebraic equations that can be solved using any standard rule, such as Gaussian elimination. We obtain the numerical solution in the form of a polynomial after obtaining the unknowns. The Hyers–Ulam and Hyers–Ulam–Rassias stability analyses are provided to demonstrate that the proposed technique is stable under certain conditions. The results of numerical experiments using the proposed technique are plotted in figures to demonstrate the accuracy of the specified approach. The results show that the suggested Bernstein approximation method for any interval is quick and effective.

1. Introduction and Fundamental Concepts

Ordinary differential equations (ODEs) are the most powerful mathematical tools found in the physical sciences. Many topics in mathematics, physics, and engineering are related to linear and nonlinear ODEs or systems of ODEs. The series of problems in ODEs with initial and boundary conditions arises in experimental physics, numerical simulation approaches, and many other scientific fields, including engineering, signals processing, economics, and acoustics. Many researchers have used various techniques to find approximate solutions to systems of ODEs. For illustration, Mesady et al. [1] presented a Jafari transformation for a system of linear ODEs with medical applications. By using the Jafari transform, ODEs can be transformed into a series of algebraic equations. Higazy and Aggarwal [2] presented the Sawi transformation for systems of ODEs. Sawi transformations have been used to obtain the concentration of chemical reactions for different examples from physical chemistry. In [3], Ricky et al. proposed the neural ordinary differential equations. This technique is based on parameterizing the derivative of the hidden state rather than specifying a deep series of concealed layers. In [4], Brown and Biggs provided methods for solving the ODEs for unconstrained optimization, while Zadunaisky [5] proposed a numerical technique to find the errors in the numerical solutions for the system of ODEs. In order to solve stiff and non-stiff systems of ODEs, in [6] Linda presented automatic selection of methods, which provides a way of determining whether or not a problem can be addressed by a set of strategies suitable for stiff or non-stiff problems. In recent years, many methods based on numerical investigation of nonlinear moving boundary problems, temperature-dependent numerical studies, the solutions of the nonlinear regularized long wave equation, and the nonlinear sinh-Gordon model have been proposed to find numerical solutions of physical models [7,8,9,10,11]. In [12], Neuberger presented the steepest descent in Hilbert spaces for a general system of linear differential equations. To forecast traffic on a limited time scale based on evolutionary algorithms, Chen et al. [13] explained time series forecasting using a system of ODEs. Farshid [14] proposed a differential transform method for systems of ODEs. The main aim of this effort was to catch the exact solution when the solution is known in terms of series expansion. Biazer et al. [15] presented the solution of a system of ODEs using the Adomian decomposition method by converting the given system into a system of first-order ODEs. Kurnaz and Galip [16] presented a comprehensive discussion of solving a system of ODEs by adjusting the step size. Their method allows for control of the truncation error used in numerical methods. In view of the above literature, a very small number of attempts have been made to find the numerical solution of a system of ODEs, and there are no works as yet addressing ways to find the numerical solution of a fourth-order ODE on any finite interval [ a , b ] . In this paper, a new method for finding a numerical solution to a system of fourth-order ODEs on any given interval is developed. It is demonstrated that the suggested method works well and is suitable for solving a linear system of fourth-order ODEs. Additionally, while the Hyers–Ulam stability of ODEs has recently been studied [17,18,19,20], no stability analysis of a system of ODEs for any given arbitrary interval has yet been carried out. In this paper, a new method is presented for solving systems of linear and nonlinear fourth-order ODEs using operational matrices of generalised symmetric Bernstein polynomials. Furthermore, Hyers–Ulam and Hyers–Ulam–Rassias stability analyses are provided for a system of fourth-order ODEs on any finite interval.

2. Preliminaries

The system of fourth-order ODEs has the layout
D x x x x v 1 ( x ) = f ( x , v q , D x v q , D x x v q , D x x x v q ) , D x x x x v 2 ( x ) = f ( x , v q , D x v q , D x x v q , D x x x v q ) , D x x x x v k ( x ) = f ( x , v q , D x v q , D x x v q , D x x x v q ) ,
for a x b and q = 1 , 2 , , k , with the following initial conditions:
v 1 ( a ) = a 1 , v 1 ( a ) = b 1 , v 1 ( a ) = c 1 , v 1 ( a ) = d 1 , v 2 ( a ) = a 2 , v 2 ( a ) = b 2 , v 2 ( a ) = c 2 , v 2 ( a ) = d 2 , v k ( a ) = a k , v k ( a ) = b k , v k ( a ) = c k , v k ( a ) = d k .
where a q , b q , c q , and d q are the given constants for q = 1 , 2 , , k and D x x x x shows the fourth-order derivative. The goal is to use all of these initial conditions to find functions v 1 ( x ) , v 2 ( x ) ,…, v k ( x ) that satisfy the differential equations in (1). The Bernstein polynomial approximation of the function v q ( x ) : [ a , b ] R is defined as
B p ( v q ; x ) = i = 0 p α q , i B i , p ( x ) , q = 1 , 2 , , k ,
and B i , p ( x ) is defined as
B i , p ( v ) = p C i ( v a ) i ( b v ) p i ( b a ) p , i = 0 , , p ,
where p C i = p ! i ! ( p i ) ! , p is the degree, and i is an index of the Bernstein basis polynomial. Each Bernstein basis function satisfies the symmetry condition as follows:
B i , p ( v ) = B i p , i ( a + b v ) .
The Bernstein basis function is non-negative, forms the partition of unity, and can be written in terms of power basis. The following result shows the convergence of Bernstein`s approximation for the system of ODEs.
Theorem 1.
Let B p ( v q ; x ) be the Bernstein polynomial approximation defined in Equation (3) and the sequence { B p ( v q ; x ) } for q = 1 , 2 , , k converge to v q ( x ) for any v q ( x ) C [ a , b ] .
Proof. 
See [21]. □
From Theorem 1, we is apparent that, for uniform continuity, for any v q ( x ) C [ a , b ] and for ϵ > 0 there is an approximate solution B p ( v q ; x ) such that
B p ( v q ; x ) v q ( x ) < ϵ .
In [22], it is proven that the degree of the polynomial p must satisfy the inequality p > Ł δ 2 ϵ ( a + b ) ( 2 a b ) , where Ł = v q .

3. The Numerical Scheme

In this section, the discretization technique for finding the solution of the system of fourth-order ODEs is presented. Here, the unknown function α q , i is replaced by the Bernstein basis function from Equation (3), then the given system of ordinary differential equations is discretized using Chebyshev nodes.

Decomposition of Fourth-Order Ordinary Differential Equations

Linear combination of Bernstein basis functions B p , i ( x ) , i = 0 , 1 , 2 , , p of degree p yields the approximate polynomial v [ A ] ( x ) of the exact function v ( x ) , as follows:
v [ A ] ( x ) = i = 0 p α i B i , p ( x ) = D ψ ( x ) ,
where D = [ α 0 , α 1 , . . . , α p ] and ψ ( x ) = [ B 0 , p , B 1 , p , . . . , B p , p ] T , such as
ψ ( x ) = p 0 ( x a ) 0 ( b x ) p 0 p 1 ( x a ) 1 ( b x ) p 1 p p ( x a ) p ( b x ) 0 .
Now, the vector ψ ( x ) can be expressed as a product of the matrix ( p + 1 ) × ( p + 1 ) and a vector ( p + 1 ) × 1 as
ψ ( x ) = B X ,
where
B = 1 ( b a ) p a 00 a 01 a 02 a 0 p a 10 a 11 a 12 a 1 p a p 0 a p 1 a p 2 a p p , X = 1 x x 2 x p .
Equation (7) implies that
X = B 1 ψ ( x ) ,
which can be solved using any standard rule (LU-method). The method for obtaining an approximate approach to the system of fourth-order ordinary differential equations is now presented. Consider v q [ A ] ( x ) as an approximation of the qth function, provided as
v q [ A ] ( x ) = i = 0 p α q , i B i , p ( x ) = D q ψ ( x ) , q = 1 , 2 , , k ,
where D q = [ α q , 0 , α q , 1 , , α q , p ] . The derivative of ψ ( x ) is denoted by E , which is the ( p + 1 ) × ( p + 1 ) operational matrix, and X = 1 , x , x 2 , , x p T . This can be represented as follows:
D x ψ ( x ) = E ψ ( x ) = B D x X = B D x 1 x x 2 x p
= B 0 1 2 x p x p 1 = B 0 0 0 0 1 0 0 0 0 2 0 0 0 0 p 0 1 x x 2 x p 1 x p = B 0 0 0 0 1 0 0 0 0 2 0 0 0 0 p 0 X .
From (9), we have
D x ψ ( x ) = B 0 0 0 0 1 0 0 0 0 2 0 0 0 0 p 0 B 1 ψ ( x ) .
This implies that
D x = B 0 0 0 0 1 0 0 0 0 2 0 0 0 0 p 0 B 1 , = B σ B 1
where
σ = 0 0 0 0 1 0 0 0 0 2 0 0 0 0 p 0 .
To solve the system of ODEs (1) using operational matrices, v q ( x ) can be approximated by the symmetric Bernstein polynomial as v q ( x ) = D q ψ ( x ) along with
D x v q ( x ) = D q E ψ ( x ) , D x x v q ( x ) = D q ( E ) 2 ψ ( x ) , D x x x v q ( x ) = D q ( E ) 3 ψ ( x ) .
As a result, system (1) can be written as
D 1 ( E ) 4 ψ ( x ) = f ( x , D q ψ ( x ) , D q E ψ ( x ) , D q ( E ) 2 ψ ( x ) , D q ( E ) 3 ψ ( x ) ) , D 2 ( E ) 4 ψ ( x ) = f ( x , D q ψ ( x ) , D q E ψ ( x ) , D q ( E ) 2 ψ ( x ) , D q ( E ) 3 ψ ( x ) ) , D k ( E ) 4 ψ ( x ) = f ( x , D q ψ ( x ) , D q E ψ ( x ) , D q ( E ) 2 ψ ( x ) , D q ( E ) 3 ψ ( x ) ) ,
with the initial conditions
D 1 ψ ( a ) = a 1 , D 1 E ψ ( a ) = b 1 , D 1 ( E ) 2 ψ ( a ) = c 1 , D 1 ( E ) 3 ψ ( a ) = d 1 , D 2 ψ ( a ) = a 2 , D 2 E ψ ( a ) = b 2 , D 2 ( E ) 2 ψ ( a ) = c 2 , D 2 ( E ) 3 ψ ( a ) = d 2 , D k ψ ( a ) = a k , D k E ψ ( a ) = b k , D k ( E ) 2 ψ ( a ) = c k , D k ( E ) 3 ψ ( a ) = d k .
Now, we substitute the collocation points in Equations (12) and (13) to obtain an algebraic equation system that can be solved easily. The collocation points in this case are the Chebyshev nodes, which are provided as
x β = 1 2 ( a + b ) + 1 2 ( b a ) cos ( 2 β 1 ) π 2 p , β = 1 , 2 , , p 3 .
The Chebyshev nodes are very important in approximation, as they form a particularly good set of nodes for polynomial approximation.

4. Hyers–Ulam Stability

The study of stability problems for various functional equations is very important and has attracted the attention of many researchers. The definition of Hyers–Ulam stability has applicable significance, as it means that when studying the Hyers–Ulam stable system one does not have to reach the exact solution, which is often quite difficult or time-consuming. This is useful in many applications, for example fluid dynamics, numerical analysis, optimization, economics, etc. In this section, we examine the Hyers–Ulam stability and Hyers–Ulam–Rassias stability of the system of linear differential equations of order four, for which we consider the qth differential of system (1)
v q ( i v ) ( x ) r = 0 3 a r ( x ) v q ( r ) ( x ) = 0 , q = 1 , 2 , , k ,
with the initial conditions
v q ( a ) = a q , v q ( a ) = b q , v q ( a ) = c q a n d v q ( a ) = d q
for all x I , v q ( x ) C 4 ( I ) , and I = [ a , b ] R .
First, we provide definitions of Hyers–Ulam–Rassias stability and Hyers–Ulam stability for the qth differential (Equation (15)) with initial conditions (Equation (16)).
Definition 1
([23]). The system of fourth-order differential equations (Equation (15)) possesses Hyers–Ulam stability with the starting conditions in Equation (16) if a constant M > 0 exists such that x C 4 ( I ) for each ϵ > 0 if
| v q ( i v ) ( x ) r = 0 3 a r ( x ) v q ( r ) ( x ) | ϵ , x I
with initial conditions
v q ( a ) = a q , v q ( a ) = b q , v q ( a ) = c q a n d v q ( a ) = d q .
Then, for each approximate solution at the pth degree B p ( v q ; x ) C 4 ( I ) satisfying the differential equation (Equation (15)) with Equation (16), we have
| v q ( x ) B p ( v q ; x ) | M ϵ x I ,
where a real number (M) is the Hyers–Ulam stability constant for the differential equation (Equation (15)) with Equation (16).
Definition 2
([23]). The system of the differential equation in Equation (15) possesses Hyers–Ulam–Rassias stability with the starting conditions in Equation (16) if ϕ ( x ) C ( I , R + ) such that x C 4 ( I ) for every ϵ > 0 if
| v q ( i v ) ( x ) r = 0 3 a r ( x ) v q ( r ) ( x ) | ϵ ϕ ( x ) , x I
with initial conditions
v q ( a ) = a q , v q ( a ) = b q , v q ( a ) = c q a n d v q ( a ) = d q .
Then, for each approximate solution at the pth degree B p ( v q ; x ) C 4 ( I ) satisfying the differential equation (Equation (15)) with Equation (16), we have
| v q ( x ) B p ( v q ; x ) | M ϵ ϕ ( x ) x I ,
where a real number (M) is the Hyers–Ulam–Rassias stability constant for the differential equation (Equation (15)) with Equation (16).
The following main results provide the Hyers–Ulam–Rassias stability and Hyers–Ulam stability of the system of differential equations ( Equation (15)) with Equation (16).
Theorem 2.
If r = 0 3 m a x | a r ( x ) | < 4 ! ( b a ) 4 , then the system of differential equations (Equation (1)) has Hyers–Ulam stability with the initial conditions in Equation (2).
Proof. 
For every ϵ > 0 , consider
| v q ( i v ) ( x ) r = 0 3 a r ( x ) v q ( r ) ( x ) | < ϵ
with initial conditions
v q ( a ) = a q , v q ( a ) = b q , v q ( a ) = c q a n d v q ( a ) = d q .
Here, v q ( x ) C 4 ( [ a , b ] ) satisfies | v q ( x ) | < | v q ( x ) | < | v q ( x ) | < | v q ( x ) | ; thus, per Taylor’s formula, we have
v q ( x ) = v q ( a ) + v q ( a ) ( x a ) + v q ( a ) ( x a ) 2 2 ! + v q ( a ) ( x a ) 3 3 ! + v q ( i v ) ( ξ ) ( x a ) 4 4 ! .
Thus,
| v q ( x ) | = | v q ( a ) + v q ( a ) ( x a ) + v q ( a ) ( x a ) 2 2 ! + v q ( a ) ( x a ) 3 3 ! + v q ( i v ) ( ξ ) ( x a ) 4 4 ! | ,
which implies that
max | v q ( x ) | 1 + ( b a ) + ( b a ) 2 2 ! + ( b a ) 3 3 ! max | v q ( x ) | + ( b a ) 4 4 ! max | v q ( i v ) ( x ) | ,
in turn implying that
max | v q ( x ) | = η max | v q ( x ) | + ( b a ) 4 4 ! max | v q ( i v ) ( x ) r = 0 3 a r ( x ) v q ( r ) ( x ) + r = 0 3 a r ( x ) v q ( r ) ( x ) | ,
where η = 1 + ( b a ) + ( b a ) 2 2 ! + ( b a ) 3 3 ! . From Equation (21), we have
max | v q ( x ) | η max | v q ( x ) | + ( b a ) 4 4 ! ϵ + r = 0 3 a r ( x ) max | v q ( x ) | = η + r = 0 3 a r ( x ) ( b a ) 4 4 ! max | v q ( x ) | + ( b a ) 4 4 ! ϵ = δ max | v q ( x ) | + ( b a ) 4 4 ! ϵ
Finally,
max | v q ( x ) | ( b a ) 4 4 ! ( 1 δ ) ϵ .
Choosing M as ( b a ) 4 4 ! ( 1 δ ) , where δ = η + r = 0 3 a r ( x ) ( b a ) 4 4 ! , we have
max | v q ( x ) | M ϵ .
Obviously, if v q [ A ] ( x ) is the approximate solution of Equation (15) with Equation (16), then
| v q ( x ) v q [ A ] ( x ) | M ϵ .
As a result of Definition 1, the system of the fourth-order differential equation in Equation (1) is Hyers–Ulam stable.
Theorem 3.
If r = 0 3 m a x | a r ( x ) | < 4 ! ( b a ) 4 , then the system of the ordinary differential equation in Equation (1) has Hyers–Ulam–Rassias stability with initial conditions Equation (2).
Proof. 
For every ϵ > 0 and φ : [ a , b ] [ 0 , ) , consider
| v q ( i v ) ( x ) r = 0 3 a r ( x ) v q ( r ) ( x ) | < φ ( x ) ϵ
with initial conditions
v q ( a ) = a q , v q ( a ) = b q , v q ( a ) = c q a n d v q ( a ) = d q .
Here, v q ( x ) C 4 ( [ a , b ] ) satisfies | v q ( x ) | < | v q ( x ) | < | v q ( x ) | < | v q ( x ) | ; thus, per Taylor’s formula we have
v q ( x ) = v q ( a ) + v q ( a ) ( x a ) + v q ( a ) ( x a ) 2 2 ! + v q ( a ) ( x a ) 3 3 ! + v q ( i v ) ( ξ ) ( x a ) 4 4 ! ,
which implies that
max | v q ( x ) | 1 + ( b a ) + ( b a ) 2 2 ! + ( b a ) 3 3 ! max | v q ( x ) | + ( b a ) 4 4 ! max | v q ( i v ) ( x ) | .
Thus,
max | v q ( x ) | η max | v q ( x ) | + ( b a ) 4 4 ! φ ( x ) ϵ + r = 0 3 a r ( x ) max | v q ( x ) | .
Finally,
max | v q ( x ) | ( b a ) 4 4 ! ( 1 δ ) φ ( x ) ϵ .
Obviously, if v q [ A ] ( x ) is the approximate solution of Equation (15) with Equation (16), then
| v q ( x ) v q [ A ] ( x ) | M φ ( x ) ϵ ,
where M = ( b a ) 4 4 ! ( 1 δ ) , δ = η + r = 0 3 a r ( x ) ( b a ) 4 4 ! , and η = 1 + ( b a ) + ( b a ) 2 2 ! + ( b a ) 3 3 ! . □
Thus, the system of fourth-order differential equations in Equation (1) exhibits Hyers–Ulam–Rassias stability as a result of Definition 2.

5. Numerical Problems

This portion considers a number of mathematical problems that can be solved utilizing numerical techniques depending on the symmetry of Bernstein’s polynomials. Table 1, Table 2, Table 3, Table 4, Table 5 and Table 6 demonstrate that the presented approach is both effective and accurate. Graphs are shown in Figure 1, Figure 2 and Figure 3 to demonstrate the actual and numerical solutions and compare their accuracy. The following formula is used to calculate the error at the pth degree:
E q , p = | v q ( x ) v q [ A ] ( x ) | ,
where E q , p is the error of the qth function, v q ( x ) is the exact solution, and v q [ A ] ( x ) is the approximate solution. MAPLE and MATLAB computer languages were used on R2018a running on an Intel(R) Core(TM) i7-1165G7 CPU 1.30 GHz processor with 8.00 GB RAM to execute the numerical computations for the solution of numerical problems.
Problem 1.
Consider the following system of fourth-order ordinary differential equations arising in beam-column theory:
e 2 x d 4 v 1 d x 4 + d v 2 d x = 2 x + e x , sin x 2 d 4 v 2 d x 4 + d v 1 d x = e x ,
subject to initial conditions
v 1 ( 0 ) = 1 , v 1 ( 0 ) = 1 , v 1 ( 0 ) = 1 , v 1 ( 0 ) = 1 , v 2 ( 0 ) = 0 , v 2 ( 0 ) = 0 , v 2 ( 0 ) = 2 , v 2 ( 0 ) = 0 ,
and with exact solutions
v 1 ( x ) = e x , v 2 ( x ) = x 2 .
Using the proposed technique described in Section 3, we obtain the following approximate solutions for degree p = 6 :
v 1 , 6 [ A ] ( x ) 0.000752383749 x 6 0.00732584334 x 5 + 0.0407880099 x 4 0.1666 x 3 + 0.5 x 2 x + 1 . v 2 , 6 [ A ] ( x ) 0.000000368277456 x 6 + 0.0000104695635 x 5 + 0.0000153956102 x 4 + x 2 .
Similarly, the approximate solutions for degree p = 10 are as follows:
v 1 , 10 [ A ] ( x ) 0.000000157599924 x 10 0.000000246680218 x 9 + 0.000024375981 x 8 0.00019800899 x 7 + 0.00138863514 x 6 0.00833322734 x 5 + 0.0416666378 x 4 0.1666666667 x 3 + 0.5 x 2 x + 1 . v 2 , 10 [ A ] ( x ) 0.00000000001 x 10 + 0.0000000007 x 9 0.000000001 x 8 + 0.000000002 x 7 0.000000002 x 6 + 0.000000001 x 5 0.00000000001 x 4 + x 2 .
The respective absolute errors of these solutions for degree p = 6 and p = 10 are provided in Table 1. A comparison of these numerical results shows the effectiveness and fast convergence of the proposed technique.
Problem 2.
Consider the system of fourth-order ordinary differential equations:
d 4 v 1 d x 4 + d 2 v 2 d x 2 = x e x 4 e x sin x , d 4 v 2 d x 4 + d 2 v 1 d x 2 = sin x + x e x 2 e x ,
subject to initial conditions
v 1 ( 0 ) = 0 , v 1 ( 0 ) = 1 , v 1 ( 0 ) = 3 , v 1 ( 0 ) = 4 , v 2 ( 0 ) = 0 , v 2 ( 0 ) = 1 , v 2 ( 0 ) = 0 , v 2 ( 0 ) = 1 ,
and with exact solutions
v 1 ( x ) = x e x , v 2 ( x ) = sin x .
Using the proposed technique described in Section 3, we obtain the following approximate solutions for degree p = 6 :
v 1 , 6 [ A ] ( x ) 0.00404440557 x 6 + 0.0350125944 x 5 0.160895381 x 4 + 0.5 x 3 x 2 + x . v 2 , 6 [ A ] ( x ) 0.00080338704 x 6 + 0.00955722545 x 5 0.0011631606 x 4 0.16 x 3 + x .
Similarly, the approximate solutions for degree p = 10 are as follows:
v 1 , 10 [ A ] ( x ) 0.0000014885022 x 10 + 0.0000217261914 x 9 0.000193901265 x 8 + 0.00138461759 x 7 0.00833065013 x 6 + 0.0416655459 x 5 0.166666361 x 4 + 0.50 x 3 x 2 + x . v 2 , 10 [ A ] ( x ) 1.0000000463173 x 10 + 0.000003139404 x 9 5.00000008280 x 8 0.000197850635 x 7 3.00000056860 x 6 + 0.00833348079 x 5 4.0000000006391 x 4 0.1666666667 x 3 + x .
Table 3 shows the respective absolute errors of these solutions for degrees p = 6 and p = 10 . Comparing these numerical results demonstrates the proposed technique’s effectiveness and rapid convergence.
Problem 3.
Consider the following nonlinear system of fourth-order ordinary differential equations:
x 2 d 4 v 1 d x 4 + sin 2 x d 3 v 2 d x 3 + d v 2 d x = 4 x 2 sin x + x 3 cos x + 2 x + 1 . 0 x 1 , cos 1 x d 4 v 2 d x 4 + x 2 d 3 v 2 d x 3 + d v 1 d x = cos x x sin x .
subject to initial conditions
v 1 ( 0 ) = 0 , v 1 ( 0 ) = 1 , v 1 ( 0 ) = 0 , v 1 ( 0 ) = 3 . v 2 ( 0 ) = 0 , v 2 ( 0 ) = 1 , v 2 ( 0 ) = 2 , v 2 ( 0 ) = 0 .
and with exact solutions
v 1 ( x ) = x cos x , v 2 ( x ) = x 2 + x .
Using the proposed technique described in Section 3, we obtain the following approximate solutions for degree p = 6 :
v 1 , 6 [ A ] ( x ) 0.0055108317 x 6 + 0.0509503842 x 5 0.00834103905 x 4 0.5 x 3 + x . v 2 , 6 [ A ] ( x ) 0.0000989355488 x 6 0.000205561024 x 5 + 0.000209351273 x 4 + x 2 + x .
Similarly, the approximate solutions for degree p = 10 are as follows:
v 1 , 10 [ A ] ( x ) 0.00000159494 x 10 + 0.000028967068 x 9 0.0000063157 x 8 0.0013828025 x 7 0.000003863542 x 6 + 0.0416682915 x 5 4.0000000450508 x 4 0.5 x 3 + x . v 2 , 10 [ A ] ( x ) 0.00000000476 x 10 0.000000023321 x 9 + 0.000000049280 x 8 0.000000058259 x 7 + 0.00000004288 x 6 0.000000019687 x 5 + 0.0000000005707 x 4 + x 2 + x .
The absolute error of these solutions for degrees p = 6 and p = 10 is shown in Table 5. Comparison of these numerical results again demonstrates the effectiveness and rapid convergence of the proposed technique.

6. Conclusions

Among the many applications of ODE systems, the simulation of dynamic processes subject to invariants is one that naturally arises in a variety of fields, for example electrical circuits, mechanical systems, and chemical reactions. These processes are frequently modeled by systems made up of differential operations. In this paper, the approximate solution of a system of fourth-order ordinary differential equations on any interval is obtained using the symmetric form of the Bernstein basis functions. The general method involves discretizing a fourth-order system of ODEs into an algebraic system of equations, solving it for unknowns using Gaussian elimination and initial conditions, then generating approximations for the system’s unknown functions. A few examples have been provided to demonstrate the effectiveness of the given technique. It can be concluded that this technique is direct, has low computational cost, and produces the best results even when the degree p is small. With a few modifications, the presented technique can be extended to systems of higher-order linear and nonlinear ordinary and partial differential equations. Furthermore, the proposed method can be used to estimate derivatives in a variety of image processing problems as well as in scattered data approximation and interpolation [24,25]. Another possible extension is to implement the cubic B-spline interpolation for systems of ordinary differential equations and [26] for systems of integral equations.

Author Contributions

Conceptualization, F.K.; Formal analysis, M.B., R.M., S.A.A.K. and F.K.; Funding acquisition, S.A.A.K.; Investigation, M.B. and K.S.; Methodology, M.B., K.S., R.M., S.A.A.K., and F.K.; Project administration, S.A.A.K.; Resources, S.A.A.K.; Supervision, F.K.; Validation, M.B., S.A.A.K. and F.K.; Visualization, F.K.; Writing—original draft, M.B., K.S. and R.M.; Writing—review and editing, S.A.A.K. and F.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data used to support the findings of the study are included within this paper.

Acknowledgments

The fourth author was fully supported by Universiti Malaysia Sabah in Malaysia. The authors are especially thankful to the Faculty of Computing and Informatics at Universiti Malaysia Sabah for providing computing facilities and support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ahmed, I.E.; Yasser, S.H.; Abdullah, M.A. Jafari transformation for solving a system of ordinary differential equations with medical application. Fractal Fract 2021, 5, 130. [Google Scholar]
  2. Higazy, M.; Aggarwal, S. Sawi Transformation for System of Ordinary Differential Equations with Application. Ain Shams Eng. J. 2021, 12, 3173–3182. [Google Scholar] [CrossRef]
  3. Chen, R.; Rubanova, Y.; Bettencourt, J.; David, D. Neural Ordinary Differential Equations; University of Toronto, Vector Institute: Toronto, ON, Canada, 2018. [Google Scholar]
  4. Brown, A.A.; Biggs, M.C. Some effective methods for unconstrained optimization based on the solution of systems of ordinary differential equations. J. Optim. Theory Appl. 1989, 62, 2. [Google Scholar] [CrossRef]
  5. Zadunaisky, P.E. A method for estimation of errors propogated in the numerical solution of a system of ordinary equations. Univ. Buenos Airs 1966, 25, 281–282. [Google Scholar]
  6. Petzold, L. Automatic selection of methods for solving stiff and nonstiff systems of ordinary differential equations. Soc. Ind. Appl. Math. 1983, 4, 196–5204. [Google Scholar] [CrossRef]
  7. Chaurasiya, V.; Jitendra, S. Numerical investigation of a non-linear moving boundary problem describing solidification of a phase change material with temperature dependent thermal conductivity and convection. J. Therm. Stress. 2023, 46, 799–822. [Google Scholar] [CrossRef]
  8. Chaurasiya, V.; Upadhyay, S.; Rai, K.N.; Jitendra, S. A temperature-dependent numerical study of a moving boundary problem with variable thermal conductivity and convection. Waves Random Complex Media 2023, 1–25. [Google Scholar] [CrossRef]
  9. Rashidinia, J.; Rasoulizadeh, M.N. Radial Basis Function Generated Finite Difference Method For the Solution of Sinh-Gordan Equation. J. Appl. Eng. Math 2021, 11, 893–905. [Google Scholar]
  10. Rasoulizadeh, M.N.; Nikan, O.; Avazzadeh, Z. The impact of LRBF-FD on the solutions of the nonlinear regularized long wave equation. Math. Sci. 2021, 15, 365–376. [Google Scholar] [CrossRef]
  11. Nikan, O.; Avazzadeh, Z.; Rasoulizadeh, M.N. Soliton solutions of the nonlinear sine-Gordon model with Neumann boundary conditions arising in crystal dislocation theory. Nonlinear Dyn. 2021, 106, 783–813. [Google Scholar] [CrossRef]
  12. Neuberger, J.W. Steepest Descent for General Systems of Linear Differential Equations in Hilbert Space; North Texas State University Denton: Denton, TX, USA, 2006. [Google Scholar]
  13. Chen, Y.; Yang, B.; Meng, Q.; Zhao, Y.; Abraham, A. Time-series forecasting using a system of ordinary differential equations. Inf. Sci. 2011, 181, 106114. [Google Scholar] [CrossRef]
  14. Mirzaee, F. Differential transform method for solving linear and nonlinear systems of ordinary differential equations. Appl. Math. Sci. 2011, 70, 3465–3472. [Google Scholar]
  15. Biazar, J.; Babolian, E.; Islam, R. Solution of the system of ordinary differential equations by Adomian decomposition method. Appl. Math. Comput. 2004, 147, 713–719. [Google Scholar] [CrossRef]
  16. Kurnaz, A.; Oturanc, G. The differential transform approximation for the system of ordinary differential equations. Int. J. Comput. Math. 2005, 82, 709–719. [Google Scholar] [CrossRef]
  17. Li, Y. Hyers-Ulam stability of linear differential equations y′′=λ2y. Thai J. Math 2010, 8, 215–219. [Google Scholar]
  18. Modebei, M.I.; Otaide, I.; Olaiya, O.O. Generalized Hyers-Ulam stability of second order linear ordinary differential equation with initial condition. Adv. Inequal. Appl. 2014, 36, 2050–7461. [Google Scholar]
  19. Corduneanu, C. Principles of Differential and Integral Equations; Chelsea Publication Company: New York, NY, USA, 1971. [Google Scholar]
  20. Gavruta, P.; Jung, S.; Li, Y. Hyers-Ulam Stability for Second Order Linear Differential Equation with Boundary Conditions. Electron. J. Differ. Equ. 2011, 80, 1–5. [Google Scholar]
  21. Powel, M.J.D. Approximation Theory and Methods; Cambridge University Press: Cambridge, UK, 1981. [Google Scholar]
  22. Basit, M.; Faheem, K. An Effective Approach to Solving the System of Fredholm Integral Equations Based on Bernstein Polynomial on Any Finite Interval. Alex. Eng. J. 2021, 61, 26112623. [Google Scholar] [CrossRef]
  23. Murali, R.; Selvan, A. Hyers-Ulam stability of nth order linear differential equation. Proyecciones 2019, 38, 222–224. [Google Scholar]
  24. Karim, S.A.A. Rational bi-quartic spline with six parameters for surface interpolation with application in image enlargement. IEEE Access 2020, 8, 115621–115633. [Google Scholar] [CrossRef]
  25. Karim, S.A.A.; Saaban, A.; Nguyen, V.T. Scattered data interpolation using quartic triangular patch for shape-preserving interpolation and comparison with mesh-free methods. Symmetry 2020, 12, 1071. [Google Scholar] [CrossRef]
  26. Abdul Karim, S.A.; Khan, F.; Basit, M. Symmetric Bernstein Polynomial Approach for the System of Volterra Integral Equations on Arbitrary Interval and Its Convergence Analysis. Symmetry 2022, 14, 1343. [Google Scholar] [CrossRef]
Figure 1. Comparison of numerical and exact results of v 1 ( x ) in (a) and v 2 ( x ) in (b) for degree p = 10 at distinct node points x [ 0 , 1 ] of Problem 1.
Figure 1. Comparison of numerical and exact results of v 1 ( x ) in (a) and v 2 ( x ) in (b) for degree p = 10 at distinct node points x [ 0 , 1 ] of Problem 1.
Mathematics 11 03076 g001
Figure 2. Comparison of numerical and exact results of v 1 ( x ) in (a) and v 2 ( x ) in (b) for degree p = 10 at distinct node points x [ 0 , 1 ] of Problem 2.
Figure 2. Comparison of numerical and exact results of v 1 ( x ) in (a) and v 2 ( x ) in (b) for degree p = 10 at distinct node points x [ 0 , 1 ] of Problem 2.
Mathematics 11 03076 g002
Figure 3. Comparison of numerical and exact results of v 1 ( x ) in (a) and v 2 ( x ) in (b) for degree p = 10 at distinct node points x [ 0 , 1 ] of Problem 3.
Figure 3. Comparison of numerical and exact results of v 1 ( x ) in (a) and v 2 ( x ) in (b) for degree p = 10 at distinct node points x [ 0 , 1 ] of Problem 3.
Mathematics 11 03076 g003
Table 1. Absolute error of the system in Equation (24) at distinct nodes between 0 and 1 for degree p = 6 and p = 10 of Bernstein basis polynomials.
Table 1. Absolute error of the system in Equation (24) at distinct nodes between 0 and 1 for degree p = 6 and p = 10 of Bernstein basis polynomials.
x E 1 , 6 ( x ) E 1 , 10 ( x ) E 2 , 6 ( x ) E 2 , 10 ( x )
0.0 0.0 0.0 0.0 0.0
0.1 7.84 × 10 8 2.04 × 10 12 1.64 × 10 9 1.53 × 10 16
0.2 1.12 × 10 8 2.42 × 10 11 2.79 × 10 7 1.34 × 10 13
0.3 5.09 × 10 6 9.58 × 10 11 1.49 × 10 7 1.18 × 10 12
0.4 1.44 × 10 6 2.46 × 10 10 4.99 × 10 7 4.96 × 10 12
0.5 3.19 × 10 5 5.04 × 10 10 1.28 × 10 6 1.45 × 10 11
0.6 6.0 × 10 5 9.00 × 10 10 2.79 × 10 6 3.41 × 10 11
0.7 1.01 × 10 5 1.46 × 10 9 5.41 × 10 6 6.96 × 10 11
0.8 1.58 × 10 4 2.22 × 10 9 9.64 × 10 6 1.28 × 10 10
0.9 2.34 × 10 4 3.20 × 10 9 1.60 × 10 5 2.18 × 10 10
1 3.31 × 10 4 4.44 × 10 9 2.54 × 10 5 3.15 × 10 10
Table 2. Relative error of the system in Equation (24) at distinct nodes between 0 and 1 for degree p = 6 and p = 10 of Bernstein basis polynomials.
Table 2. Relative error of the system in Equation (24) at distinct nodes between 0 and 1 for degree p = 6 and p = 10 of Bernstein basis polynomials.
x E 1 , 6 ( x ) E 1 , 10 ( x ) E 2 , 6 ( x ) E 2 , 10 ( x )
0.0 0.0 0.0 0.0 0.0
0.1 8.66 × 10 8 2.04 × 10 10 1.81 × 10 9 1.53 × 10 14
0.2 1.36 × 10 8 6.05 × 10 11 3.40 × 10 7 3.35 × 10 12
0.3 6.87 × 10 6 1.06 × 10 9 1.10 × 10 7 1.31 × 10 11
0.4 2.14 × 10 6 1.53 × 10 9 3.34 × 10 7 3.11 × 10 11
0.5 5.25 × 10 6 2.01 × 10 9 7.76 × 10 7 5.81 × 10 11
0.6 1.09 × 10 4 2.50 × 10 9 1.53 × 10 6 9.47 × 10 11
0.7 2.03 × 10 5 2.97 × 10 9 2.68 × 10 6 1.42 × 10 10
0.8 3.51 × 10 4 3.46 × 10 9 4.33 × 10 6 2.00 × 10 10
0.9 5.75 × 10 4 3.95 × 10 9 6.50 × 10 6 2.69 × 10 10
1 8.99 × 10 4 4.44 × 10 9 9.34 × 10 6 3.15 × 10 10
Table 3. Absolute error of the system in Equation (25) at distinct nodes between 0 and 1 for degree p = 6 and p = 10 of Bernstein basis polynomials.
Table 3. Absolute error of the system in Equation (25) at distinct nodes between 0 and 1 for degree p = 6 and p = 10 of Bernstein basis polynomials.
x E 1 , 6 ( x ) E 1 , 10 ( x ) E 2 , 6 ( x ) E 2 , 10 ( x )
0.0 0.0 0.0 0.0 0.0
0.1 5.14 × 10 7 2.16 × 10 11 1.04 × 10 8 2.89 × 10 12
0.2 7.36 × 10 6 2.57 × 10 10 1.51 × 10 6 3.47 × 10 11
0.3 3.34 × 10 5 1.01 × 10 9 6.98 × 10 5 1.39 × 10 10
0.4 9.50 × 10 5 2.61 × 10 9 2.20 × 10 5 3.66 × 10 10
0.5 2.09 × 10 4 5.35 × 10 9 4.54 × 10 5 7.74 × 10 10
0.6 3.94 × 10 4 9.56 × 10 9 8.75 × 10 5 1.43 × 10 9
0.7 6.68 × 10 4 1.55 × 10 8 1.51 × 10 4 2.42 × 10 9
0.8 1.04 × 10 3 2.36 × 10 8 2.44 × 10 4 3.85 × 10 9
0.9 1.54 × 10 3 3.41 × 10 8 3.73 × 10 4 5.83 × 10 9
1 2.19 × 10 3 4.75 × 10 8 5.46 × 10 4 8.53 × 10 9
Table 4. Relative error of the system in Equation (25) at distinct nodes between 0 and 1 for degree p = 6 and p = 10 of Bernstein basis polynomials.
Table 4. Relative error of the system in Equation (25) at distinct nodes between 0 and 1 for degree p = 6 and p = 10 of Bernstein basis polynomials.
x E 1 , 6 ( x ) E 1 , 10 ( x ) E 2 , 6 ( x ) E 2 , 10 ( x )
0.0 0.0 0.0 0.0 0.0
0.1 4.49 × 10 6 1.56 × 10 10 4.32 × 10 6 9.94 × 10 9
0.2 1.50 × 10 5 4.54 × 10 9 1.33 × 10 4 2.65 × 10 9
0.3 3.54 × 10 4 9.73 × 10 9 3.15 × 10 2 5.24 × 10 8
0.4 3.54 × 10 4 9.73 × 10 9 3.15 × 10 3 5.24 × 10 8
0.5 6.89 × 10 4 1.76 × 10 8 5.20 × 10 3 8.86 × 10 8
0.6 1.19 × 10 3 2.90 × 10 8 8.35 × 10 3 1.36 × 10 7
0.7 1.92 × 10 3 4.45 × 10 8 1.23 × 10 2 1.98 × 10 7
0.8 2.89 × 10 3 6.56 × 10 8 1.74 × 10 2 2.75 × 10 7
0.9 4.20 × 10 3 9.31 × 10 8 2.37 × 10 2 3.71 × 10 7
1 2.19 × 10 3 4.75 × 10 8 3.12 × 10 2 4.88 × 10 7
Table 5. Absolute error of the system in Equation (26) at distinct nodes between 0 and 1 for degree p = 6 and p = 10 of Bernstein basis polynomials.
Table 5. Absolute error of the system in Equation (26) at distinct nodes between 0 and 1 for degree p = 6 and p = 10 of Bernstein basis polynomials.
x E 1 , 6 ( x ) E 1 , 10 ( x ) E 2 , 6 ( x ) E 2 , 10 ( x )
0.0 0.0 0.0 0.0 0.0
0.1 7.46 × 10 7 3.15 × 10 11 1.89 × 10 8 4.11 × 10 13
0.2 1.07 × 10 5 3.75 × 10 10 2.75 × 10 7 4.94 × 10 12
0.3 4.87 × 10 5 1.48 × 10 9 1.26 × 10 6 1.97 × 10 11
0.4 1.38 × 10 4 3.81 × 10 9 3.65 × 10 6 5.14 × 10 11
0.5 3.06 × 10 4 7.81 × 10 9 8.20 × 10 6 1.08 × 10 10
0.6 5.77 × 10 4 1.39 × 10 9 1.57 × 10 5 2.00 × 10 10
0.7 9.77 × 10 4 2.27 × 10 9 2.73 × 10 5 3.43 × 10 10
0.8 1.53 × 10 3 3.45 × 10 8 4.43 × 10 5 5.06 × 10 10
0.9 2.26 × 10 3 4.98 × 10 8 6.85 × 10 5 8.85 × 10 10
1 3.20 × 10 3 6.91 × 10 8 1.02 × 10 4 1.36 × 10 9
Table 6. Relative error of the system in Equation (26) at distinct nodes between 0 and 1 for degree p = 6 and p = 10 of Bernstein basis polynomials.
Table 6. Relative error of the system in Equation (26) at distinct nodes between 0 and 1 for degree p = 6 and p = 10 of Bernstein basis polynomials.
x E 1 , 6 ( x ) E 1 , 10 ( x ) E 2 , 6 ( x ) E 2 , 10 ( x )
0.0 0.0 0.0 0.0 0.0
0.1 7.46 × 10 6 3.15 × 10 10 1.71 × 10 7 3.73 × 10 12
0.2 5.35 × 10 5 1.87 × 10 9 1.14 × 10 6 2.05 × 10 11
0.3 1.62 × 10 4 3.7 × 10 9 3.23 × 10 6 5.05 × 10 11
0.4 3.45 × 10 4 9.52 × 10 9 3.65 × 10 6 9.17 × 10 11
0.5 6.12 × 10 4 1.56 × 10 9 1.09 × 10 5 1.44 × 10 10
0.6 9.61 × 10 4 2.31 × 10 9 1.63 × 10 5 2.08 × 10 10
0.7 1.39 × 10 3 3.24 × 10 9 2.29 × 10 5 2.88 × 10 10
0.8 1.91 × 10 3 4.31 × 10 8 3.07 × 10 5 3.51 × 10 10
0.9 2.51 × 10 3 5.53 × 10 8 4.00 × 10 5 5.17 × 10 10
1 3.20 × 10 3 6.91 × 10 8 5.10 × 10 4 6.8 × 10 10
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Basit, M.; Shahnaz, K.; Malik, R.; Karim, S.A.A.; Khan, F. An Effective Approach Based on Generalized Bernstein Basis Functions for the System of Fourth-Order Initial Value Problems for an Arbitrary Interval. Mathematics 2023, 11, 3076. https://doi.org/10.3390/math11143076

AMA Style

Basit M, Shahnaz K, Malik R, Karim SAA, Khan F. An Effective Approach Based on Generalized Bernstein Basis Functions for the System of Fourth-Order Initial Value Problems for an Arbitrary Interval. Mathematics. 2023; 11(14):3076. https://doi.org/10.3390/math11143076

Chicago/Turabian Style

Basit, Muhammad, Komal Shahnaz, Rida Malik, Samsul Ariffin Abdul Karim, and Faheem Khan. 2023. "An Effective Approach Based on Generalized Bernstein Basis Functions for the System of Fourth-Order Initial Value Problems for an Arbitrary Interval" Mathematics 11, no. 14: 3076. https://doi.org/10.3390/math11143076

APA Style

Basit, M., Shahnaz, K., Malik, R., Karim, S. A. A., & Khan, F. (2023). An Effective Approach Based on Generalized Bernstein Basis Functions for the System of Fourth-Order Initial Value Problems for an Arbitrary Interval. Mathematics, 11(14), 3076. https://doi.org/10.3390/math11143076

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop