Approximate solutions and Hyers–Ulam stability for a system of the coupled fractional thermostat control model via the generalized differential transform

In this paper, we consider a new coupled system of fractional boundary value problems based on the thermostat control model. With the help of fixed point theory, we investigate the existence criterion of the solution to the given coupled system. This property is proved by using the Krasnoselskii’s fixed point theorem and its uniqueness is proved via the Banach principle for contractions. Further, the Hyers–Ulam stability of solutions is investigated. Then, we find the approximate solution of the coupled fractional thermostat control system by using a numerical technique called the generalized differential transform method. To show the consistency and validity of our theoretical results, we provide two illustrative examples.


Introduction
Fractional calculus is the most important field of applied mathematics that extends all existing integer-order operators to arbitrary-order ones. This importance is due to the high accuracy of fractional operators in modeling of many real-world phenomena in the context of different fractional boundary value problems. Some instances regarding applications of fractional differential equations (FDEqs) can be seen in various fields of science such as aerodynamics, physics, bioengineering, image processing, biochemistry viscoelasticity, biophysics, electrochemistry, mathematical biology, etc. (see [1][2][3]). In the last few decades, an extensive area of studies in relation to the existence/nonexistence theory of solutions to the aforementioned FDEqs has got much attention from mathematicians. Diverse results and theoretical findings can be followed in the literature regarding the existence theory for different structures of FDEqs; for more analysis and review, see . Because of the modeling of most of the applied processes in the framework of coupled systems of arbitrary order FDEqs, many researchers have been focused on the establishment of the existence theory for such systems. In such a direction, plenty of research manuscripts can be observed in the literature, see [34][35][36][37][38]. Recently, Shah, Wang, Khalil, and Khan [39] devoted their focus to establishing some results on the following category of coupled systems of integral fractional boundary value problems (FBVPs) for FDEqs with movable boundary conditions (BCs) given as in which f , g ∈ C(I × [0, ∞)), 1 < σ , δ ≤ 2, and 0 < a, b < 1. More recently, Alrabaiah et al. [40] implemented a qualitative study of a nonlinear system of coupled integral pantograph delay FDEqs as follows: where CC D q,η j t 0 denotes the Caputo conformable derivatives of order η j ∈ (1, 2) with q ∈ (0, 1] for j = 1, 2, RC I q,θ t 0 is the RL-conformable integral of order θ > 0, c, d ∈ (t 0 , T), a 1 , a 2 , b 1 , b 2 , m 1 , m 2 ∈ R, 0 < λ < 1, and f j : [t 0 , T] × R × R → R are continuous for j = 1, 2.
Other aspects in relation to the stability analysis and numerical methods to FDEqs have been investigated by many mathematicians in recent years. Most of the time, finding exact solutions to nonlinear BVPs is challenging and sometimes it is a difficult task. Hence, these issues motivated researchers to search for the best approximate solutions for given nonlinear BVPs. To achieve such an aim, mathematicians utilized different techniques and procedures like homotopy methods [42], q-homotopy analysis transform method [43], decomposition methods [44,45], integral transform techniques [46], etc.
One of the strongest approaches to find numerical solutions to linear and nonlinear BVPs of FDEqs is the generalized differential transform method (GDT-method). In various papers, this transform has been applied to solve and analyze nonlinear BVPs of FDEqs for approximate solutions; see [47,48]. Note that there is no classical standard method to analyze nonlinear BVPs of FDEqs for getting solutions explicitly. This is due to the complexity of fractional calculus involved in the supposed BVPs. Accordingly, a reliable method is needed to search for approximate solutions in the context of series in relation to the given FBVPs. On the other hand, it is also applicable if, along with the approximate solutions, their stability is investigated. That's why analysis of the stability of existing FDEqs has received great attention. In other words, it gives a better result in different applications like optimization, economics, numerical analysis, physics, where obtaining the exact solution is a quite difficult task; see [49][50][51].
Amongst the important models of physio-electrical type, we choose and study a fractional model of thermostat control. A thermostat is a measurement tool that measures and regulates the temperature of an arbitrary physical system and takes actions subject to its temperature which is maintained near an appropriate and desired degree. This instrument is applied in any controlling devices and industrial systems, including building central heating, air conditioners, water heaters, ovens, refrigerators, car engines, and even medical incubators, which raise or decrease the temperature.
The study of the mathematical model of a thermostat control was implemented by Infante and Webb [52] in 2006, which is insulated at t = 0 via the controller at the time t = 1 and is given as where b ∈ I is a real constant and k > 0. The performance of the thermostat is interpreted by this 2nd-order mathematical model in such a way that at the moment t = b, the system's heat is added or discharged based on the temperature detected by the existing sensors. Infante et al. analyzed some results on the solution's existence for the suggested model by utilizing the fixed point index theory. In the next years, Nieto and Pimentel [53] transformed the above ordinary thermostat model to a new model of the fractional type as follows: x(t)) = 0 (∀t ∈ I := [0, 1]), and extended the obtained results in the paper published by Infante et al. to the more general and accurate findings. Note that in that model, c D γ is the Caputo derivative of order γ ∈ (1, 2]. In this paper, we focus our intention on some qualitative aspects of possible solutions for a system of the coupled fractional thermostat model. In more precise words, we consider the following construction of a coupled thermostat model inspired by the model (1) as follows:  in which σ ∈ (1, 2], δ ∈ (1, 2], 0 < σ -1 ≤ 1, 0 < δ -1 ≤ 1, b, c ∈ (0, 1), p, q > 0, and c D 1 = d dt . Along with these, the existing mappings K, M : I × R ≥0 → R ≥0 are continuous and c D ς denotes the derivation operator of order ς ∈ {1, σ , σ -1, δ, δ -1} in the sense of Caputo.
To prove the existence and uniqueness of solutions to the above coupled fractional sys-

