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

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.


Introduction
e 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. e 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][9][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. e main aim of this work is to introduce new useful algorithms depending on the DT method to solve systems of linear differential equations.

Analysis of Differential Transform
e basic definition and the fundamental theorems of the DTM and its applicability for various kinds of differential equations are given in [13][14][15][16][17][18][19][20]. For the convenience of the reader, we will present a review of the DTM. To do this, we assume that f(x) ∈ C ∞ (I); then, for any point x 0 ∈ I, the Taylor series of f(x) about x 0 can be represented by Definition 1. Let f (x) be an analytic function about x 0 ; then, the k th order differential transform of f (x) is defined as From equations (1) and (2), we get which 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 f(x) is expressed by a finite series, and equation (3) can be written as where n is decided by the convergence of natural frequency. e 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 v (x), respectively; then, we have the following properties.
Let U (k) and V (k) be the differential transformations of the functions u (x) and v (x), respectively; then, we have for the constants α and β the following properties: e proof of this theorem can be found in [2-5, 8, 10, 12-15, 20].

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. e systems are assumed be autonomous, which means that the independent variable t is not present explicitly. Such a system has the appearance X′ � AX, where A is an n × n matrix given by A � a 11 a 12 · · · a 1n a 21 a 22 · · · a 2n · · · · · · · · · · · · a n1 a n2 · · · a nn Theorem 2. Let X(k) be the k th differential transform of X (t) and A k . be the k th power of the matrix A. en, the solution of the system given by (5) can be expressed as Proof. By using the differential transform property (3) of eorem 1, we get From this equation, we have the following recurrence relation: en, for k � 0, 1, 2, . . . , we get for k � j, we get where A j is the power the matrix A; then, the solution of the system is given by Now, the proof is completed.
□ e principles that have been established can now be applied to the nonhomogeneous problem where W can be a vector function of t.
Theorem 3. Let X(k) and W (k) be the k th differential transforms of X (t) and W (t), respectively, and A k be the k th power of the matrix A. en, the solution of the system given by (13) can be expressed as Proof. By taking the k th differential transform of both sides of (13), we get According to eorem 1 and the operations of differential transform, we have the following recurrence relation: Consequently, . (17) erefore, this leads to the general form of X(k), where en, the general solution of the system is given by Now, the proof is completed.

Applications and Numerical Results
In order to illustrate the advantages and the accuracy of the results given by eorems 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. (20) is system can be written in the matrix form as or Applying the result of eorem 2, we get en, the solutions are given by x 2 (t) � 7, which is the exact solution of the given problem.
Example 2. Consider the system Applying the result of eorem 2, consequently, the solutions of the given system are given as en, the solutions of the given system are (27) e solutions can be written in matrix form as which 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: (29) e exact solutions are obtained by [17] as follows To find the solutions of the given system applying the results in eorem 3, the approximate solutions up to k � 6 are given by (31) Journal of Mathematics erefore, en, the solutions are which is the exact solution of the given problem.
Example 4. Consider the following system: e exact solution of this system is given by [17] as follows: To find the approximate solutions of the given system, applying the results of eorem 3, we have 6 Journal of Mathematics erefore,