About this Journal Submit a Manuscript Table of Contents
Abstract and Applied Analysis
Volume 2012 (2012), Article ID 270106, 16 pages
doi:10.1155/2012/270106
Research Article

On Sumudu Transform Method in Discrete Fractional Calculus

Department of Mathematics and Computer Science, Cankaya University, Eskisehir Yolu 29 km, 06810 Ankara, Turkey

Received 16 April 2012; Revised 4 June 2012; Accepted 5 June 2012

Academic Editor: Paul Eloe

Copyright © 2012 Fahd Jarad and Kenan Taş. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

In this paper, starting from the definition of the Sumudu transform on a general time scale, we define the generalized discrete Sumudu transform and present some of its basic properties. We obtain the discrete Sumudu transform of Taylor monomials, fractional sums, and fractional differences. We apply this transform to solve some fractional difference initial value problems.

1. Introduction

The fractional calculus, which is as old as the usual calculus, deals with the generalization of the integration and differentiation of integer order to arbitrary order. It has recently received a lot of attention because of its interesting applications in various fields of science, such as, viscoelasticity, diffusion, neurology, control theory, and statistics, see [16].

The analogous theory for discrete fractional calculus was initiated by Miller and Ross [7], where basic approaches, definitions, and properties of the theory of fractional sums and differences were reported. Recently, a series of papers continuing this research has appeared. We refer the reader to the papers [812] and the references cited therein.

In the early 1990's, Watugala [13, 14] introduced the Sumudu transform and applied it to solve ordinary differential equations. The fundamental properties of this transform, which are thought to be an alternative to the Laplace transform were then established in many articles [1519].

The Sumudu transform is defined over the set of functions 𝐴 = 𝑓 ( 𝑡 ) 𝑀 , 𝜏 1 , 𝜏 2 | | | | > 0 , 𝑓 ( 𝑡 ) < 𝑀 𝑒 | 𝑡 | / 𝜏 𝑗 , i f 𝑡 ( 1 ) 𝑗 × [ 0 , ) ( 1 . 1 ) by 1 𝐹 ( 𝑢 ) = 𝕊 { 𝑓 } ( 𝑢 ) = 𝑢 0 𝑓 ( 𝑡 ) 𝑒 ( 𝑡 / 𝑢 ) 𝑑 𝑡 , 𝑢 𝜏 1 , 𝜏 2 . ( 1 . 2 )

Although the Sumudu transform of a function has a deep connection to its Laplace transform, the main advantage of the Sumudu transform is the fact that it may be used to solve problems without resorting to a new frequency domain because it preserves scales and unit properties. By these properties, the Sumudu transform may be used to solve intricate problems in engineering and applied sciences that can hardly be solved when the Laplace transform is used. Moreover, some properties of the Sumudu transform make it more advantageous than the Laplace transform. Some of these properties are(i)The Sumudu transform of a Heaviside step function is also a Heaviside step function in the transformed domain.(ii) 𝕊 { 𝑡 𝑛 } ( 𝑢 ) = 𝑛 ! 𝑢 𝑛 .(iii) l i m 𝑢 𝜏 1 𝐹 ( 𝑢 ) = l i m 𝑡 𝑓 ( 𝑡 ) .(iv) l i m 𝑢 𝜏 2 𝐹 ( 𝑢 ) = l i m 𝑡 𝑓 ( 𝑡 ) .(v) l i m 𝑡 0 𝑓 ( 𝑡 ) = l i m 𝑢 0 𝐹 ( 𝑢 ) .(vi) F o r a n y r e a l o r c o m p l e x n u m b e r 𝑐 , 𝕊 { 𝑓 ( 𝑐 𝑡 ) } ( 𝑢 ) = 𝐹 ( 𝑐 𝑢 ) .

In particular, since constants are fixed by the Sumudu transform, choosing 𝑐 = 0 , it gives 𝐹 ( 0 ) = 𝑓 ( 0 ) .

In dealing with physical applications, this aspect becomes a major advantage, especially in instances where keeping track of units, and dimensional factor groups of constants, is relevant. This means that in problem solving, u and G(u) can be treated as replicas of t and f(t), respectively [20].

Recently, an application of the Sumudu and Double Sumudu transforms to Caputo-fractional differential equations is given in [21]. In [22], the authors applied the Sumudu transform to fractional differential equations.

Starting with a general definition of the Laplace transform on an arbitrary time scale, the concepts of the h-Laplace and consequently the discrete Laplace transform were specified in [23]. The theory of time scales was initiated by Hilger [24]. This theory is a tool that unifies the theories of continuous and discrete time systems. It is a subject of recent studies in many different fields in which dynamic process can be described with discrete or continuous models.

In this paper, starting from the definition of the Sumudu transform on a general time scale, we define the discrete Sumudu transform and present some of its basic properties.

The paper is organized as follows: in Sections 2 and 3, we introduce some basic concepts concerning the calculus of time scales and discrete fractional calculus, respectively. In Section 4, we define the discrete Sumudu transform and present some of its basic properties. Section 5 is devoted to an application.

2. Preliminaries on Time Scales

A time scale 𝕋 is an arbitrary nonempty closed subset of the real numbers . The most well-known examples are 𝕋 = , 𝕋 = , and 𝕋 = 𝑞 = { 𝑞 𝑛 𝑛 } { 0 } , where 𝑞 > 1 . The forward and backward jump operators are defined by 𝜎 ( 𝑡 ) = i n f { 𝑠 𝕋 𝑠 > 𝑡 } , 𝜌 ( 𝑡 ) = s u p { 𝑠 𝕋 𝑠 < 𝑡 } , ( 2 . 1 ) respectively, where i n f = s u p 𝕋 and s u p = i n f 𝕋 . A point 𝑡 𝕋 is said to be left-dense if 𝑡 > i n f 𝕋 and 𝜌 ( 𝑡 ) = 𝑡 , right-dense if 𝑡 < s u p 𝕋 and 𝜎 ( 𝑡 ) = 𝑡 , left-scattered if 𝜌 ( 𝑡 ) < 𝑡 , and right-scattered if 𝜎 ( 𝑡 ) > 𝑡 . The graininess function 𝜇 𝕋 [ 0 , ) is defined by 𝜇 ( 𝑡 ) = 𝜎 ( 𝑡 ) 𝑡 . For details, see the monographs [25, 26].

The following two concepts are introduced in order to describe classes of functions that are integrable.

Definition 2.1 (see [25]). A function 𝑓 𝕋 is called regulated if its right-sided limits exist at all right-dense points in 𝕋 and its left-sided limits exist at all left-dense points in 𝕋 .

Definition 2.2 (see [25]). A function 𝑓 𝕋 is called rd-continuous if it is continuous at right-dense points in 𝕋 and its left-sided limits exist at left-dense points in 𝕋 .

The set 𝕋 𝜅 is derived from the time scale 𝕋 as follows: if 𝕋 has a left-scattered maximum 𝑚 , then 𝕋 𝜅 = 𝕋 { 𝑚 } . Otherwise, 𝕋 𝜅 = 𝕋 .

