A General Scheme for Solving Systems of Linear First-Order Differential Equations Based on the Differential Transform Method

Abstract

In this study, we develop the differential transform method in a new scheme to solve systems of first-order differential equations. The differential transform method is a procedure to obtain the coefficients of the Taylor series of the solution of differential and integral equations. So, one can obtain the Taylor series of the solution of an arbitrary order, and hence, the solution of the given equation can be obtained with required accuracy. Here, we first give some basic definitions and properties of the differential transform method, and then, we prove some theorems for solving the linear systems of first order. Then, these theorems of our system are converted to a system of linear algebraic equations whose unknowns are the coefficients of the Taylor series of the solution. Finally, we give some examples to show the accuracy and efficiency of the presented method.

1. Introduction

The differential transform was first introduced by Zhou [1], and up to now, the DT method has been developed for solving various kinds of differential and integral equations in many literatures. For example, Ali [2] has developed the DT method for solving partial differential equations and Ayaz [3, 4] has applied this method to differential algebraic equations. Arikoglu and Ozkol [5] have solved the integrodifferential equations with boundary value conditions by the DT method. Odibat [6] has used the DT method for solving Volterra integral equations with separable kernels. Tari and Ziyaee [7] have solved the system of two-dimensional nonlinear Volterra integrodifferential equations by the DT method. The systems of integral and integrodifferential equations, the multiorder fractional differential equations, the system of fractional differential equations, the singularly perturbed Volterra integral equations, and the time-fractional diffusion equation have been solved by the DT method in [2, 6, 8–10]. Also, the DT method has been applied to nonlinear parabolic-hyperbolic partial differential equations, and a modified approach of DT has been developed to nonlinear partial differential equations in studies by Biazarand Abdul Halim-Haasan [8]. Patil and Kembayat [11] have solved the two-dimensional Fredholm integral equations. Abdewahid [12] introduced a new basic formula s for the one-dimensional differential transform. The main aim of this work is to introduce new useful algorithms depending on the DT method to solve systems of linear differential equations.

2. Analysis of Differential Transform

The basic definition and the fundamental theorems of the DTM and its applicability for various kinds of differential equations are given in [13–20]. For the convenience of the reader, we will present a review of the DTM. To do this, we assume that ; then, for any point , the Taylor series of about can be represented by

Definition 1. Let f (x) be an analytic function about ; then, the order differential transform of f (x) is defined asFrom equations (1) and (2), we getwhich implies that the concept of differential transform is derived from Taylor series expansion, but the method does not evaluate the derivatives symbolically. However, the relative derivatives are calculated by an iterative way which are described by the transformed equations of the original function. In real applications, the function is expressed by a finite series, and equation (3) can be written aswhere n is decided by the convergence of natural frequency. The fundamental operations performed by differential transform can be readily obtained and are given.

Theorem 1. Let U (k) and V (k) be the differential transformations of the functions u (x) and (x), respectively; then, we have the following properties.

Let U (k) and V (k) be the differential transformations of the functions u (x) and (x), respectively; then, we have for the constants the following properties:(1)If , then (2)If , then (3)If , then (4)If , then (5)If , then (6)If , then (7)If , then (8)If , then (9)If , then (10)If , then (11)If , then (12)If , then (13)If , then (14)If , then

The proof of this theorem can be found in [2–5, 8, 10, 12–15, 20].

3. Description of the Method

In the first part of this section, we introduce a general algorithm depending on the differential transform method to solve systems of n linear differential equations with constant coefficients. The systems are assumed be autonomous, which means that the independent variable t is not present explicitly. Such a system has the appearancewhere A is an

Theorem 2. Let be the differential transform of X (t) and . be the power of the matrix A. Then, the solution of the system given by (5) can be expressed as

Proof. By using the differential transform property (3) of Theorem 1, we getFrom this equation, we have the following recurrence relation:Then, for for , we getwhere is the power the matrix ; then, the solution of the system is given byNow, the proof is completed.

The principles that have been established can now be applied to the nonhomogeneous problemwhere W can be a vector function of t.

Theorem 3. Let and W (k) be the differential transforms of X (t) and W (t), respectively, and be the power of the matrix A. Then, the solution of the system given by (13) can be expressed as

Proof. By taking the differential transform of both sides of (13), we getAccording to Theorem 1 and the operations of differential transform, we have the following recurrence relation:Consequently,Therefore, this leads to the general form of , whereThen, the general solution of the system is given byNow, the proof is completed.

4. Applications and Numerical Results

In order to illustrate the advantages and the accuracy of the results given by Theorems 2 and 3 for solving homogeneous and nonhomogeneous problems, we have applied the method to the following examples.

Example 1. Consider the following system of homogeneous differential equations.This system can be written in the matrix form asorApplying the result of Theorem 2, we getThen, the solutions are given bywhich is the exact solution of the given problem.

Example 2. Consider the systemApplying the result of Theorem 2, consequently, the solutions of the given system are given asThen, the solutions of the given system areThe solutions can be written in matrix form aswhich is the exact solution of the given problem.
Now, we will give some examples of nonhomogeneous initial value problems.

Example 3. Consider the following system:The exact solutions are obtained by [17] as followsTo find the solutions of the given system applying the results in Theorem 3, the approximate solutions up to k = 6 are given byTherefore,Then, the solutions arewhich is the exact solution of the given problem.

Example 4. Consider the following system:The exact solution of this system is given by [17] as follows:To find the approximate solutions of the given system, applying the results of Theorem 3, we haveTherefore,Then, the solution of the given system iswhich are the exact solutions of the given system.

5. Conclusions

In this work, we successfully apply the differential transform method to find useful formulas for solving homogeneous or nonhomogeneous systems of n first-order differential equations. The present methods reduce the computational difficulties of the traditional methods, and all the calculations can be made simple manipulations. Several examples were tested by applying the results of Theorems 2 and 3, and the results have shown a remarkable performance.

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that they have no conflicts of interest.