Next Article in Journal
Inverse Problem for a Partial Differential Equation with Gerasimov–Caputo-Type Operator and Degeneration
Next Article in Special Issue
Adsorption on Fractal Surfaces: A Non Linear Modeling Approach of a Fractional Behavior
Previous Article in Journal
Numerical Simulation for the Treatment of Nonlinear Predator–Prey Equations by Using the Finite Element Optimization Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Solutions of Bernoulli Equations in the Fractional Setting

Dipartimento di Scienze di Base e Applicate per l’Ingegneria, Sapienza Università di Roma, 00161 Rome, Italy
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Submission received: 5 May 2021 / Revised: 9 June 2021 / Accepted: 13 June 2021 / Published: 17 June 2021
(This article belongs to the Special Issue Fractional Behavior in Nature 2021)

Abstract

:
We present a general series representation formula for the local solution of the Bernoulli equation with Caputo fractional derivatives. We then focus on a generalization of the fractional logistic equation and present some related numerical simulations.

1. Introduction

Interest in time fractional evolutive systems has progressively grown in recent years. Models arising with aspects related to non-local behavior need to be studied in the fractional setting; see [1,2] for an overview and [2,3] for fractional growth models for social and biological dynamics. The fractional derivatives are indeed non-local operators, that is, convolution-type operators. In the applied sciences, the main interest in fractional models is due to the fact that such models introduce the so-called memory effect. This effect is mainly justified by the non-locality of the time-fractional derivative, and it seems to be relevant in the characterization of many applied models. A second reading is given in terms of the delaying effect. Indeed, the time-fractional derivative introduces a different clock for the underlying model, as in the case of the relaxation equation. When the order of the fractional derivative is one, the underlying model emerges.
Here, we locally solve the following Cauchy system, involving a fractional Bernoulli equation of the form
C D t β u + a 0 u = a 1 u p + 1 , u ( 0 ) = u 0
where p N , with p 1 , a 0 and a 1 are real numbers and C D t β denotes the Caputo derivative. If β = 1 then C D t β u = u and (1) is the Bernoulli equation, studied by Jacob Bernoulli (1695). We remark that Bernoulli equations (with β = 1 ) arise in non-linear models of production and capital accumulation, in particular when polynomial production functions are considered, see ([4] Chapter 6.3). As an example of an application to mathematical bio-medicine, we refer to ([5] Chapter 4), where the parameter p is associated with the exponential growth rate of a homogeneous avascular tumor.
As a particular case, the exact solution in the case a 0 = a 1 = 1 and β = 1 is given by
u ( t ) = e t c 0 + e p t p ; where c 0 : = u 0 p 1 .
Moving to the fractional setting, we have similar approaches that cannot be followed. As it is well known, the solution of the fractional logistic equation—corresponding to p = 1 and a 0 = a 1 = 1 in (1)—was an open problem, and in [6], the first and the third author were able to solve the fractional logistic equation by series representation, giving a detailed formula involving Euler numbers for u 0 = 1 / 2 . The use of a logistic equation to model population dynamics dates back to the mid-19th century when Pierre François Verhulst devoted a series of papers [7,8,9] to the description of a model of population growth. Since then, the logistic equation is largely used in bio-mathematical modeling, for instance, to describe cancer growth models and epidemic models. In the particular case of social dynamics, the need to model memory effects and non-local behavior lead to a decline of the above models in a fractional setting. At this early stage, the problem of addressing non-linear dynamics and series expansions proved, at least in a local setting, to be a viable approach. This methodology was applied to SIS epidemic models in [10] (a compartmental epidemic model in which the population is split into susceptible and infective individuals) and also further investigated in [11]. The present study extends the result in [6] to general initial data and to Bernoulli equations of general degree p + 1 : we present a recursive formula for the coefficients of the solutions and explicit closed formulas for the first terms. Note that the relation with Euler numbers for general initial data, even in the logistic case p = 1 , appears to be lost, but the general recursive formula preserves its structure, based on generalized binomial coefficients that were introduced in [6] and further investigated in [12]. Then the proposed method is applied to the particular case a 0 = a 1 = ± 1 , related to the fractional logistic equation, and we present a qualitative analysis of the solutions based on numerical simulations—see Figure 1.