Definition 2.3 (see [25]). A function 𝑓 𝕋 is said to be delta differentiable at a point 𝑡 𝕋 𝜅 if there exists a number 𝑓 Δ ( 𝑡 ) with the property that given any 𝜀 > 0 , there exists a neighborhood 𝑈 of 𝑡 such that | | [ ] 𝑓 ( 𝜎 ( 𝑡 ) ) 𝑓 ( 𝑠 ) 𝑓 Δ [ ] | | | | | | ( 𝑡 ) 𝜎 ( 𝑡 ) 𝑠 𝜀 𝜎 ( 𝑡 ) 𝑠 𝑠 𝑈 . ( 2 . 2 )

We will also need the following definition in order to define the exponential function on an arbitrary time scale.

Definition 2.4 (see [25]). A function 𝑝 𝕋 is called regressive provided 1 + 𝜇 ( 𝑡 ) 𝑝 ( 𝑡 ) 0 for all 𝑡 𝕋 𝜅 .

The set of all regressive and rd-continuous functions forms an Abelian group under the “circle plus” addition defined by ( 𝑝 𝑞 ) ( 𝑡 ) = 𝑝 ( 𝑡 ) + 𝑞 ( 𝑡 ) + 𝜇 ( 𝑡 ) 𝑝 ( 𝑡 ) 𝑞 ( 𝑡 ) 𝑡 𝕋 𝜅 . ( 2 . 3 )

The additive inverse 𝑝 of 𝑝 is defined by ( 𝑝 ) ( 𝑡 ) = 𝑝 ( 𝑡 ) 1 + 𝜇 ( 𝑡 ) 𝑝 ( 𝑡 ) 𝑡 𝕋 𝜅 . ( 2 . 4 )

Theorem 2.5 (see [25]). Let 𝑝 and 𝑡 0 𝕋 be a fixed point. Then the exponential function 𝑒 𝑝 ( , 𝑡 0 ) is the unique solution of the initial value problem 𝑦 Δ 𝑡 = 𝑝 ( 𝑡 ) 𝑦 , 𝑦 0 = 1 . ( 2 . 5 )

3. An Introduction to Discrete Fractional Calculus

In this section, we introduce some basic definitions and a theorem concerning the discrete fractional calculus.

Throughout, we consider the discrete set 𝑎 = { 𝑎 , 𝑎 + 1 , 𝑎 + 2 , } , w h e r e 𝑎 i s x e d . ( 3 . 1 )

Definition 3.1 (see [27]). Let 𝑓 𝑎 and 𝜈 > 0 be given. Then the 𝜈 t h -order fractional sum of 𝑓 is given by Δ 𝑎 𝜈 1 𝑓 ( 𝑡 ) = Γ ( 𝜈 ) 𝑡 𝜈 𝑠 = 𝑎 ( 𝑡 𝜎 ( 𝑠 ) ) 𝜈 1 𝑓 ( 𝑠 ) f o r 𝑡 𝑎 + 𝜈 . ( 3 . 2 ) Also, we define the trivial sum by Δ 𝑎 0 𝑓 ( 𝑡 ) = 𝑓 ( 𝑡 ) f o r 𝑡 𝑁 𝑎 . ( 3 . 3 ) Note that the fractional sum operator Δ 𝑎 𝜈 maps functions defined on 𝑎 to functions defined on 𝑎 + 𝜈 .

In the above equation the term ( 𝑡 𝜎 ( 𝑠 ) ) 𝜈 1 is the generalized falling function defined by 𝑡 𝜈 = Γ ( 𝑡 + 1 ) Γ ( 𝑡 + 1 𝜈 ) ( 3 . 4 ) for any 𝑡 , 𝜈 for which the right-hand side is well defined. As usual, we use the convention that division by a pole yields zero.

Definition 3.2 (see [27]). Let 𝑓 𝑎 and 𝜈 0 be given, and let 𝑁 be chosen such that 𝑁 1 < 𝜈 𝑁 . Then the 𝜈 t h -order Riemann-Liouville fractional difference of 𝑓 is given by Δ 𝜈 𝑎 𝑓 ( 𝑡 ) = Δ 𝑁 Δ 𝑎 ( 𝑁 𝜈 ) 𝑓 ( 𝑡 ) f o r 𝑡 𝑎 + 𝑁 𝜈 . ( 3 . 5 )

It is clear that, the fractional difference operator Δ 𝜈 𝑎 maps functions defined on 𝑎 to functions defined on 𝑎 + 𝑁 𝜈 .

As stated in the following theorem, the composition of fractional operators behaves well if the inner operator is a fractional difference.

Theorem 3.3 (see [27]). Let 𝑓 𝑎 be given and suppose 𝜈 , 𝜇 > 0 with 𝑁 1 < 𝜈 𝑁 . Then Δ 𝜈 𝑎 + 𝜇 Δ 𝑎 𝜇 𝑓 ( 𝑡 ) = Δ 𝑎 𝜈 𝜇 𝑓 ( 𝑡 ) f o r 𝑡 N 𝑎 + 𝜇 + 𝑁 𝜈 . ( 3 . 6 )

A disadvantage of the Riemann-Liouville fractional difference operator is that when applied to a constant 𝑐 , it does not yield 0. For example, for 0 < 𝑣 < 1 , we have Δ 𝜈 𝑎 𝑐 = 𝑐 ( 𝑡 𝑎 ) 𝜈 . Γ ( 1 𝜈 ) ( 3 . 7 )

In order to overcome this and to make the fractional difference behave like the usual difference, the Caputo fractional difference was introduced in [12].

Definition 3.4 (see [12]). Let 𝑓 𝑎 and 𝜈 0 be given, and let 𝑁 be chosen such that 𝑁 1 < 𝜈 𝑁 . Then the 𝜈 t h -order Caputo fractional difference of 𝑓 is given by 𝐶 Δ 𝜈 𝑎 𝑓 ( 𝑡 ) = Δ 𝑎 ( 𝑁 𝜈 ) Δ 𝑁 𝑓 ( 𝑡 ) f o r 𝑡 𝑎 + 𝑁 𝜈 . ( 3 . 8 )

It is clear that the Caputo fractional difference operator 𝐶 Δ 𝜈 𝑎 maps functions defined on 𝑎 to functions defined on 𝑎 + 𝑁 𝜈 as well. And it follows from the definition of the Caputo fractional difference operator that 𝐶 Δ 𝜈 𝑎 𝑐 = 0 . ( 3 . 9 )

4. The Discrete Sumudu Transform

The following definition is a slight generalization of the one introduced by Jarad et al. [28].

