Transform mathematics has traditionally been utilized for obtaining solutions of differential and integro-differential equations (DEs and IDEs) which arise in control theory, engineering, and related areas such as pharmacologic and mathematical modeling   . In general, transform mathematics allows for the conversion of differentiation and integration into algebraic processes which yield a preliminary solution that is expressed within the transform domain. Subsequent inversion of this transform function then produces the actual solution of the original DE or IDE which is represented in terms of moments.
The Laplace transform is classically utilized for this purpose  . Typically, the inversion of a Laplace transform is accomplished using pre-existing tables which facilitate this process. Partial fraction expansion may be helpful, in simplifying an expression within the Laplace domain, but is not always reliable in generating a readily invertible result.
Therefore, when a tabular result is unavailable, inversion of the Laplace transform may be difficult; requiring integration within the complex plane. Approximate Laplace transform inversion techniques exist which are based upon numerical methods , Fourier analysis , or repetitive symbolic differentiation .
Relatively recently, alternative mathematical transforms have been developed. Specifically, the Sumudu transform is defined as :
It should be noted that the Sumudu transform of a real function yields results, within the u-domain, which are also real. Additionally, the Sumudu transform “preserves dimensions.” Thus, whatever physical dimension has will also have   . This is particularly helpful when checking for algebraic accuracy.
For simplicity, the Sumudu transform process will be referred to using the S operator:
Whereas the inversion process will be referred to as:
Note that the traditional Laplace transform is defined using the L operator:
A “duality” between the Laplace and Sumudu transforms therefore exists in terms of similar mathematical properties regarding linearity, convolution, differentiation, and integration. These topics have been explored and discussed previously  . Furthermore, tables of Sumudu transforms are available which allow for straightforward conversions, between the u-domain and t-domains, for many commonly used functions  .
If t in the above equation has units of time, then s consequently has units of complex frequency, whereas u would also have the dimension of time. As previously stated, the Sumudu transform “preserves units.”
The purpose of this paper is to demonstrate an approximate inversion process of the Sumudu transform utilizing a geometric power series technique which uses non-negative integer values for n:
Furthermore, this approximate inversion process may only apply within a narrow numerical range owing to the limitations of a truncated geometric power series as a means of representing a function.
Note that in (6) would also be continuously differentiable or “smooth.” In addition, is a constant coefficient. Moreover, this geometric power series method can also be combined, utilizing superposition, with single or multiple known “pre-existing” Sumudu transform or transforms; such as those available from a table:
Multiple geometric power series may also be combined in an additive or multiplicative manner.
The Sumudu Transform of a Power Series
The Sumudu transform of a power series in the t-domain is a “factorial-based” amplified power series in the u-domain (See Appendix A) . Therefore:
Thus, for integer values of :
A summation of multiple power terms, in the time domain, could then be expressed as a power series:
The corresponding Sumudu transform of the above equation would therefore be:
Use of factorial-based coefficient diminution (FBCD) subsequently yields the inversion of (15) “back to” (14):
Note that t has to be substituted for u on the RHS of the above equation. Therefore, functions which can be approximated with a geometric power series, expressed with the form of (6) in the u-domain, may be readily inverted, back to the t-domain, using the aforementioned technique. Thus, approximate solutions of certain types of both linear differential equations, as well as linear integro-differential equations, can be generated. However, these approximate solutions may also be range-limited.
In addition, a Sumudu transform is frequently a rational function expressed within the u-domain:
However, may sometimes be expressed utilizing transcendental functions. Nonetheless, approximating as a geometric power series in the form of (6) can often be accomplished using commonly-known mathematical techniques and algorithms. Computer-based symbolic processors can also be utilized. As previously stated, geometric power series representations of functions may have limitations which restrict the acceptable accuracy of this approximation to that of a relatively narrow numerical range.
Lastly, a Sumudu transform may sometimes be expressed as a product of two or more rational functions:
Each rational function would then be approximated as a unique geometric power series. Consequently, would be represented as the product of two or more infinite series:
Note that: , , and all represent series-specific constant coefficients. Furthermore, the above product of multiple series can be combined using discrete convolution (the Cauchy product)  . However, each series must be truncated.
As will be shown later, this technique can also be utilized as an approximation-based means to invert a partial fraction expansion from the u-domain back to the t-domain.
Numerical analysis and the conversion of rational functions into geometric power series were accomplished using Mathcad (PTC Corporation, MA, USA). Graphs were prepared using Excel (Microsoft Corporation, WA, USA).
3. Preliminary Examples of the Exponential, Sine, and Cosine Functions
3.1. Exponential Function
The Sumudu transform of the exponential function is :
Therefore, the Sumudu transform of the exponential function subsequently expressed as a power series is:
Use of FBCD with substitution of t for u then yields the well-known Taylor’s series of an exponential function. Thus, the above Sumudu transform is approximately inverted without the need for integration within the complex plane:
3.2. Sine Function
Using (1) the Sumudu transform of the sine function is :
Expressing the above as a power series:
By using FBCD, the above Sumudu transform, expressed as a power series, is approximately inverted. This yields the Taylor series of a sine function. Note that u has been replaced by t:
3.3. Cosine Function
The Sumudu transform of the cosine function is found in a manner similar to that of the sine function :
The Taylor’s series for the cosine function is then approximated using FBCD along with substitution of t for u within the above equation:
4. Differentiation and Integration within the Sumudu Domain
Differentiation and integration, with the Sumudu transform, have an “inverse resemblance” to that of Laplace transforms  . The first derivative is:
Whereas the second derivative is:
Higher-order derivatives can also be determined:
Furthermore, integration within the Sumudu transform domain can be summarized as:
Sumudu Shift Theorem
The Sumudu transform of a function can be expressed using the Sumudu shift theorem and :
Making use of the product rule yields:
5. Approximate Solutions to Certain Linear Differential Equationsand Integro-Differential Equations
5.1. Example 1
Consider the fourth-order linear inhomogeneous differential equation with all initial conditions equal to zero:
Using the Sumudu transform technique:
Which is equivalent to:
Expressing the above equation as a power series yields:
The approximate solution is readily obtained by means of FBCD and substitution of t for u:
The fourth derivative of the above is:
Inspection of (48) and (50) demonstrates the solution:
The numerical results of this are graphically illustrated in Figure 1.
5.2. Example 2
The FBCD technique can also be used to determine the approximate solution to this inhomogeneous first order integro-differential equation (IDE):
Note that a straightforward ordinary differential equation (ODE) results:
Assume that . Solving for :
Expressing the above as a power series:
Representing the above using series notation:
By applying FBCD and substituting t for u results in the inversion of to :
Figure 1. The numerical solution to is displayed.
5.3. Example 3
Consider the following integro-differential equation (IDE) which incorporates the sine integral (Si) function:
where: . Note that that . The Sumudu transform of the above IDE is:
where . This is derived in Appendix B. Rearranging (63):
Figure 2. An approximate solution to: is illustrated.
Figure 3. The above graph displays , , , and from Example 2.
The above can be expressed as:
Note that :
Thus, the product of the above two power series, (68) and (69), yields an approximation for :
Use of the Cauchy product, or discrete convolution with two truncated series, results in a double or “nested” series:
This allows for a single expression for u as a power function and ultimately a power series. Thus, the two series, which have undergone convolution, can also be algebraically combined:
Use of FBCD and substitution of t for u, yields the approximate inversion of :
In addition, a double series which resulted from the discrete convolution of two series, is readily integrated and differentiated utilizing the same “term-by-term” methodology as a single series:
Using a value of results in range-limited approximation of the Si function as shown in Figure 4. Furthermore, and are graphically displayed in Figure 5. It should be noted that the value of was utilized owing to the upper limits of the Mathcad factorial function.
5.4. Example 4
Consider the following third-order differential equation in which all initial conditions are equal to zero:
Figure 4. The sine integral function is approximated by from Example 3.
Figure 5. The above graph illustrates and from Example 3.
It is expressed within the Sumudu transform domain as:
Note that the following relationship: is from a pre-existing Sumudu transform table  .
Solving for :
5.4.1. Solution Method 1: Series Representation of the Partial Fraction Expansion
Use of the partial fraction expansion method results in the following expression:
The above equation is then represented as a sum of five series:
Use of FBCD and substitution of t for u yields the following five power series:
5.4.2. Solution Method 2: Use of Discrete Convolution (The Cauchy Product)
Reiterating the original Sumudu transform equation:
Algebraic rearrangement subsequently yields the product of two rational functions:
Each of these rational functions can be approximated using infinite series representations:
Making use of the distributive property:
Applying discrete convolution (the Cauchy product) with two truncated series, each of m terms, yields:
Utilizing further algebraic simplification:
The use of FBCD and substitution of t for u yields the approximate solution:
“Term-by-term” differentiation results in the following:
The results of this are graphically illustrated in Figure 6 and Figure 7. Note that a value of has been utilized owing to the upper limits of the Mathcad factorial function. Although not illustrated, both method 1 and method 2 yielded numerical results with similar range-limited accuracy.
5.5. Example 5
The following double integral equation can be expressed as an ordinary differential equation (ODE) using a Sumudu transform and the shift theorem:
Separating variables and simplifying:
Figure 6. The function and its first, second, and third derivatives are illustrated from Example 4.
Figure 7. The approximate and exact solutions for Example 4 are illustrated.
Integrating both sides:
This subsequently yields:
Exponentiation results in:
Allowing the constant of integration, C, to equal zero yields:
The above equation can subsequently be expressed as an infinite series:
Distribution of the term yields:
The above Sumudu transform can be inverted by separating it into two distinct functions within the u-domain. Firstly:
where is the Dirac delta function.
Secondly, combining the Sumudu transform along with FBCD for the remaining terms of (116) yields the following:
Additionally, t has replaced u in the above equation as well. Subsequent integration yields:
where is the Heaviside step function which is the indefinite integral of the Dirac delta function. The double integral of (118) is:
Note that the variable t is the indefinite integral of the Heaviside step function. Furthermore, t is also the double indefinite integral of the Dirac delta function. Lastly, the constants of integration for (119) and (120) are equal to zero. Figure 8 and Figure 9 illustrate the solution of this example.
5.6. Example 6
The aforementioned techniques can also be applied to certain partial differential equations (PDEs). As an example:
The above equation can be readily converted to an ODE using the Sumudu transform:
Figure 8. The function , from Example 5, is graphically displayed.
Figure 9. Based on Example 5, and are additive inverses.
An integrating factor is then utilized to take advantage of the product rule:
Expressing the above equation as two separate indefinite integrals:
This yields an expression such that the integrating factor now cancels from both sides of the resulting equation:
The Sumudu transform is then simplified:
Use of FBCD and substitution of t for u yields:
The partial derivatives of (128) demonstrate the solution:
5.7. Example 7
The Sumudu transform, with FBCD, and can also be utilized with this mixed partial integro-differential equation (IDE):
Expressing the above results in an ODE within the Sumudu domain:
Separating variables yields:
Integrating both sides:
This results in an expression which utilizes the natural log function:
Exponentiation then yields :
With the above equation can then be approximated as a power series:
Utilizing FBCD and substitution of t for u results in the following expression:
Note that both infinite series, in (140) and (141), use indices starting at to avoid a negative factorial. Furthermore:
5.8. Example 8
Consider a partial IDE:
Use of a Sumudu transform leads to an ODE:
Integrating both sides:
Allowing C to equal zero:
Note that cannot be readily represented as a power series with the form of (6). However, can be approximated as the product of two power series each of which uses the format of (6):
Figure 10. Examination of the above graph demonstrates how and are additive inverses from Example 7. This graph is based upon and .
Figure 11. Inspection of the above graph, from Example 7, illustrates how and are again additive inverses. In this case, and .
Discrete convolution (the Cauchy product) can then be utilized to combine the product of the above two series into a truncated double series:
Algebraically combining terms yields:
Use of FBCD, with substitution of t for u, produces the solution:
Figure 12. Using Example 8, sum of , , and is equal to zero. Note that: and .
Figure 13. Again, using Example 8, the sum of , , and is also equal to zero. Note that: and .
By taking advantage of the relationship between a geometric power series in the u-domain and its inversion back to the t-domain utilizing FBCD, an approximate range-limited solution to certain differential equations as well as integro-differential equations may be obtained. This process may also be facilitated with the utilization of convolution. Further research and applications of this technique, particularly with non-linear and fractional differential equations, may be warranted.
Appendix A: Determining the Sumudu Transform of tn Using Its Laplace Transform
The Laplace transform of is:
Using (5) the corresponding Sumudu transform is therefore:
Simplifying the above:
For integer values of . Therefore, the use of FBCD with substitution of t for u, leads to inversion of the above Sumudu transform:
Appendix B: The Sumudu Transform of the Sine Integral (Si) Function
Note that the Si function is based upon the following:
Using (31), is approximated as:
Use of “factorial-based” coefficient amplification yields the following Sumudu transform:
Note that the factorial-based coefficient amplification results in a cancellation, in both the numerator and denominator, of the above expression:
It is well established that :
 Atlas, G., Li, K.-J.J., Amin, S. and Hahn, R.G. (2017) Development and Retrospective Clinical Assessment of a Patient-Specific Closed-Form Integro-Differential Equation Model of Plasma Dilution. Biomedical Engineering and Computational Biology, 8, 1-20.
 Atlas, G. (2008) Development and Application of a Logistic-Based Systolic Model for Hemodynamic Measurements Using the Esophageal Doppler Monitor. Cardiovascular Engineering, 8, 159-173.
 Hassanzadeh, H. and Pooladi-Darvish, M. (2007) Comparison of Different Numerical Laplace Inversion Methods for Engineering Applications. Applied Mathematics and Computation, 189, 1966-1981.
 Crump, K.S. (1976) Numerical Inversion of Laplace Transforms Using a Fourier Series Approximation. Journal of the Association for Computing Machinery, 23, 89-96.
 Watugala, G.K. (1993) Sumudu Transform: A New Integral Transform to Solve Differential Equations and Control Engineering Problems. International Journal of Mathematical Education in Science and Technology, 24, 35-43.
 Belgacem, F.B.M., Karaballi, A.A. and Kalla, S.L. (2003) Analytical Investigations of the Sumudu Transform and Applications to Integral Production Equations. Mathematical Problems in Engineering, 3, 103-118.
 Belgacem, F.B.M. and Karaballi, A.A. (2006) Sumudu Transform Fundamental Properties Investigations and Applications. Journal of Applied Mathematics and Stochastic Analysis, 1-23.