On Sumudu Transform and System of Differential Equations

and Applied Analysis 3 For each complex number x,ΨP x define a linear mapping of C|N| into C. If anyNj is zero, Ψ i,j P is the empty matrix for all i and the corresponding column of matrices in ΨP is absent. If Nj 0 for all j,ΨP x is defined to be the unique linear mapping of {0} C0 into C; its matrix representation is then the empty matrix. In particular consider ΨP x ( x4 2x2 1 − 2x 3x4 2x 4x3 )


Introduction
A differential equation by itself is inherently underconstrained in the absence of initial values as well as boundary conditions.It is also well known that a differential equation along with the initial values or boundary conditions can be represented by an integral equation by using this integral representation, it becomes possible to solve the problem.However one of the most important achievements, and applications of integral transform methods is solving the partial differential equations PDEs of second order.For this purpose recently a new integral transform, which is called Sumudu transform, was introduced by Watugala 1, 2 and used by Weerakoon 3 for partial derivatives of Sumudu transform, provided the complex inversion formula in order to solve the differential equations in different applications of system engineering, control theory and applied physics.The convolution theorem of Sumudu transform was proved by Asiru in 4 .This new transform was applied to the solution of ordinary differential equations and control engineering problems; see 1, 5 .In 6 , some fundamental properties of the Sumudu transform were established.In 7 , this new transform was applied to the one-dimensional neutron transport equation.In fact the relationship between double Sumudu and double Laplace transforms were studied in 8, 9 .Furthermore in 10 , the Sumudu transform was extended to the distributions and some of their properties were also studied.Thus, there have been several works on Sumudu transform and applied to different kind of problems.
In this paper, we prove Sumudu transform of convolution for the matrices and use to solve the regular system of differential equations.
Throughout the paper we use a square matrix, P P ij of regular system having size n × n of polynomials and the associated determinant, det P .If det P is not the zero polynomial which we write as det P / 0 , we have deg det P ≤ |N P | where N P is the degree of the polynomials in the regular matrix P. The case of equality is so important that we make the following statement.We say that P is regular if det P / 0 and the condition where N j P considered as the highest power of the variable term that occurs in the jth column of matrix P , that is, Next we extend the result given in 4 as follows.For each i and j we define Ψ i,j P x to be the 1 × N j matrix of polynomials given by the matrix product where the a k are the coefficients of P ij and P ij N j k 0 a k /x k .In terms of 7.1 in 4 we have where the number of zero indicated is N j −deg P ij .We define Ψ P to be matrix of polynomials and having size m × |N| defined in terms of the array of matrices: For each complex number x, Ψ P x define a linear mapping of C |N| into C m .If any N j is zero, Ψ i,j P is the empty matrix for all i and the corresponding column of matrices in Ψ P is absent.If N j 0 for all j, Ψ P x is defined to be the unique linear mapping of {0} C 0 into C m ; its matrix representation is then the empty matrix.In particular consider Then we have N 1 4, N 2 3, and N 4, 3 .Thus Ψ P is the 3 × 7 matrix computed as In general, if f f 1 , f 2 , . . ., f p is a sequence of functions on a, b , h h 1 , h 2 , . . ., h p with each h i being an integer ≥ 0 and f i being h i − 1 times differentiable on a, b , we will write, using the notation in 4 , whenever the limits exist.If any h i is zero, the corresponding string is absent.If h i 0 for all i we define Φ f, a; h f, b; h 0 ∈ C 0 .The following proposition was proved in 4 .
Proposition 1.1 Sumudu transform of higher derivatives .Let f be n times differentiable on 0, ∞ and let f t and, for any polynomial P of degree n, S P Ḋ f P u S f u − M P u φ f, n 1.9 for u ∈ dom Sf .In particular (with φ f; n here written as a column vector).For n 2 we have Now we want to extend the above proposition to the system of differential equations as follows.
Proposition 1.2.Let f = f 1 , f 2 , f 3 , . . ., f n be functions and let h = h 1 , h 2 , . . ., h n be a sequence of integers ≥ 0 such that m × n, a matrix of polynomial P, satisfies N P h.Furthermore for each j we let f j be h j times differentiable on 0, ∞ and we let f t 0 for t < 0. Suppose that and, we have, In the following theorem we discuss the Sumudu transform of convolution of matrices.Similarly, in case n 3 it is also true that Assuming that n k is true for the case n k 1, we have the following matrices: having Sumudu transforms given by then we have where

1.23
Thus, The proof completes.
The inverse P −1 u of P u will exist provided that u is not a root of the equation det P u 0; hence let P denote the adjugate matrix of P ; by elementary matrix theory we have Proposition 1.4 Solution of homogeneous equation of regular system .Let P be regular and let y be N P times differentiable on 0, ∞ and zero on −∞, 0 and suppose that then y is given (except at 0) by the formula y S −1 P u −1 Ψ P u φ y, N P .

1.27
Proof.By using 1.12 and assuming that y is Sumudu transformable, we have S P Ḋ y u P u S y u − Ψ P u φ y, h .1.28 Equation 1.26 becomes P u S y u Ψ P u φ y, h , 1.29 using 1.25 S y u P u −1 Ψ P u φ y, h .

1.30
Finally by taking the inverse Sumudu transform of the above equation we have where we assume that the inverse transform exists.
The next proposition was proved in 4 for the single differential equation, and we extend it to the regular system of differential equations.Proposition 1.5.Let P be regular and let b be the greatest of the real parts of the roots of the equation Let f f i be continuous on 0, ∞ , and zero on −∞, 0 , locally integrable, and Sumudu transformable and furthermore suppose that P D y f, 1.33 then we have Note that the most general system of equations can be written in the matrix form as 1.35 which denotes the system n j 1 Here the P ij are polynomials and if we let P P ij , f f ij , respectively, and the vectors y y i , g g i with n components, then the above equation can be written in the form of under the initial condition y 0 y 0 , y 0 y 1 , . . ., y n−1 0 y n−1 .

1.38
Since y k is locally integrable, thus, Sumudu transformable for 0 ≤ k ≤ n and for every such k, then Sumudu transform of 1.37 is given by P u S y u uS f u S g u Ψ P u φ y, N P , 1.39 where P u is a matrix defined by and Ψ P u by 1.41 .In order to find the solution of 1.37 , first of all we multiply 1.39 by the inverse matrix P −1 u , then we get S y u P u −1 uS f u S g u Ψ P u φ y, N P .

1.41
Now, by taking the inverse Sumudu transform for both sides of 1.41 y t S −1 P u −1 uS f u S g u S −1 P −1 u Ψ P u φ y, N P 1.42 provided that the inverse exists for each term in the right-hand side of 1.42 .To illustrate our method, we give the following example.
Example 1.6.Solve for t > 0 the system of two equations

1.43
The matrix and we have det P u 1/u 4 4 which has degree 4 N P .Thus P, regular.Now by applying Sumudu transform to the above system we have where Ψ P u φ y, N P is given by On using 1.25 we obtain Thus based on the above discussions we note that the Sumudu transform can be applied for system of differential equations thus can be used in many engineering problems.