Preliminaries on Fractional Calculus

Let us consider the set A C ( [ a , b ] ) of continuous functions with a derivative in L 1 ( [ a , b ] ) . Thus, v A C ( [ a , b ] ) is continuous and such that v = g L 1 ( [ a , b ] ) , that is, v has the representation
v ( t ) = v ( a ) + a t g ( s ) d s , t [ a , b ] .
We notice that the space A C ( [ a , b ] ) coincides with the Sobolev space
W 1 , 1 ( [ a , b ] ) = { v L 1 ( [ a , b ] ) : v L 1 ( [ a , b ] ) }
endowed with the norm
v W 1 , 1 = v L 1 + v L 1 .
For v A C ( [ a , b ] ) and β ( 0 , 1 ) , we introduce the Riemann–Liouville derivative of v,
D t β v ( t ) : = 1 Γ ( 1 β ) d d t a t v ( s ) ( t s ) β d s , t [ a , b ]
and the Caputo derivative of v,
C D t β v ( t ) : = 1 Γ ( 1 β ) a t v ( s ) ( t s ) β d s , t [ a , b ] .
Further on, we use the following relation between derivatives
D t β v ( t ) = 1 Γ ( 1 β ) v ( a ) ( t a ) β + C D t β v ( t ) .
The relation (5), together with the existence of the derivatives (3) and (4), hold a.e. on [ a , b ] and D t β v L q ( [ a , b ] ) with 1 q 1 / β (see for example ([13] page 28)).
In this study, we consider the fractional equations on [ 0 , b ) . Let us emphasize that if v ( 0 ) = 0 , then Formula (5) gives the equivalence
D t β v = C D t β v .

2. Fractional Bernoulli Equations

