Existence, uniqueness, and approximate solutions for the general nonlinear distributed-order fractional differential equations in a Banach space

The purpose of this paper is to provide sufficient conditions for the local and global existence of solutions for the general nonlinear distributed-order fractional differential equations in the time domain. Also, we provide sufficient conditions for the uniqueness of the solutions. Furthermore, we use operational matrices for the fractional integral operator of the second kind Chebyshev wavelets and shifted fractional-order Jacobi polynomials via Gauss–Legendre quadrature formula and collocation methods to reduce the proposed equations into systems of nonlinear equations. Also, error bounds and convergence of the presented methods are investigated. In addition, the presented methods are implemented for two test problems and some famous distributed-order models, such as the model that describes the motion of the oscillator, the distributed-order fractional relaxation equation, and the Bagley–Torvik equation, to demonstrate the desired efficiency and accuracy of the proposed approaches. Comparisons between the methods proposed in this paper and the existing methods are given, which show that our numerical schemes exhibit better performances than the existing ones.

In [5], Atanacković et al. studied the existence and uniqueness of solutions for DOFDEs of the form in L 1 loc (R) ∩ C 1 ([0, ∞)). Such equations arise in distributed derivative models of system identification theory and visco-elasticity.
As the realm of DOFDEs describing the real-life response of physical systems grows, the demand for numerical solutions to analyze the behavior of these equations becomes more pronounced in order to overcome the mathematical complexity of analytical solutions. Therefore, the development of effective and easy-to-use numerical schemes for solving such equations acquires an increasing interest. While several numerical techniques have been proposed to solve many different problems (see, for instance, [1-3, 13, 15, 25, 33-35, 43, 45, 48] and the references therein), there have been few research studies that developed numerical methods to solve general DOFDEs (see [19,21,36,38,47,50,53]). The development, however, for efficient numerical methods to solve DOFDEs is still an important issue [21].
The aim of this paper is to provide sufficient conditions for the existence and uniqueness of solutions for Eqs. (1) and (2). Also, we are going to approximate solutions for the mentioned equations with high precision. To do this, in Sect. 2, we present a review of fractional calculus, an introduction of the second kind Chebyshev wavelets (SKCWs), shifted fractional-order Jacobi polynomials (SFOJPs), function approximations, and operational matrices for the Riemann-Liouville fractional integral operator. In Sect. 3, we provide sufficient conditions for the existence and uniqueness of solutions for general DOFDEs. In Sect. 4, by using operational matrices, mentioned in Sect. 2, we approximate the solution of Eqs. (1) and (2). In Sect. 5, we obtain the error bounds for the approximations. In Sect. 6, we solve two test problems and some famous distributed-order models, such as the model that describes the motion of the oscillator, the distributed-order fractional relaxation equation, and the Bagley-Torvik equation, to show that our approaches will increase the accuracy of the methods used for such operational matrices. Finally, a conclusion is given in Sect. 7.

Fractional calculus
The Riemann-Liouville integral of fractional order ι > 0 is defined as follows: where * and (·) are the convolution product and gamma function, respectively.

Definition 2.2 ([41])
The Caputo derivative of fractional order ι > 0 is defined as follows: with the following properties:

SKCWs and function approximation
SKCWs are as follows: where i = 1, 2, . . . , 2 k , j = 0, 1, . . . , M -1. Here, T * j (t) is the shifted Chebyshev polynomial of the second kind of degree j ≥ 0, defined on the interval [0, 1] by can be expanded by SKCWs as follows: wherê and ω(t) = w( 2 k t f ti + 1). We truncate the infinite series given in Eq. (6), and then we approximate a function f (t) in the following form: wherê

SFOJPs and function approximation
SFOJPs of order i are defined on the interval [0, t f ] by the following formula [22]: where θ , ϑ ∈ R and 0 < λ < 1.

Operational matrices of the Riemann-Liouville fractional integral operator
Following [23,51], we can obtain the operational matrix of the Riemann-Liouville fractional integral operator based on SKCWs for t ∈ [0, t f ) in the following theorem.

Theorem 2.1 Let (t) be the vector of SKCWs. Then
where Also, is the operational matrix of fractional integral operator for the block-pulse functions, where According to our previous work [37], we have the following theorem.

Theorem 2.2 Let (t) be the vector of SFOJPs. Then
where I t f ,q is the operational matrix of the Riemann-Liouville integral operator of fractional order q with entries Here, B(., .) is a beta function and k, l = 0, 1, . . . , N .

