arxiv: v2 [math.ca] 13 Oct 2015

Similar documents
SMOOTHNESS PROPERTIES OF SOLUTIONS OF CAPUTO- TYPE FRACTIONAL DIFFERENTIAL EQUATIONS. Kai Diethelm. Abstract

This work has been submitted to ChesterRep the University of Chester s online research repository.

Multi-Term Linear Fractional Nabla Difference Equations with Constant Coefficients

A Survey of Numerical Methods in Fractional Calculus

ANALYSIS AND APPLICATION OF DIFFUSION EQUATIONS INVOLVING A NEW FRACTIONAL DERIVATIVE WITHOUT SINGULAR KERNEL

Maximum principle for the fractional diusion equations and its applications

Math 104: Homework 7 solutions

DIfferential equations of fractional order have been the

arxiv: v2 [math.ca] 8 Nov 2014

arxiv: v2 [math.ap] 20 Nov 2017

arxiv: v1 [math.co] 28 Oct 2016

Fractional differential equations with integral boundary conditions

ON FRACTIONAL RELAXATION

Math 141: Section 4.1 Extreme Values of Functions - Notes

Higher order numerical methods for solving fractional differential equations

Existence of global solutions of some ordinary differential equations

Numerical solution of the Bagley Torvik equation. Kai Diethelm & Neville J. Ford

Caculus 221. Possible questions for Exam II. March 19, 2002

A Mean Value Theorem for the Conformable Fractional Calculus on Arbitrary Time Scales

6.2 Deeper Properties of Continuous Functions

AN EFFICIENT PARALLEL ALGORITHM FOR THE NUMERICAL SOLUTION OF FRACTIONAL DIFFERENTIAL EQUATIONS. Kai Diethelm 1,2

Upper and lower solutions method and a fractional differential equation boundary value problem.

An Analytical Scheme for Multi-order Fractional Differential Equations

Estimates for probabilities of independent events and infinite series

Oscillation theorems for nonlinear fractional difference equations

5.5 Deeper Properties of Continuous Functions

Bounded Tiling-Harmonic Functions on the Integer Lattice

The Lusin Theorem and Horizontal Graphs in the Heisenberg Group

EXISTENCE AND UNIQUENESS OF SOLUTIONS FOR A SYSTEM OF FRACTIONAL DIFFERENTIAL EQUATIONS 1. Yong Zhou. Abstract

A computationally effective predictor-corrector method for simulating fractional order dynamical control system

A CONSTRUCTION OF ARITHMETIC PROGRESSION-FREE SEQUENCES AND ITS ANALYSIS

Chapter 8. Infinite Series

Constructive proof of deficiency theorem of (g, f)-factor

4 3A : Increasing and Decreasing Functions and the First Derivative. Increasing and Decreasing. then

1 Definition of the Riemann integral

Lesson 59 Rolle s Theorem and the Mean Value Theorem

Math 421, Homework #6 Solutions. (1) Let E R n Show that = (E c ) o, i.e. the complement of the closure is the interior of the complement.

MATH 409 Advanced Calculus I Lecture 16: Mean value theorem. Taylor s formula.

Applied Mathematics Letters

ELEMENTARY NUMBER THEORY AND METHODS OF PROOF

THE CONLEY ATTRACTORS OF AN ITERATED FUNCTION SYSTEM

8. Limit Laws. lim(f g)(x) = lim f(x) lim g(x), (x) = lim x a f(x) g lim x a g(x)

A truncation regularization method for a time fractional diffusion equation with an in-homogeneous source

( f ^ M _ M 0 )dµ (5.1)

On intermediate value theorem in ordered Banach spaces for noncompact and discontinuous mappings

Nonlinear Analysis 71 (2009) Contents lists available at ScienceDirect. Nonlinear Analysis. journal homepage:

New computational method for solving fractional Riccati equation

Notes 6 : First and second moment methods

Equilibrium points, stability and numerical solutions of fractional-order predator prey and rabies models

MTH4100 Calculus I. Week 8 (Thomas Calculus Sections 4.1 to 4.4) Rainer Klages. School of Mathematical Sciences Queen Mary, University of London

Large topological cliques in graphs without a 4-cycle

The Fractional-order SIR and SIRS Epidemic Models with Variable Population Size

MATH 140B - HW 5 SOLUTIONS

arxiv: v1 [math.co] 4 Jan 2018

DETERMINATION OF AN UNKNOWN SOURCE TERM IN A SPACE-TIME FRACTIONAL DIFFUSION EQUATION

Premiliminary Examination, Part I & II

arxiv: v1 [math.mg] 15 Jul 2013

Numerical Fractional Calculus Using Methods Based on Non-Uniform Step Sizes. Kai Diethelm. International Symposium on Fractional PDEs June 3 5, 2013