Let us introduce
u ( t ) = n 0 c n ( 1 ) t β n Γ ( β n + 1 ) , t ( 0 , r )
where r > 0 is the radius of convergence and c n ( 1 ) are real coefficients.
Next, we introduce the generalized binomial coefficient
n k β = Γ ( n β + 1 ) Γ ( k β + 1 ) Γ ( ( n k ) β + 1 ) , 0 k n , k , n N , β ( 0 , 1 ) .
Theorem 1.
The unique continuous solution on [ 0 , b ) [ 0 , r ) to
C D t β u + a 0 u = a 1 u p + 1 , u ( 0 ) = u 0 , p N , a 0 , a 1 R , β ( 0 , 1 )
has the series representation (6) on [ 0 , b ) where
c 0 ( 1 ) = u 0 c n + 1 ( 1 ) = a 0 c n ( 1 ) + a 1 k = 0 n n k β c k ( p ) c n k ( 1 ) f o r n 0 .
Proof. 
We consider
C D t β u + a 0 u = a 1 u p + 1 , u ( 0 ) = u 0
where C D t β denotes the Caputo derivative. Let us denote the power of u as follows
u ( t ) = n 0 c n ( 1 ) t β n Γ ( β n + 1 ) , u 2 ( t ) = n 0 c n ( 2 ) t β n Γ ( β n + 1 ) ,
u 3 ( t ) = n 0 k = 0 n n k β c k ( 2 ) c n k ( 1 ) t β n Γ ( β n + 1 ) = n 0 c n ( 3 ) t β n Γ ( β n + 1 ) ,
where
c k ( 2 ) = s = 0 k k s β c s ( 1 ) c k s ( 1 ) c n ( 3 ) = k = 0 n n k β c k ( 2 ) c n k ( 1 )
and, by further iterations,
u p + 1 ( t ) = n 0 c n ( p + 1 ) t β n Γ ( β n + 1 )
where
c n ( p + 1 ) = k = 0 n n k β c k ( p ) c n k ( 1 ) .
Note that for p N
c 0 ( p ) = c 0 ( 1 ) p .
Therefore,
a 0 u a 1 u p + 1 = u ( a 0 a 1 u p ) = i c i ( 1 ) t β i Γ ( β i + 1 ) · a 0 a 1 j c j ( p ) t β j Γ ( β j + 1 ) = a 0 i c i ( 1 ) t β i Γ ( β i + 1 ) a 1 i c i ( 1 ) t β i Γ ( β i + 1 ) j c j ( p ) t β j Γ ( β j + 1 ) = a 0 j c j ( 1 ) t β j Γ ( β j + 1 ) a 1 n c n ( p + 1 ) t β n Γ ( β n + 1 ) = n a 0 c n ( 1 ) a 1 c n ( p + 1 ) t β n Γ ( β n + 1 ) = n a 0 c n ( 1 ) a 1 k = 0 n n k β c k ( p ) c n k ( 1 ) t β n Γ ( β n + 1 ) .
On the other hand, from (5), we have that
C D t β u ( t ) = u 0 t β Γ ( 1 β ) + D t β n 0 c n ( 1 ) t β n Γ ( β n + 1 )
where, after some calculation, from (3),
D t β t β n Γ ( β n + 1 ) = t β ( n 1 ) Γ ( β ( n 1 ) + 1 ) , t [ 0 , b ) .
Thus, we obtain
C D t β u ( t ) = u 0 t β Γ ( 1 β ) + n 0 c n ( 1 ) t β ( n 1 ) Γ ( β ( n 1 ) + 1 ) = u 0 t β Γ ( 1 β ) + c 0 ( 1 ) t β Γ ( β + 1 ) + n 1 c n ( 1 ) t β ( n 1 ) Γ ( β ( n 1 ) + 1 ) .
From the fact that u 0 = c 0 ( 1 ) by construction, we write
C D t β u ( t ) = n 0 c n + 1 ( 1 ) t β n Γ ( β n + 1 ) , t [ 0 , b ] .
Then the solution to
C D t β u = a 0 u + a 1 u p + 1 = u ( a 0 a 1 u p )
can be written in terms of the coefficients c n ( 1 ) , n N given by
c n + 1 ( 1 ) = a 0 c n ( 1 ) + a 1 c n ( p + 1 ) = c n ( 1 ) a 0 a 1 c 0 ( p ) + a 1 k = 1 n n k β c k ( p ) c n k ( 1 ) , = a 0 c n ( 1 ) + a 1 k = 0 n n k β c k ( p ) c n k ( 1 ) for n 0 .

Some Closed Formulas

The first few elements of the sequence c n ( p ) , p N (over the index p) are
c n ( 1 ) , n N c n ( 2 ) = s = 0 n n s β c s ( 1 ) c n s ( 1 ) c n ( 3 ) = k = 0 n n k β s = 0 k k s β c s ( 1 ) c k s ( 1 ) c n k ( 1 ) .
We compute the first terms of c n ( 1 ) . Fix p 1 and note that c 0 ( 1 ) = u 0 and if n = 1 then
c 1 ( 1 ) = u 0 ( a 0 a 1 u 0 p ) .
To compute c 2 ( 1 ) we need c 1 ( h ) for h = 1 , , p . By (9) one can prove by induction that for h N
c 1 ( h ) = h u 0 h ( a 0 a 1 u 0 p ) .
Hence, we apply (7) and we get
c 2 ( 1 ) = c 1 ( 1 ) a 0 a 1 c 0 ( p ) + a 1 1 1 β c 1 ( p ) c 0 ( 1 )
from which we deduce by a direct computation
c 2 ( 1 ) = u 0 ( a 0 a 1 u 0 p ) ( a 0 a 1 ( p + 1 ) u 0 p ) .
To compute c 3 ( 1 ) , we also need c 2 ( h ) for h = 1 , , p . Using (9), we obtain
c 2 ( h ) = k = 0 2 2 k β c k ( h 1 ) c 2 k ( 1 )
and, by an inductive argument, using also the fact that c 0 ( p ) = u 0 p for all p 1 , one can prove the closed formula for all h 1
c 2 ( h ) = u 0 h ( a 0 a 1 u 0 p ) h ( a 0 a 1 ( p + 1 ) u 0 p ) + h ( h 1 ) 2 2 1 β ( a 0 a 1 u 0 p ) .
Then
c 3 ( 1 ) = a 0 c 2 ( 1 ) + a 1 c 2 ( p + 1 ) = a 0 u 0 ( a 0 a 1 u 0 p ) ( a 0 a 1 ( p + 1 ) u 0 p ) + a 1 u 0 p + 1 ( a 0 a 1 u 0 p ) ( ( p + 1 ) ( a 0 a 1 ( p + 1 ) u 0 p ) + p ( p + 1 ) 2 2 1 β ( a 0 a 1 u 0 p ) ) = u 0 ( a 0 a 1 u 0 p ) ( a 0 a 1 ( p + 1 ) u 0 p ) 2 + p ( p + 1 ) 2 2 1 β a 1 u 0 p + 1 ( a 0 a 1 u 0 p ) 2 = u 0 ( a 0 a 1 u 0 p ) ( ( a 0 a 1 ( p + 1 ) u 0 p ) 2 p ( p + 1 ) 2 2 1 β a 1 u 0 p ( a 0 a 1 u 0 p ) ) .