Definition 4.1. The Sumudu transform of a regulated function 𝑓 𝕋 𝑎 is given by 𝕊 𝑎 1 { 𝑓 } ( 𝑢 ) = 𝑢 𝑎 𝑒 ( 1 / 𝑢 ) ( 𝜎 ( 𝑡 ) , 𝑎 ) 𝑓 ( 𝑡 ) Δ 𝑡 𝑢 𝒟 { 𝑓 } , ( 4 . 1 ) where 𝑎 is fixed, 𝕋 𝑎 is an unbounded time scale with infimum 𝑎 and 𝒟 { 𝑓 } is the set of all nonzero complex constants 𝑢 for which 1 / 𝑢 is regressive and the integral converges.

In the special case, when 𝕋 𝑎 = 𝑎 , every function 𝑓 𝑎 is regulated and its discrete Sumudu transform can be written as 𝕊 𝑎 { 1 𝑓 } ( 𝑢 ) = 𝑢 𝑘 = 0 𝑢 𝑢 + 1 𝑘 + 1 𝑓 ( 𝑘 + 𝑎 ) ( 4 . 2 ) for each 𝑢 { 1 , 0 } for which the series converges. For the convergence of the Sumudu transform, we need the following definition.

Definition 4.2 (see [27]). A function 𝑓 𝑎 is of exponential order 𝑟 ( 𝑟 > 0 ) if there exists a constant 𝐴 > 0 such that | | | | 𝑓 ( 𝑡 ) 𝐴 𝑟 𝑡 f o r s u c i e n t l y l a r g e 𝑡 . ( 4 . 3 )

The following lemma can be proved similarly as in Lemma 12 in [27].

Lemma 4.3. Suppose 𝑓 𝑎 is of exponential order 𝑟 > 0 . Then 𝕊 𝑎 | | | { 𝑓 } ( 𝑢 ) e x i s t s f o r a l l 𝑢 { 1 , 0 } s u c h t h a t 𝑢 + 1 𝑢 | | | > 𝑟 . ( 4 . 4 )

The following lemma relates the shifted Sumudu transform to the original.

Lemma 4.4. Let 𝑚 0 and 𝑓 𝑎 𝑚 and 𝑔 𝑎 are of exponential order 𝑟 > 0 . Then for all 𝑢 { 1 , 0 } such that | ( 𝑢 + 1 ) / 𝑢 | > 𝑟 , 𝕊 𝑎 𝑚 𝑢 { 𝑓 } ( 𝑢 ) = 𝑢 + 1 𝑚 𝕊 𝑎 1 { 𝑓 } ( 𝑢 ) + 𝑢 𝑚 1 𝑘 = 0 𝑢 𝑢 + 1 𝑘 + 1 𝕊 𝑓 ( 𝑘 + 𝑎 𝑚 ) , ( 4 . 5 ) 𝑎 + 𝑚 { 𝑔 } ( 𝑢 ) = 𝑢 + 1 𝑢 𝑚 𝕊 𝑎 1 { 𝑔 } ( 𝑢 ) 𝑢 𝑚 1 𝑘 = 0 𝑢 + 1 𝑢 𝑚 1 𝑘 𝑔 ( 𝑘 + 𝑎 ) . ( 4 . 6 )

Proof. For all 𝑢 { 1 , 0 } such that | ( 𝑢 + 1 ) / 𝑢 | > 𝑟 , we have 𝕊 𝑎 𝑚 { 1 𝑓 } ( 𝑢 ) = 𝑢 𝑘 = 0 𝑢 𝑢 + 1 𝑘 + 1 = 1 𝑓 ( 𝑘 + 𝑎 𝑚 ) 𝑢 𝑘 = 𝑚 𝑢 𝑢 + 1 𝑘 + 1 1 𝑓 ( 𝑘 + 𝑎 𝑚 ) + 𝑢 𝑚 1 𝑘 = 0 𝑢 𝑢 + 1 𝑘 + 1 = 1 𝑓 ( 𝑘 + 𝑎 𝑚 ) 𝑢 𝑘 = 0 𝑢 𝑢 + 1 𝑘 + 𝑚 + 1 1 𝑓 ( 𝑘 + 𝑎 ) + 𝑢 𝑚 1 𝑘 = 0 𝑢 𝑢 + 1 𝑘 + 1 = 𝑢 𝑓 ( 𝑘 + 𝑎 𝑚 ) 𝑢 + 1 𝑚 𝕊 𝑎 1 { 𝑓 } ( 𝑢 ) + 𝑢 𝑚 1 𝑘 = 0 𝑢 𝑢 + 1 𝑘 + 1 𝕊 𝑓 ( 𝑘 + 𝑎 𝑚 ) , 𝑎 + 𝑚 1 { 𝑔 } ( 𝑢 ) = 𝑢 𝑘 = 0 𝑢 𝑢 + 1 𝑘 + 1 = 1 𝑔 ( 𝑘 + 𝑎 + 𝑚 ) 𝑢 𝑘 = 𝑚 𝑢 𝑢 + 1 𝑘 𝑚 + 1 = 1 𝑔 ( 𝑘 + 𝑎 ) 𝑢 𝑘 = 0 𝑢 𝑢 + 1 𝑘 𝑚 + 1 1 𝑔 ( 𝑘 + 𝑎 ) 𝑢 𝑚 1 𝑘 = 0 𝑢 𝑢 + 1 𝑘 𝑚 + 1 = 𝑔 ( 𝑘 + 𝑎 ) 𝑢 + 1 𝑢 𝑚 𝕊 𝑎 1 { 𝑔 } ( 𝑢 ) 𝑢 𝑚 1 𝑘 = 0 𝑢 + 1 𝑢 𝑚 1 𝑘 𝑔 ( 𝑘 + 𝑎 ) . ( 4 . 7 )

Taylor monomials are very useful for applying the Sumudu transform in discrete fractional calculus.

Definition 4.5 (see [27]). For each 𝜇 ( ) , define the 𝜇 t h -Taylor monomial to be 𝜇 ( 𝑡 , 𝑎 ) = ( 𝑡 𝑎 ) 𝜇 Γ ( 𝜇 + 1 ) f o r 𝑡 𝑎 . ( 4 . 8 )

Lemma 4.6. Let 𝜇 ( ) and 𝑎 , 𝑏 such that 𝑏 𝑎 = 𝜇 . Then for all 𝑢 { 1 , 0 } such that | ( 𝑢 + 1 ) / 𝑢 | > 1 , one has 𝕊 𝑏 𝜇 ( , 𝑎 ) ( 𝑢 ) = ( 𝑢 + 1 ) 𝜇 . ( 4 . 9 )

