arxiv: v1 [math.na] 15 Nov 2013
|
|
- Ami Cross
- 6 years ago
- Views:
Transcription
1 NUMERICAL APPROXIMATIONS FOR FRACTIONAL DIFFERENTIAL EQUATIONS arxiv: v1 [math.na] 15 Nov 013 Yuri Dimitrov Department of Applied Mathematics and Statistics University of Rousse 8 Studentsa str. Rousse 7017, Bulgaria ymdimitrov@uni-ruse.bg November 18, 013 Abstract The Grünwald and shifted Grünwald formulas for the function yx yb are first order approximations for the Caputo fractional derivative of the function yx with lower limit at the point b. We obtain second and third order approximations for the Grünwald and shifted Grünwald formulas with weighted averages of Caputo derivatives when sufficient number of derivatives of the function yx are equal to zero at b, using the estimate for the error of the shifted Grünwald formulas. We use the approximations to determine implicit difference approximations for the sub-diffusion equation which have second order accuracy with respect to the space and time variables, and second and third order numerical approximations for ordinary fractional differential equations. 010 Mathematics Subject Classification: 6A33, 34A08, 65M1 Key Words and Phrases: fractional differential equation, implicit difference approximation, Grünwald formula, stability, convergence 1 Introduction Fractional derivatives are an effective tool for modeling diffusion processes in complex systems. Mathematical models with partial fractional differential equations have been used to describe complex processes in physics, biology, chemistry and economics [10-0]. The time fractional diffusion equation is 1
2 a parabolic partial fractional differential equation obtained from the heatdiffusion equation by replacing the time derivative with a fractional derivative of order α. α ux,t = ux,t +Gx,t. 1 t α x When 0 < α < 1 the equation is called fractional sub-diffusion equation and it is a model of a slow diffusion process. When the order of the fractional deriative is between one and two the equation is called fractional superdiffusion equation. There is a growing need to design efficient algorithms for numerical solution of partial fractional differential equations. Finitedifference approximations for the heat-diffusion α = and time-fractional diffusion equations have been studied [38-50] for their importance in practical applications as well as for evaluation of their performance. A common way to approximate the Caputo derivative of order α, when 0 < α < 1, is approximation 6. The finite difference approximation for the fractional sub-diffusion equation which uses approximation 6 for the time fractional derivative and central difference approximation for the second derivative with respect to x has accuracy Oτ α +h, where h and τ are the step sizes of the discretizations with respect to the space and time variables x and t. Tadjeran et al. [3] use the estimate for the error for the Grünwald formula 13 to design an algorithm for a second order numerical approximation of the solution of the space fractional diffusion equation of order α, when 1 < α <. ux, t t = dx α ux,t x α +Gx,t. The algorithm uses a Cran-Nicholson approximation with respect to the time variable t and an extrapolation with respect to the space variable x. Ding and Li [4] compute a numerical solution of the fractional diffusionwave equation with reaction term α ux,t t α = K α ux,t x C α ux,t+gx,t, where K α > 0 and C α > 0 are the diffusion and reaction coefficients, using a compact difference approximation with an accuracy Oτ +h 4. Gorenflo [5] showed that the shifted Grünwald formula h α α h,α/yx is a second order approximation for the fractional derivative y α x = h α α h,α/ yx+o h
3 when the transition of yx to zero is sufficiently smooth at the lower limit of fractional differentiation. The smoothness condition requires that yb = 0. This approximation is a special case of Theorem 1i, when p = α/. There is a significant interest in designing efficient numerical solutions for ordinary and partial fractional differential equations, which stems from the possibility to use fractional derivatives to explain complex processes in nature and social sciences and their relation to integer order differential equations. While approximations of fractional derivatives with accuracy Oh α have been studied extensively, new algorithms with second and higher order accuracy [7-37] have been proposed and successfully applied for numerical solution of ordinary and partial fractional differential equations. In this paper we construct a second order implicit difference approximation for the fractional sub-diffusion equation 1 and second and third order approximations for the ordinary fractional differential equation y α x+yx = fx when equations1 and have sufficiently smooth solutions. In section 3 we use the estimate for the error of the Grünwald and shifted Grünwald formulas 13 to obtain second and third order approximations for the Grünwald and shifted Grünwald formulas with weighted averages of Caputo derivatives on a uniform grid, when sufficient number of derivatives of the function yx are equal to zero at the point b. In section 4 we derive recurrence relations 5, 6 and 33 for second and third order approximations to the solution of ordinary fractional differential equation and we determine estimates for the Grünwald weights. In section 5 we use approximations 9 and 10 to construct implicit difference approximations 56 and 57 for the solution of the sub-diffusion equation 1 and we show that they have second order accuracy Oτ +h with respect to the space and time variables. Preliminaries The fractional derivatives are generalizations of the integer order derivatives. Let yx be a real-valued function defined for x b. The Riemann-Liouville and Caputo fractional derivatives of order α, when 0 < α < 1 are defined as D α RLyx = 1 d Γ1 α dx x b yξ x ξ αdξ, 3 Dx α yx = yα x = dα dx αyx = 1 x y ξ Γ1 α b x ξ αdξ. 4 3
4 The Caputo fractional derivative of the constant function 1 is zero and the Riemann-Liouvillederivativeof1isx b α /Γ1 α. TheCaputoderivative of the function yx satisfies D α x yx = Dα x yx yb. If a function yx is Caputo differentiable of order α then it is differentiable in the sense of the definition of Riemann-Liouville derivative. The classes Caputo differentiable functions C α, when 0 < α < 1, include the C 1 functions and are suitable for numerical computations for fractional differentiable equations. For a strict definition of a fractional derivative we assume Lebesgue integration in the definitions of Caputo and Riemann-Liouville fractional derivatives. The Caputo and Riemann-Liouville fractional derivatives satisfy [1, p. 53] D α RLyx = D α Cyx+ yb Γ1 αx b α. 5 The Caputo and Riemann-Liouville fractional derivatives of the function yx are equal when yb = 0. The Miller-Ross sequential fractional derivative of order α 1 +α for the Caputo derivative is defined as y α 1+α x = D α 1 x D α x yx. In the special cases α 1 = 1,α 1 = and α = α, when 0 < α < 1 y 1+α x = d dx yα x, y +α x = d dx yα x. The local behavior of a differentiable function is described with its Taylor series and Taylor polynomials. The properties of a function yx close to the lower limit b can be described with its Caputo and Miller-Ross derivatives at the point b. The following theorem is a generalization of the Mean-Value Theorem for differentiable functions [57]. Theorem. Generalized Mean-Value Theorem Let y C α [b,x]. Then yx = yb+ x bα Γα+1 yα ξ x b ξ x x. Fractional Taylor series for Caputo and Miller-Ross derivatives [57, 58] can be defined using an approach similar to the derivation of the classical Taylor series which involves only integer order derivatives. While the local 4
5 properties of a function at the lower limit b can be explained with its fractionalderivativestheirvaluesatanyotherpointxdependonthevaluesofthe function on the interval [b,x]. We can observe a similarity between the Taylor series of the function yx and relation 13 for its fractional derivatives and the shifted Grünwald formulas. Two important special functions in fractional calculus are the gamma and Mittag-Leffler functions. The gamma function has properties Γ0 = 1, Γz +1 = zγz. When n is a positive integer Γn = n 1!. The one-parameter and twoparameter Mittag-Leffler functions are defined for α > 0 as E α z = n=0 z n Γαn+1, E α,βz = n=0 z n Γαn+β. Some special cases of the one-parameter Mittag-Leffler function E 1 z = e z, E z = cosz, E1z = e z erfc z, where ercfz is the complimentary error function ercfz = e t dt. π The Mittag-Leffler functions appear in the solutions of ordinary and fractional differential equations [, chap. 1]. The ordinary differential equation y x = yx,y0 = 1,y 0 = 1 hassolutionyx = E,1 x +xe, x = coshx+sinhx = e x. Thefractional differential equation y α x = λyx,y0 = 1 has solution yx = E α λx α. We can determine the analytical solutions of linear ordinary and partial fractional differential equations using integral transforms. The following formulas for the Laplace transform of the derivatives of the Mittag-Leffler functions are often used to determine the analytical solutions of linear fractional differential equations []. L{t α+β 1 E α,β ±atα }s =!sα β s α a +1. z 5
6 Analytical solutions of ordinary and partial fractional differential equations canbefoundonlyforspecial cases oftheequationsandtheinitialandboundary conditions. We can determine numerical approximations for the solutions of a much larger class of equations which include nonlinear fractional differential equations. Approximations for the Caputo and Riemann-Liouville derivatives are obtained from the Grünwald-Letniov fractional derivative and by approximating the fractional integral in the definition. Let x n = b+nh be a uniform grid on the x-axis starting from the point b, and y n = yx n = yb + nh, where h > 0 is a small number. The following approximation of the Caputo derivative is derived from quadrature approximations of the fractional integral for y x in the definition of Caputo fractionalderivative[39]ontheinterval[b,x n ], byapproximatingtheintegrals on all subintervals of length h. where c α 0 = 1/Γ α and c α y n α = 1 n 1 c α h α y n +Oh α, 6 =0 = +11 α 1 α α. Γ α When the function yx has continuous second derivative, approximation 6 has accuracy Oh α. The weights c α satisfy c α 0 > 0, c α 1 < c α < < c α < < 0, c α = 0. The Grünwald-Letniov fractional derivative is closely related to Caputo and Riemann-Liouville derivatives DGL α yx = lim x 0 1 x α [ x b x] =0 α 1 n yx n x. n n=0 The two most often used values of the lower limit of the fractional differentiation b are zero and. When the lower limit b = the upper limit of the sum in the definition of Grünwald-Letniov derivative is. The fractional binomial coefficients are defined similarly to the integer binomial coefficients with the gamma function α Γα+1 = n Γn+1Γα n+1 = αα 1 α n+1. n! 6
7 The Riemann-Liouville, Caputo and Grünwald-Letniov derivatives are equal when yb = 0 and y C 1 [b,x] [1, p. 43]. Let s denote by α hyx and α h,pyx the Grünwald difference operator and shifted Grünwald difference operators for the function yx. N x,h α h yx = α 1 n yx nh, n n=0 N x,h α h,p yx = α 1 n yx n ph. n n=0 where N x,h = [ ] x b h, and h > 0 is a small number. The Grünwald operator is a special case of the shifted Grünwald operator when the shift value p = 0. When yb = 0, we derive approximations for the Riemann-Liouville and Caputo derivatives from the definition of Grünwald-Letniov derivative. D α x yx = Dα RL yx = Dα GL yx h α α h yx h α α h,p yx. Wewillcall theapproximationsh α α h yxandh α α h,p yxofthefractional derivative of order α - Grünwald formula and shifted Grünwald formulas of the function yx. Let w n α be the weights of the Grünwald formulas. α w n α = 1 n. n When yx is a continuously-differentiable function the shifted Grünwald formulas for the function yx yb are first-order approximations for the Caputo derivative of yx. y α x = 1 h α N x,h n=0 w α n yx n ph yb+oh. 7 In Theorem 1i we show that the shifted Grünwald formulas of the function yx yb are second-order approximations for the Caputo fractional derivative at the point x+p α/h, when the function yx is sufficiently differentiableontheinterval[b,x]andthevaluesofitsfirstandsecondderivatives are equal to zero at the point b. y α x+ p α h = 1 N x,h w α h α n yx n ph yb+o h. 8 n=0 7
8 WederivetheaboveformulainTheorem1ifromrelation13fortheshifted Grünwald formulas and the fractional derivatives of the function yx. When the conditions of Theorem 1 are satisfied the Riemann-Liouville and Caputo derivatives are equal. The values of a function and the shifted Grünwald formulas satisfy β 1 yx 1 +β yx = yβ 1 x 1 +β x +O h, β 1 α h,p yx 1+β α h,q yx = α h,β 1 p+β q yβ 1x 1 +β x +O h. when β 1 +β = 1 and yx is a sufficiently smooth function. An alternative way to obain 8 is to apply the approximation for average value of shifted Grunwald formulas to approximation.15 in Tian et al. [30]. An important special case of the above formula is when yb = 0 and p = 0. y α x αh = 1 N x,h w α h α n yx nh+o h. 9 n=0 This approximation is closely related to formula.9 in [5] and is suitable for constructing second order weighted numerical approximations for fractional differential equations on a uniform grid. As a direct consequence of 9 we obtain a second order approximation for the Grünwald formula using average values of Caputo derivatives on consecutive nodes of a uniform grid. N 1 x,h α w α h α n yx nh = y α n α y n α +O h. 10 n=0 In Corollary 5 we determine a third order approximation for the Grünwald formula using a weighted average of three consecutive values of the Caputo derivative on a uniform grid, for sufficiently differentiable functions yx which satisfy the conditions of Theorem 1ii. N 1 x,h h α n=0 a w n α yx nh = 8 5a a 4 + a 8 y α n + 11a y α n a 4 y n α +O h Approximations 9,10 and 11 are suitable for numerical computations for fractional differential equations on a uniform grid, when the solutions are sufficiently differentiable functions. In section 4 we determine second and third order approximations for ordinary differential equation, and in 8
9 section 5 we construct stable difference approximations for the fractional subdiffusion equation which have second order accuracy with respect to the space and time variables. The Grünwald weights w n α are computed recursively with w α 0 = 1,w α 1 = α and w n α = 1 α+1 w α n n 1. The numbers w n α are the coefficients of the binomial series 1 z α = w n α z n. When α is a positive integer the sum is finite, and the binomial series converges at the point z = 1 when 0 < α < 1. The weights w n α have the following properties. n=0 w α 0 > 0, w α 1 < w α < < w α n < < 0, n=0 w α n = 0. 1 When the lower limit of fractional differentiation b, the upper limit of the sum is finite and N x,h n=0 w n α > 0. The shifted Grünwald formulas for yx are first order approximations for the Riemann-Liouville derivative of the function yx. The approximation error can be represented as a sum of higher order Riemann-Liouville fractional derivatives 13. This estimate is obtained in Tadjeran et al. [3] when the order α of the Riemann-Liouville derivative is between one and two using Fourier transform of the shifted Grünwald formulas. The estimate for the error of the shifted Grünwald formulas is generalized in Hejazi et al. [6] for arbitrary positive α and p. Theorem. Let α and pbe positivenumbers, and supposethaty C [α]+n+ R and all derivatives of y up to order [α] + n + belong to L 1 R. Then if b =, there exist constants c l independent of h,y,x such that n 1 h α α h,pyx = DRLyx+ α l=1 c l h l D α+l RL yx+ohn. 13 The numbers c l are the coefficients of the series expansion of the function 1 e z α ω α,p z = e pz. z If a function yx is defined on a finite interval I we can extend it to the real line by setting yx = 0 when x / I. In this way the Caputo derivative of 9
10 the extended function is equal to zero when x / I. In the next section we use 13 to detrmine second and third order approximations for the Grünwald and shifted Grünwald formulas with weighted averages of Caputo derivatives on consecutive points of a uniform grid. 3 Approximations for Grünwald and shifted Grünwald formulas The definitions of Riemann-Liouville and Caputo derivatives are the two most commonly used definitions for fractional derivatives. The Caputo and Riemann-Liouville fractional derivatives of order n + α, where n is a positive integer and 0 < α < 1 are defined as D n+α x yx = 1 Γ1 α D n+α RL yx = 1 Γ1 α x b d n+1 x dx n+1 b y n+1 ξ x ξ αdξ, yξ x ξ αdξ. When yx is sufficiently differentiable function on the interval [b, x], the Riemann-Liouville and Caputo derivatives are related as [1, p. 53] D n+α RL yx = Dn+α x yx + n =0 y b Γ α 1 x b α n. 14 We can compute the value of the Riemann-Liouville derivative of the function yx from the value of the Caputo derivative and its integer order derivatives at the lower limit b. One disadvantage of of the Riemann-Liouville derivative is that it has a singularity at the lower limit b. When α > 1 the singularity is non-integrable. The class of Caputo differentiable functions of order n + α includes the functions with n + 1 continuous derivatives. In section we discussed properties of the Caputo derivatives at the lower limit b, which are related to the properties of the integer order derivatives. We often prefer to study fractional differential equations with the Caputo derivative because its properties mae it an attractive fractional derivative for numerical solution of fractional differential equations. Let h be the step size of a uniform grid on the x-axis staring from the lower limit of fractional differentiation b and x n = b+nh, y n = yx n = yb+nh. InTheorem1weusetheestimatefortheerroroftheshiftedGrünwaldformulas 13 and relation 14 for the Riemann-Liouville and Caputo derivatives to obtain second and third order approximations for the the shifted Grünwald formulas using fractional order Caputo and Miller-Ross derivatives. 10
11 Theorem 1. Approximations for the shifted Grünwald formulas i Let yb = y b = y b = 0 and y C 4 [b,x]. Then h α α h,p x+ yx = yα p α h +O h ; 15 ii Let y0 = y 0 = y 0 = y b = 0 and y C 5 [b,x]. Then h α α h,pyx = y x+ α p α h + α 4 h y +α x+o h 3. Proof. The function ω α,p z has power series expansion 1 e z α ω α,p z = e pz = c 0 +c 1 z +c z +, z where c 0 = 1,c 1 = p α/ and [3] c = 1 4 1p 1αp+α+3α = α 4 +1 p αp+ α = α p α. From 13 and n = we obtain h α α h,p yx = Dα RL yx+ p α hd α+1 RL yx+o h. Let 0 < β < 1. The Caputo and Riemann-Liouville fractional derivatives of order 3 + β satisfy 14 with n = 3. The Riemann-Liouville derivative D 3+β RL yx has non-integrable singularities of orders 1+β,+β and 3+β at thelowerlimitb, withcoefficientsyb,y bandy b. TheCaputoderivative of yx of order 3 + β and the function x b β are integrable on a finite interval. Therefore D 3+β RL yx L1 [b,x] when yb = y b = y b = 0. Similarly, D 4+β RL yx L1 [b,x] when yb = y b = y b = y 0 = 0. The Riemann-Liouville and Caputo derivatives of order α for the functions yx which satisfy the conditions of Theorem 1 are equal D α RL yx = yα x and we can represent D α+1 RL yx with the Caputo derivative yα x as D α+1 RL yx = 1 Γ1 α d dx x b yξ x ξ αdξ = d dx Dα RLyx = d dx yα x. We obtain the following relation for the shifted Grünwald formulas and the Caputo derivative of the function yx. h α α h,pyx = y α x+ p α h d dx yα x+o h. 11
12 From the mean value theorem for the function y α x we have that y α x+ p α h = y α x+ p α h d dx yα x+o h. Therefore y x+ α p α h = h α α h,p yx+o h. Now we prove ii. From 13 and n = 3 we obtain h α α h,pyx = DRLyx+ α p α hd α+1 RL yx+ α p α h D +α RL yx+o h 3. We have that D +α RL yx = 1 Γ1 α d 3 dx 3 x Then h α α h,pyx =y α x+ p α 1 p α b yξ d d x ξ αdξ = dx Dα RLyx = x. dx yα h d dx yα x+ h d x+ α d dx yα h x+o h 3. 4 dx yα From the Mean-Value Theorem for y α x we obtain h α α h,p x+ yx = yα p α h + α 4 h y +α +O h 3. The most useful special case of 15 is when p = 0. Corollary. Let yb = y b = y b = 0 and y C 4 [b,x]. Then h α α h yx = yα x αh +O h. 16 When α = 1 and α = approximation 16 becomes a central difference approximation for the first and second derivatives of the function yx y x h = yx yx h +O h, h 1
13 y x h = yx yx h+yx h h +O h. When α = n, where n is a positive integer, the weights w α = 0 for > n. This case is discussed in [5]. Experimental results suggest that when the conditions of Theorem 1 are not satisfied, the order of approximation 15 fluctuates. While in some cases the order may still be two, in many cases it is lower, even the order may be one and lower than one Section Let yx be a sufficiently differentiable function. Claim 3. Approximations for values of a function on a uniform grid y n β = βy n 1 +1 βy n +O h, 17 y n β = 1 ββ 1y n +β βy n β 1β y n +O h Proof. From the Mean Value Theorem there exist numbers θ 1 and θ, such that y n = y n β +βhy n β + β h y n θ 1, Hence, y n 1 = y n β 1 βhy n β + 1 β h y n θ. βy n 1 +1 βy n = y n β + 1 ββ h y n θ 1 + β1 β h y n θ. Let D be an upper bound for the second derivative. Then βy n 1 +1 βy n y n β 1 ββd h. The proof for approximation 18 uses third order expansions and is similar to the proof of 17. Lemma 4. Approximations for the shifted Grünwald formulas i Let yb = y b = y b = 0 and y C 4 [b,x]. Then α h α α h,p y n = p y α n p α y n α +O h ; 19 ii Let yb = y b = y b = y b = 0 and y C 5 [b,x]. Then h α α h,p y n = β 1 y α n +β y α n 1 +β 3y α n +O h 3, 0 where β 1 = p + p 5α 4 pα + α 8 and β = p p + 11α 1 +pα α 4, β 3 = 1+ 3p + p 17α 4 pα + α 8. 13
14 Proof. From Theorem 1i with x = x n we have h α α h,py n = y x α n + p α h +O h. From 17 with β = α p we obtain α h α α h,p y n = p y α n p α y n α +O h. Now we use the formula from Theorem 1ii to determine a third order approximation for the shifted Grünwald formulas. h α α h,p y n = y x α n + p α h + α 4 h y n +α +O h 3. The central difference approximation for y n +α with nodes {x n,x n 1,x n } has order Oh Then y +α n h α α h,p y n = y α x n + = yα n y α n 1 +yα n +Oh, h p α h + αh y α 4 n y α n 1 +y α n h +Oh +O h 3. From 18 with β = α p we obtain approximation 0. We obtain second and third order approximations for the Grünwald formula from 19 and 0 with p = 0. Corollary 5. Approximations for the Grünwald formula i Let yb = y b = y b = 0 and y C 4 [b,x]. Then α h α α hy n = y α n α y n α +O h ; ii Let yb = y b = y b = y b = 0 and y C 5 [b,x]. Then a h α α h y n = 8 5a y α 11a 4 n + 1 a y α 4 n a 4 + a 8 y α n +O h 3. In Corollary and Corollary 5 we determined second and third order approximations for the Grünwald formula using values of Caputo derivatives. The three approximations are suitable for algorithms for numerical solution of fractional differential equations on a uniform grid. In the next lemma, we discuss a property of the Caputo derivative of a continuously differentiable function in a neighborhood of the lower limit b. 14
15 Lemma 6. Let 0 < α < 1 and y C 1 [b,b+ǫ], where ǫ > 0. Then y α b = 0. Proof. The function y is bounded on the interval [b,b+ǫ]. Let When b < x < b+ǫ we have y α x 1 Γ1 α y α x x From the squeeze low of limits δ = max b x b+ǫ y x. y ξ b x ξ αdξ δ Γ1 α δ x ξ1 α x Γ1 α 1 α 0 lim x b y α x lim x b δx b 1 α Γ α b x b x ξ α dξ, δx b1 α Γ α. 1 = 0. Hence, y α b = lim x b y α x = 0. In the next two sections we compute approximations for the solutions of equations 1 and when b = 0 and the solutions are sufficiently differentiable functions. We use the result from Lemma 6 to determine the values of the derivatives of the solution of equation at the lower limit x = 0, and the partial derivatives u t x,0 and u tt x,0 of the solution of the fractional sub-diffusion equation 1 when t = 0. 4 Numerical solution of ordinary fractional differential equations The field of numerical computations for fractional differential equations has been rapidly gaining popularity for the last several decades. The fractional differential equations have diverse algorithms for computation of numerical solutions and a potential for practical applications. The analytical solutions of linear fractional differential equations with constant coefficients are determined with the integral transforms method []. The algorithms for numerical solution can be used for a much larger class of equations including fractional differential equations with non-constant coefficients. Dithelm et al. [60] proposed a prediction-correction algorithm for numerical approximation 15
16 for ordinary fractional differential equations with accuracy O h min{,1+α}. Deng [63] presented an improved prediction-correction algorithm with accuracy O h min{,1+α}. Higher order prediction-correction algorithms are discussed in [37]. While the number of computations for numerical solution of ordinary fractional differential equations is much smaller than the number of computations for partial fractional differential equations it is greater than the number of computations for ordinary differential equations. An acceptable approximation for the solution of equation is obtained from the Grünwald formula approximation for the Caputo derivative. It converges to the exact solution with accuracy Oh. An improved approximation 3 with accuracy Oh α is obtained when we use approximation 6 instead of the Grünwald formula. Numerical experiments for equation 4 and approximations and 3 on the interval [0,1] are given in Table 1. In Figure 1 the two approximations are compared with the exact solution and the second order approximation 5, when h = 0.1. In the present section we use approximations 10 and 11 to obtain recurrence relations for second and third order approximations to the solution of ordinary fractional differential equation on the interval [0, 1], y α x+yx = fx and we give a proof for the convergence of the algorithm. We can assume that equation has initial condition y0 = 0 because the function yx = yx y0 is a solution of equation with right-hand side fx = fx y0. The exact solution of equation is obtained with the Laplace transform method [] yx = x 0 ξ α 1 E α,α ξ α fx ξdξ. An alternative approach to nonzero initial condition is to use approximations 9,10 and 11 for the function yx y0 as in formulas 7 and 8. Let h = T/N where T > 0 and N is a positive integer. By approximating the fractional derivative at the point x n = nh using the Grünwald formula we obtain y n + 1 h α n =0 n =1 w α y n +y n f n, w α y n +h α y n h α f n. 16
17 Table 1: Maximum error and order of approximations and 3 for equation 4 on the interval [0,1] when α = /3. h Error Order h Error Order The truncation errors of the two approximations are Oh and Oh 1+α. We compute an approximation ỹ n to the exact solution of at the point x n with ỹ 0 = 0 and the recurrence relations ỹ n = 1 1+h α h α f n n =1 w α ỹn. Approximation has accuracy Oh. By approximating at the Caputo derivative at the points x n with 6 we obtain similar recurrence relations 1 n ỹ n = h α f h α +c α n c α ỹn. 3 0 The accuracy of approximation 3 is Oh α. We evaluate numerically the performance of approximations and 3 for the equation =1 y α x+yx = x +α +Γ3+αx. 4 Equation 4 has solution yx = x +α. When h = 0.1 and α = /3 the maximum error of approximations and 3 are and In Table 1 we compute the maximum errors and the order of approximations and 3 for α = /3 and different values of h. The graphs of approximations and 3 filled squares and empty circles and the solution of 4 on the interval [0,1] are given in Figure Second-order approximations In Corollary 5 we determined a second order approximation 10 for the Grünwald formula with average values of Caputo derivatives. We use the approximation and 17 to obtain second order numerical solutions 5 and 6 for equation. α h α α hy n = y α n α y n α +O h. 17
18 From equation we have Then α h y n = h α α y α n 1 = f n 1 y n 1, y α n = f n y n. f n 1 y n 1 +h α 1 α f n y n +O h +α. The Grünwald formula for y n is defined as α h y n = n =0 w α y n = y n + n =1 w α y n. Let γ = 1 1+h α 1 α. The exact solution of equation satisfies α y n = γ h f α n 1 y n 1 +h α 1 α f n n =1 w α y n +O h +α. We can obtain a more convenient form of the above formula using 17 α f n α = f n α f n +O h, 1 1 α h α f n α + α h α hα y n 1 w α y n +O h +α. = We compute second order approximations ỹ n to the exact solution y n of y n = 1+ equation with ỹ 0 = 0 and the recurrence relations 1 ỹ n = 1+ 1 α h α f n α + α h α hα ỹ n 1 ỹ n = γ h α 1 α n n = f n + α f n 1 + α hα ỹ n 1 w α ỹ n, 5 n = w α ỹn. 6 18
19 Theorem 7. Suppose that the solution yx of equation is sufficiently differentiable function on an interval [0, T] and y0 = y 0 = y 0 = 0. Then approximations 5 and 6 converge to the solution of equation with accuracy Oh. The proof of Theorem 7 is similar to the proof of Theorem 14. The two numerical solutions 5 and 6 have second order accuracy Oh. Experimental results suggest that in many cases the error of approximation 5 is smaller than the error of approximation 6 because it has a smaller truncation error. In Table we compute the maximum error and order of approximation 5 for equation 4 and different values of h. When h = 0.1 and α = /3 the error of approximation 5 is In Figure we compare approximations,3 and 5 with the exact solution of equation 4. Figure 1: Graph of the exact solution of equation 4 and approximations, 3 and 5 for h = 0.1 and α = / Suppose that the solution yx of equation is sufficiently differentiable function. Denote L 1 = y 0, L = y 0, L 3 = y 0. In the next lemma we determine L 1 and L from the function fx. Lemma 8. Suppose that yx is a sufficiently differentiable solution of. f0 = 0, L 1 = limf 1 α x, x 0 d L = lim x 0 dx f1 α x L 1x α 1. Γα 19
20 Table : Maximum error and order of approximation 5 for equation 4 on the interval [0,1] and α = /3. h Error Ratio log Ratio Proof. From Lemma 6 we have that Then y α 0 = 0, and y 1 α 0 = 0. 0 = y0 = f0 y α 0 = f0. By applying fractional differentiation of order 1 α to both sides of Hence, y x+y 1 α x = f 1 α x. L 1 = y 0 = f 1 α 0 y 1 α 0 = f 1 α 0 = lim x 0 f 1 α x. By differentiating equation we obtain y x+y 1+1 α x = f 1+1 α x. 7 We determine the value of the Miller-Ross derivative y 1+1 α 0 using differentiation and integration by parts. y 1+1 α x = d dx y1 α x = 1 d x y ξ Γα dx 0 x ξ 1 αdξ, Γαy 1+1 α x = d x y ξd x ξα, dx 0 α Γ1+αy 1+1 α x = d x y ξx ξ α ] x0 dx y ξx ξ α dξ, Γ1+αy 1+1 α x = d dx x y 0x α y ξx ξ α dξ 0 0 0,
21 From 7, We have that D α x y 1+1 α x = y 0x α 1 Γα y 1+1 α x = y 0x α x y ξ Γα 0 x ξ 1 αdξ, 8 Γα y x = f 1+1 α x L 1x α 1 Γα +Dx α yx. Dx α yx. y0 = 0, when y C [0,1]. Hence, L = y 0 = lim f 1+1 α x L 1x α 1 x 0 Γα. Let zx = yx L 1 x L x. 9 Lemma 9. The function zx satisfies z0 = z 0 = z 0 = 0 and is a solution of the ordinary fractional differential equation where Fx = fx L 1 x L x z α x+zx = Fx, 30 L 1 Γ α x1 α Proof. z0 = y0 = 0. By differentiating 9 L Γ3 α x α. z 0 = y 0 L 1 = 0, z 0 = y 0 L = 0. The function zx has fractional derivative of order α Then z α x = y α x L 1 L Γ α x1 α Γ3 α x α. z α x+zx = y α x+yx L 1 x L x L 1 Γ α x1 α L Γ3 α x α, z α x+zx = fx L 1 x L x L 1 Γ α x1 α L Γ3 α x α. Approximations 5 and 6 are second order numerical solutions of equation 30, because the function zx satisfies the conditions of Lemma 7. 1
22 4. Estimates for Grünwald weights In section 4.1 we determined recurrence relations 5 and 6 for second order numerical solutions of equation. In section 4.3 and section 5 we obtain a third order numerical solution of equation and second order difference approximations for the fractional sub-diffusion equation. In Theorem 7, Theorem 14 and Theorem 9 we discuss the convergence properties of the approximations. The proofs rely on the lower bound 3 for the tail of the sum of Grunwald weights 31. The Grünwald weights w n α are computed recursively as where { w α n } n=1 w α 0 = 1, w α 1 = α, w α = αα 1, α w n α = 1 n = 1 α+1 w α n n n 1, is an increasing sequence of negative numbers with sum n=1 w α n = The Grünwald weights converge to zero with an asymptotic rate [4] w α n α 1 as n. Γ1 α n 1+α We determine bounds for the Grünwald weights using the properties of the exponential function. Lemma 10. Inequalities for the exponential function i 1 x < e x for 0 < x < 1; ii 1 x > e x x for 0 < x < /3. Proof. Let hx = e x +x 1. The function hx has a positive derivative. Therefore, h x = 1 e x > 1 e 0 = 0. hx = e x +x 1 > h0 = 0. Now we prove ii. By taing logarithm from both sides ln1 x > x x.
23 The function ln1 x has Maclaurin series ln1 x = n=1 x n n = x x n=3 x n n. Inequality ii is equivalent to We have that 0 < x < /3. Then x 3 + x 4 + x3 5 + n=6 x > n=3 x n n, 1 > x 3 + x 4 + x3 5 + x n < x 3 + x 4 + x3 5 + n=6 n=6 x n n. n < n 3 n= n 6 < x n n < = < 1. n 1 n 3 n=6 n 6, 3 In the next two lemmas we determine upper and lower bounds for and w α. =n Lemma 11. Estimates for Grünwald weights π e α α1 α α n α+1 < w α α α+1 n < n+1 α+1. n=6 w n α Proof. We determine bounds for the Grünwald weights from the recursive formula and the inequalities for the exponential function. w α n = 1 α+1 w α n 1 < e α+1 n w α n 1, n w α n < e α+1 n w α n 1 < e α+1 n e α+1 n 1 w α n, w α n < e α+1 n e α+1 n 1 e α+1 w α 1 = αe α+1 n = 1. 3
24 The function 1/x is decreasing for x 0. Then and n = 1 n+1 > 1 dx = lnn+1 ln x w α n < αe α+1 n = 1 < αe α+1lnn+1 ln = αα+1 n+1 α+1. Now we use the inequality for the exponential function from Lemma 10ii to determine a lower bound for the Grünwald weights. w α n = 1 α+1 w α n 1 > e α+1 n α+1 n α w n 1, n w α n > e α+1 n α+1 n α w n 1 > e α+1 n α+1 n e α+1 n 1 α+1 n 1 α w n, w α n > e α+1 n α+1 n e α+1 n 1 α+1 α+1 n 1 e 3 α+1 3 α w, w α α1 α n > w α α1 α n > We have that [5] =3 e α+1 n =3 1 α+1 n =3 1, e α+1 n =3 1 e α+1 = = = π =3 =1 Then w α α1 α n > e α+1 n =3 1 α+1 e π The function 1/x is decreasing for x 0. Then n 1 n < 1 dx = lnn ln. x Hence, w α α1 α n > w α α+1 n > e π π e α e α+1lnn ln, α1 α α n α+1. 4
25 We use the estimates for the Grünwald weights to determine bounds for the tail of 31. Lemma 1. Bounds for sums of Grünwald weights 1 α 5 α < n Proof. From Lemma 11 we have =n w α =n π e α α1 α α < α+1 < w α =n α. 3 n < =n The function 1/x α+1 is a decreasing for x 0 and we have: Hence, =n 1 +1 < α+1 n =n 1 x α+1dx < =n 1 +1 < 1 α+1 αn < 1 α α+1. π e α α α < n α =n The number α is between zero and one. Then π e α w α =n > e 4 π > α+1, < α+1 n. α α α+1 +1 α Third-order approximation In the present section we use approximation 11 to determine recurrence relations for a third order approximation for equation, when the solution is sufficiently differentiable function. In section 4. we determined an estimate for sums of Grünwald weights. In Lemma 13 and Theorem 14 we use the lower bound to prove that the approximation converges to the solution with accuracy Oh 3, when the solution satisfies the conditions of Corollary 5. a h α α hy n = 8 5a 4 y α n + 11a 1 a 4 5 y α n a 4 + a y n α +O h 3, 8
26 a α hy n = h α 8 5a 11a f n y n +h α 4 1 a 4 We have that Let α h y n = n =0 h α 1 17a 4 + a 8 w α y n = y n + n =1 γ = 1+h α 1 17α 4 + a. 8 w α y n. f n 1 y n 1 + f n y n +O h 3+α. The solution of equation satisfies [ 11a y n = γ 1 h α 1 a a f n 1 y n 1 +h α 4 8 5a f n y n 4 + h α 1 17a ] n 4 + a f n w α 8 y n +O h 3+α. We compute a numerical solution ỹ n with ỹ 0 = ỹ 1 = 0 and [ 11a ỹ n = γ 1 h f α 1 a a n 1 ỹ n 1 +h α 4 8 5a f n ỹ n 4 + h α 1 17a ] n 4 + a f n w α 8 ỹn. 33 Let e n = y n ỹ n be the error of approximation 33. The numbers e n satisfy [ 11a e n = γ 1 h e α 1 a a n 1 +h α 4 8 5a ] n e n + w α 4 e n +A n, where A n is the truncation error of approximation 33. The numbers A n satisfy A n < Ah 3+α, where A is a constant such that =1 A > 1 αα D 3+α, 10Γ4+α =1 =1 and D 3+α = max y 3+α x. 0 x 1 6
27 We can represent the recurrence relations for the errors e n as [ 11a e n = γ 1 α h α 1 a α1 α 5a e n 1 + +h α 4 4 a e n 8 ] n w α e n +A n, =3 where 11a γ 1 = α h α 1 a 4 e n = γ 1 γ 1 e n 1 +γ e n + n γ e n +A n, 34 =3,γ = α1 α 5a +h α 4 a 8 The numbers γ n are positive, because h α < 1 and 0 < α < 1.,γ = w α 3. Lemma 13. Suppose that equation has sufficiently differentiable solution on the interval [0,T] and y0 = y 0 = y 0 = y 0 = 0. Then 10A e n < n α h 3+α α α Proof. We prove 35 by induction on n. e 0 = y 0 ỹ 0 = 0. From the Generalized Mean Value Theorem e 1 < e 1 = y 1 ỹ 1 = yh = y3+α ξ Γ4+α h3+α, y 3+α ξ Γ4+α h3+α D 3+α Γ4+α h3+α < 10A 1 α αh3+α. Suppose that 35 holds for all n n 1. n 1 e n γ γ 1 1 e n 1 +γ e n + γ e n + A n. By the induction hypothesis e n < 10A 1 α αn α h 3+α < = α αnα h 3+α 7 = 1,,n 1.
28 Then =1 e n < 10Aγ 1 n 1 h 3+α γ 1 α αnα +Ah 3+α. 36 We have that n 1 11a γ = h α 1 a 4 Hence, n 1 a γ = h α 8 17a + 4 =1 n 1 a γ = 1+h α 8 17a 4 =1 From 36 and 37, n 1 γ < γ =1 e n < 10Aγ 1 1 α αnα h 3+α =1 5a +h α 4 a n w α =1 w α =n w α =n γ w α =1 w α =n = γ h α w α =n The number γ = 1+h 1 α 17α + a satisfies 4 8 because h α < 1 and Then and e n 1 17α 4 + a α 4 + a 8 10A 1 α αnα h 3+α 1 < γ <, > > 7 4 > 0, < 1 17α a 4 < 1. 1 < γ 1 < 1 1 γ 1 8,, w α =n.. 37 w α =n +Ah 3+α. +Ah 3+α,
29 e n < 10A h 3+α 5A h 3+α 1 α αnα 1 α αnα In Lemma 1 we showed that Hence, e n < e n < w α =n 10A 1 α αnα h 3+α > 1 α 5 5A 1 α 1 α α 5 α. n 10A 1 α αnα h 3+α Ah 3+α +Ah 3+α = w α =n +Ah 3+α. α n α h 3+α +Ah 3+α, n 10A 1 α αnα h 3+α. In the next theorem we determine the order of approximation 33 on the interval [0,T]. Theorem 14. Suppose that the solution of equation satisfies y0 = y 0 = y 0 = y 0 = 0. Then approximation 33 converges to the solution with accuracy Oh 3. Proof. The point x n = nh is in the interval [0,T] when because T = Nh. From Lemma 13, e n < n N = T h 10A 1 α αnα h 3+α 10AT α e n < h 3. 1 α α 10A 1 α α T α h αh3+α, Lemma 15. Suppose that equation has a sufficiently smooth solution. d L 3 = lim x L 1x α L x α 1. x 0 dx f1 α Γα Γα 1 9
30 Proof. By differentiating equation we obtain From 8, y +1 α x = d dx y1+1 α x = d dx Using integration by parts we obtain Then y x+y 1+1 α x = f 1+1 α x, y x+y +1 α x = f +1 α x. y 0x α 1 Γα y +1 α x = y 0x α Γα 1 + y 0x α Γα Γα ThevalueofD 3 α x + 1 x y ξ Γα 0 x ξ 1 αdξ. x y +1 α x = L x α 1 + L 1x α Γα Γα 1 +D3 α x yx. 0 y ξ x ξ 1 αdξ, y x = f +1 α x L x α 1 L 1x α Γα Γα 1 D3 α x yx. y0 is zero, when yx is sufficiently differentiable function. d x L x α 1 L 1x α. dx f1 α Γα Γα 1 y 0 = lim x 0 Let zx = yx L 1 x L x L 3 6 x3. Lemma 16. The function zx satisfies z0 = z 0 = z 0 = z 0 = 0 and is a solution of the ordinary fractional differential equation where z α x+zx = Fx, 38 Fx = fx L 1 x L x L 3 6 x3 L 1x 1 α Γ α L x α Γ3 α L 3x 3 α Γ4 α. From Theorem 14, we can compute a third order numerical solution of equation 38 with recurrence relations
31 4.4 Numerical examples In Lemma 8, we showed that when equation has a sufficiently differentiable solution the function fx satisfies the following conditions. C1 f0 = 0; C The following limits exist limf 1 α x = L 1, lim x 0 x 0 d dx f1 α x L 1x α 1 = L. Γα When L 1 = L = 0 the solution yx of satisfies y 0 = y 0 = 0. The value of L 3 = y 0 is determined from the following limit. d lim x L 1x α L x α 1 = L x 0 dx f1 α 3. Γα Γα 1 InExample1weusethealgorithmsfromLemma9andLemma16tocompute second and third order numerical solutions for equation with right-hand side 39. In Example we consider equations 4 and 43 for which the conditions C1 and C are not satisfied Example 1 Consider the ordinary fractional differential equation { y α x+yx = fx, with right-hand side y0 = 0. fx = x+3x +4x 3 +6x 3+α + x1 α Γ α + 6x α Γ3 α + 39 The function fx satisfies f0 = 0, f 1 α x = xα Γ1+α + 6x1+α Γ+α + 4x+α Γ3+α + Γ4+αx +α Γ3+α 4x 3 α Γ4 α +Γ4+αx3. ++6x+1x +63+αx +α. 31
32 Then L 1 = f 1 α 0 =. Now we compute the value of L = y 0. Let d dx f1 α x = xα 1 Γα + 6xα Γ1+α + 4x1+α Γ+α + Γ4+αx1+α Γ+α +6+4x+6+α3+αx 1+α d L = lim x 0 dx f1 α x xα 1 = 6. Γα F 1 x = fx L 1 x L x L 1 Γ α x1 α L Γ3 α x α, F 1 x = 4x 3 +6x 3+α + 4x3 α Γ4 α +Γ4+αx3. The function z 1 x = yx x 3x is a solution of z α 1 x+z 1 x = 4x 3 +6x 3+α + 4x3 α Γ4 α +Γ4+αx3. 40 The solution of equation 40 is z 1 x = 4x 3 + 6x 3+α. From Theorem 7, we can compute second order numerical solutions of equation 40 with approximations 5 and 6. Experimental results for aprroximation 6 and α = 0.5 are given in Table 3 and Figure. Let d dx f1 α x = xα 1 Γα + 6xα 1 Γα + 4xα Γ1+α + Γ4+αxα Γ1+α + 4+6α+1α+α+3x α, d L 3 = lim x xα x 0 dx f1 α Γα 6xα 1 = 4. Γα 1 F x = fx L 0 L 1 x L x L 3 6 x3 L 1x 1 α Γ α L x α Γ3 α L 3x 3 α Γ4 α, F x = 6x 3+α +Γ4+αx 3. The function z x = yx 1 x 3x 4x 3 is a solution of the equation z α x+z x = 6x 3+α +Γ4+αx Equation 41 has solution z x = 6x 3+α. Experimental results for a third order numerical solution of equation 41 using recurrence relations 33 are given in Table 4 and Figure 3. 3
33 Figure : Graph of the solution of equation40 and second order approximation 6 for h = 0.05 and α = 0.5. Figure 3: Graph of the solution of equation 41 and third order approximation 33 for h = 0.05 and α = Table 3: Maximum error and order of approximation 6 for equation 40 with α = 1/4. h Error Ratio Order Table 4: Maximum error and order of approximation 33 for equation 41 with α = 3/4. h Error Ratio Order
34 Table 5: Maximum error and order of approximation 5 for equation 4. h Error Ratio log Ratio Example We compute numerical solutions for two ordinary fractional differential equations for which the conditions C1 and C for differentiable solution are not satisfied. y 0.5 x+yx = x 0.5 +Γ The solution of equation 4 is yx = x 0.5. The solution is not a continuously differentiable function, because condition C1 is not satisfied f0 = Γ We compute a numerical solution of equation 4 with recurrence relations 5. When h = 0.05 the error is The approximation converges to the solution with a very slow rate and the approximation order is smaller than one. Experimental results are given in Table 5. The following equation has solution yx = x 1.5. y 0.5 x+yx = Γ1.5x+x The solution is equation 43 has better differentiability properties than the solution of equation 4. We can expect that approximation 5 has higher accuracy for equation 43. fx = Γ1.5x+x 1.5, f0 = 0. f 0.75 x = x Γ.5 Γ1.5 x0.5, L 1 = lim x 0 f 0.75 x = 0. d dx f0.75 x = 0.5 Γ.5 + x0.5 Γ0.5x 0.5, L d = lim x 0 dx f0.75 x =. The solution of equation 43 doesn t have a continuous second derivative on the interval [0,1], because condition C is not satisfied. When h = 0.05 the error is The accuracy of approximation 5 for equation 43 is around Oh 1.31, when h > Table 6. 34
35 Table 6: Maximum error and order of approximation 5 for equation 4. h Error Ratio Order Second-order implicit difference approximations for the time fractional sub-diffusion equation In the present section we determine implicit difference approximations for the fractional sub-diffusion equation. We show that when the solution of the sub-diffusion equation is a sufficiently differentiable function the difference approximations have second order accuracy Oτ +h. The analytic solution of the fractional sub-diffusion equation can be determined using Laplace-Fourier transform [] or separation of variables for special cases of the boundary conditions and the function Gx, t. The fractional diffusion equation α ux,t = ux,t, t α x 44 u0,t = u1,t = 0,ux,0 = gx. has an analytical solution [6]. ux,t = c n E α n π t α sinnπx 45 n=1 on the domain {0 x 1,t 0},where c n = 1 0 gξ sinnπξdξ and E α is the one-parameter Mittag-Leffler function. Each term E α n π t α sinnπx is a solution of 1 and its coefficient c n is the coefficient of the Fourier sine series of the function gx. The graph of the analytical solution of 44 when gx = x x 1 is given in Figure 4. 35
36 Figure 4: Analytical solution of 44 for α = 1/, gx = x x 1 and 0 t Second order implicit difference approximation The fractional sub-diffusion equation is an important equation in fractional calculus. The implicit difference approximation which uses approximation 6 for the fractional derivative and central difference approximation for the second derivative with respect to x, has accuracy Oτ α +h [41]. Finite difference approximations are convenient way to approximate the solution of partial fractional differential equations. They combine simple description with stability and high accuracy. Even when the exact solution of the timefractional diffusion equation is available, the finite difference approximations may have higher accuracy than approximations using the exact solution. The approximation error of a numerical solution of 44 computed by truncating 45 includes errors from the truncation of the Fourier series at the endpoints and the approximations of the coefficients c n and the Mittag-Leffler functions. In this section we determine second order difference approximations for equation 1 on the domain D = [0,1] [0,T]. We can assume that equation 1 has homogeneous initial and boundary conditions α ux,t = ux,t +Gx,t, t α x 46 u0,t = u1,t = 0,ux,0 = 0. 36
37 If equation 1 is given with non-homogeneous initial or boundary conditions the substitution ux,t = ux,t ux,0 1 xu0,t u0,0 xu1,t u1,0 converts the equation to equation which has the same form and homogenous initial and boundary conditions. In Lemma 6, we showed that when the solution ux, t is sufficiently differentiable function with respect to the time variable t, the fractional derivative α ux,0 is zero when t = 0. When the t α solution ux,t has continuous second derivative u xx x,t with respect to t, the function Gx,t satisfies the condition Gx,0 = 0. Gx,0 = α ux,t t α ux,t t=0 x = 0. t=0 This compatibility condition corresponds to the condition f0 = 0 for differentiable solution of ordinary fractional differential equation, with initial condition y0 = 0. In order to construct second order difference approximations for equation 44 using approximations 9 and 10 for the Caputo derivative, the first step is to ensure that the first and second derivatives of the solution u t x,t and u tt x,t with respect to the time variable t are equal to zero when t = 0. Let L 1 x = ux,t t, L x = ux,0 t=0 t, Lx = 3 ux,0 t=0 t x. t=0 By applying time fractional derivative of order 1 α to equation1 we obtain ux, t t = 3 α ux,t t 1 α x + 1 α Gx,t t 1 α. 47 When u xx x,t is bounded in D, we have that 3 α ux,t = 0. Then t 1 α x L 1 x = 1 α Gx,t t 1 α. 48 t=0 By differentiating 47 with respect to t we obtain ux,t t = 3 α ux,t + 1 α Gx,t. 49 t t 1 α x t t 1 α The Caputo derivative of order 1 α of the function u xx x,t with respect to t is defined as 3 α ux,t t 1 α x = 1 Γα t ux,ξ t x t ξ α 1 dξ.
38 After integration by parts and differentiation with respect to t we obtain 3 α ux,t = tα 1 3 ux,t t t 1 α x Γα t x + 1 t 4 ux,ξ t=0 Γα t x t ξα 1 dξ. From the definition of Caputo derivative of order α with respect to t 0 So, 1 Γα t 0 4 ux,ξ t x t ξα 1 dξ = Dt α ux,t. x We have that 3 ux,t t x 3 α ux,t = tα 1 t t 1 α x = t α ux,t t α Γα Lx+D α t ux,t x. 50 Gx,t = α ux,t t t α t Gx,t. By integrating by parts and differentiating with respect to t Then α ux,t t t α 3 ux,t t x = ux,0 t = ux,0 t 3 ux,t t x t α Γ1 α + 1 t α Γ1 α + 1 Γ1 α The function Lx is computed as Lx = 3 ux,t t x = lim t 0 t=0 Γ1 α t 0 t 0 ux,ξ t ξ α dξ. t ux,ξ t ξ α dξ t t Gx,t, = L 1x Γ1 αt α +D1+α t ux,t t Gx,t. L1 x Γ1 αt α t Gx,t 51 because the Caputo derivative Dt 1+α ux,0 is zero when u tt x,t is bounded. The function L x is computed from 49 and 50 ux,t t ux,t Γα Lx+D α t x = tα α Gx,t, t t 1 α L x = ux,0 t Lx = lim t 0 Γαt + 1 α t 1 α Gx,t. 5 t 1 α 38
39 Let vx,t = ux,t L 1 xt L x t. The partial derivatives v t x,t and v tt x,t of the function vx,t are equal to zero when t = 0. α vx,t = α ux,t L 1xt 1 α t α t α Γ α L xt α Γ3 α, vx,t = ux,t L L x x 1 xt x t. The function vx, t is solution of the fractional sub-diffusion equation α vx,t = vx,t +Hx,t, t α x v0,t = v1,t = 0,vx,0 = 0, where 53 Hx,t = Gx,t L 1xt 1 α Γ α L xt α Γ3 α +L 1xt+ L x t. 54 In section 4 we used approximation 10 to determine second order numerical solutions 5 and 6 for equation. Now we use 10 and a central difference approximation for u xx x,t to construct implicit difference approximations 56 and 57 for equation 53. In Theorem 9 we show that the difference approximations are unconditionally stable and have second order accuracy Oτ +h. Let h = 1/N and τ = T/M where M and N are positive integers, and x n = nh, t m = mτ, v m n = vx n,t m, H m n = Hx n,t m. From approximation 10 and equation 53 α α α vx n,t m 1 hvx n,t m = + 1 α α vx n,t m +O τ t α t α α vx n,t m 1 = + 1 α vx n,t m x x α + Hn m α Hn m +O τ. By approximating the second derivatives u xx x n,t m 1 and u xx x n,t m with central difference formulas we obtain 1 m τ α =0 w α vm n +O τ +h = + α v m 1 n 1 vn m 1 +v m 1 n+1 1 α v m n 1 vn m +vm n+1 h h + α 39 H m 1 n + 1 α Hn m.
arxiv: v1 [math.na] 9 Jun 2018
Computational and Applied Mathematics manuscript No. will be inserted by the editor Asymptotic expansions and approximations for the Caputo derivative Yuri Dimitrov Radan Miryanov Venelin Todorov arxiv:806.0342v
More informationNumerical solutions of ordinary fractional differential equations with singularities
arxiv:86.3269v [mat.na] 8 Jun 28 Numerical solutions of ordinary fractional differential equations wit singularities Yuri Dimitrov, Ivan Dimov, Venelin Todorov Abstract Te solutions of fractional differential
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Implicit Schemes for the Model Problem The Crank-Nicolson scheme and θ-scheme
More informationA generalized Gronwall inequality and its application to fractional differential equations with Hadamard derivatives
A generalized Gronwall inequality and its application to fractional differential equations with Hadamard derivatives Deliang Qian Ziqing Gong Changpin Li Department of Mathematics, Shanghai University,
More information3 Applications of partial differentiation
Advanced Calculus Chapter 3 Applications of partial differentiation 37 3 Applications of partial differentiation 3.1 Stationary points Higher derivatives Let U R 2 and f : U R. The partial derivatives
More informationApplied Mathematics Letters
Applied Mathematics Letters 24 (211) 219 223 Contents lists available at ScienceDirect Applied Mathematics Letters journal homepage: www.elsevier.com/locate/aml Laplace transform and fractional differential
More informationSection Taylor and Maclaurin Series
Section.0 Taylor and Maclaurin Series Ruipeng Shen Feb 5 Taylor and Maclaurin Series Main Goal: How to find a power series representation for a smooth function us assume that a smooth function has a power
More informationDifferential equations, comprehensive exam topics and sample questions
Differential equations, comprehensive exam topics and sample questions Topics covered ODE s: Chapters -5, 7, from Elementary Differential Equations by Edwards and Penney, 6th edition.. Exact solutions
More informationarxiv: v2 [math.ap] 20 Nov 2017
ON THE MAXIMUM PRINCIPLE FOR A TIME-FRACTIONAL DIFFUSION EQUATION YURI LUCHKO AND MASAHIRO YAMAMOTO arxiv:172.7591v2 [math.ap] 2 Nov 217 Abstract. In this paper, we discuss the maximum principle for a
More informationDegree Master of Science in Mathematical Modelling and Scientific Computing Mathematical Methods I Thursday, 12th January 2012, 9:30 a.m.- 11:30 a.m.
Degree Master of Science in Mathematical Modelling and Scientific Computing Mathematical Methods I Thursday, 12th January 2012, 9:30 a.m.- 11:30 a.m. Candidates should submit answers to a maximum of four
More informationHigher order numerical methods for solving fractional differential equations
Noname manuscript No. will be inserted by the editor Higher order numerical methods for solving fractional differential equations Yubin Yan Kamal Pal Neville J Ford Received: date / Accepted: date Abstract
More informationarxiv: v1 [math.oc] 28 Mar 2011
Fractional variational calculus for nondifferentiable functions arxiv:3.546v [math.oc] 28 Mar 2 Ricardo Almeida ricardo.almeida@ua.pt Delfim F. M. Torres delfim@ua.pt Department of Mathematics, University
More informationElena Gogovcheva, Lyubomir Boyadjiev 1 Dedicated to Professor H.M. Srivastava, on the occasion of his 65th Birth Anniversary Abstract
FRACTIONAL EXTENSIONS OF JACOBI POLYNOMIALS AND GAUSS HYPERGEOMETRIC FUNCTION Elena Gogovcheva, Lyubomir Boyadjiev 1 Dedicated to Professor H.M. Srivastava, on the occasion of his 65th Birth Anniversary
More informationA FEW REMARKS ON THE SUPREMUM OF STABLE PROCESSES P. PATIE
A FEW REMARKS ON THE SUPREMUM OF STABLE PROCESSES P. PATIE Abstract. In [1], Bernyk et al. offer a power series and an integral representation for the density of S 1, the maximum up to time 1, of a regular
More informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
More informationA Survey of Numerical Methods in Fractional Calculus
A Survey of Numerical Methods in Fractional Calculus Fractional Derivatives in Mechanics: State of the Art CNAM Paris, 17 November 2006 Kai Diethelm diethelm@gns-mbh.com, k.diethelm@tu-bs.de Gesellschaft
More information) k ( 1 λ ) n k. ) n n e. k!(n k)! n
k ) λ ) k λ ) λk k! e λ ) π/!. e α + α) /α e k ) λ ) k λ )! λ k! k)! ) λ k λ k! λk e λ k! λk e λ. k! ) k λ ) k k + k k k ) k ) k e k λ e k ) k EX EX V arx) X Nα, σ ) Bp) Eα) Πλ) U, θ) X Nα, σ ) E ) X α
More informationApplied Mathematics Masters Examination Fall 2016, August 18, 1 4 pm.
Applied Mathematics Masters Examination Fall 16, August 18, 1 4 pm. Each of the fifteen numbered questions is worth points. All questions will be graded, but your score for the examination will be the
More informationMath 4263 Homework Set 1
Homework Set 1 1. Solve the following PDE/BVP 2. Solve the following PDE/BVP 2u t + 3u x = 0 u (x, 0) = sin (x) u x + e x u y = 0 u (0, y) = y 2 3. (a) Find the curves γ : t (x (t), y (t)) such that that
More informationENGI 9420 Lecture Notes 8 - PDEs Page 8.01
ENGI 940 Lecture Notes 8 - PDEs Page 8.01 8. Partial Differential Equations Partial differential equations (PDEs) are equations involving functions of more than one variable and their partial derivatives
More informationLEGENDRE POLYNOMIALS AND APPLICATIONS. We construct Legendre polynomials and apply them to solve Dirichlet problems in spherical coordinates.
LEGENDRE POLYNOMIALS AND APPLICATIONS We construct Legendre polynomials and apply them to solve Dirichlet problems in spherical coordinates.. Legendre equation: series solutions The Legendre equation is
More informationAdditional Homework Problems
Additional Homework Problems These problems supplement the ones assigned from the text. Use complete sentences whenever appropriate. Use mathematical terms appropriately. 1. What is the order of a differential
More informationAnswers to Problem Set Number MIT (Fall 2005).
Answers to Problem Set Number 5. 18.305 MIT Fall 2005). D. Margetis and R. Rosales MIT, Math. Dept., Cambridge, MA 02139). November 23, 2005. Course TA: Nikos Savva, MIT, Dept. of Mathematics, Cambridge,
More informationReview For the Final: Problem 1 Find the general solutions of the following DEs. a) x 2 y xy y 2 = 0 solution: = 0 : homogeneous equation.
Review For the Final: Problem 1 Find the general solutions of the following DEs. a) x 2 y xy y 2 = 0 solution: y y x y2 = 0 : homogeneous equation. x2 v = y dy, y = vx, and x v + x dv dx = v + v2. dx =
More informationAP Calculus Testbank (Chapter 9) (Mr. Surowski)
AP Calculus Testbank (Chapter 9) (Mr. Surowski) Part I. Multiple-Choice Questions n 1 1. The series will converge, provided that n 1+p + n + 1 (A) p > 1 (B) p > 2 (C) p >.5 (D) p 0 2. The series
More informationA finite element solution for the fractional equation
A finite element solution for the fractional equation Petra Nováčková, Tomáš Kisela Brno University of Technology, Brno, Czech Republic Abstract This contribution presents a numerical method for solving
More informationTwo special equations: Bessel s and Legendre s equations. p Fourier-Bessel and Fourier-Legendre series. p
LECTURE 1 Table of Contents Two special equations: Bessel s and Legendre s equations. p. 259-268. Fourier-Bessel and Fourier-Legendre series. p. 453-460. Boundary value problems in other coordinate system.
More informationReflected Brownian Motion
Chapter 6 Reflected Brownian Motion Often we encounter Diffusions in regions with boundary. If the process can reach the boundary from the interior in finite time with positive probability we need to decide
More informationMulti-Term Linear Fractional Nabla Difference Equations with Constant Coefficients
International Journal of Difference Equations ISSN 0973-6069, Volume 0, Number, pp. 9 06 205 http://campus.mst.edu/ijde Multi-Term Linear Fractional Nabla Difference Equations with Constant Coefficients
More informationAn Efficient Numerical Method for Solving. the Fractional Diffusion Equation
Journal of Applied Mathematics & Bioinformatics, vol.1, no.2, 2011, 1-12 ISSN: 1792-6602 (print), 1792-6939 (online) International Scientific Press, 2011 An Efficient Numerical Method for Solving the Fractional
More informationSeries Solution of Linear Ordinary Differential Equations
Series Solution of Linear Ordinary Differential Equations Department of Mathematics IIT Guwahati Aim: To study methods for determining series expansions for solutions to linear ODE with variable coefficients.
More informationHOMOTOPY PERTURBATION METHOD TO FRACTIONAL BIOLOGICAL POPULATION EQUATION. 1. Introduction
Fractional Differential Calculus Volume 1, Number 1 (211), 117 124 HOMOTOPY PERTURBATION METHOD TO FRACTIONAL BIOLOGICAL POPULATION EQUATION YANQIN LIU, ZHAOLI LI AND YUEYUN ZHANG Abstract In this paper,
More informationChapter 2 Analytical Approximation Methods
Chapter 2 Analytical Approximation Methods 2.1 Introduction As we mentioned in the previous chapter, most of the nonlinear ODEs have no explicit solutions, i.e., solutions, which are expressible in finite
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University A Model Problem and Its Difference Approximations 1-D Initial Boundary Value
More informationPower series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) a n z n. n=0
Lecture 22 Power series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) Recall a few facts about power series: a n z n This series in z is centered at z 0. Here z can
More informationDIfferential equations of fractional order have been the
Multistage Telescoping Decomposition Method for Solving Fractional Differential Equations Abdelkader Bouhassoun Abstract The application of telescoping decomposition method, developed for ordinary differential
More informationCALCULUS JIA-MING (FRANK) LIOU
CALCULUS JIA-MING (FRANK) LIOU Abstract. Contents. Power Series.. Polynomials and Formal Power Series.2. Radius of Convergence 2.3. Derivative and Antiderivative of Power Series 4.4. Power Series Expansion
More informationEquivalence of the Initialized Riemann-Liouville Derivatives and the Initialized Caputo Derivatives arxiv: v1 [math.
Equivalence of the Initialized Riemann-Liouville Derivatives and the Initialized Caputo Derivatives arxiv:1811.11537v1 [math.gm] 14 Nov 218 Jian Yuan 1, College of Mathematic and Information Science, Shandong
More informationSEVERAL RESULTS OF FRACTIONAL DERIVATIVES IN D (R + ) Chenkuan Li. Author's Copy
RESEARCH PAPER SEVERAL RESULTS OF FRACTIONAL DERIVATIVES IN D (R ) Chenkuan Li Abstract In this paper, we define fractional derivative of arbitrary complex order of the distributions concentrated on R,
More informationChapter 3a Topics in differentiation. Problems in differentiation. Problems in differentiation. LC Abueg: mathematical economics
Chapter 3a Topics in differentiation Lectures in Mathematical Economics L Cagandahan Abueg De La Salle University School of Economics Problems in differentiation Problems in differentiation Problem 1.
More informationExtended Adomian s polynomials for solving. non-linear fractional differential equations
Theoretical Mathematics & Applications, vol.5, no.2, 25, 89-4 ISSN: 792-9687 (print), 792-979 (online) Scienpress Ltd, 25 Extended Adomian s polynomials for solving non-linear fractional differential equations
More informationFractional Diffusion Theory and Applications Part II
Fractional Diffusion Theory and Applications Part II p. 1/2 Fractional Diffusion Theory and Applications Part II 22nd Canberra International Physics Summer School 28 Bruce Henry (Trevor Langlands, Peter
More informationFinite Difference Method of Fractional Parabolic Partial Differential Equations with Variable Coefficients
International Journal of Contemporary Mathematical Sciences Vol. 9, 014, no. 16, 767-776 HIKARI Ltd, www.m-hiari.com http://dx.doi.org/10.1988/ijcms.014.411118 Finite Difference Method of Fractional Parabolic
More informationA Novel Approach with Time-Splitting Spectral Technique for the Coupled Schrödinger Boussinesq Equations Involving Riesz Fractional Derivative
Commun. Theor. Phys. 68 2017 301 308 Vol. 68, No. 3, September 1, 2017 A Novel Approach with Time-Splitting Spectral Technique for the Coupled Schrödinger Boussinesq Equations Involving Riesz Fractional
More informationON THE SOLUTIONS OF NON-LINEAR TIME-FRACTIONAL GAS DYNAMIC EQUATIONS: AN ANALYTICAL APPROACH
International Journal of Pure and Applied Mathematics Volume 98 No. 4 2015, 491-502 ISSN: 1311-8080 (printed version); ISSN: 1314-3395 (on-line version) url: http://www.ijpam.eu doi: http://dx.doi.org/10.12732/ijpam.v98i4.8
More informationRepresentations of solutions of general nonlinear ODEs in singular regions
Representations of solutions of general nonlinear ODEs in singular regions Rodica D. Costin The Ohio State University O. Costin, RDC: Invent. 2001, Kyoto 2000, [...] 1 Most analytic differential equations
More informationAnalysis/Calculus Review Day 3
Analysis/Calculus Review Day 3 Arvind Saibaba arvindks@stanford.edu Institute of Computational and Mathematical Engineering Stanford University September 15, 2010 Big- Oh and Little- Oh Notation We write
More information8.7 Taylor s Inequality Math 2300 Section 005 Calculus II. f(x) = ln(1 + x) f(0) = 0
8.7 Taylor s Inequality Math 00 Section 005 Calculus II Name: ANSWER KEY Taylor s Inequality: If f (n+) is continuous and f (n+) < M between the center a and some point x, then f(x) T n (x) M x a n+ (n
More informationExercises for Chap. 2
1 Exercises for Chap..1 For the Riemann-Liouville fractional integral of the first kind of order α, 1,(a,x) f, evaluate the fractional integral if (i): f(t)= tν, (ii): f(t)= (t c) ν for some constant c,
More informationMarch Algebra 2 Question 1. March Algebra 2 Question 1
March Algebra 2 Question 1 If the statement is always true for the domain, assign that part a 3. If it is sometimes true, assign it a 2. If it is never true, assign it a 1. Your answer for this question
More informationProblem List MATH 5173 Spring, 2014
Problem List MATH 5173 Spring, 2014 The notation p/n means the problem with number n on page p of Perko. 1. 5/3 [Due Wednesday, January 15] 2. 6/5 and describe the relationship of the phase portraits [Due
More informationSome Results Based on Generalized Mittag-Leffler Function
Int. Journal of Math. Analysis, Vol. 6, 2012, no. 11, 503-508 Some Results Based on Generalized Mittag-Leffler Function Pratik V. Shah Department of Mathematics C. K. Pithawalla College of Engineering
More informationThis ODE arises in many physical systems that we shall investigate. + ( + 1)u = 0. (λ + s)x λ + s + ( + 1) a λ. (s + 1)(s + 2) a 0
Legendre equation This ODE arises in many physical systems that we shall investigate We choose We then have Substitution gives ( x 2 ) d 2 u du 2x 2 dx dx + ( + )u u x s a λ x λ a du dx λ a λ (λ + s)x
More information1 + lim. n n+1. f(x) = x + 1, x 1. and we check that f is increasing, instead. Using the quotient rule, we easily find that. 1 (x + 1) 1 x (x + 1) 2 =
Chapter 5 Sequences and series 5. Sequences Definition 5. (Sequence). A sequence is a function which is defined on the set N of natural numbers. Since such a function is uniquely determined by its values
More informationTMA4120, Matematikk 4K, Fall Date Section Topic HW Textbook problems Suppl. Answers. Sept 12 Aug 31/
TMA420, Matematikk 4K, Fall 206 LECTURE SCHEDULE AND ASSIGNMENTS Date Section Topic HW Textbook problems Suppl Answers Aug 22 6 Laplace transform 6:,7,2,2,22,23,25,26,4 A Sept 5 Aug 24/25 62-3 ODE, Heaviside
More informationTaylor and Maclaurin Series. Copyright Cengage Learning. All rights reserved.
11.10 Taylor and Maclaurin Series Copyright Cengage Learning. All rights reserved. We start by supposing that f is any function that can be represented by a power series f(x)= c 0 +c 1 (x a)+c 2 (x a)
More informationFRACTIONAL DIFFERENTIAL EQUATIONS
FRACTIONAL DIFFERENTIAL EQUATIONS An Introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of their Solution and some of their Applications by Igor Podlubny Technical University
More informationTAYLOR AND MACLAURIN SERIES
TAYLOR AND MACLAURIN SERIES. Introduction Last time, we were able to represent a certain restricted class of functions as power series. This leads us to the question: can we represent more general functions
More informationHigh Order Numerical Methods for the Riesz Derivatives and the Space Riesz Fractional Differential Equation
International Symposium on Fractional PDEs: Theory, Numerics and Applications June 3-5, 013, Salve Regina University High Order Numerical Methods for the Riesz Derivatives and the Space Riesz Fractional
More informationOn boundary value problems for fractional integro-differential equations in Banach spaces
Malaya J. Mat. 3425 54 553 On boundary value problems for fractional integro-differential equations in Banach spaces Sabri T. M. Thabet a, and Machindra B. Dhakne b a,b Department of Mathematics, Dr. Babasaheb
More informationswapneel/207
Partial differential equations Swapneel Mahajan www.math.iitb.ac.in/ swapneel/207 1 1 Power series For a real number x 0 and a sequence (a n ) of real numbers, consider the expression a n (x x 0 ) n =
More informationOrdinary Differential Equation Theory
Part I Ordinary Differential Equation Theory 1 Introductory Theory An n th order ODE for y = y(t) has the form Usually it can be written F (t, y, y,.., y (n) ) = y (n) = f(t, y, y,.., y (n 1) ) (Implicit
More informationPartial Differential Equations for Engineering Math 312, Fall 2012
Partial Differential Equations for Engineering Math 312, Fall 2012 Jens Lorenz July 17, 2012 Contents Department of Mathematics and Statistics, UNM, Albuquerque, NM 87131 1 Second Order ODEs with Constant
More informationMATH 251 Final Examination December 16, 2015 FORM A. Name: Student Number: Section:
MATH 5 Final Examination December 6, 5 FORM A Name: Student Number: Section: This exam has 7 questions for a total of 5 points. In order to obtain full credit for partial credit problems, all work must
More informationLecture 4: Ito s Stochastic Calculus and SDE. Seung Yeal Ha Dept of Mathematical Sciences Seoul National University
Lecture 4: Ito s Stochastic Calculus and SDE Seung Yeal Ha Dept of Mathematical Sciences Seoul National University 1 Preliminaries What is Calculus? Integral, Differentiation. Differentiation 2 Integral
More informationMath 113 (Calculus 2) Exam 4
Math 3 (Calculus ) Exam 4 November 0 November, 009 Sections 0, 3 7 Name Student ID Section Instructor In some cases a series may be seen to converge or diverge for more than one reason. For such problems
More information7.3 Singular points and the method of Frobenius
284 CHAPTER 7. POWER SERIES METHODS 7.3 Singular points and the method of Frobenius Note: or.5 lectures, 8.4 and 8.5 in [EP], 5.4 5.7 in [BD] While behaviour of ODEs at singular points is more complicated,
More informationCompletion Date: Monday February 11, 2008
MATH 4 (R) Winter 8 Intermediate Calculus I Solutions to Problem Set #4 Completion Date: Monday February, 8 Department of Mathematical and Statistical Sciences University of Alberta Question. [Sec..9,
More informationYunhi Cho and Young-One Kim
Bull. Korean Math. Soc. 41 (2004), No. 1, pp. 27 43 ANALYTIC PROPERTIES OF THE LIMITS OF THE EVEN AND ODD HYPERPOWER SEQUENCES Yunhi Cho Young-One Kim Dedicated to the memory of the late professor Eulyong
More informationMATH 205C: STATIONARY PHASE LEMMA
MATH 205C: STATIONARY PHASE LEMMA For ω, consider an integral of the form I(ω) = e iωf(x) u(x) dx, where u Cc (R n ) complex valued, with support in a compact set K, and f C (R n ) real valued. Thus, I(ω)
More informationLECTURE 14: REGULAR SINGULAR POINTS, EULER EQUATIONS
LECTURE 14: REGULAR SINGULAR POINTS, EULER EQUATIONS 1. Regular Singular Points During the past few lectures, we have been focusing on second order linear ODEs of the form y + p(x)y + q(x)y = g(x). Particularly,
More informationVibrating-string problem
EE-2020, Spring 2009 p. 1/30 Vibrating-string problem Newton s equation of motion, m u tt = applied forces to the segment (x, x, + x), Net force due to the tension of the string, T Sinθ 2 T Sinθ 1 T[u
More informationMathematical Economics. Lecture Notes (in extracts)
Prof. Dr. Frank Werner Faculty of Mathematics Institute of Mathematical Optimization (IMO) http://math.uni-magdeburg.de/ werner/math-ec-new.html Mathematical Economics Lecture Notes (in extracts) Winter
More informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More informationarxiv:cs/ v1 [cs.na] 23 Aug 2004
arxiv:cs/0408053v1 cs.na 23 Aug 2004 WEIGHTED AVERAGE FINITE DIFFERENCE METHODS FOR FRACTIONAL DIFFUSION EQUATIONS Santos B. Yuste,1 Departamento de Física, Universidad de Extremadura, E-0071 Badaoz, Spain
More informationSummer 2017 MATH Solution to Exercise 5
Summer 07 MATH00 Solution to Exercise 5. Find the partial derivatives of the following functions: (a (xy 5z/( + x, (b x/ x + y, (c arctan y/x, (d log((t + 3 + ts, (e sin(xy z 3, (f x α, x = (x,, x n. (a
More informationTaylor and Maclaurin Series
Taylor and Maclaurin Series MATH 211, Calculus II J. Robert Buchanan Department of Mathematics Spring 2018 Background We have seen that some power series converge. When they do, we can think of them as
More informationFinite Differences for Differential Equations 28 PART II. Finite Difference Methods for Differential Equations
Finite Differences for Differential Equations 28 PART II Finite Difference Methods for Differential Equations Finite Differences for Differential Equations 29 BOUNDARY VALUE PROBLEMS (I) Solving a TWO
More informationHomework #3 Solutions
Homework #3 Solutions Math 82, Spring 204 Instructor: Dr. Doreen De eon Exercises 2.2: 2, 3 2. Write down the heat equation (homogeneous) which corresponds to the given data. (Throughout, heat is measured
More informationLecture notes: Introduction to Partial Differential Equations
Lecture notes: Introduction to Partial Differential Equations Sergei V. Shabanov Department of Mathematics, University of Florida, Gainesville, FL 32611 USA CHAPTER 1 Classification of Partial Differential
More informationu n 2 4 u n 36 u n 1, n 1.
Exercise 1 Let (u n ) be the sequence defined by Set v n = u n 1 x+ u n and f (x) = 4 x. 1. Solve the equations f (x) = 1 and f (x) =. u 0 = 0, n Z +, u n+1 = u n + 4 u n.. Prove that if u n < 1, then
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University A Model Problem in a 2D Box Region Let us consider a model problem of parabolic
More informationRiemann Integrable Functions
Riemann Integrable Functions MATH 464/506, Real Analysis J. Robert Buchanan Department of Mathematics Summer 2007 Cauchy Criterion Theorem (Cauchy Criterion) A function f : [a, b] R belongs to R[a, b]
More informationLecture 13 - Wednesday April 29th
Lecture 13 - Wednesday April 29th jacques@ucsdedu Key words: Systems of equations, Implicit differentiation Know how to do implicit differentiation, how to use implicit and inverse function theorems 131
More informationResearch Article Almost Periodic Viscosity Solutions of Nonlinear Parabolic Equations
Hindawi Publishing Corporation Boundary Value Problems Volume 29, Article ID 873526, 15 pages doi:1.1155/29/873526 Research Article Almost Periodic Viscosity Solutions of Nonlinear Parabolic Equations
More information1 Directional Derivatives and Differentiability
Wednesday, January 18, 2012 1 Directional Derivatives and Differentiability Let E R N, let f : E R and let x 0 E. Given a direction v R N, let L be the line through x 0 in the direction v, that is, L :=
More informationLecture 32: Taylor Series and McLaurin series We saw last day that some functions are equal to a power series on part of their domain.
Lecture 32: Taylor Series and McLaurin series We saw last day that some functions are equal to a power series on part of their domain. For example f(x) = 1 1 x = 1 + x + x2 + x 3 + = ln(1 + x) = x x2 2
More informationẋ = f(x, y), ẏ = g(x, y), (x, y) D, can only have periodic solutions if (f,g) changes sign in D or if (f,g)=0in D.
4 Periodic Solutions We have shown that in the case of an autonomous equation the periodic solutions correspond with closed orbits in phase-space. Autonomous two-dimensional systems with phase-space R
More informationENGI 9420 Lecture Notes 8 - PDEs Page 8.01
ENGI 940 ecture Notes 8 - PDEs Page 8.0 8. Partial Differential Equations Partial differential equations (PDEs) are equations involving functions of more than one variable and their partial derivatives
More informationStarting from Heat Equation
Department of Applied Mathematics National Chiao Tung University Hsin-Chu 30010, TAIWAN 20th August 2009 Analytical Theory of Heat The differential equations of the propagation of heat express the most
More informationComputers and Mathematics with Applications. Fractional variational calculus for nondifferentiable functions
Computers and Mathematics with Applications 6 (2) 397 34 Contents lists available at ScienceDirect Computers and Mathematics with Applications journal homepage: www.elsevier.com/locate/camwa Fractional
More informationINFINITE SEQUENCES AND SERIES
11 INFINITE SEQUENCES AND SERIES INFINITE SEQUENCES AND SERIES In section 11.9, we were able to find power series representations for a certain restricted class of functions. INFINITE SEQUENCES AND SERIES
More informationSubdiffusion in a nonconvex polygon
Subdiffusion in a nonconvex polygon Kim Ngan Le and William McLean The University of New South Wales Bishnu Lamichhane University of Newcastle Monash Workshop on Numerical PDEs, February 2016 Outline Time-fractional
More informationMATH 425, FINAL EXAM SOLUTIONS
MATH 425, FINAL EXAM SOLUTIONS Each exercise is worth 50 points. Exercise. a The operator L is defined on smooth functions of (x, y by: Is the operator L linear? Prove your answer. L (u := arctan(xy u
More information2tdt 1 y = t2 + C y = which implies C = 1 and the solution is y = 1
Lectures - Week 11 General First Order ODEs & Numerical Methods for IVPs In general, nonlinear problems are much more difficult to solve than linear ones. Unfortunately many phenomena exhibit nonlinear
More informationA Note on the Differential Equations with Distributional Coefficients
MATEMATIKA, 24, Jilid 2, Bil. 2, hlm. 115 124 c Jabatan Matematik, UTM. A Note on the Differential Equations with Distributional Coefficients Adem Kilicman Department of Mathematics, Institute for Mathematical
More informationThe method of lines (MOL) for the diffusion equation
Chapter 1 The method of lines (MOL) for the diffusion equation The method of lines refers to an approximation of one or more partial differential equations with ordinary differential equations in just
More informationON THE NUMERICAL SOLUTION FOR THE FRACTIONAL WAVE EQUATION USING LEGENDRE PSEUDOSPECTRAL METHOD
International Journal of Pure and Applied Mathematics Volume 84 No. 4 2013, 307-319 ISSN: 1311-8080 (printed version); ISSN: 1314-3395 (on-line version) url: http://www.ijpam.eu doi: http://dx.doi.org/10.12732/ijpam.v84i4.1
More informationFinal Exam May 4, 2016
1 Math 425 / AMCS 525 Dr. DeTurck Final Exam May 4, 2016 You may use your book and notes on this exam. Show your work in the exam book. Work only the problems that correspond to the section that you prepared.
More informationHW2 Solutions. MATH 20D Fall 2013 Prof: Sun Hui TA: Zezhou Zhang (David) October 14, Checklist: Section 2.6: 1, 3, 6, 8, 10, 15, [20, 22]
HW2 Solutions MATH 20D Fall 2013 Prof: Sun Hui TA: Zezhou Zhang (David) October 14, 2013 Checklist: Section 2.6: 1, 3, 6, 8, 10, 15, [20, 22] Section 3.1: 1, 2, 3, 9, 16, 18, 20, 23 Section 3.2: 1, 2,
More information