3. Fractional Logistic Equations

Here, we extend some of the results established in [6] for the fractional logistic equation with initial datum u 0 = 1 / 2 to the case of general initial data u 0 ( 0 , 1 ) . Applying the above method, if p = 1 and a 0 = a 1 = 1 , then the solution of
C D t β u = u u 2 ; u ( 0 ) = u 0
can be represented in series form
u ( t ) = n 0 c n ( 1 ) t β n Γ ( β n + 1 )
where
c 0 ( 1 ) = u 0 ; c n + 1 ( 1 ) = c n ( 1 ) k = 0 n n k β c k ( 1 ) c n k ( 1 ) .
Note that, as shown in [6], when u 0 = 1 / 2 , the above formula reduces to
c 0 ( 1 ) = 1 2 ; c 1 ( 1 ) = 1 4 ; c 2 n ( 1 ) = 0 ; c 2 n + 1 ( 1 ) = k = 0 n 2 n 2 k + 1 β c 2 k + 1 ( 1 ) c 2 ( n k ) 1 ( 1 ) for n 1 .
Keeping a 0 = a 1 = 1 and considering the equation with generic p, we have that the solution of
C D t β u = u u p + 1 ; u ( 0 ) = u 0
can be represented in series form
u ( t ) = n 0 c n ( p ) t β n Γ ( β n + 1 )
where
c 0 ( 1 ) = u 0 c 1 ( 1 ) = u 0 ( 1 u 0 p ) ; c n + 1 ( 1 ) = c n ( 1 ) k = 0 n n k β c k ( p ) c n k ( 1 ) .
Remark 1
(On the null coefficients in the general case). Note that if u 0 = 1 / ( p + 1 ) p then, in view of (11), c 2 ( 1 ) = 0 , in agreement with the case p = 1 . However, it is not possible to deduce that c 2 n ( 1 ) 0 as in the case p = 1 , because even for p = 2 , choosing u 0 = 1 / 3 one can numerically verify that c 2 ( 1 ) = 0 and c 4 ( 1 ) 0 . Then, one may look for some other generalization, for instance imposing c p + 1 ( 1 ) = 0 and guess whether c n ( p + 1 ) ( 1 ) = 0 for some n 1 . Furthermore, in this case, the answer is negative. Indeed, using the last equality (12) to solve the equation c 3 ( 1 ) = 0 with respect to the initial datum u 0 , we get by a direct computation that if u 0 = 12 + 6 π + 144 + 96 π 24 + 18 π and if p = 2 then c 3 ( 1 ) = 0 , but symbolic numerical computations yield c n ( 1 ) 0 for n = 3 , , 30 .
Consider now the case a 0 = a 1 = 1 . Then the solution of
C D t β u = u u 2 ; u ( 0 ) = u 0
can be represented as
u ( t ) = n 0 c ¯ n ( p ) t β n Γ ( β n + 1 )
where
c 0 ( 1 ) = u 0 c 1 ( 1 ) = u 0 ( 1 u 0 p ) ,
c n + 1 ( 1 ) = c n ( 1 ) + k = 0 n n k β c k ( p ) c n k ( 1 ) .
We compare the coefficients c n ( 1 ) and c ¯ n ( 1 ) . We have c ¯ 0 ( 1 ) = c 0 ( 1 ) , c ¯ 1 ( 1 ) = c 1 ( 1 ) and, by induction, c ¯ 1 ( p ) = c 1 ( p ) . Moreover, c ¯ 2 ( 1 ) = c 2 ( 1 ) . However, this symmetry breaks as soon as we consider c ¯ 2 ( 2 ) ; indeed, we have
c ¯ 2 ( 2 ) = c 2 ( 2 ) + 2 2 1 β ( c 1 ( 1 ) ) 2 .