Existence and uniqueness of solutions
In the following theorem, by using Schauder's fixed point theorem [57], we prove the local existence of solutions for general DOFDEs in a Banach space.
Then general DOFDEs have at least one solution on .
It is clear that the set 0 is convex, closed, and bounded.
Define the operator for any f ∈ 0 . Clearly, we have Therefore, Tf ≤ b, and we can deduce T( 0 ) ⊂ 0 . Furthermore, for any t 1 , t 2 ∈ 0 such that t 2 > t 1 , we have Since G 1 is Lipschitz with the constant ς , we get

Now we can write
Note that if t 2 → t 1 , then the right-hand side of (10) tends to zero. Therefore, T : 0 → 0 is equicontinuous, and consequently, from the Arzela-Ascoli theorem [12], the closure of where v ∈ 0 . We need to show that T is continuous. Clearly, we have Suppose that, for any ε > 0, there exists δ > 0 such that |f (t)v(t)| < δ. Assume that assumptions (C1)-(C4) hold, therefore and the proof is completed. Now, by using Tychonoff 's fixed point theorem [57], we are going to discuss a global existence result for general DOFDEs.
has a solution u(t) which exists for every t ≥ 0, where κ 0 ≤ |κ 0 | and |κ i | ≤ κ i for i = 1, 2, . . . , r. Also, then for every Proof Let V be a real space of all continuous functions from (0, ∞) into R n . The topology on V being that induced by the family of pseudo-norms Under this topology, V is a linear space, locally convex and complete. Now consider where u(t) is a solution of (11). Clearly, in the topology of V, V 0 is bounded, convex, and closed. Consider (11) whose fixed point corresponds to a solution of (1). Evidently, in the topology of V, the map T is compact. Hence, in view of the boundedness of V 0 , the closure of T(V 0 ) is compact.
Considering assumptions (D1)-(D4) yields Since u(t) is a solution of (11), using the definition of V 0 gives |Tf (t)| ≤ u(t). Thus, T(V 0 ) ⊂ V 0 , and from Tychonoff 's fixed point theorem [57], T has a fixed point in V 0 . Therefore, the proof of this theorem is completed.
Now we are going to prove the uniqueness of the solution for general DOFDEs.
Proof Let Then we have for any t ∈ and f , f 1 ∈ C( , R n ). Therefore, Consequently, it has a unique fixed point, and therefore the general DOFDE has a unique solution f ∈ C( , R n ).

The methods of solution 4.1 Explanation of the SKCWs method
In this section, without loss of generality we suppose that β ≥ α r . Now we approximate C 0 D β t f (t) by the SKCWs as follows: By applying Eqs. (3) and (9), one can obtain Now we take the operator C 0 D q t of Eq. (12). So, we have Note that from properties (4) and (5), C 0 D q t (t j ) can be determined. Also, we take the operator C 0 D α i t , i = 1, 2, . . . , r, of Eq. (12). Using the obtained results in (1) gives Using the Gauss-Legendre formula and collocating the obtained equation at where w n and τ n are weights and nods of Gauss-Legendre quadrature rule [9], respectively. By using the "fsolve" command of Maple 2018, we solve the arising system, and then we determineF. Finally, from Eq. (7), an approximate solution for Eqs. (1) and (2) can be obtained.

Explanation of the SFOJPs method
Now similar to Sect. 4.1, by using SFOJPs, we convert Eqs. (1) and (2) to a system of equations as follows: where t m , m = 0, 1, . . . , N , are roots of SFOJPs. Also, w n and τ n are weights and nods of Gauss-Legendre quadrature rule [9], respectively. By the "fsolve" command of Maple 2018, we solve the above system, and then the unknown vectorF can be determined. Finally, from Eq. (8), we obtain an approximate solution for Eqs. (1) and (2).

Error bounds for the SKCWs method
In this subsection, we present error bounds for the SKCWs method. To do this, we define which are inner product and norm on the space L 2 ω ([0, t f )), respectively. Now we recall the following theorems from our previous work [47].
converges uniformly to f (t) and the coefficients in (6) explicitly satisfy Then we have

Error bounds for the SFOJPs method
Here, we discuss error bounds for the SFOJPs method. To do this, first, we define the following inner product and norm on the weighted space L 2 be the fractional-polynomial space of finite dimension. where Proof Since f N (t) is the best approximation to f (t) from N , defined in (13), we have Considering the generalized Taylors formula Taking L 2 w (λ,θ ,ϑ) t f -norm in both sides of inequality (15) leads to By substituting the above relation into (17) and from (16), we can write By taking the square root of (19), we obtain inequality (14).