Proof. By the general binomial formula ( 𝑥 + 𝑦 ) 𝜈 = 𝑘 = 0 𝑣 𝑘 𝑥 𝑘 𝑦 𝑣 𝑘 ( 4 . 1 0 ) for 𝜈 , 𝑥 , 𝑦 such that | 𝑥 | < | 𝑦 | , where 𝑣 𝑘 𝜈 = 𝑘 𝑘 ! , ( 4 . 1 1 ) as in [27], it follows from (4.10) and 𝑘 𝑣 = ( 1 ) 𝑘 𝑘 + 𝑣 1 𝑣 1 , ( 4 . 1 2 ) where 𝑘 0 that 1 ( 1 𝑦 ) 𝜈 = ( ( 𝑦 ) + 1 ) 𝜈 = 𝑘 = 0 𝑦 𝑘 + 𝑣 1 𝑣 1 𝑘 ( 4 . 1 3 ) for 𝜈 and | 𝑦 | < 1 .
And since 𝑏 𝑎 = 𝜇 , we have for all 𝑢 { 1 , 0 } such that | ( 𝑢 + 1 ) / 𝑢 | > 1 , ( 𝑢 + 1 ) 𝜇 = 1 1 𝑢 + 1 ( 1 ( 𝑢 / ( 𝑢 + 1 ) ) ) 𝜇 + 1 = 1 𝑢 + 1 𝑘 = 0 𝜇 𝑢 𝑘 + 𝜇 𝑢 + 1 𝑘 = 1 𝑢 𝑘 = 0 𝜇 𝑢 𝑘 + 𝜇 𝑢 + 1 𝑘 + 1 = 1 𝑢 𝑘 = 0 ( 𝑘 + 𝜇 ) 𝜇 𝑢 Γ ( 𝜇 + 1 ) 𝑢 + 1 𝑘 + 1 = 1 𝑢 𝑘 = 0 𝜇 𝑢 ( 𝑘 + 𝑏 , 𝑎 ) 𝑢 + 1 𝑘 + 1 = 𝕊 𝑏 𝜇 ( , 𝑎 ) ( 𝑢 ) . ( 4 . 1 4 )

Definition 4.7 (see [27]). Define the convolution of two functions 𝑓 , 𝑔 𝑎 by ( 𝑓 𝑔 ) ( 𝑡 ) = 𝑡 𝑟 = 𝑎 𝑓 ( 𝑟 ) 𝑔 ( 𝑡 𝑟 + 𝑎 ) f o r 𝑡 𝑎 . ( 4 . 1 5 )

Lemma 4.8. Let 𝑓 , 𝑔 𝑎 be of exponential order 𝑟 > 0 . Then for all 𝑢 { 1 , 0 } such that | ( 𝑢 + 1 ) / 𝑢 | > 𝑟 , 𝕊 𝑎 { 𝑓 𝑔 } ( 𝑢 ) = ( 𝑢 + 1 ) 𝕊 𝑎 { 𝑓 } ( 𝑢 ) 𝕊 𝑎 { 𝑔 } ( 𝑢 ) . ( 4 . 1 6 )

Proof. Since 𝕊 𝑎 { 1 𝑓 𝑔 } ( 𝑢 ) = 𝑢 𝑘 = 0 𝑢 𝑢 + 1 𝑘 + 1 ( = 1 𝑓 𝑔 ) ( 𝑘 + 𝑎 ) 𝑢 𝑘 = 0 𝑢 𝑢 + 1 𝑘 + 1 𝑘 + 𝑎 𝑟 = 𝑎 = 1 𝑓 ( 𝑟 ) 𝑔 ( ( 𝑘 + 𝑎 ) 𝑟 + 𝑎 ) 𝑢 𝑘 𝑘 = 0 𝑟 = 0 𝑢 𝑢 + 1 𝑘 + 1 𝑓 ( 𝑟 + 𝑎 ) 𝑔 ( 𝑘 𝑟 + 𝑎 ) , ( 4 . 1 7 ) the substitution 𝜏 = 𝑘 𝑟 yields 𝕊 𝑎 { 1 𝑓 𝑔 } ( 𝑢 ) = 𝑢 𝜏 = 0 𝑟 = 0 𝑢 𝑢 + 1 𝜏 + 𝑟 + 1 1 𝑓 ( 𝑟 + 𝑎 ) 𝑔 ( 𝜏 + 𝑎 ) = ( 𝑢 + 1 ) 𝑢 𝑟 = 0 𝑢 𝑢 + 1 𝑟 + 1 1 𝑓 ( 𝑟 + 𝑎 ) 𝑢 𝜏 = 0 𝑢 𝑢 + 1 𝜏 + 1 𝑔 ( 𝜏 + 𝑎 ) = ( 𝑢 + 1 ) 𝕊 𝑎 { 𝑓 } ( 𝑢 ) 𝕊 𝑎 { 𝑔 } ( 𝑢 ) ( 4 . 1 8 ) for all 𝑢 { 1 , 0 } such that | ( 𝑢 + 1 ) / 𝑢 | > 𝑟 .

Theorem 4.9. Suppose 𝑓 𝑎 is of exponential order 𝑟 1 and let 𝜈 > 0 with 𝑁 1 < 𝜈 𝑁 . Then for all 𝑢 { 1 , 0 } such that | ( 𝑢 + 1 ) / 𝑢 | > 𝑟 , 𝕊 𝑎 + 𝜈 Δ 𝑎 𝜈 𝑓 ( 𝑢 ) = ( 𝑢 + 1 ) 𝜈 𝕊 𝑎 𝕊 { 𝑓 } ( 𝑢 ) , ( 4 . 1 9 ) 𝑎 + 𝜈 𝑁 Δ 𝑎 𝜈 𝑓 𝑢 ( 𝑢 ) = 𝑁 ( 𝑢 + 1 ) 𝑁 𝜈 𝕊 𝑎 { 𝑓 } ( 𝑢 ) . ( 4 . 2 0 )