Numerical Simulations

Our tests rely on a direct application of the above explicit formulas. Coefficients are calculated via Wolfram Mathematica software. In our tests, we focused on the logistic case a 0 = a 1 = 1 and on the case a 0 = a 1 = 1 . We computed the coefficients c n ( 1 ) using the recursive Formulas (9) and (10), and we approximated the solution of (1) with the partial sum
u ( t ) = n = 0 N c n ( 1 ) t β n Γ ( β n + 1 )
with N = 200 , the parameter N was tuned so that no appreciable difference can be noted with respect to higher-order approximations. The method was validated by comparison with the exact solutions of (1), which can be explicitly computed in the ordinary case β = 1 .
In Figure 2, we evaluated u ( t ) with fixed initial datum u 0 = 1 / 2 and p = 1 and several orders of fractional derivatives. We may note the expected damping effect of fractional derivation.
Figure 3 compares the solutions u ( t ) with different initial data, setting β = 1 / 2 and p = 1 . The resulting set of ordered curves suggests the local uniqueness of the solutions; however, the investigation of this is beyond the purpose of the present paper.
We then investigated higher degree fractional Euler equations, setting β = u 0 = 1 / 2 and letting p vary between 1 and 3, see Figure 4. At least near 0, from a qualitative point of view, the solutions display a similar behavior, and no intersections between the solutions are detected.
For the sake of comparison, we collected some of the above results in Figure 5, showing the combined effect of varying initial data u 0 and degrees p.
Finally we propose some numerical estimations for the radius of convergence of the series (8), by computing the sequence
r n : = Γ ( β n + 1 ) | c n ( 1 ) | 1 β n .
In Figure 6 and Figure 7, we plotted the first 300 terms of r n with varying degrees p and orders of derivation β . The asymptotic behavior of r n suggests an exponential increase for the series coefficients c n ( 1 ) in all the cases under exam. Furthermore, their comparison shows the radius of convergence r : = lim r n to be decreasing with respect to both the degree p (Figure 6) and order of derivation β (Figure 7). Finally, no substantial difference between the case a 0 = a 1 = 1 and the case a 0 = a 1 = 1 emerged.

4. Conclusions

Although the fractional differential calculus is as old as the ordinary one (the first introduction of a fractional derivative dates back to Leibniz in 1695), the complexity of the computations involved still prevents a full exploitation of its applicative potential. Here, we consider the Caputo fractional derivative. The presence of a singular kernel introduces a memory, and it gives a non-local character to the dynamics. These non-local aspects are realistic in the applications. By introducing fractional derivatives, we obtain a global effect, which may slow down the underlying dynamics. In particular, for non-linear equations, we still observe the classical delaying effect introduced by the fractional derivative for linear equations. This motivates our analysis. Further directions include estimations of the radius of convergence, global existence results, and possible generalization to other non-linear fractional dynamics.

Author Contributions

