# On Solutions of Linear Fractional Differential Equations with Uncertainty

Authors: T. Allahviranloo; S. Abbasbandy; M. R. Balooch Shahryari; S. Salahshour; D. Baleanu
Journal: Abstract and Applied Analysis
 (2013)

---

## Abstract

The solutions of linear fuzzy fractional differential equations (FFDEs) under the Caputo differentiability have been investigated. To this end, the fuzzy Laplace transform was used to obtain the solutions of FFDEs. Then, some new results regarding the relation between some types of differentiability have been obtained. Finally, some applicable examples are solved in order to show the ability of the proposed method.


---

## Body

## 1. Introduction

The topic of fuzzy differential equations (FDEs) has been rapidly growing in recent years (see, e.g., [1–8] and the references therein). Recently, applying fractional differential equations has grown increasingly and have excited as a considerable interest both in mathematics and applications, such that they were used in modeling of many physical and chemical processes and in engineering [9–15].In addition, many research papers are published to consider solutions of fractional differential equations (see, e.g., [16–18] and the references therein).Recently, Agarwal et al. [19] proposed the concept of solutions for fractional differential equations with uncertainty. They considered the Riemann-Liouville differentiability to solve FFDEs which is a combination of Hukuhara difference and Riemann-Liouville derivative. There are some basic papers which are written by Bede and Gal [20] who discussed on shortcomings of applications of Hukuhara difference. So, we adopt a generalization of  a strongly generalized differentiability to a fractional case.In the following, we propose the Riemann-Liouville differentiability by using Hukuhara difference which is named as the Riemann-Liouville H-differentiability. Therefore, a direct procedure is adopted to derive such concept which is constructed based on the combination of strongly generalized differentiability [21] and the Riemann-Liouville derivative [14, 15]. Moreover, we suggest the concept of fractional derivatives under Caputo’s differentiability by applying the Hukuhara difference which is named as Caputo’s H-differentiability. Similar to the deterministic cases, construction of Caputo’s derivatives is based on the definitions of the Riemann-Liouville derivatives in fuzzy cases.Consequently, we intend to propose an analytical method to solve FFDEs. Since, considering solutions of FFDEs is a new subject, presets the fractional Green’s functions for fuzzy fractional differential equations is considered and, as particular cases, we obtain the classical harmonic oscillator, the damped harmonic oscillator, relaxation equation, all in the fuzzy fractional versions by using fuzzy laplace transforms; so we should first implement analytical method to solve it; then numerical methods can be applied. To this end, we adopt fractional Green’s functions to solve FFDEs by using fuzzy Laplace transforms method. One can see some useful papers about fuzzy Laplace transforms in [2, 22, 23].Recently, Salahshour et al. proposed some new results toward existence and uniqueness of the  solutions of FFDE [24]. Then, Mazandarani and Kamyad numerically solved the FFDE using the Euler method [25]. Also, Agarwal et al. [26] investigated fuzzy fractional integral equation under compactness type condition.This paper is structured as follows. In Section1, we recall some well-known definitions of fuzzy numbers and present some needed concepts. In Section 2, Caputo’s H-differentiability is introduced, and the relation between the Riemann-Liouville differentiability and Caputo’s H-differentiability and some of their properties is considered. Consequently, the fuzzy Laplace transforms are considered for fuzzy-valued function, and an essential theorem for the Laplace transform of f under Caputo’s H-derivative is given in Section 3. The solutions of FFDEs are investigated by using the fuzzy Laplace transforms and their inverses in Section 4. In Section 5, some examples are solved to illustrate the method. Finally, conclusion is drawn in Section 6.


## 2. Preliminaries