Proof. First note that the shift formula (4.5) implies that for all 𝑢 { 1 , 0 } such that | ( 𝑢 + 1 ) / 𝑢 | > 𝑟 , 𝕊 𝑎 + 𝜈 𝑁 Δ 𝑎 𝜈 𝑓 ( 1 𝑢 ) = 𝑢 𝑘 = 0 𝑢 𝑢 + 1 𝑘 + 1 Δ 𝑎 𝜈 = 𝑢 𝑓 ( 𝑘 + 𝑎 + 𝜈 𝑁 ) 𝑢 + 1 𝑁 𝕊 𝑎 + 𝜈 Δ 𝑎 𝜈 𝑓 1 ( 𝑢 ) + 𝑢 𝑁 1 𝑘 = 0 𝑢 𝑢 + 1 𝑘 + 1 Δ 𝑎 𝜈 = 𝑢 𝑓 ( 𝑘 + 𝑎 + 𝜈 𝑁 ) 𝑢 + 1 𝑁 𝕊 𝑎 + 𝜈 Δ 𝑎 𝜈 𝑓 ( 𝑢 ) , ( 4 . 2 1 ) taking 𝑁 zeros of Δ 𝑎 𝜈 𝑓 into account. Furthermore, by (4.9), (4.15), and (4.16), 𝕊 𝑎 + 𝜈 Δ 𝑎 𝜈 𝑓 ( 1 𝑢 ) = 𝑢 𝑘 = 0 𝑢 𝑢 + 1 𝑘 + 1 Δ 𝑎 𝜈 = 1 𝑓 ( 𝑘 + 𝑎 + 𝜈 ) 𝑢 𝑘 = 0 𝑢 𝑢 + 1 𝑘 + 1 𝑘 + 𝑎 𝑟 = 𝑎 ( 𝑘 + 𝑎 + 𝜈 𝜎 ( 𝑟 ) ) 𝜈 1 = 1 Γ ( 𝜈 ) 𝑓 ( 𝑟 ) 𝑢 𝑘 = 0 𝑢 𝑢 + 1 𝑘 + 1 𝑘 + 𝑎 𝑟 = 𝑎 𝑓 ( 𝑟 ) 𝜈 1 = 1 ( ( 𝑘 + 𝑎 ) 𝑟 + 𝑎 , 𝑎 ( 𝜈 1 ) ) 𝑢 𝑘 = 0 𝑢 𝑢 + 1 𝑘 + 1 𝑓 𝜈 1 ( , 𝑎 ( 𝜈 1 ) ) ( 𝑘 + 𝑎 ) = 𝕊 𝑎 𝑓 𝜈 1 ( , 𝑎 ( 𝜈 1 ) ) ( 𝑢 ) = ( 𝑢 + 1 ) 𝕊 𝑎 { 𝑓 } ( 𝑢 ) 𝕊 𝑎 𝜈 1 ( , 𝑎 ( 𝜈 1 ) ) = ( 𝑢 + 1 ) ( 𝑢 + 1 ) 𝜈 1 𝕊 𝑎 { 𝑓 } ( 𝑢 ) = ( 𝑢 + 1 ) 𝜈 𝕊 𝑎 { 𝑓 } ( 𝑢 ) . ( 4 . 2 2 ) Then we obtain 𝕊 𝑎 + 𝜈 𝑁 Δ 𝑎 𝜈 𝑓 𝑢 ( 𝑢 ) = 𝑢 + 1 𝑁 𝕊 𝑎 + 𝜈 Δ 𝑎 𝜈 𝑓 = 𝑢 ( 𝑢 ) 𝑁 ( 𝑢 + 1 ) 𝑁 𝜈 𝕊 𝑎 { 𝑓 } ( 𝑢 ) . ( 4 . 2 3 )

Theorem 4.10. Suppose 𝑓 𝑎 is of exponential order 𝑟 1 and let 𝜈 > 0 with 𝑁 1 < 𝜈 𝑁 . Then for all 𝑢 { 1 , 0 } such that | ( 𝑢 + 1 ) / 𝑢 | > 𝑟 , 𝕊 𝑎 + 𝑁 𝜈 Δ 𝜈 𝑎 𝑓 ( 𝑢 ) = ( 𝑢 + 1 ) 𝑁 𝜈 𝑢 𝑁 𝕊 𝑎 { 𝑓 } ( 𝑢 ) 𝑁 1 𝑘 = 0 𝑢 𝑘 𝑁 Δ 𝑎 𝜈 𝑁 + 𝑘 𝑓 ( 𝑎 + 𝑁 𝜈 ) . ( 4 . 2 4 )

Proof. Let 𝑓 , 𝑟 , 𝜈 , and 𝑁 be as in the statement of the theorem. We already know from Theorem 3.8 in [28] that (4.24) holds when 𝜈 = 𝑁 , that is, 𝕊 𝑎 Δ 𝑁 𝑓 1 ( 𝑢 ) = 𝑢 𝑁 𝕊 𝑎 { 𝑓 } ( 𝑢 ) 𝑁 1 𝑘 = 0 𝑢 𝑘 𝑁 Δ 𝑘 𝑓 ( 𝑎 ) . ( 4 . 2 5 ) If 𝑁 1 < 𝜈 < 𝑁 , then 0 < 𝑁 𝜈 < 1 and hence it follows from (3.6), (4.19), and (4.25) that 𝕊 𝑎 + 𝑁 𝜈 Δ 𝜈 𝑎 𝑓 ( 𝑢 ) = 𝕊 𝑎 + 𝑁 𝜈 Δ 𝑁 Δ 𝑎 ( 𝑁 𝜈 ) 𝑓 = 1 ( 𝑢 ) 𝑢 𝑁 𝕊 𝑎 + 𝑁 𝜈 Δ 𝑎 ( 𝑁 𝜈 ) 𝑓 ( 𝑢 ) 𝑁 1 𝑘 = 0 𝑢 𝑘 𝑁 Δ 𝑘 Δ 𝑎 ( 𝑁 𝜈 ) = 𝑓 ( 𝑎 + 𝑁 𝜈 ) ( 𝑢 + 1 ) 𝑁 𝜈 𝑢 𝑁 𝕊 𝑎 { 𝑓 } ( 𝑢 ) 𝑁 1 𝑘 = 0 𝑢 𝑘 𝑁 Δ 𝑎 𝜈 𝑁 + 𝑘 𝑓 ( 𝑎 + 𝑁 𝜈 ) . ( 4 . 2 6 )

In the following theorem the Sumudu transform of the Caputo fractional difference operator is presented.

Theorem 4.11. Suppose 𝑓 𝑎 is of exponential order 𝑟 1 and let 𝜈 > 0 with 𝑁 1 < 𝜈 𝑁 . Then for all 𝑢 { 1 , 0 } such that | ( 𝑢 + 1 ) / 𝑢 | > 𝑟 , 𝕊 𝑎 + 𝑁 𝜈 𝐶 Δ 𝜈 𝑎 𝑓 ( 𝑢 ) = ( 𝑢 + 1 ) 𝑁 𝜈 𝑢 𝑁 𝕊 𝑎 { 𝑓 } ( 𝑢 ) 𝑁 1 𝑘 = 0 𝑢 𝑘 Δ 𝑘 . 𝑓 ( 𝑎 ) ( 4 . 2 7 )

Proof. Let 𝑓 , 𝑟 , 𝜈 , and 𝑁 be as in the statement of the theorem. We already know from (4.25) that 𝑣 = 𝑁 , (4.27) holds. If 𝑁 1 < 𝜈 < 𝑁 , then 0 < 𝑁 𝜈 < 1 and hence it follows from (4.19) and (4.25) that 𝕊 𝑎 + 𝑁 𝜈 𝐶 Δ 𝜈 𝑎 𝑓 ( 𝑢 ) = 𝕊 𝑎 + 𝑁 𝜈 Δ 𝑎 ( 𝑁 𝜈 ) Δ 𝑁 𝑓 = ( 𝑢 ) ( 𝑢 + 1 ) 𝑁 𝜈 𝕊 𝑎 Δ 𝑁 𝑓 = ( 𝑢 ) ( 𝑢 + 1 ) 𝑁 𝜈 𝑢 𝑁 𝕊 𝑎 { 𝑓 } ( 𝑢 ) 𝑁 1 𝑘 = 0 𝑢 𝑘 Δ 𝑘 . 𝑓 ( 𝑎 ) ( 4 . 2 8 )