Conceptualization; methodology; software; validation; formal analysis; investigation;writing—original draft preparation; writing—review and editing: M.D., A.C.L. and P.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Sapienza Università di Roma Project “Fractional derivatives in Science and Engineering” (2018).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Podlubny, I. Fractional Differential Equations; Mathematics in Science and Engineering; Academic Press: Cambridge, MA, USA, 1999; Volume 198. [Google Scholar]
  2. Valentim, C.A., Jr.; Oliveira, N.A.; Rabi, J.A.; David, S.A. Can fractional calculus help improve tumor growth models? J. Comput. Appl. Math. 2020, 379, 112964. [Google Scholar] [CrossRef]
  3. Yang, X.J.; Tenreiro Machado, J. A new insight into complexity from the local fractional calculus view point: Modelling growths of populations. Math. Mod. Meth. Appl. Sci. 2017, 40, 6070–6075. [Google Scholar] [CrossRef] [Green Version]
  4. Haavelmo, T. A Study in the Theory of Economic Evolution; North-Holland: Amsterdam, The Netherlands, 1964. [Google Scholar]
  5. Helen, M. Byrne, Modelling Avascular Tumour Growth. In Cancer Modelling and Simulation; Preziosi, L., Ed.; Chapman & Hall/CRC: Boca Raton, FL, USA, 2003. [Google Scholar]
  6. D’Ovidio, M.; Loreti, P. Solutions of fractional logistic equations by Euler’s numbers. Phys. A Stat. Mech. Appl. 2018, 506, 1081–1092. [Google Scholar] [CrossRef] [Green Version]
  7. Verhulst, P.-F. Notice sur la loi que la population poursuit dans son accroissement. Corresp. Math. Phys. 1838, 10, 113–121. [Google Scholar]
  8. Verhulst, P.-F. Recherches mathématiques sur la loi d’accroissement de la population. Nouv. Mémoires L’Académie R. Sci. Belles Lett. Brux. 1845, 18, 8. [Google Scholar]
  9. Verhulst, P.-F. Deuxième mémoire sur la loi d’accroissement de la population. Mémoires L’Académie R. Sci. Lett. Beaux-Arts Belg. 1847, 20, 1–32. [Google Scholar]
  10. Balzotti, C.; D’Ovidio, M.; Loreti, P. Fractional SIS Epidemic Models. Fractal Fract. 2020, 4, 44. [Google Scholar] [CrossRef]
  11. Area, I.; Nieto, J.J. Power series solution of the fractional logistic equation. Phys. A Stat. Mech. Appl. 2021, 573, 125947. [Google Scholar] [CrossRef]
  12. D’Ovidio, M.; Lai, A.C.; Loreti, P. Generalized binomials in fractional calculus. arXiv 2020, arXiv:2010.05610. [Google Scholar]
  13. Diethelm, K. The Analysis of Fractional Differential Equations. In Lecture Notes in Mathematics; Springer: Berlin/Heidelberg, Germany, 2010. [Google Scholar]