The basic definition of fuzzy numbers can be seen in [27].We denote byℝ the set of all real numbers and the set of all fuzzy numbers on ℝ is indicated by 𝔼. A fuzzy number is a mapping u:ℝ→[0,1] with the following properties.(a)
uis upper semicontinuous.(b)
u is fuzzy convex; that is, u(λx+(1-λ)y)≥min{u(x),u(y)} for all x,y∈ℝ, λ∈[0,1].(c)
u is normal; that is; ∃x0∈ℝ for which u(x0)=1.(d)
suppu={x∈ℝ∣u(x)>0} is the support of the u, and its closure cl(suppu) is compact.An equivalent parametric definition is also given in [28–30] as follows.Definition 1.
A fuzzy numberu in parametric form is a pair (u_,u¯) of functions u_(r), u¯(r), 0≤r≤1, which satisfy the following requirements:  (1)
u_(r) is a bounded nondecreasing left continuous function in (0,1]; and right continuous at 0;(2)
u¯(r) is a bounded nonincreasing left continuous function in (0,1]; and right continuous at 0;(3)
u_(r)≤u¯(r),0≤r≤1.Moreover, we also can present ther-cut representation of fuzzy number as [u]r=[u_(r),u¯(r)] for all 0≤r≤1.According to Zadeh’s extension principle, operation of addition on𝔼 is defined by
(1)(u⊕v)(x)=supy∈ℝmin{u(y),v(x-y)},x∈ℝ,
and scalar multiplication of a fuzzy number is given by
(2)(k⊙u)(x)={u(xk),k>0,0~,k=0,
where 0~∈𝔼.The Hausdorff distance between fuzzy numbers is given byd:𝔼×𝔼→ℝ+∪{0},
(3)d(u,v)=supr∈[0,1]max{|u_(r)-v_(r)|,|u¯(r)-v¯(α)|},
where u=(u_(r),u¯(r)), v=(v_(r),v¯(r))⊂ℝ is utilized in [20]. Then, it is easy to see that d is a metric in 𝔼 and has the following properties (see [31]):(1)
d(u+w,v+w)=d(u,v), for all u,v,w∈𝔼,(2)
d(ku,kv)=|k|d(u,v), for all k∈ℝ, u,v∈𝔼,(3)
d(u+v,w+e)≤d(u,w)+d(v,e), for all u,v,w,e∈𝔼,(4)
(d,𝔼) is a complete metric space.Theorem 2 (see [32]).
Letf(x) be a fuzzy-valued function on [a,∞), and it is represented by (f_(x;r),f¯(x;r)). For any fixed r∈[0,1], assume that f_(x;r) and f¯(x;r) are Riemann-integrable on [a,b] for every b≥a, and assume that there are two positive M_(r) and M¯(r) such that ∫ab‍|f_(x;r)|dx≤M_(r) and ∫ab‍|f¯(x;r)|dx≤M¯(r) for every b≥a. Then f(x) is improper fuzzy Riemann-integrable on [a,∞), and the improper fuzzy Riemann-integral is a fuzzy number. Furthermore, one has:
(4)∫a∞f(x)dx=(∫a∞f_(x;r)dx,∫a∞f¯(x;r)dx).Definition 3.
Letx,y∈𝔼. If there exists z∈𝔼 such that x=y+z, then z is called the H-difference of x and y, and it is denoted by x⊖y.In this paper, the sign “⊖” always stands for H-difference, and also note that x⊖y≠x+(-1)y.


### 2.1. Caputo’s H-Differentiability

In this section, the concept of the fuzzy Caputo derivatives has been reviewed [33–35]. Also, we denote by C𝔽[a,b] a space of all fuzzy-valued functions which are continuous on [a,b]. Also, we denote the space of all Lebesgue integrable fuzzy-value functions on the bounded interval [a,b]⊂ℝ by L𝔽[a,b], and we denote the space of fuzzy-value functions f(x) which have continuous H-derivative up to order n-1 on [a,b] such that f(n-1)(x)∈AC𝔽([a,b]) by AC(n)𝔽([a,b]).Definition 4.
Letf∈C𝔽[a,b]∩L𝔽[a,b]; the fuzzy Riemann-Liouville integral of fuzzy-valued function f is defined as follows:
(5)(Ia+βf)(x)=1Γ(β)∫ax‍f(t)dt(x-t)1-β,x>a,0<β<1.Sincef(x;r)=[f_(x;r),f¯(x;r)], for all 0≤r≤1, then one can indicate the fuzzy Riemann-Liouville integral of fuzzy-valued function f based on the lower and upper functions as follows.Theorem 5.
Letf(x)∈C𝔽[a,b]∩L𝔽[a,b]; the fuzzy Riemann-Liouville integral of fuzzy-valued function f is defined by
(6)(Ia+βf)(x;r)=[(Ia+βf_)(x;r),(Ia+βf¯)(x;r)],0≤r≤1,
where
(7)(Ia+βf_)(x;r)=1Γ(β)∫ax‍f_(t;r)dt(x-t)1-β,0≤r≤1,(Ia+βf¯)(x;r)=1Γ(β)∫ax‍f¯(t;r)dt(x-t)1-β,0≤r≤1.Now, we define the fuzzy Riemann-Liouville fractional derivatives of ordern-1<β<n for fuzzy-valued function f (which is a direct extension of strongly generalized H-differentiability [20] and lateral type of H-differentiability [36] in the fractional literature) as follows.Definition 6.
Letfn∈C(n)𝔽[a,b]∩L1𝔽[a,b], let x0∈(a,b), and let Φ(x)=(1/Γ(n-β))∫ax‍(f(t)dt/(x-t)β-n+1) (n=[β]+1;x>a). We say that f is fuzzy Riemann-Liouville fractional differentiable of order β, at x0, if there exists an element (Da+βf)(x0)∈𝔼, such that for all h>0 sufficiently small(i)(8)(Da+βf)(x0)=limh→0Φ(n-1)(x0+h)⊖Φ(n-1)(x0)h=limh→0Φ(n-1)(x0)⊖Φ(n-1)(x0-h)h
or(ii)(9)(Da+βf)(x0)=limh→0Φ(n-1)(x0)⊖Φ(n-1)(x0+h)-h=limh→0Φ(n-1)(x0-h)⊖Φ(n-1)(x0)-h.For sake of simplicity, we say that a fuzzy-valued functionf is RL[(i)-β]-differentiable if it is differentiable as in Definition 6 case (i) and is RL[(ii)-β]-differentiable if it is differentiable as in Definition 6 case (ii).Theorem 7.
Letf(x)∈C(n)𝔽[a,b]∩L𝔽[a,b], x0∈(a,b), (β>0), and (n=[β]+1) such that for all 0≤r≤1, then
(i) iff(x) is a RL[(i)-β]-differentiable fuzzy-valued function, then
(10)(Da+βRLf)(x0;r)=[Da+βRLf_(x0,r),RLDa+βf¯(x0,r)];
(ii) iff(x) is a RL[(ii)-β]-differentiable fuzzy-valued function, then
(11)(DRLa+βf)(x0;r)=[Da+βf-RL(x0;r),RLDa+βf_(x0,r)],
where
(12)(Da+βf_RL)(x0;r)=[1Γ(n-β)dndxn∫ax‍f_(t;r)dt(x-t)β-n+1]x=x0,(Da+βf¯RL)(x0;r)=[1Γ(n-β)dndxn∫ax‍f¯(t,r)dt(x-t)β-n+1]x=x0.Definition 8.
Letf∈C(n)𝔽[a,b]∩L1𝔽[a,b], and f(n) is integrable; then the right fuzzy Caputo derivative of f for n-1<β<n and x∈[a,b] is denoted by (CDa+βf)(x)∈𝔼 and defined by
(13)(DCa+βf)(x)=1Γ(n-β)⊙∫ax‍(x-t)-β+n-1⊙f(n)(t)dt.Theorem 9.
Letf(x)∈C(n)𝔽[a,b]∩L𝔽[a,b], x∈(a,b), and n-1<β<n, such that for all 0≤r≤1, then
(i) iff is a C[(i)-β]-differentiable fuzzy-valued function, then
(14)(Da+βCf)(x;r)=[(Da+βCf_)(x;r),(Da+βCf¯)(x;r)];
(ii) iff(x) is a C[(ii)-β]-differentiable fuzzy-valued function, then
(15)(Da+βCf)(x;r)=[(Da+βCf¯)(x;r),(Da+βCf_)(x;r)],
where
(16)(Da+βCf_)(x;r)=Da+βRL[f_(t;r)-∑k=0n-1xkk!f_(k)(a;r)](x;r)=1Γ(n-β)dndxn×∫ax‍f_(t;r)-∑k=0n-1(tk/k!)f_(k)(a;r)(x-t)β-n+1dt,(Da+βf¯C)(x;r)=Da+βRL[f¯(t;r)-∑k=0n-1‍xkk!f¯(k)(a;r)](x;r)=1Γ(n-β)dndxn×∫ax‍f¯(t,r)-∑k=0n-1(tk/k!)f¯(k)(a;r)(x-t)β-n+1dt.Proof.
We prove the case ofC[(i)-β]-differentiability, (β)>0,  (n=[β]+1), and for the case of C[(ii)-β]-differentiability the proof is completely similar to previous one, and hence it is omitted.
Let us consider thatf is C[(i)-β]-differentiability and that x∈(a,b); then we have the following:
(17)[Φ(n-1)(x+h)⊖Φ(n-1)(x)]r=[-Φ_(n-1)(x;r),Φ¯(n-1)(x+h;r)-Φ¯(n-1)(x;r)Φ_(n-1)(x+h;r)-Φ_(n-1)(x;r),Φ¯(n-1)(x+h;r)-Φ¯(n-1)(x;r)],
and multiplying by 1/h, we get the following:
(18)[Φ(n-1)(x+h)⊖Φ(n-1)(x)]rh=[Φ_(n-1)(x+h;r)-Φ_(n-1)(x;r)h,Φ¯(n-1)(x+h;r)-Φ¯(n-1)(x;r)h].
Similarly, we get the following:(19)[Φ(n-1)(x)⊖Φ(n-1)(x-h)]rh=[Φ_(n-1)(x;r)-Φ_(n-1)(x-h;r)h,Φ¯(n-1)(x;r)-Φ¯(n-1)(x-h;r)h].
Passing to the limit we obtain(20)(DCa+βf)(x;r)=[DCa+βf_(x;r),DCa+βf¯(x;r)],
which proves the theorem.


## 2.1. Caputo’s H-Differentiability

In this section, the concept of the fuzzy Caputo derivatives has been reviewed [33–35]. Also, we denote by C𝔽[a,b] a space of all fuzzy-valued functions which are continuous on [a,b]. Also, we denote the space of all Lebesgue integrable fuzzy-value functions on the bounded interval [a,b]⊂ℝ by L𝔽[a,b], and we denote the space of fuzzy-value functions f(x) which have continuous H-derivative up to order n-1 on [a,b] such that f(n-1)(x)∈AC𝔽([a,b]) by AC(n)𝔽([a,b]).Definition 4.
Letf∈C𝔽[a,b]∩L𝔽[a,b]; the fuzzy Riemann-Liouville integral of fuzzy-valued function f is defined as follows:
(5)(Ia+βf)(x)=1Γ(β)∫ax‍f(t)dt(x-t)1-β,x>a,0<β<1.Sincef(x;r)=[f_(x;r),f¯(x;r)], for all 0≤r≤1, then one can indicate the fuzzy Riemann-Liouville integral of fuzzy-valued function f based on the lower and upper functions as follows.Theorem 5.
Letf(x)∈C𝔽[a,b]∩L𝔽[a,b]; the fuzzy Riemann-Liouville integral of fuzzy-valued function f is defined by
(6)(Ia+βf)(x;r)=[(Ia+βf_)(x;r),(Ia+βf¯)(x;r)],0≤r≤1,
where
(7)(Ia+βf_)(x;r)=1Γ(β)∫ax‍f_(t;r)dt(x-t)1-β,0≤r≤1,(Ia+βf¯)(x;r)=1Γ(β)∫ax‍f¯(t;r)dt(x-t)1-β,0≤r≤1.Now, we define the fuzzy Riemann-Liouville fractional derivatives of ordern-1<β<n for fuzzy-valued function f (which is a direct extension of strongly generalized H-differentiability [20] and lateral type of H-differentiability [36] in the fractional literature) as follows.Definition 6.
Letfn∈C(n)𝔽[a,b]∩L1𝔽[a,b], let x0∈(a,b), and let Φ(x)=(1/Γ(n-β))∫ax‍(f(t)dt/(x-t)β-n+1) (n=[β]+1;x>a). We say that f is fuzzy Riemann-Liouville fractional differentiable of order β, at x0, if there exists an element (Da+βf)(x0)∈𝔼, such that for all h>0 sufficiently small(i)(8)(Da+βf)(x0)=limh→0Φ(n-1)(x0+h)⊖Φ(n-1)(x0)h=limh→0Φ(n-1)(x0)⊖Φ(n-1)(x0-h)h
or(ii)(9)(Da+βf)(x0)=limh→0Φ(n-1)(x0)⊖Φ(n-1)(x0+h)-h=limh→0Φ(n-1)(x0-h)⊖Φ(n-1)(x0)-h.For sake of simplicity, we say that a fuzzy-valued functionf is RL[(i)-β]-differentiable if it is differentiable as in Definition 6 case (i) and is RL[(ii)-β]-differentiable if it is differentiable as in Definition 6 case (ii).Theorem 7.
Letf(x)∈C(n)𝔽[a,b]∩L𝔽[a,b], x0∈(a,b), (β>0), and (n=[β]+1) such that for all 0≤r≤1, then
(i) iff(x) is a RL[(i)-β]-differentiable fuzzy-valued function, then
(10)(Da+βRLf)(x0;r)=[Da+βRLf_(x0,r),RLDa+βf¯(x0,r)];
(ii) iff(x) is a RL[(ii)-β]-differentiable fuzzy-valued function, then
(11)(DRLa+βf)(x0;r)=[Da+βf-RL(x0;r),RLDa+βf_(x0,r)],
where
(12)(Da+βf_RL)(x0;r)=[1Γ(n-β)dndxn∫ax‍f_(t;r)dt(x-t)β-n+1]x=x0,(Da+βf¯RL)(x0;r)=[1Γ(n-β)dndxn∫ax‍f¯(t,r)dt(x-t)β-n+1]x=x0.Definition 8.
Letf∈C(n)𝔽[a,b]∩L1𝔽[a,b], and f(n) is integrable; then the right fuzzy Caputo derivative of f for n-1<β<n and x∈[a,b] is denoted by (CDa+βf)(x)∈𝔼 and defined by
(13)(DCa+βf)(x)=1Γ(n-β)⊙∫ax‍(x-t)-β+n-1⊙f(n)(t)dt.Theorem 9.
Letf(x)∈C(n)𝔽[a,b]∩L𝔽[a,b], x∈(a,b), and n-1<β<n, such that for all 0≤r≤1, then
(i) iff is a C[(i)-β]-differentiable fuzzy-valued function, then
(14)(Da+βCf)(x;r)=[(Da+βCf_)(x;r),(Da+βCf¯)(x;r)];
(ii) iff(x) is a C[(ii)-β]-differentiable fuzzy-valued function, then
(15)(Da+βCf)(x;r)=[(Da+βCf¯)(x;r),(Da+βCf_)(x;r)],
where
(16)(Da+βCf_)(x;r)=Da+βRL[f_(t;r)-∑k=0n-1xkk!f_(k)(a;r)](x;r)=1Γ(n-β)dndxn×∫ax‍f_(t;r)-∑k=0n-1(tk/k!)f_(k)(a;r)(x-t)β-n+1dt,(Da+βf¯C)(x;r)=Da+βRL[f¯(t;r)-∑k=0n-1‍xkk!f¯(k)(a;r)](x;r)=1Γ(n-β)dndxn×∫ax‍f¯(t,r)-∑k=0n-1(tk/k!)f¯(k)(a;r)(x-t)β-n+1dt.Proof.
We prove the case ofC[(i)-β]-differentiability, (β)>0,  (n=[β]+1), and for the case of C[(ii)-β]-differentiability the proof is completely similar to previous one, and hence it is omitted.
Let us consider thatf is C[(i)-β]-differentiability and that x∈(a,b); then we have the following:
(17)[Φ(n-1)(x+h)⊖Φ(n-1)(x)]r=[-Φ_(n-1)(x;r),Φ¯(n-1)(x+h;r)-Φ¯(n-1)(x;r)Φ_(n-1)(x+h;r)-Φ_(n-1)(x;r),Φ¯(n-1)(x+h;r)-Φ¯(n-1)(x;r)],
and multiplying by 1/h, we get the following:
(18)[Φ(n-1)(x+h)⊖Φ(n-1)(x)]rh=[Φ_(n-1)(x+h;r)-Φ_(n-1)(x;r)h,Φ¯(n-1)(x+h;r)-Φ¯(n-1)(x;r)h].
Similarly, we get the following:(19)[Φ(n-1)(x)⊖Φ(n-1)(x-h)]rh=[Φ_(n-1)(x;r)-Φ_(n-1)(x-h;r)h,Φ¯(n-1)(x;r)-Φ¯(n-1)(x-h;r)h].
Passing to the limit we obtain(20)(DCa+βf)(x;r)=[DCa+βf_(x;r),DCa+βf¯(x;r)],
which proves the theorem.


## 3. The Fuzzy Laplace Transforms

In this section, we consider the fuzzy Laplace transform for fuzzy-valued function; then derivative theorem is given which is essential to determine solutions of FFDEs.In this way, Allahviranloo and Ahmadi [2] suggested the concept of Laplace transforms for fuzzy-valued function as follows.Definition 10 (see [2]).
Letf(x) be continuous fuzzy-value function. Suppose that f(x)⊙e-px is improper fuzzy Riemann-integrable on [0,∞); then ∫0∞‍f(x)⊙e-pxdx is called fuzzy Laplace transforms and is denoted by
(21)L[f(x)]=∫0∞‍f(x)⊙e-pxdx(p>0and integer).Remark 11.
In [8], the authors have investigated under what conditions the fuzzy-valued functions can possess the fuzzy Laplace transform. So, we suppose that the given fuzzy-valued functions have mentioned conditions throughout the paper.From Theorem2, we have
(22)∫0∞‍f(x)⊙e-pxdx=(∫0∞‍f_(x;r)e-pxdx,∫0∞‍f¯(x;r)e-pxdx).
Also by using the definition of classical Laplace transform one has
(23)ℓ[f_(x;r)]=∫0∞‍f_(x;r)e-pxdx,ℓ[f¯(x,r)]=∫0∞‍f¯(x;r)e-pxdx;
then we get
(24)L[f(x)]=(ℓ[f_(x;r)],ℓ[f¯(x;r)]).Theorem 12 (see [2]).
Letf(x), g(x) be continuous-fuzzy-valued functions; suppose that c1, c2 are constant; then
(25)L[(c1⊙f(x))⊕(c2⊙g(x))]=(c1⊙L[f(x)])⊕(c2⊙L[g(x)]).Lemma 13 (see [2]).
Letf(x) be continuous fuzzy-value function on [0,∞) and let λ∈ℝ; then
(26)L[λ⊙f(x)]=λ⊙L[f(x)].Lemma 14 (see [2]).
Letf(x) be continuous fuzzy-value function and g(x)≥0. Suppose that (f(x)⊙g(x))⊙e-px is improper fuzzy Riemann-integrable on [0,∞); then
(27)∫0∞‍(f(x)⊙g(x))⊙e-pxdx=(∫0∞‍g(x)f_(x,r)e-pxdx,∫0∞‍g(x)f¯(x,r)e-pxdx).Theorem 15 (first translation theorem (see [2])).
Letf be continuous fuzzy value function and let L[f(x)]=F(p). Then
(28)L[eax⊙f(x)]=F(p-a),
where eax is real value function.In order to solve FFDEs, it is necessary to know the fuzzy Laplace transform of Caputo’s H-derivative off, (DCa+βf)(x). The virtue of L[(DCa+βf)(x)] is that it can be written in terms of L[f(x)].Theorem 16 (derivative theorem).
Suppose thatf(x)∈C(n)𝔽[0,∞)∩L𝔽[0,∞); then(29)L[(DCa+βf)(x)]=pβL[f(t)]⊖(∑k=0n-1‍pβ-k-1f(k))(0)
if D(k-1)f is C[(i)-β]-differentiable, and also
(30)L[(DCa+βf)(x)]=-(∑k=0n-1‍pβ-k-1f(k))(0)⊖(-pβL[f(t)])
if D(k-1)f is C[(ii)-β]-differentiable.Proof.
For arbitrary fixedr∈[0,1] we have
(31)pβL[f(t)]⊖pβ-1f(0)=[(∑k=0n-1‍pβ-k-1f¯(k))pβℓ[f_(x;r)]-(∑k=0n-1‍pβ-k-1f_(k))(0;r),pβℓ[f¯(x;r)]-(∑k=0n-1‍pβ-k-1f¯(k))(0;r)].
Since D(k-1)f isC[(i)-β], we get the following:
(32)(Da+βCf)(x;r)=[(Da+βCf_)(x;r),(Da+βCf¯)(x;r)]=[(Da+βCf_)(x;r),(Da+βf¯C)(x;r)].
Hence, we get
(33)ℓ[(Da+βCf¯)(x;r)]=ℓ[(Da+βCf¯)(x;r)]=pβℓ[f¯(x;r)]-(∑k=0n-1pβ-k-1f¯(k))(0;r),ℓ[(Da+βCf_)(x;r)]=ℓ[(Da+βCf_)(x;r)]=pβℓ[f_(x;r)]-(∑k=0n-1pβ-k-1f_(k))(0;r).
Then, we obtained the following:
(34)pβL[f(x)]⊖(∑k=0n-1‍pβ-k-1f(k))(0)=(ℓ[(DCa+βf_)(x;r)],ℓ[(DCa+βf¯)(x;r)]).
Then, we have
(35)pβL[f(x)]⊖(∑k=0n-1‍pβ-k-1f(k))(0)=L[(Da+βCf_)(x;r),(Da+βCf¯)(x;r)].
Using (32) leads to obtain
(36)pβL[f(x)]⊖(∑k=0n-1‍pβ-k-1f(k))(0)=L[(DCa+βf)(x)].
Now, we assume that D(k-1)f is the C[(ii)-β]; then for arbitrary fixed r∈[0,1] we have
(37)-(∑k=0n-1‍pβ-k-1f(k))(0)⊖(-pβL[f(t)])=(-f2(t;r),-f1(t;r)),-f2(t;r)=-(∑k=0n-1‍pβ-k-1f¯(k))(0;r)+pβℓ[f¯(x;r)],-f1(t;r)=-(∑k=0n-1‍pβ-k-1f_(k))(0;r)+pβℓ[f_(x;r)].
Since D(k-1)f is C[(ii)-β], we get
(38)(DCa+βf)(x;r)=[(DCa+βf_)(x;r),(DCa+βf¯)(x;r)]=[(DCa+βf¯)(x;r),(DCa+βf_)(x;r)],
and, hence, we get
(39)-(∑k=0n-1‍pβ-k-1f(k))(0)⊖(-pβL[f(t)])=pβℓ[f¯(x;r)]-(∑k=0n-1‍pβ-k-1f¯(k))(0;r),pβℓ[f_(x;r)]-(∑k=0n-1‍pβ-k-1f_(k))(0;r).
So, we get
(40)-(∑k=0n-1‍pβ-k-1f(k))(0)⊖(-pβL[f(t)])=[ℓ[(DCa+βf¯)(x;r),ℓ(DCa+βf_)(x;r)]].
Thus, we obtain the following:
(41)-(∑k=0n-1‍pβ-k-1f(k))(0)⊖(-pβL[f(t)])=L[(DCa+βf)(x)],
which completes the proof.


## 4. Fractional Green’s Functions

In the following, we present the so-called fundamental solutions of fractional Green’s functions for fuzzy fractional differential equations under Caputo’s H-differentiability. First we consider the nonhomogeneous fuzzy fractional differential equation of order0<β<1 together with nonhomogeneous initial condition
(42)(DC0+2βΦ)(x)⊕(DC0+2β)(x)=φ[x;Φ(x)],0<β<1,Φ(0)=f(x)∈𝔼,
where t>0. Fractional Green’s functions for case so-called the classical Green’s functions, that is, considering φ(x)=δ(x) and Φ(0)=0~, are discussed.


### 4.1. Determining Algebraic Solutions

The algebraic solutions of fuzzy fractional differential equations have been investigated under Caputo’s H-differentiability. We provided the fuzzy Laplace transform to obtain the solutions of FFDE (42). By taking Laplace transform on the both sides of (42), we get
(43)L[(D0+2βCΦ)(x)⊕(D0+βCΦ)(x)]=L[φ(x,Φ(x))].
After that, based on the types of Caputo’s H-differentiability we have the following cases.Case I.
Let us consider thatDΦ(x) is a C[(i)-β]-differentiable function; then (43) is extended based on the its lower and upper functions as follows:
(44)ℓ[φ_(x,Φ(x);r)]=p2βℓ[Φ_(x,r)]-p2β-1Φ_(0,r)+pβℓ[Φ_(x,r)]-pβ-1Φ_(0,r),ℓ[φ¯(x,Φ(x);r)]=p2βℓ[Φ¯(x,r)]-p2β-1Φ¯(0,r)+pβℓ[Φ¯(x,r)]-pβ-1Φ¯(0,r),
where
(45)φ_(x,Φ(x);r)=min{[Φ_(x;r),Φ¯(x;r)]φ(x,u)∣u∈[Φ_(x;r),Φ¯(x;r)]},0≤r≤1,φ¯(x,Φ(x);r)=max{[Φ_(x;r),Φ¯(x;r)]φ(x,u)∣u∈[Φ_(x;r),Φ¯(x;r)]},0≤r≤1.
In order to solve the linear system (44), for simplify we assume that
(46)ℓ[Φ_(x;r)]=H1(p;r),0≤r≤1,ℓ[Φ¯(x;r)]=K1(p;r),0≤r≤1,
where H1(p;r) and K1(p;r) are solutions of system (44). Then, by using the inverse Laplace transform, Φ_(x;r) and Φ¯(x;r) are computed as follows:
(47)Φ_(x,r)=ℓ-1[H1(p;r)]=Φ_(0;r)G1(t)+∫0t‍G1(t;τ)φ_(t-τ;r)dτ,0≤r≤1,Φ¯(x;r)=ℓ-1[K1(p;r)]=Φ¯(0;r)G1(t)+∫0t‍G1(t;τ)φ¯(t-τ;r)dτ,0≤r≤1.
To explicit the corresponding Green’s function, which is solution of the above relation satisfying the condition Φ(0)=0~ and taking φ(x)=δ(x), we use the linearity of the inverse Laplace transform to obtain
(48)G1(t∣0)=ℓ-1[Φ(x;r)].Case II.
Let us consider thatDΦ(x) is C[(ii)-β]-differentiable; then (43) can be written as follows:
(49)ℓ[φ_(x,Φ(x);r)]=p2βℓ[Φ_(x,r)]-p2β-1Φ_(0,r)+pβℓ[Φ_(x,r)]-pβ-1Φ_(0,r),ℓ[φ¯(x,Φ(x);r)]=p2βℓ[Φ¯(x,r)]-p2β-1Φ¯(0,r)+pβℓ[Φ¯(x,r)]-pβ-1Φ¯(0,r),
where
(50)φ_(x,Φ(x);r)=min{φ(x,u)∣u∈[Φ_(x;r),Φ¯(x;r)]},0≤r≤1,φ¯(x,Φ(x);r)=max{φ(x,u)∣u∈[Φ_(x;r),Φ¯(x;r)]},0≤r≤1.
Consequently, to solve the linear system (49), we set for simplify
(51)ℓ[Φ_(x;r)]=H2(p;r),ℓ[Φ¯(x;r)]=K2(p;r),
where H2(p;r) and K2(p;r) are solutions of system (49). Then by making use of the inverse Laplace transform, Φ_(x;r) and Φ¯(x;r) are computed as follows:
(52)Φ_(x;r)=ℓ-1[H2(p;r)]=Φ_(0;r)G2(t)+∫0t‍G2(t;τ)φ_(t-τ;r)dτ,0≤r≤1,Φ¯(x;r)=ℓ-1[K2(p;r)]=Φ¯(0;r)G2(t)+∫0t‍G2(t;τ)φ¯(t-τ;r)dτ,0≤r≤1.
To show the corresponding Green’s function, which is solution of (42) satisfying the condition Φ(0)=0~ and taking φ(x)=δ(x), we use the linearity of the inverse Laplace transform to obtain
(53)G2(t∣0)=ℓ-1[Φ(x;r)].


## 4.1. Determining Algebraic Solutions

The algebraic solutions of fuzzy fractional differential equations have been investigated under Caputo’s H-differentiability. We provided the fuzzy Laplace transform to obtain the solutions of FFDE (42). By taking Laplace transform on the both sides of (42), we get
(43)L[(D0+2βCΦ)(x)⊕(D0+βCΦ)(x)]=L[φ(x,Φ(x))].
After that, based on the types of Caputo’s H-differentiability we have the following cases.Case I.
Let us consider thatDΦ(x) is a C[(i)-β]-differentiable function; then (43) is extended based on the its lower and upper functions as follows:
(44)ℓ[φ_(x,Φ(x);r)]=p2βℓ[Φ_(x,r)]-p2β-1Φ_(0,r)+pβℓ[Φ_(x,r)]-pβ-1Φ_(0,r),ℓ[φ¯(x,Φ(x);r)]=p2βℓ[Φ¯(x,r)]-p2β-1Φ¯(0,r)+pβℓ[Φ¯(x,r)]-pβ-1Φ¯(0,r),
where
(45)φ_(x,Φ(x);r)=min{[Φ_(x;r),Φ¯(x;r)]φ(x,u)∣u∈[Φ_(x;r),Φ¯(x;r)]},0≤r≤1,φ¯(x,Φ(x);r)=max{[Φ_(x;r),Φ¯(x;r)]φ(x,u)∣u∈[Φ_(x;r),Φ¯(x;r)]},0≤r≤1.
In order to solve the linear system (44), for simplify we assume that
(46)ℓ[Φ_(x;r)]=H1(p;r),0≤r≤1,ℓ[Φ¯(x;r)]=K1(p;r),0≤r≤1,
where H1(p;r) and K1(p;r) are solutions of system (44). Then, by using the inverse Laplace transform, Φ_(x;r) and Φ¯(x;r) are computed as follows:
(47)Φ_(x,r)=ℓ-1[H1(p;r)]=Φ_(0;r)G1(t)+∫0t‍G1(t;τ)φ_(t-τ;r)dτ,0≤r≤1,Φ¯(x;r)=ℓ-1[K1(p;r)]=Φ¯(0;r)G1(t)+∫0t‍G1(t;τ)φ¯(t-τ;r)dτ,0≤r≤1.
To explicit the corresponding Green’s function, which is solution of the above relation satisfying the condition Φ(0)=0~ and taking φ(x)=δ(x), we use the linearity of the inverse Laplace transform to obtain
(48)G1(t∣0)=ℓ-1[Φ(x;r)].Case II.
Let us consider thatDΦ(x) is C[(ii)-β]-differentiable; then (43) can be written as follows:
(49)ℓ[φ_(x,Φ(x);r)]=p2βℓ[Φ_(x,r)]-p2β-1Φ_(0,r)+pβℓ[Φ_(x,r)]-pβ-1Φ_(0,r),ℓ[φ¯(x,Φ(x);r)]=p2βℓ[Φ¯(x,r)]-p2β-1Φ¯(0,r)+pβℓ[Φ¯(x,r)]-pβ-1Φ¯(0,r),
where
(50)φ_(x,Φ(x);r)=min{φ(x,u)∣u∈[Φ_(x;r),Φ¯(x;r)]},0≤r≤1,φ¯(x,Φ(x);r)=max{φ(x,u)∣u∈[Φ_(x;r),Φ¯(x;r)]},0≤r≤1.
Consequently, to solve the linear system (49), we set for simplify
(51)ℓ[Φ_(x;r)]=H2(p;r),ℓ[Φ¯(x;r)]=K2(p;r),
where H2(p;r) and K2(p;r) are solutions of system (49). Then by making use of the inverse Laplace transform, Φ_(x;r) and Φ¯(x;r) are computed as follows:
(52)Φ_(x;r)=ℓ-1[H2(p;r)]=Φ_(0;r)G2(t)+∫0t‍G2(t;τ)φ_(t-τ;r)dτ,0≤r≤1,Φ¯(x;r)=ℓ-1[K2(p;r)]=Φ¯(0;r)G2(t)+∫0t‍G2(t;τ)φ¯(t-τ;r)dτ,0≤r≤1.
To show the corresponding Green’s function, which is solution of (42) satisfying the condition Φ(0)=0~ and taking φ(x)=δ(x), we use the linearity of the inverse Laplace transform to obtain
(53)G2(t∣0)=ℓ-1[Φ(x;r)].


## 5. Examples

In this section, we consider several examples, and we find explicit expressions for fractional Green’s functions, including some versions of the fuzzy harmonic oscillator equations and fuzzy relaxation equation under Caputo’s H-differentiability.Example 1 (Green’s function for fuzzy damped harmonic oscillator).
Let us examine the following FFDE:(54)a⊙(D0+2βCΦ)(x)+b⊙(D0+βCΦ)(x)=λw2⊙Φ(x)+φ(x),0<β<1,Φ(0)∈𝔼,
where a,b are real positive constants and w is the frequency of the harmonic oscillator. We will solve this example according to the two following cases for λ∈ℝ.Case I.
Suppose thatλ∈ℝ+=(0,∞); then taking Laplace transform on the both sides of above equation, we have
(55)L[a⊙(D0+2βCΦ)(x)+b⊙(D0+βCΦ)(x)]=L[λw2⊙Φ(x)+φ(x)].
Using C[(i)-β]-differentiability, we get
(56)λw2ℓ[Φ_(x;r)]+ℓ[φ_(x;r)]=ap2βℓ[Φ_(x,r)]-ap2β-1Φ_(0,r)+bpβℓ[Φ_(x,r)]-bpβ-1Φ_(0,r),λw2ℓ[Φ¯(x;r)]+ℓ[φ¯(x;r)]=ap2βℓ[Φ¯(x,r)]-ap2β-1Φ¯(0,r)+bpβℓ[Φ¯(x,r)]-bpβ-1Φ¯(0,r).
Then, after some calculations, we obtain the following:
(57)(ap2β+bpβ-λw2)ℓ[Φ_(x;r)]=(ap2β-1+bpβ-1)Φ_(0;r)+ℓ[φ_(x;r)],(ap2β+bpβ-λw2)ℓ[Φ¯(x;r)]=(ap2β-1+bpβ-1)Φ¯(0;r)+ℓ[φ¯(x;r)].
Consequently, taking the inverse of Laplace on both sides of (57) we have
(58)Φ_(x;r)=Φ_(0;r)ℓ-1×[(ap2β-1+bpβ-1)Φ_(0;r)+ℓ[φ_(x;r)]ap2β+bpβ-λw2],0≤r≤1,Φ¯(x;r)=Φ¯(0;r)ℓ-1×[(ap2β-1+bpβ-1)Φ¯(0;r)+ℓ[φ¯(x;r)]ap2β+bpβ-λw2],0≤r≤1.
Finally, we determine the solution of FFDE as follows:
(59)Φ_(x;r)=Φ_(0;r)∑j=0∞‍(-ba)jxβjE2β,βj+1(j+1)(λw2ax2β)+baxβE2β,βj+β+1(j+1)(λw2ax2β)+1a∑j=0∞‍(-ba)j∫0x‍φ_(t-y;r)yβj+2β-1E2β,βj+2β(j+1)×(λw2ay2β)dy,0≤r≤1,Φ¯(x;r)=Φ¯(0;r)∑j=0∞‍(-ba)jxβjE2β,βj+1(j+1)(λw2ax2β)+baxβE2β,βj+β+1(j+1)(λw2ax2β)+1a∑j=0∞‍(-ba)j∫0x‍φ¯(t-y;r)yβj+2β-1E2β,βj+2β(j+1)×(λw2ay2β)dy,0≤r≤1.
To explicit the corresponding Green’s function under C[(i)-β]-differentiability, which is solution of (54) satisfying the condition Φ(0)=0~ and taking φ(x)=δ(x), we use the linearity of the inverse Laplace transform to obtain
(60)G1(x∣0)≡ℓ-1[ℓ[Φ(x;r)]]=1a∑j=0∞‍(-ba)jtβj+2β-1E2β,βj+2β(j+1)(λw2at2β).
On the other hand, taking φ(x)=0~, we have a homogeneous fuzzy fractional differential equation whose solution under C[(i)-β]-H-differentiability is given by
(61)Φ_(x;r)=Φ_(0;r)∑j=0∞‍(-ba)jxβjE2β,βj+1(j+1)(λw2ax2β)+baxβE2β,βj+β+1(j+1)(λw2ax2β),0≤r≤1,Φ¯(x;r)=Φ¯(0;r)∑j=0∞‍(-ba)jxβjE2β,βj+1(j+1)(λw2ax2β)+baxβE2β,βj+β+1(j+1)(λw2ax2β),0≤r≤1.Case II.
Ifλ∈ℝ-=(-∞,0), then by using C[(ii)-β]-differentiability and Theorem 16, and their inverses, fractional Green’s function will be obtained similar to (60).Example 2 (Green’s function for fuzzy the driven harmonic oscillator).
Let us investigate the following FFDEs; namely,(62)(D0+2βCΦ)(x)=λw2β⊙Φ(x)+φ(x),0<β<1Φ(0)∈𝔼,withb=0 and 1/2<β<1; the fuzzy fractional differential equation associated with the driven harmonic oscillator is obtained. We obtain the same result taking b=0, a=1 and solve this example according to the two following cases for λ∈ℝ.Case I.
Suppose thatλ∈ℝ+=(0,+∞); then taking Laplace transform on the both sides of above equation, we have the following:
(63)L[(Da+2βCΦ)(x)]=L(λw2β⊙Φ(x)+φ(x)).
Using C[(i)-β]-differentiability, we obtain
(64)λw2βℓ[Φ_(x;r)]+ℓ[φ_(x;r)]=p2βℓ[Φ_(x,r)]-p2β-1Φ_(0,r),λw2βℓ[Φ¯(x;r)]+ℓ[φ¯(x;r)]=p2βℓ[Φ¯(x,r)]-p2β-1Φ_(0,r).
Then, after some manipulations the following result is reported:
(65)(p2β-λw2β)ℓ[Φ_(x;r)]=p2β-1Φ_(0,r)+ℓ[φ_(x;r)],(p2β-λw2β)ℓ[Φ¯(x;r)]=p2β-1Φ¯(0,r)+ℓ[φ¯(x;r)].
By taking inverse of Laplace on the both sides of (75) we have
(66)Φ_(x;r)=Φ_(0;r)ℓ-1[p2β-1Φ_(0,r)+ℓ[φ_(x;r)]p2β-λw2β],0≤r≤1,Φ¯(x;r)=Φ¯(0;r)ℓ-1[p2β-1Φ_(0,r)+ℓ[φ¯(x;r)]p2β-λw2β],0≤r≤1.
Finally, we determine the solution of FFDE as follows:
(67)Φ_(x;r)=Φ_(0;r)E2β(λw2βx2β)+∫0x‍φ_(t-y;r)y2β-1E2β,2β(λw2βy2β)dy,0≤r≤1,Φ¯(x;r)=Φ¯(0;r)E2β(λw2βx2β)+∫0x‍φ¯(t-y;r)y2β-1E2β,2β(λw2βy2β)dy,0≤r≤1.
To show the corresponding Green’s function under C[(i)-β]-differentiability, which is solution of (62) satisfying the condition Φ(0)=0~ and taking φ(x)=δ(x), we use the linearity of the inverse Laplace transform  to obtain
(68)G1(x∣0)≡x2β-1E2β,2β(λw2βx2β).
In the case β=1, we report
(69)G1(x∣0)=xE2,2(λw2x2)
which is the classic Green’s function associated with fuzzy the driven harmonic oscillator.Takingφ(x)=0~, we have a homogeneous fuzzy fractional differential equation whose solution under C[(i)-β]-H-differentiability is given by
(70)Φ_(x;r)=Φ_(0;r)E2β(λw2βx2β),0≤r≤1,Φ¯(x;r)=Φ¯(0;r)E2β(λw2βx2β),0≤r≤1.
For β=1, we get
(71)Φ_(x;r)=Φ_(0;r)coswt,0≤r≤1,Φ¯(x;r)=Φ¯(0;r)coswt,0≤r≤1.Case II.
Let us suppose thatλ∈ℝ-=(-∞,0); then using C[(ii)-β]-differentiability and Theorem 16 and their inverses, fractional Green’s function will be obtained similar to (68).Example 3 (Green’s functions for fuzzy fractional relaxation equation).
The next step is to consider the following FFDE:(72)b⊙(DC0+βΦ)(x)=λw2⊙Φ(x)+φ(x),0<β<1,Φ(0)∈𝔼,
where a=0 in (54) and 0<β<1. We solve this example according to the two following cases corresponding to λ∈ℝ.Case I.
Suppose thatλ∈ℝ+=(0,+∞); then taking Laplace transform on the both sides of above equation, we have
(73)L[b⊙(Da+βCΦ)(x)]=L[λw2⊙Φ(x)+φ(x)].
Using C[(i)-β]-differentiability, we get
(74)λwβℓ[Φ_(x;r)]+ℓ[φ_(x;r)]=bpβℓ[Φ_(x,r)]-bpβ-1Φ_(0,r),λwβℓ[Φ¯(x;r)]+ℓ[φ¯(x;r)]=bpβℓ[Φ¯(x,r)]-bpβ-1Φ_(0,r).
Then, after some manipulations we get the following:
(75)(bpβ-λw2)ℓ[Φ_(x;r)]=bpβ-1Φ_(0,r)+ℓ[φ_(x;r)],(bpβ-λw2)ℓ[Φ¯(x;r)]=bpβ-1Φ¯(0,r)+ℓ[φ¯(x;r)].
By taking inverse of Laplace on the both sides of (75) we have
(76)Φ_(x;r)=Φ_(0;r)ℓ-1[bpβ-1Φ_(0,r)+ℓ[φ_(x;r)]bpβ-λw2],0≤r≤1,Φ¯(x;r)=Φ¯(0;r)ℓ-1[bpβ-1Φ_(0,r)+ℓ[φ¯(x;r)]bpβ-λw2],0≤r≤1.
Thus, we determine the solution of FFDE as follows:
(77)Φ_(x;r)=Φ_(0;r)Eβ(λw2bxβ)+1b∫0x‍φ_(t-y;r)yβ-1Eβ,β(λw2byβ)dy,0≤r≤1,Φ¯(x;r)=Φ¯(0;r)Eβ(λw2bxβ)+1b∫0x‍φ¯(t-y;r)yβ-1Eβ,β(λw2byβ)dy,0≤r≤1.
Here, we obtain explicitly Green’s functions under C[(i)-β]-differentiability, introducing φ(x)=δ(x) and Φ(0)=0~; that is,
(78)G1(x∣0)≡1bxβ-1Eβ,β(λwβbxβ).
When β=1, we get the classic Green’s function under strongly generalized differentiability
(79)G1(x∣0)=1bexp(λw2bx).
On the other hand, taking φ(x)=0~, we have a homogeneous fuzzy fractional differential equation whose solution under C[(i)-β]-H-differentiability is given by
(80)Φ_(x;r)=Φ_(0;r)Eβ(λw2bxβ),0≤r≤1,Φ¯(x;r)=Φ¯(0;r)Eβ(λw2bxβ),0≤r≤1.
Finally, in the case β=1 we get
(81)Φ_(x;r)=Φ_(0;r)exp(λw2bx),0≤r≤1,Φ¯(x;r)=Φ¯(0;r)exp(λw2bx),0≤r≤1.Case II.
Suppose thatλ∈ℝ-=(-∞,0); then using C[(ii)-β]-differentiability and Theorem 16, and their inverses, fractional Green’s function will be obtained similar to (78).Now, we provide some examples such that the obtained solutions coincide with the previously reported solution in integer case [37].Example 4.
For a special case, let us consider the following FFDE:(82)(D0+βCΦ)(x)=λ⊙Φ(x),0<β,x<1,Φ(0)∈𝔼,
where y represents the number of radionuclides present in a given radioactive and λ is a decay constant.Case I.
Suppose thatλ∈ℝ+=(0,+∞); then using C[(i)-β]-differentiability and Theorem 16, the solution is derived as follows:
(83)Φ_(x;r)=Φ_(0,r)Eβ,β[λxβ],0≤r≤1,Φ¯(x;r)=Φ¯(0,r)Eβ,β[λxβ],0≤r≤1.Case II.
Suppose thatλ∈ℝ-=(-∞,0); then using C[(ii)-β]-differentiability and Theorem 16, the solution will be obtained similar to (83). For more details see Figures 1, 2, 3, and 4.Figure 1
Solution of Example4, Case I, β=1.Figure 2
Solution of Example4, Case I, β=0.5.Figure 3
Solution of Example4, Case II, β=1.Figure 4
Solution of Example4, Case II, β=0.5.Example 5.
Let us discuss the following FFDE:(84)(DC0+βΦ)(x)=(-1)⊙Φ(x)+x+1,0<β,x<1,Φ(0)∈𝔼.
Then, using C[(ii)-β]-differentiability and Theorem 16, the solution becomes as follows:
(85)y_(x;r)=Φ_(0,r)Eβ,β[λxβ]+∫0x‍(x-t)β-1Eβ,β[λ(x-t)β](t+1)dt,0≤r≤1,y¯(x;r)=Φ¯(0,r)Eβ,β[λxβ]+∫0x‍(x-t)β-1Eβ,β[λ(x-t)β](t+1)dt,0≤r≤1.
For more details see Figures 5 and 6.Figure 5
Solution of Example5, β=1.Figure 6
Solution of Example5, β=0.5.


## 6. Conclusion

In this paper, the fractional Green’s functions for the solutions of FFDEs by using the fuzzy Laplace transforms under Caputo’s differentiability were applied. Some particular cases were also presented, involving some fuzzy versions of the harmonic oscillator. To this end, some needed concepts are given like the Caputo derivatives, fuzzy Laplace transforms, and their inverses, and so forth. For future research, we extend the boundary value problem for FDE [38] to the fractional case.



---

*Source: 178378-2013-11-21.xml*