Illustrative examples
In this section, we present five problems which are tested by Maple 2018. Also, we obtain the absolute errors by Note that, in all tables,m denotes the numbers of bases.
Example 6.1 Consider the following nonlinear DOFDE: , with the exact solution f (t) = t 3 √ t. In the above problem, the distributed-order term is discretized with the seven-point Gauss-Legendre quadrature rule. The numerical results at the interval [0, 1), obtained by the SKCWs method, are reported in Table 1 by selectingm = 2 k M = 6, 12. Graphs of the exact and approximate solutions and also absolute errors, obtained by the mentioned method with k = 2, M = 3, are plotted in Fig. 1. This figure and Table 1 illustrate the efficiency and accuracy of the method. Example 6.2 Consider the DOFDE in the following form [14,19,21,39]: ln t , with the exact solution f (t) = t 5 .    In the above problem, we discretize the distributed-order term with the eight-point Gauss-Legendre quadrature rule. In Table 2, we report the absolute errors at the interval [0, 1) for the SKCWs method by selectingm = N + 1 = 9, 12 with λ = 1 2 , θ = 3 2 , ϑ = 3; for the SFOJPs method by selectingm = 2 k M = 16, 20; for the hybrid of Legendre polynomials and block-pulse functions (HLBP) method [39] by selectingm = NM = 24; and for the block-pulse functions (BPFs) method [39] by selectingm = N = 32. Graphs of the exact and approximate solutions and also absolute errors, obtained by the SKCWs and SFOJPs methods with k = 1, M = 10, N = 11, λ = 1 2 , θ = 3 2 , ϑ = 3, are plotted in Fig. 2. This figure and Table 2 show the efficiency and accuracy of the new methods in comparison with the other methods reported in [39]. Example 6.3 Consider the DOFDE in the following form [19-21, 38, 50, 53]: The above equations describe the motion of the oscillator, where γ , a, b are constants; ω is the eigen frequency of the undamped system; g(t) is the external forcing function; f (t)  and σ (t) are the displacement and the dissipation force, respectively. In this problem, the forced vibrations of the distributed-order oscillator subjected to the harmonic excitation g(t) = g 0 sin( t) are studied. The solution of this problem is obtained with g 0 = 1, = 1.2ω, ω = 3, and γ = 1. If a = b, the solution is identical to the elastic with ω el = √ 1 + ω 2 = √ 10, and the exact solution is In Table 3, we report the absolute errors at the interval [0, 10) for the SFOJPs method by selectingm = N + 1 = 39, 49 with λ = 1 2 , θ = ϑ = 0; for the second kind Chebyshev wavelets (SKCWs) method [47] by selectingm = N(M + 1) = 52; and for the Müntz-Legendre wavelets (MLWs) method [36] by selectingm = N(M + 1) = 60. We emphasize that the results reported in [36,47] have been compared with the results reported in [19,38,39,53], and it was concluded that the methods in [36,47] are more accurate than the other methods. Therefore, here, we just compare our method with the methods of [36,47]. Graphs of the exact and approximate solutions and also absolute errors, obtained by the SKCWs and SFOJPs with k = 4, M = 5, N = 49, λ = 1 2 , θ = ϑ = 0, are plotted in Fig. 3. This figure and Table 3 show the efficiency and accuracy of the SFOJPs method in comparison with the methods reported in [36,47].   (s) = ln(s)(6s + 6) -12s + 12 (ln(s)) 3 .
In the above problem, we discretize the distributed-order term with the three-point Gauss-Legendre quadrature rule. In Table 4, we report the absolute errors at the interval [0, 10) for the SFOJPs method by selectingm = N + 1 = 4, 9 with λ = 1 2 , θ = 3 2 , ϑ = 3; and for the Müntz-Legendre wavelets (MLWs) method [36] withm = 2 k-1 M = 16, 32. Figure 4 shows that by using 8 number of SFOJPs the obtained results are better than the results of [50] that obtained by using 1000 BPFs for solving this problem. Example 6.5 Consider the following Bagley-Torvik equation [7,8], where the damping term is expressed in terms of distributed-order derivatives [46]: This equation is called fractional oscillator equation, when the order of damping term is constant.
In the above problem, the distributed-order term is discretized with the three-point

Conclusion
In this research paper, based on Schauder's and Tychonoff 's fixed point theorems, sufficient conditions for the local and global existence of solutions were provided for general DOFDEs. Also, sufficient conditions were provided for the uniqueness of the solutions. Furthermore, we proposed new methods to solve DOFDEs of the general form in the time domain. By using these methods, the mentioned equations were reduced to systems of algebraic equations. We solved these systems by using the "fsolve" command of Maple 2018. The error bounds of the methods have been discussed. In addition, the presented methods were implemented for two test problems and some famous distributed-order models, such as the model that describes the motion of the oscillator, the distributed-order fractional relaxation equation, and the Bagley-Torvik equation. It showed that by applying the SKCWs and SFOJPs methods, the obtained results are better than the other existing methods. We deduce that the proposed methods are efficient numerical tools for solving DOFDEs.