MULTISTAGE HOMOTOPY ANALYSIS METHOD FOR SOLVING NON- LINEAR RICCATI DIFFERENTIAL EQUATIONS

Vertex-Coloring Edge-Weighting of Bipartite Graphs with Two Edge Weights

NON-MONOTONICITY HEIGHT OF PM FUNCTIONS ON INTERVAL. 1. Introduction

arxiv: v1 [math.ca] 28 Feb 2014

Received: 11 January 2019; Accepted: 25 January 2019; Published: 4 February 2019

Oscillatory Solutions of Nonlinear Fractional Difference Equations

Polynomial functions on subsets of non-commutative rings a link between ringsets and null-ideal sets

SOME STRUCTURE THEOREMS FOR INVERSE LIMITS WITH SET-VALUED FUNCTIONS

Math 321 Final Examination April 1995 Notation used in this exam: N. (1) S N (f,x) = f(t)e int dt e inx.

ζ(u) z du du Since γ does not pass through z, f is defined and continuous on [a, b]. Furthermore, for all t such that dζ

arxiv:math/ v1 [math.co] 6 Dec 2005

Optimal primitive sets with restricted primes

Unmixed Graphs that are Domains

THE HOT SPOTS CONJECTURE FOR NEARLY CIRCULAR PLANAR CONVEX DOMAINS

Stabilization of fractional positive continuous-time linear systems with delays in sectors of left half complex plane by state-feedbacks

Algebraic methods toward higher-order probability inequalities

ON THE OSCILLATION OF THE SOLUTIONS TO LINEAR DIFFERENCE EQUATIONS WITH VARIABLE DELAY

SOME RESULTS FOR BOUNDARY VALUE PROBLEM OF AN INTEGRO DIFFERENTIAL EQUATIONS WITH FRACTIONAL ORDER

Tangent spaces, normals and extrema

Institut für Mathematik

Picard s Iterative Method for Caputo Fractional Differential Equations with Numerical Results

MATH 117 LECTURE NOTES

Bull. Math. Soc. Sci. Math. Roumanie Tome 60 (108) No. 1, 2017, 3 18

Dynamic diffusion-type equations on discrete-space domains

arxiv: v3 [math.ra] 10 Jun 2016

Sequences of Real Numbers

A generalized Gronwall inequality and its application to fractional differential equations with Hadamard derivatives

Analytic solution of fractional integro-differential equations

The Lost Cousin of the Fundamental Theorem of Algebra

#A69 INTEGERS 13 (2013) OPTIMAL PRIMITIVE SETS WITH RESTRICTED PRIMES

Calculus (Math 1A) Lecture 4

Local strong convexity and local Lipschitz continuity of the gradient of convex functions

Calculus (Math 1A) Lecture 4

Leighton Coles Wintner Type Oscillation Criteria for Half-Linear Impulsive Differential Equations

EXTENDED LAGUERRE INEQUALITIES AND A CRITERION FOR REAL ZEROS

Numerical Solution of Multiterm Fractional Differential Equations Using the Matrix Mittag Leffler Functions

18.785: Analytic Number Theory, MIT, spring 2006 (K.S. Kedlaya) Dirichlet series and arithmetic functions

First In-Class Exam Solutions Math 410, Professor David Levermore Monday, 1 October 2018

NORMS ON SPACE OF MATRICES

Cauchy Integral Formula Consequences

APPLICATIONS OF DIFFERENTIATION

Solutions Final Exam May. 14, 2014

Transcription:

arxiv:1510.03285v2 [math.ca] 13 Oct 2015 MONOTONICITY OF FUNCTIONS AND SIGN CHANGES OF THEIR CAPUTO DERIVATIVES Kai Diethelm 1,2 Abstract It is well known that a continuously differentiable function is monotone in an interval [a,b] if and only if its first derivative does not change its sign there. We prove that this is equivalent to requiring that the Caputo derivatives of all orders α 0,1) with starting point a of this function do not have a change of sign there. In contrast to what is occasionally conjectured, it not sufficient if the Caputo derivatives have a constant sign for a few values of α 0,1) only. MSC 2010: Primary 26A33 Key Words and Phrases: fractional calculus, Caputo derivative, monotone function, sign change. 1. Introduction Classical integer order calculus provides numerous results that link the shape properties of a function s graph such as, e. g., monotonicity or the location of extrema) with the signs of its first or second derivative. The following theorems contain some specific examples; since they are so well known, we omit the proofs. Theorem 1.1. A function f C 1 [a,b] is monotone on [a,b] if and only if its first derivative f does not change its sign in this interval. October 12, 2015