Figure 1. On the (left), solutions of the logistic equation C D β u u = u 2 (on the left) with β = 1 / 2 and u 0 = 1 / 2 , 1 / 3 , 1 / 4 , 1 / 5 ; darker lines correspond to bigger u 0 s. On the (right), solutions of the logistic equation C D β u u = u p + 1 with β = u 0 = 1 / 2 and p = 1 , 2 , 3 ; darker lines correspond to bigger ps.
Figure 1. On the (left), solutions of the logistic equation C D β u u = u 2 (on the left) with β = 1 / 2 and u 0 = 1 / 2 , 1 / 3 , 1 / 4 , 1 / 5 ; darker lines correspond to bigger u 0 s. On the (right), solutions of the logistic equation C D β u u = u p + 1 with β = u 0 = 1 / 2 and p = 1 , 2 , 3 ; darker lines correspond to bigger ps.
Fractalfract 05 00057 g001
Figure 2. Numerical solutions of C D β u u = u 2 (on the left) and of C D β u = u + u 2 (on the right) with u 0 = 1 / 2 and β = 1 , 1 / 2 , 1 / 3 , 1 / 4 . Darker lines correspond to bigger β ’s.
Figure 2. Numerical solutions of C D β u u = u 2 (on the left) and of C D β u = u + u 2 (on the right) with u 0 = 1 / 2 and β = 1 , 1 / 2 , 1 / 3 , 1 / 4 . Darker lines correspond to bigger β ’s.
Fractalfract 05 00057 g002
Figure 3. Solutions of C D β u u = u 2 (on the left) and of C D β u + u = u 2 (on the right) with β = 1 / 2 and u 0 = 1 / 2 , 1 / 3 , 1 / 4 , 1 / 5 . Darker lines correspond to bigger u 0 ’s.
Figure 3. Solutions of C D β u u = u 2 (on the left) and of C D β u + u = u 2 (on the right) with β = 1 / 2 and u 0 = 1 / 2 , 1 / 3 , 1 / 4 , 1 / 5 . Darker lines correspond to bigger u 0 ’s.
Fractalfract 05 00057 g003
Figure 4. Solutions of C D β u u = u p + 1 (on the left) and of D β u + u = u p + 1 (on the right) with β = u 0 = 1 / 2 and p = 1 , 2 , 3 . Darker lines correspond to bigger p’s.
Figure 4. Solutions of C D β u u = u p + 1 (on the left) and of D β u + u = u p + 1 (on the right) with β = u 0 = 1 / 2 and p = 1 , 2 , 3 . Darker lines correspond to bigger p’s.
Fractalfract 05 00057 g004
Figure 5. Solutions of C D β u + u = u p + 1 (on the left) and of C D β u = u u ( p + 1 ) (on the right) with β = 1 / 2 and u 0 = 1 / 2 , 1 / 3 , 1 / 4 , 1 / 5 and p = 1 , 2 . Darker lines correspond to bigger u 0 ’s; continuous lines correspond to p = 2 and dashed lines correspond to p = 1 .
Figure 5. Solutions of C D β u + u = u p + 1 (on the left) and of C D β u = u u ( p + 1 ) (on the right) with β = 1 / 2 and u 0 = 1 / 2 , 1 / 3 , 1 / 4 , 1 / 5 and p = 1 , 2 . Darker lines correspond to bigger u 0 ’s; continuous lines correspond to p = 2 and dashed lines correspond to p = 1 .
Fractalfract 05 00057 g005
Figure 6. Asymptotic behavior of the sequence approximating the radius of convergence r n for a 0 = a 1 = 1 (on the left) and a 0 = a 1 = 1 (on the right) with β = 1 / 2 , u 0 = 1 / 3 and p = 1 , 2 , 3 in blue, orange, and green, respectively.
Figure 6. Asymptotic behavior of the sequence approximating the radius of convergence r n for a 0 = a 1 = 1 (on the left) and a 0 = a 1 = 1 (on the right) with β = 1 / 2 , u 0 = 1 / 3 and p = 1 , 2 , 3 in blue, orange, and green, respectively.
Fractalfract 05 00057 g006
Figure 7. Asymptotic behavior of the sequence approximating the radius of convergence r n for a 0 = a 1 = 1 (on the left) and a 0 = a 1 = 1 (on the right) with u 0 = 1 / 3 , p = 1 , and β = 1 , 1 / 2 , 1 / 3 , 1 / 4 in blue, orange, green, and red, respectively.
Figure 7. Asymptotic behavior of the sequence approximating the radius of convergence r n for a 0 = a 1 = 1 (on the left) and a 0 = a 1 = 1 (on the right) with u 0 = 1 / 3 , p = 1 , and β = 1 , 1 / 2 , 1 / 3 , 1 / 4 in blue, orange, green, and red, respectively.
Fractalfract 05 00057 g007
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

D’Ovidio, M.; Lai, A.C.; Loreti, P. Solutions of Bernoulli Equations in the Fractional Setting. Fractal Fract. 2021, 5, 57. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract5020057

AMA Style

D’Ovidio M, Lai AC, Loreti P. Solutions of Bernoulli Equations in the Fractional Setting. Fractal and Fractional. 2021; 5(2):57. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract5020057

Chicago/Turabian Style

D’Ovidio, Mirko, Anna Chiara Lai, and Paola Loreti. 2021. "Solutions of Bernoulli Equations in the Fractional Setting" Fractal and Fractional 5, no. 2: 57. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract5020057

Article Metrics

Back to TopTop