Lemma 4.12. Let 𝑓 𝑎 be given. For any 𝑝 0 and 𝜈 > 0 with 𝑁 1 < 𝜈 𝑁 , one has 𝐶 Δ 𝑎 𝜈 + 𝑝 𝑓 ( 𝑡 ) = 𝐶 Δ 𝜈 𝑎 Δ 𝑝 𝑓 ( 𝑡 ) 𝑓 𝑜 𝑟 𝑡 𝑎 + 𝑁 𝑣 . ( 4 . 2 9 )

Proof. Let f, v, N, and p be given as in the statement of the lemma. Then 𝐶 Δ 𝑎 𝜈 + 𝑝 𝑓 ( 𝑡 ) = Δ 𝑎 ( 𝑁 + 𝑝 𝜈 𝑝 ) Δ 𝑁 + 𝑝 𝑓 ( 𝑡 ) = Δ 𝑎 ( 𝑁 𝜈 ) Δ 𝑁 Δ 𝑝 = 𝑓 ( 𝑡 ) 𝐶 Δ 𝜈 𝑎 Δ 𝑝 𝑓 ( 𝑡 ) . ( 4 . 3 0 )

Corollary 4.13. Suppose 𝑓 𝑎 is of exponential order 𝑟 1 , 𝜈 > 0 with 𝑁 1 < 𝜈 𝑁 and 𝑝 0 . Then for all 𝑢 { 1 , 0 } such that | ( 𝑢 + 1 ) / 𝑢 | > 𝑟 , 𝕊 𝑎 + 𝑁 𝜈 𝐶 Δ 𝑎 𝜈 + 𝑝 𝑓 ( 𝑢 ) = ( 𝑢 + 1 ) 𝑁 𝜈 𝑢 𝑁 + 𝑝 𝕊 𝑎 { 𝑓 } ( 𝑢 ) 𝑁 + 𝑝 1 𝑘 = 0 𝑢 𝑘 Δ 𝑘 . 𝑓 ( 𝑎 ) ( 4 . 3 1 )

Proof. The proof follows from (4.25), (4.27), and (4.29).

5. Applications

In this section, we will illustrate the possible use of the discrete Sumudu transform by applying it to solve some initial value problems. The following initial value problem was solved in Theorem 23 in [27] by using the Laplace transforms.

Example 5.1. Suppose 𝑓 𝑎 is of exponential order 𝑟 1 and let 𝜈 > 0 with 𝑁 1 < 𝜈 𝑁 . The unique solution to the fractional initial value problem Δ 𝜈 𝑎 + 𝜈 𝑁 𝑦 ( 𝑡 ) = 𝑓 ( 𝑡 ) , 𝑡 𝑎 Δ 𝑘 𝑦 ( 𝑎 + 𝜈 𝑁 ) = 𝐴 𝑘 , 𝑘 { 0 , 1 , , 𝑁 1 } , 𝐴 𝑘 ( 5 . 1 ) is given by 𝑦 ( 𝑡 ) = 𝑁 1 𝑘 = 0 𝛼 𝑘 ( 𝑡 𝑎 ) 𝜈 + 𝑘 𝑁 + Δ 𝑎 𝜈 𝑓 ( 𝑡 ) , 𝑡 𝑎 + 𝜈 𝑁 , ( 5 . 2 ) where 𝛼 𝑘 = Δ 𝜈 𝑁 + 𝑘 𝑎 + 𝜈 𝑁 𝑦 ( 𝑎 ) = Γ ( 𝜈 + 𝑘 𝑁 + 1 ) 𝑘 𝑝 = 0 𝑘 𝑝 𝑗 = 0 ( 1 ) 𝑗 𝑘 ! ( 𝑘 𝑗 ) 𝑁 𝜈 𝑘 𝑝 𝑗 𝐴 𝑘 𝑝 𝑝 ( 5 . 3 ) for 𝑘 { 0 , 1 , , 𝑁 1 } .