2 K. Diethelm Theorem 1.2. If a function f C 1 [a,b] attains a global maximum or a global minimum at some point t a,b) then f t ) = 0. If a function f C 2 [a,b] attains a global maximum at some point t a,b) then f t ) = 0 and f t ) 0; if t a,b) is the location of a global minimum of f then f t ) = 0 and f t ) 0. It is quite natural to ask for generalizations of such statements that involve fractional, rather than integer order, derivatives. In this paper we concentrate on fractional derivatives of Caputo type [2, Chapter 3], i. e. D a α ft) := J α α a D α ft) 1.1) for α > 0 and α / N where denotes the ceiling function, D m is the standard mth order differential operator for m N and Ja β is the Riemann- Liouville integral operator of order β > 0 with starting point a, viz. Ja β 1 gt) := Γβ) t a t s) β 1 gs)ds. 1.2) In this case, the following fractional extension of the first part of Theorem 1.2 has been shown by Luchko [4, Theorem 1]: Theorem 1.3. If a function f C 1 [a,b] attains a global maximum at some point t a,b) then D α a ft ) 0 for all α 0,1). If t a,b) is thelocation of theglobal minimum of f then D α a ft ) 0for all α 0,1). In fact, Luchko s result is more general than our formulation here since he admits functions f from a slightly larger class of functions; however, for the sake of simplicity, we have chosen this marginally restricted version. Some authors have believed that this result can be extended to fractional derivatives of order α 1, 2); specifically they claimed that, under the additional condition f C 2 [a,b], the inequality D α aft ) 0 holds for these α if t is a global minimum point of f. While Luchko s result Theorem 1.3) extends the first part of the classical Theorem 1.2, such aresultwould providean extension of thesecond part of Theorem 1.2. However, recently Al-Refai [1] who also provides some references to claims of such results in the literature) explicitly constructed a counterexample and hence demonstrated that the second part of Theorem 1.2 cannot be extended in this manner. Our goal here is to perform a similar investigation for possible fractional versions of Theorem 1.1. In a certain sense, it thus continues the author s earlier work [3] that transferred other classical calculus results,

MONOTONICITY AND SIGN OF DERIVATIVE... 3 including the mean value theorems of differential and integral calculus, to the fractional setting. 2. Caputo derivatives of monotone functions We begin our investigation of a generalization of Theorem 1.1 to the fractional case with a very simple but nevertheless useful observation. Theorem 2.1. Let the function f C 1 [a,b] be monotone on [a,b]. Then, for all α 0,1) and all t [a,b], we have that { D a α ft) 0 if f is decreasing, 0 if f is increasing. In this formulation, we use the terms increasing and decreasing in the non-strict sense. P r o o f. Assume that f is increasing. Then, f is nonnegative in [a,b], and hence, by the definition of the Caputo operator given in eq. 1.1) and eq. 1.2), D aft) α 1 t = t s) α f s)ds 0. Γ1 α) a In the case of a decreasing function f, we argue in an analog way. We obtain an immediate consequence that can be interpreted as the fractional generalization of the forward direction of Theorem 1.1. Corollary 2.1. Let the function f C 1 [a,b] be monotone on [a,b]. Then, for all α 0,1), the Caputo derivative D a α f does not have a change of sign in [a,b]. But we can also obtain a fractional generalization of the backward direction of Theorem 1.1: Theorem 2.2. Assume that f C 1 [a,b] is such that D α aft) 0 for all t [a,b] and all α α 0,1) with some α 0 0,1). Then, f is monotone increasing. Similarly, if D α aft) 0 for all t and α mentioned above, then f is monotone decreasing. P r o o f. Our observation in the first case is based on the definition of the Caputo derivative and certain well known properties of the Riemann- Liouville integral. Specifically, consider first the case that t a, b]. Then,

4 K. Diethelm in view of eq. 1.1) and our assumption, 0 D α aft) = J 1 α a f t). Since f is assumed to be continuous, we may invoke [2, Theorem 2.10] to determine the limit α 1 of the right-hand side of this inequality; this yields 0 lim α 1 Dα aft) = lim α 1 J1 α a f t) = f t). We thus conclude that f is nonnegative on the half-open interval a,b]. The continuity of f on the closed interval [a,b] then implies that f a) 0 as well; hence f is nonnegative on the entire interval [a,b], and it follows that f is increasing. The second claim follows immediately if we apply the result that we have just shown to f instead of f. Theorem 2.2 has an interesting consequence that we explicitly note. Theorem 2.3. Assume that f C 1 [a,b] is such that D a α ft) 0 for all t [a,b] and all α α 0,1) with some α 0 0,1). Then, D a α ft) 0 for all t [a,b] and all α 0,1). The analog statement holds when the sign of the derivatives is reversed. Thus, if we can prove that the Caputo derivatives of f are of a constant signfor all α sufficiently close to1, thenit already follows that thisproperty even holds for all α from the complete interval 0,1) and, by Theorem 2.2, also for the classical case α = 1. P r o o f. If the Caputo derivatives D aft) α are nonnegative for the indicated values of t and α then Theorem 2.2 tells us that f is increasing; the claim then follows from Theorem 2.1. The statement for the case of opposite signs of D a α ft) follows in an analog manner. It thus turns out that, contrary to an occasionally stated conjecture, the lack of a change of sign of D α af for some, but not all, α 0,1) is not sufficient to imply the monotonicity of f. One may ask how a non-monotone functions which has a fractional derivative without a change of sign may look like, and in particular whether such a function must necessarily be of a very exotic or even pathological nature. The following example demonstrates that this is not the case; indeed a simple polynomial of a quite low degree can already have this property. An interesting aspect in this connection is the fact that our counterexample is the same one that Al-Refai used in his investigations from [1, Section 2] regarding the fractional version of the second part of Theorem 1.2 that we mentioned above.