Primitive notions
As the notions of the Riemann-Liouville integral and Caputo derivative have a key task in our research, we provide them at this moment.
if the integral's value is finite.
For a continuous mapping x : R ≥0 → R, the σ th Riemann-Liouville derivation operator is defined by if the integral's value is finite.
For an absolutely continuous mapping x on R ≥0 , the σ th Caputo derivation operator is introduced as if the integral's value is finite.

Results of existence
We here compute and derive the equivalent system of coupled integral equations corresponding to the coupled thermostat control model (3) and, further, we present required criteria for the existence of solutions to the mentioned FBVP system (3).
is given by the integral equation Proof If x * satisfies the fractional thermostat linear equation (4), then c D σ x * (t) + g(t) = 0 holds. By integrating the latter equation and by virtue of 1 < σ ≤ 2, we get in which it is necessary that we find coefficients m 0 , m 1 ∈ R. On the other hand, the properties of the Caputo derivative yield and Using the condition c D 1 x(0) = 0 and (7), we figure out that m 1 = 0. Moreover, the equations (6), (8), and the condition and thus we reach At last, we substitute the obtained coefficients m 0 and m 1 into (6), and so (6) becomes thus, this argument is finished.
With due attention to the above proposition, we can formulate here an equivalent structure of coupled integral equations to the investigated system of coupled BVPs for the fractional thermostat model (3) in the next proposition.
Then an equivalent configuration of the system of coupled BVPs of the fractional thermostat model is provided by the coupled integral equations for any t ∈ I.
Thanks to the above proposition and for the sake of our subsequent arguments, we aim to introduce two operators g 1 : B → B and g 2 : B → B which take forms and So, we find the system of coupled operator equations in a format as follows: Hence, on the product space, we define g : In the sequel, we deal with the fixed points of g . For this we exploit Banach's and Krasnoselskii's fixed point theorems.