Proof. Since 𝑓 is of exponential order 𝑟 , then 𝕊 𝑎 { 𝑓 } ( 𝑢 ) exists for all 𝑢 { 1 , 0 } such that | ( 𝑢 + 1 ) / 𝑢 | > 𝑟 . So, applying the Sumudu transform to both sides of the fractional difference equation in (5.1), we have for all 𝑢 { 1 , 0 } such that | ( 𝑢 + 1 ) / 𝑢 | > 𝑟 , 𝕊 𝑎 Δ 𝜈 𝑎 + 𝜈 𝑁 𝑦 ( 𝑢 ) = 𝕊 𝑎 { 𝑓 } ( 𝑢 ) . ( 5 . 4 ) Then from (4.24), it follows ( 𝑢 + 1 ) 𝑁 𝜈 𝑢 𝑁 𝕊 𝑎 + 𝜈 𝑁 { 𝑦 } ( 𝑢 ) 𝑁 1 𝑘 = 0 𝑢 𝑘 𝑁 Δ 𝜈 𝑁 + 𝑘 𝑎 + 𝜈 𝑁 𝑦 ( 𝑎 ) = 𝕊 𝑎 { 𝑓 } ( 𝑢 ) ( 5 . 5 ) and hence 𝕊 𝑎 + 𝜈 𝑁 𝑢 { 𝑦 } ( 𝑢 ) = 𝑁 ( 𝑢 + 1 ) 𝑁 𝜈 𝕊 𝑎 { 𝑓 } ( 𝑢 ) + 𝑁 1 𝑘 = 0 𝑢 𝑘 ( 𝑢 + 1 ) 𝑁 𝜈 Δ 𝜈 𝑁 + 𝑘 𝑎 + 𝜈 𝑁 𝑦 ( 𝑎 ) . ( 5 . 6 ) By (4.20), we have 𝑢 𝑁 ( 𝑢 + 1 ) 𝑁 𝜈 𝕊 𝑎 { 𝑓 } ( 𝑢 ) = 𝕊 𝑎 + 𝜈 𝑁 Δ 𝑎 𝜈 𝑓 ( 𝑢 ) . ( 5 . 7 ) Considering the terms in the summation, by using the shifting formula (4.5), we see that for each 𝑘 { 0 , 1 , , 𝑁 1 } , 𝑢 𝑘 ( 𝑢 + 1 ) 𝑁 𝜈 = 𝑢 𝑢 + 1 𝑘 ( 𝑢 + 1 ) 𝜈 + 𝑘 𝑁 = 𝑢 𝑢 + 1 𝑘 𝕊 𝑎 + 𝜈 + 𝑘 𝑁 𝜈 + 𝑘 𝑁 ( , 𝑎 ) ( 𝑢 ) = 𝕊 𝑎 + 𝜈 𝑁 𝜈 + 𝑘 𝑁 ( ( 1 , 𝑎 ) 𝑢 ) 𝑢 𝑘 1 𝑖 = 0 𝑢 𝑢 + 1 𝑖 + 1 𝜈 + 𝑘 𝑁 ( 𝑖 + 𝑎 + 𝜈 𝑁 , 𝑎 ) = 𝕊 𝑎 + 𝜈 𝑁 𝜈 + 𝑘 𝑁 ( , 𝑎 ) ( 𝑢 ) ( 5 . 8 ) since 𝜈 + 𝑘 𝑁 ( 𝑖 + 𝑎 + 𝜈 𝑁 , 𝑎 ) = ( 𝑖 + 𝜈 𝑁 ) 𝜈 + 𝑘 𝑁 Γ = ( 𝜈 + 𝑘 𝑁 + 1 ) Γ ( 𝑖 + 𝜈 𝑁 + 1 ) Γ ( 𝑖 𝑘 + 1 ) Γ ( 𝜈 + 𝑘 𝑁 + 1 ) = 0 ( 5 . 9 ) for 𝑖 { 0 , 𝑘 1 } .
Consequently, we have 𝕊 𝑎 + 𝜈 𝑁 { 𝑦 } ( 𝑢 ) = 𝕊 𝑎 + 𝜈 𝑁 Δ 𝑎 𝜈 𝑓 ( 𝑢 ) + 𝑁 1 𝑘 = 0 Δ 𝜈 𝑁 + 𝑘 𝑎 + 𝜈 𝑁 𝑦 ( 𝑎 ) 𝕊 𝑎 + 𝜈 𝑁 𝜈 + 𝑘 𝑁 ( , 𝑎 ) ( 𝑢 ) = 𝕊 𝑎 + 𝜈 𝑁 𝑁 1 𝑘 = 0 Δ 𝜈 𝑁 + 𝑘 𝑎 + 𝜈 𝑁 𝑦 ( 𝑎 ) 𝜈 + 𝑘 𝑁 ( , 𝑎 ) + Δ 𝑎 𝜈 𝑓 ( 𝑢 ) . ( 5 . 1 0 ) Since Sumudu transform is a one-to-one operator (see [28, Theorem 3.6]), we conclude that for 𝑡 𝑎 + 𝜈 𝑁 , 𝑦 ( 𝑡 ) = 𝑁 1 𝑘 = 0 Δ 𝜈 𝑁 + 𝑘 𝑎 + 𝜈 𝑁 𝑦 ( 𝑎 ) 𝜈 + 𝑘 𝑁 ( 𝑡 , 𝑎 ) + Δ 𝑎 𝜈 = 𝑓 ( 𝑡 ) 𝑁 1 𝑘 = 0 Δ 𝜈 𝑁 + 𝑘 𝑎 + 𝜈 𝑁 𝑦 ( 𝑎 ) Γ ( 𝜈 + 𝑘 𝑁 + 1 ) ( 𝑡 𝑎 ) 𝜈 + 𝑘 𝑁 + Δ 𝑎 𝜈 𝑓 ( 𝑡 ) , ( 5 . 1 1 ) where Δ 𝜈 𝑁 + 𝑘 𝑎 + 𝜈 𝑁 𝑦 ( 𝑎 ) = Γ ( 𝜈 + 𝑘 𝑁 + 1 ) 𝑘 𝑝 = 0 𝑘 𝑝 𝑗 = 0 ( 1 ) 𝑗 𝑘 ! ( 𝑘 𝑗 ) 𝑁 𝜈 𝑘 𝑝 𝑗 Δ 𝑘 𝑝 𝑘 𝑦 ( 𝑎 + 𝜈 𝑁 ) , ( 5 . 1 2 ) (see [27, Theorem 11]).

Example 5.2. Consider the initial value problem (5.1) with the Riemann-Liouville fractional difference replaced by the Caputo fractional difference. 𝐶 Δ 𝜈 𝑎 + 𝜈 𝑁 𝑦 ( 𝑡 ) = 𝑓 ( 𝑡 ) , 𝑡 𝑎 , Δ 𝑘 𝑦 ( 𝑎 + 𝜈 𝑁 ) = 𝐴 𝑘 , 𝑘 { 0 , 1 , , 𝑁 1 } , 𝐴 𝑘 . ( 5 . 1 3 )
Applying the Sumudu transform to both sides of the difference equation, we get for all 𝑢 { 1 , 0 } such that | ( 𝑢 + 1 ) / 𝑢 | > 𝑟 , 𝕊 𝑎 𝐶 Δ 𝜈 𝑎 + 𝜈 𝑁 𝑦 ( 𝑢 ) = 𝕊 𝑎 { 𝑓 } ( 𝑢 ) . ( 5 . 1 4 )
Then from (4.27), it follows ( 𝑢 + 1 ) 𝑁 𝜈 𝑢 𝑁 𝕊 𝑎 + 𝜈 𝑁 { 𝑦 } ( 𝑢 ) 𝑁 1 𝑘 = 0 𝑢 𝑘 𝐴 𝑘 = 𝕊 𝑎 { 𝑓 } ( 𝑢 ) . ( 5 . 1 5 ) By (4.20), we have 𝕊 𝑎 + 𝜈 𝑁 { 𝑦 } ( 𝑢 ) = 𝑁 1 𝑘 = 0 𝑢 𝑘 𝐴 𝑘 + 𝑢 𝑁 ( 𝑢 + 1 ) 𝑁 𝜈 𝕊 𝑎 = { 𝑓 } ( 𝑢 ) 𝑁 1 𝑘 = 0 𝑢 𝑘 𝐴 𝑘 + 𝕊 𝑎 + 𝜈 𝑁 Δ 𝑎 𝜈 𝑓 ( 𝑢 ) . ( 5 . 1 6 ) Since from [28], we have 𝕊 0 { 𝑡 𝑛 } ( 𝑢 ) = 𝑛 ! 𝑢 𝑛 , 𝑛 0 , ( 5 . 1 7 ) hence 𝑦 ( 𝑡 ) = 𝑁 1 𝑘 = 0 𝐴 𝑘 ( 𝑡 𝑎 𝜈 + 𝑁 ) 𝑘 𝑘 ! + Δ 𝑎 𝜈 𝑓 ( 𝑡 ) . ( 5 . 1 8 )

Remark 5.3. The initial value problem (5.1) can also be solved by using Proposition 15 in [12].