MONOTONICITY AND SIGN OF DERIVATIVE... 5 Example 2.1. Define the function g by gt) := t t 1 ) t 1) = t 3 3 2 2 t2 + 1 2 t, so that and D α 0gt) = = g t) = 3t 2 3t+ 1 2. 2.1) 6 Γ4 α) t3 α 6t 1 α Γ4 α) 3 1 Γ3 α) t2 α + 2Γ2 α) t1 α t 2 1 2 3 α)t+ 1 ) 12 3 α)2 α) 2.2) for 0 < α < 1. It follows from eq. 2.2) that the only zeros of D 0 α g are located at t 1 α) = 0, t 2 α) = 3 α ) 3 α)1+α)/3 /4 and t 3 α) = 3 α+ ) 3 α)1+α)/3 /4. Clearly, t 1 α) < t 2 α) t 3 α), andhenceweconcludethatd 0 α g doesnothaveachangeofsignon[0,t 2α)]. But now, ) d dα t 2α) = 1 1 α 1 < 0 4 33 α)1+α) for 0 < α < 1. Therefore, if we denote by I α := [0,t 2 α)] the largest interval that has the left end point at 0 and that is such that D 0 α g does not have a change of sign on I α, then these intervals monotonically shrink as α increases. Specifically, 1 2 = t 20) > t 2 α 1 ) > t 2 α 2 ) > t 2 1) = 3 3 0.211325, 6 for 0 < α 1 < α 2 < 1, and hence [ 0, 1 ] [ ] = I 0 I α1 I α2 I 1 = 0, 3 3. 2 6 Therefore, if we pick some α 0,1), then we can see that D α 0 g is nonnegative throughout the interval I := I α = [0,t 2 α )]. Moreover, it follows from our monotonicity observation on t 2 that all derivatives D α 0 g with α 0,α ) are nonnegative on the interval I as well. Nevertheless, an elementary computation reveals that the only zeros of g are located at τ 1 = 3 3)/6 = t 2 1) and τ 2 = 3+ 3)/6 = 1 τ 1 0.788675. From this observation and the fact that g is a polynomial of degree 3 with a positive leading coefficient, we conclude that g has a local maximum at the point t 2 1) which is located in the interior of our interval I ; hence, in spite

6 K. Diethelm of the fact that D α 0 gt) > 0 for almost all t I with the sole exception of the point t = 0 for which D α 0 gt) = 0 holds) and all α 0,α ], the function g is not monotone on I. Remark 2.1. We point out that all our investigations are based on the fact that the left end point of the interval on which the monotonicity is discussed coincides with the starting point of the Caputo differential operator. This is a very natural requirement because of the non-locality of the fractional differential operator. Specifically, if we had a < a < b and wanted to find out whether a function f is monotonic on [a,b], then a criterion based on the behaviour of D α af would not be of any help since the values of D α af may de influenced to an arbitrarily large extent by the behaviour of f on the interval [a,a ], i. e. outside of the interval in which we are interested. Therefore, in a situation like this only integer-order derivatives can be of any use because they are local operators and hence do not depend on the behaviour of f outside of the interval of interest. References [1] M. Al-Refai, On the fractional derivatives at extreme points. Electron. J. Qual. Theory Differ. Equ. 2012), Paper No. 55. [2] K. Diethelm, Fractional Differential Equations. Springer, Berlin 2010). [3] K. Diethelm, The mean value theorems and a Nagumo-type uniqueness theorem for Caputo s fractional calculus. Fract. Calc. Appl. Anal. 15, No. 2 2012), 304 313. [4] Y. Luchko, Maximum principle for the generalized time-fractional diffusion equation. J. Math. Anal. Appl. 351 2009), 218 223. 1 AG Numerik Institut Computational Mathematics Technische Universität Braunschweig Pockelsstr. 14 38106 Braunschweig, GERMANY e-mail: k.diethelm@tu-bs.de 2 GNS Gesellschaft für numerische Simulation mbh Am Gaußberg 2 38114 Braunschweig, GERMANY e-mail: diethelm@gns-mbh.com