Theorem 3.3 ([55]
) Let X be a nonempty complete metric space and Ψ : X → X a contraction mapping. Then we can find a unique point z ∈ X with Ψ (z) = z.
(ii) g 1 is a contraction mapping, (iii) g 2 is compact and continuous. Then there exists z ∈ Σ such that z = g 1 z + g 2 z.
Now, in relation to the proofs of the main theorems, the following hypotheses are needed: (H1) K, M : I × R ≥0 → R ≥0 are continuous, for t ∈ I and x, x, y, y ∈ R.
for t ∈ I and x, x ∈ R. (H3) There exists a real positive constant C M which satisfies for t ∈ I and y, y ∈ R. The following notations are used for convenience of calculations: Now, everything is ready to start the establishment of the first main theorem. Proof First, let us choose where θ 1 = max t∈I |K(t, 0)| and θ 2 = max t∈I |M(t, 0)|, and set which is bounded, closed, and convex. Let (x, y) ∈ Ξ , then, in view of (H1) and (H2), we can write Therefore, In a similar way, we show that Thus, Hence, from (14) and (15), it follows that Hence g (Ξ ) ⊆ Ξ . Now, we show that g is a contraction operator. Let (x, y), ( x, y) ∈ Ξ and t ∈ I, then we have Then, (16) leads to By the same calculation techniques used to get (16), we find Now, by exploiting (17) and (18), together with the assumptions C K δ 1 < 1 and C M δ 2 < 1, we arrive at Hence, g is a contraction operator. Then, from Theorem 3.3, our system of coupled BVPs for the fractional thermostat model (3) admits a uniqu solution, and this ends the proof.
In addition, we suppose that the following hypothesis holds: (H4) There are four real positive constants for each t ∈ I and any x, y ∈ B.
Theorem 3.6 Let the hypotheses (H1)-(H4) be satisfied. If, in addition, the conditions hold, then the considered system of coupled BVPs for the fractional thermostat model (3) admits a solution.
Proof Since K and M are continuous functions, then we have the same for the operator g . Let D be a bounded subset of Ξ ∈ B × B. Then, based on the hypothesis (H4), we have, for all (x, y) ∈ D, which leads to A similar argument gives us In view of (21) and (22), the boundedness of (D) is ensured. Now, we show that is an equicontinuous operator. For this, let τ 1 , τ 2 ∈ I, where τ 1 < τ 2 and (x, y) ∈ B × B, then we have By the same arguments used to get (23), one can find the following result: The inequalities (23) and (24) imply that | 1 y(τ 1 )-1 y(τ 2 )| → 0 and | 2 x(τ 1 )-2 x(τ 2 )| → 0, when τ 1 → τ 2 . Then, by using the Arzela-Ascoli theorem, we conclude that the operator is continuous and compact. Now, to complete checking the hypotheses of Theorem 3.4, it remains to show that is a contraction operator. So, if we take (y, y) ∈ B × B, we get Consequently, With the same calculation method, we arrive at the following estimate: Hence, from (25) and (26), together with the assumption (20), we have that is a contraction operator. Now, all the assumptions of Theorem 3.4 are fulfilled. So, this ensures that g admits at least one fixed point which is also a solution of the considered system of coupled BVPs for the fractional thermostat model (3).

Results for Hyers-Ulam stability
Fractional differential equations have been extensively studied from different angles. Among these, stability analysis in the Hyers-Ulam sense is an important aspect that gained proper attention from researchers [56][57][58]. Based on the fundamental definition of Hyers-Ulam stability of a system, the notion was later modified to more general types, and their results were successfully applied to various problems [59][60][61]. In this section, we will adopt a number of sufficient conditions to review the Hyers-Ulam-type stability results for the considered system of coupled BVPs for the fractional thermostat model (3).
is Hyers-Ulam stable if for the given inequality (∀t ∈ I) ∃ω g > 0 so that for any y ∈ C(I, R) satisfying equation (27), one can find the solution h ∈ C(I, R) of (27) uniquely, provided that ∀t ∈ I, is Hyers-Ulam-stable if for the following systems of inequalities: one can find two positive constants ω g 1 , ω g 2 so that for each (x, y) satisfying (28), a solution ( h, h) of above system (28) exists uniquely, provided that ∀t ∈ I, For establishing the formal theorems on the Hyers-Ulam stability for the considered system of coupled BVPs of the fractional thermostat model (3), we indicate the following conditions first.