Example 5.4. Consider the initial value problem 𝐶 Δ 𝜈 + 1 𝑎 + 𝜈 1 𝑦 ( 𝑡 ) 𝐶 Δ 𝜈 𝑎 + 𝜈 1 𝑦 ( 𝑡 ) = 0 , 𝑡 𝑎 , Δ 𝑘 𝑦 ( 𝑎 + 𝑣 𝑁 ) = 𝐴 𝑘 , 𝑘 { 0 . 1 } , 𝐴 𝑘 , ( 5 . 1 9 ) where 0 < 𝜈 1 . Applying the Sumudu transform to both sides of the equation and using (4.31) and (4.27), respectively, we get ( 𝑢 + 1 ) 1 𝜈 𝑢 2 𝕊 𝑎 + 𝜈 1 { 𝑦 } ( 𝑢 ) 𝐴 0 𝑢 𝐴 1 ( 𝑢 + 1 ) 1 𝜈 𝑢 𝕊 𝑎 + 𝜈 1 { 𝑦 } ( 𝑢 ) 𝐴 0 = 0 . ( 5 . 2 0 ) Hence we get 𝕊 𝑎 + 𝜈 1 𝐴 { 𝑦 } ( 𝑢 ) = 0 𝐴 1 + 𝐴 1 . 1 𝑢 ( 5 . 2 1 ) Since from [28], we have 𝕊 0 ( 1 + 𝜆 ) 𝑡 1 ( 𝑢 ) = | | | 1 𝜆 𝑢 f o r ( 1 + 𝜆 ) 𝑢 | | | 𝑢 + 1 < 1 , ( 5 . 2 2 ) then 𝑦 𝐴 ( 𝑡 ) = 0 𝐴 1 + 𝐴 1 2 𝑡 𝑎 𝜈 + 1 . ( 5 . 2 3 )

References

  1. A. A. Kilbas, H. M. Srivastava, and J. J. Trujillo, Theory and Applications of Fractional Differential Equations, vol. 204 of North-Holland Mathematics Studies, Elsevier Science, Amsterdam, The Netherlands, 2006.
  2. S. G. Samko, A. A. Kilbas, and O. I. Marichev, Fractional Integrals and Derivatives—Theory and Applications, Gordon and Breach Science Publishers, Linghorne, Pa, USA, 1993.
  3. I. Podlubny, Fractional Differential Equations, vol. 198 of Mathematics in Science and Engineering, Academic Press, San Diego, Calif, USA, 1999.
  4. R. L. Magin, Fractional Calculus in Bioengineering, Begell House Publisher, Redding, Conn, USA, 2006.
  5. B. J. West, M. Bologna, and P. Grigolini, Physics of Fractal Operators, Institute for Nonlinear Science, Springer, New York, NY, USA, 2003. View at Publisher · View at Google Scholar
  6. N. Heymans and I. Podlubny, “Physical interpretation of initial conditions for fractional differential equations with Riemann-Liouville fractional derivatives,” Rheologica Acta, vol. 45, pp. 765–771, 2006.
  7. K. S. Miller and B. Ross, “Fractional difference calculus,” in Proceedings of the Univalent Functions, Fractional Calculus, and Their Applications, pp. 139–152, Nihon University, 1989.
  8. F. M. Atici and P. W. Eloe, “A transform method in discrete fractional calculus,” International Journal of Difference Equations, vol. 2, no. 2, pp. 165–176, 2007.
  9. F. M. Atici and P. W. Eloe, “Initial value problems in discrete fractional calculus,” Proceedings of the American Mathematical Society, vol. 137, no. 3, pp. 981–989, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH
  10. F. M. Atıcı and P. W. Eloe, “Discrete fractional calculus with the nabla operator,” Electronic Journal of Qualitative Theory of Differential Equations, no. 3, pp. 1–12, 2009. View at Zentralblatt MATH
  11. T. Abdeljawad and D. Baleanu, “Fractional differences and integration by parts,” Journal of Computational Analysis and Applications, vol. 13, no. 3, pp. 574–582, 2011. View at Zentralblatt MATH
  12. T. Abdeljawad, “On Riemann and Caputo fractional differences,” Computers & Mathematics with Applications, vol. 62, no. 3, pp. 1602–1611, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH
  13. G. K. Watugala, “Sumudu transform: a new integral transform to solve differential equations and control engineering problems,” International Journal of Mathematical Education in Science and Technology, vol. 24, no. 1, pp. 35–43, 1993. View at Publisher · View at Google Scholar · View at Zentralblatt MATH
  14. G. K. Watugala, “The Sumudu transform for functions of two variables,” Mathematical Engineering in Industry, vol. 8, no. 4, pp. 293–302, 2002. View at Zentralblatt MATH
  15. M. A. Asiru, “Sumudu transform and the solution of integral equations of convolution type,” International Journal of Mathematical Education in Science and Technology, vol. 32, no. 6, pp. 906–910, 2001. View at Publisher · View at Google Scholar · View at Zentralblatt MATH
  16. M. A. Aşiru, “Further properties of the Sumudu transform and its applications,” International Journal of Mathematical Education in Science and Technology, vol. 33, no. 3, pp. 441–449, 2002. View at Publisher · View at Google Scholar · View at Zentralblatt MATH
  17. F. B. M. Belgacem, A. A. Karaballi, and S. L. Kalla, “Analytical investigations of the Sumudu transform and applications to integral production equations,” Mathematical Problems in Engineering, no. 3-4, pp. 103–118, 2003. View at Publisher · View at Google Scholar
  18. F. B. M. Belgacem and A. A. Karballi, “Sumudu transform fundemantal properties investigations and applications,” Journal of Applied Mathematics and Stochastic Analysis, vol. 2006, Article ID 91083, 23 pages, 2006. View at Publisher · View at Google Scholar
  19. A. Kılıçman and H. Eltayeb, “On the applications of Laplace and Sumudu transforms,” Journal of the Franklin Institute, vol. 347, no. 5, pp. 848–862, 2010. View at Publisher · View at Google Scholar
  20. F. B. M. Belgacem, “Introducing and analysing deeper Sumudu properties,” Nonlinear Studies, vol. 13, no. 1, pp. 23–41, 2006. View at Zentralblatt MATH
  21. F. Jarad and K. Tas, “Application of Sumudu and double Sumudu transforms to Caputo-Fractional dierential equations,” Journal of Computational Analysis and Applications, vol. 14, no. 3, pp. 475–483, 2012.
  22. Q. D. Katatbeh and F. B. M. Belgacem, “Applications of the Sumudu transform to fractional differential equations,” Nonlinear Studies, vol. 18, no. 1, pp. 99–112, 2011. View at Zentralblatt MATH
  23. M. Bohner and G. Sh. Guseinov, “The h-Laplace and q-Laplace transforms,” Journal of Mathematical Analysis and Applications, vol. 365, no. 1, pp. 75–92, 2010. View at Publisher · View at Google Scholar
  24. S. Hilger, “Analysis on measure chains—a unified approach to continuous and discrete calculus,” Results in Mathematics, vol. 18, no. 1-2, pp. 18–56, 1990. View at Zentralblatt MATH
  25. M. Bohner and A. Peterson, Dynamic Equations on Time Scales, Birkhäuser, Boston, Mass, USA, 2001. View at Publisher · View at Google Scholar
  26. M. Bohner and A. Peterson, Advances in Dynamic equations on Time Scales, Birkhäuser, Boston, Mass, USA, 2003. View at Publisher · View at Google Scholar
  27. M. T. Holm, The theory of discrete fractional calculus: development and application [Ph.D. thesis], 2011.
  28. F. Jarad, K. Bayram, T. Abdeljawad, and D. Baleanu, “On the discrete sumudu transform,” Romanian Reports in Physics. In press.