Lemma 4.4
Assume that (x, y) ∈ (C(I, R)) 2 is a solution of the inequality system (29). Then we have the following system of inequalities: , and g 1 , g 2 are the operators defined by (11) and (12), respectively.
Proof In accordance with (2) in Remark 4.3, we have ∀t ∈ I, Thanks to the Proposition 3.2, the solution of problem (30) can be reformulated immediately as Since t ∈ I := [0, 1], from (31) we have, on the one hand, which means that On the other hand, with the same arguments, we obtain that is, and this concludes the proof.

Theorem 4.5 Assume that the hypotheses (H2), (H3) hold. If
where δ 1 δ 2 C K C M < 1, then the solution of the system of coupled BVPs of the fractional thermostat model (9) is Hyers-Ulam-stable.
Proof Taking any solution (x, y) ∈ (C(I, R)) 2 of the system of inequalities defined by and the unique solution ( h, h) belonging to (C(I, R)) 2 of the following system: and from Proposition 3.2, together with (10), the solution of (32) can be represented by In view of (33), we can write Exploiting Lemma 4.4, with some computations, we obtain With the same calculation, we get Since δ 1 δ 2 C K C M < 1, from (35), (36), and (37), we get, after simple computations, the following estimates: Hence the solution of the system of coupled BVPs of the fractional thermostat model (9) is Hyers-Ulam-stable.

GDT-method for approximation of solutions
In many situations, the explicit solution of boundary problems for classical nonlinear differential equations is difficult and sometimes impossible. Then, it would be quite difficult or impossible to find the exact solutions of FBVPs with FDEqs. Therefore, it has become necessary to think about new methods to solve such problems. In this work we are interested in a numerical method, called the differential transform method, which was presented by Zhou in [62], which was extended to its generalized form by Odibat and Momani in [63] and named the GDT-method. It is an iterative method which gives us analytical solutions in the form of the Taylor series expansion to fractional differential equations with boundary or initial conditions. Now, we will apply the GDT-method to find an approximate solution to our system of coupled BVPs of the fractional thermostat model (3). Note that this technique is based on the generalized Taylor's formula. The GDT of the sth derivative of two given functions x(t) and y(t) of one variable is defined by (see [64]) where ( c D ) s = c D · c D · · · c D s-times ; and their inverses are defined as Then the approximate solution of our problem (3) is written as finite series of the analytical polynomial forms, where is named the order of the differential transformation and must be taken such that l = 1, m = σ , and n = δ with l, m, n ∈ N, and x (s), y (s) are the GDT of x(t) and y(t), respectively, expressed by where K (s, y (s)) and M (s, x (s)) denote the th-order differential transformations of K(s, y(s)) and M(s, x(s)), respectively. Since c D 1 x(0) = 0 and c D 1 y(0) = 0, their differential transforms give x (l) = 0, y (l) = 0, and x (s) = 0, for all s that satisfy 0 < s < 1 or l < s < m, and also y (s) = 0, for all s that satisfy 0 < s < 1 or l < s < n, and finally, x (0) = c, y (0) = c, where c and c are two undetermined real constants which can be evaluated using the second initial condition of problem (3).
Thanks to the recursive relationship (40), the solution (x(t), y(t)) of the system of coupled BVPs of the fractional thermostat model (3) can be written as the following finite series: where x c, c (s), y c, c (s) are coefficients depending on c and c which can be determined using the second initial condition of problem (3). From (41), we can write Therefore (42) gives us the following system of equations: Now, we solve (43) with respect to c and c and replace their values in (41), and ultimately we find our desired approximate solution of the system of coupled BVPs of the fractional thermostat model (3).

Conclusion
In this work some existence and uniqueness results have been obtained by using the Banach and Krasnoselskii's fixed point theorems. Also some necessary conditions for Hyers-Ulam stability of solutions of a given system of coupled BVPs of the fractional thermostat model (3) have been discussed. In addition, some approximate solutions by the GDTmethod have been found, and some illustrative examples have been presented as applications of the GDT-method on some problems as in (3). These examples show that this numerical method can give us a good and accurate approximate solution for nonlinear FBVP of FDEqs. In the next works, we are going to implement such an analysis (theoretical and numerical) for other nonlinear fractional systems of FDEqs arising in different applied models featuring generalized fractional operators.