Boundary layers in a two-point boundary value problem with fractional derivatives
|
|
- Bertram Barnett
- 5 years ago
- Views:
Transcription
1 Boundary layers in a two-point boundary value problem with fractional derivatives J.L. Gracia and M. Stynes Institute of Mathematics and Applications (IUMA) and Department of Applied Mathematics, University of Zaragoza (Spain) National University of Ireland (Cork) and Beijing Computational Science Research Center (China) Third Fractional Calculus, Probability and Non-Local Operators: Applications and Recent Developments. Bilbao Nov Research partially supported by IUMA, project MTM P and Diputación General de Aragón
2 Talk overview 1 The Caputo fractional derivative problem: Preliminaries Boundary layers in the constant coefficient problem: Laplace transform Boundary layers in the variable-coefficient problem: Maximum principle/comparison principle An example with an interior layer 2 Some comments and results for the Riemann-Liouville fractional derivative problem 3 Conclusions
3 It models, for example, anomalous diffusion in some porous media flows. Preliminaries We consider the following two-point boundary value problem D δ u(x) + b(x)u (x) + c(x)u(x) = f (x) for x (0, 1), u(0) α 0 u (0) = γ 0, u(1) + α 1 u (1) = γ 1, where 1 < δ < 2 and D δ is the Caputo fractional derivative D δ g(x) := 1 Γ(2 δ) x t=0 (x t) 1 δ g (t) dt for 0 < x 1. The constants α 0, α 1, γ 0, γ 1 and functions b, c and f are given. Assume that α 1 0 and α 0 1 δ 1.
4 Existence, and bounds on derivatives of the solution [M. Stynes & JLG, IMA J. Numer. Anal. 2015] Let b, c, f be sufficiently smooth. Assume that c 0, α 1 0 and α 0 1/(δ 1). Then the boundary value problem has a unique solution u with u C 1 [0, 1] C q+1 (0, 1] (q depends on b, c, f ), and for all x (0, 1] there exists a constant C, which is independent of x but depends on the problem data including δ, such that { C if i = 0, 1, u (i) (x) Cx δ i if i = 2, 3,..., q + 1.
5 An example Set r(x) = x δ with 1 < δ < 2. It is solution of But note that D δ u(x) = Γ(δ + 1) for x (0, 1), u(0) 1 1 δ u (0) = 0, u(1) = 1. r (x) = δx δ 1, r (x) = δ(δ 1)x δ 2 and the second derivative blows up at x = 0.
6 Existence, and bounds on derivatives of the solution [M. Stynes & JLG, IMA J. Numer. Anal. 2015] Let b, c, f be sufficiently smooth. Assume that c 0, α 1 0 and α 0 1/(δ 1). Then the boundary value problem has a unique solution u with u C 1 [0, 1] C q+1 (0, 1] (q depends on b, c, f ), and for all x (0, 1] there exists a constant C, which is independent of x but depends on the problem data including δ, such that { C u (i) if i = 0, 1, (x) Cx δ i if i = 2, 3,..., q + 1.
7 Existence, and bounds on derivatives of the solution [M. Stynes & JLG, IMA J. Numer. Anal. 2015] Let b, c, f be sufficiently smooth. Assume that c 0, α 1 0 and α 0 1/(δ 1). Then the boundary value problem has a unique solution u with u C 1 [0, 1] C q+1 (0, 1] (q depends on b, c, f ), and for all x (0, 1] there exists a constant C, which is independent of x but depends on the problem data including δ, such that { C u (i) if i = 0, 1, (x) Cx δ i if i = 2, 3,..., q + 1. but, how do the constants C i depend on δ? From now on, the constants are independent of x and δ.
8 A sample problem Recall that 1 < δ < 2. Consider the particular problem D δ u(x) + 2u (x) = 1 for x (0, 1), u(0) 1 δ 1 u (0) = 0.4, u(1) = 1.7 Exact solution of this problem can be computed using the Laplace transform
9 Exact solutions for δ = 1.6, Exact solution x Exact solution
10 Exact solutions for δ = 1.3, Exact solution x Exact solution
11 Exact solution for δ = Exact solution u x
12 A heuristic argument explaining the occurrence of a layer? (imitating an argument from singularly perturbed ODEs) Because lim δ 1 + Dδ u(x) = u (x) u (0), the ODE approaches a 1st-order ODE as δ 1 +, and then there are too many boundary conditions, so a boundary layer develops...
13 1D Convection-diffusion problems Example: εu + u = x, x (0, 1), u(0) = 0, u(1) = 1. u R (x) = x, x (0, 1), u R(0) = 0, u R (x) = x 2 / Computed solution x Solution for ε = 1, 10 1, 10 2 and 10 6
14 An example with no boundary layer This argument is incorrect as we now show by another constant-coefficient example whose exact solution can be computed.
15 An example with no boundary layer This argument is incorrect as we now show by another constant-coefficient example whose exact solution can be computed. Consider the particular problem D δ u(x) u (x) = 1 for x (0, 1), u(0) 1 δ 1 u (0) = 0.4, u(1) = 1.7
16 Solutions for δ = 1.6, Exact solution u x Exact solution u x
17 Solutions for δ = 1.2, Exact solution u x 4 Exact solution u x
18 Constant coefficient problem Consider the particular BVP D δ u(x) + bu (x) = f for x (0, 1), u(0) α 0 u (0) = γ 0, u(1) + α 1 u (1) = γ 1, where 1 < δ < 2 and the constants α 0, α 1, γ 0, γ 1, b and f are given, with α 1 0 and α 0 1/(δ 1). We solve the BVP exactly by applying a Laplace transform. The argument can be extended to the case of a variable rhs f (x). Inversion of Laplace transform solution leads naturally to Mittag-Leffler functions: E α,β (z) = k=0 z k Γ(αk + β) which is well defined for α, β, z C with Re(α) > 0.
19 Constant coefficient problem Applying L to the differential equation and observing that L{D δ u} = s δ 2 [ s 2 L{u} su(0) u (0) ] we obtain f L{u} = s 2 (s δ 1 b) + u(0) s + sδ 3 u (0) s δ 1 b f = s 2 (s δ 1 b) + γ 0 + α 0 u (0) s where we used the boundary condition at x = 0. To invert the Laplace transform, use that L { x β 1 E α,β (±ax α ) } = k!sα β s α a. + sδ 3 u (0) s δ 1 b,
20 Constant coefficient problem Thus, u(x) = γ 0 + α 0 u (0) + fx [ b + u (0) f ] b xe δ 1,2 (bx δ 1 ). Differentiating this expression, we obtain that u (x) = f [ b + u (0) f ] E δ 1,1 (bx δ 1 ). b Imposing the boundary condition u(1) + α 1 u (1) = γ 1, we get u (0) = f b + γ 1 γ 0 (1 + α 0 + α 1 )f /b α 0 + E δ 1,2 (b) + α 1 E δ 1,1 (b).
21 Exact solution of particular BVP The solution of the BVP with c 0 and b and f constant is u(x) =γ 0 + (α 0 + x)f b ( ) α 0 + xe δ 1,2 (bx δ 1 ) + γ 1 γ 0 (1 + α 0 + α 1 ) f ; α 0 + E δ 1,2 (b) + α 1 E δ 1,1 (b) b and u (x) = f b + E δ 1,1 (bx δ 1 ) α 0 + E δ 1,2 (b) + α 1 E δ 1,1 (b) ( γ 1 γ 0 (1 + α 0 + α 1 ) f b ). Recall that the layer appeared at x = 1 when δ 1 +. The only possible causes of this layer are E δ 1,n (bx δ 1 ).
22 Some estimates for E δ 1,n (t), t R Recall that 1 < δ < 2 and n is a positive integer. Elementary manipulations of the infinite series show that 1 if t 0, (n 1)! 2 E δ 1,n (t) if 0 t < 1, 1 t δe if t = 1. [Factor 1/(δ 1) here is sharp] δ 1 Conclusion: as δ 1 +, the value of E δ 1,n (t) does not blow up for t < 1 but does blow up for t = 1.
23 Some estimates for E δ 1,n (t), t R Recall that 1 < δ < 2 and n is a positive integer. Elementary manipulations of the infinite series show that 1 if t 0, (n 1)! 2 E δ 1,n (t) if 0 t < 1, 1 t δe if t = 1. [Factor 1/(δ 1) here is sharp] δ 1 Conclusion: as δ 1 +, the value of E δ 1,n (t) does not blow up for t < 1 but does blow up for t = 1. In addition, E δ 1,n ( ) is an increasing function so it blows up as δ 1 + for t > 1 but much more spectacularly... E δ 1,n (t) 1 δ 1 t(1 n)/(δ 1) exp for each fixed t > 1 as δ 1 +. (t 1/(δ 1))
24 Contour for integrating M-L function [I. Podlubny, 1999] For z G + (1, ϕ), 0 < α < 2, and πα/2 < ϕ min{π, πα}, E α,β (z) = 1 α z(1 β)/α exp(z 1/α )+ 1 ζ (1 β)/α exp(ζ 1/α ) 2απi ζ z γ(1,ϕ) where γ(1, ϕ) is the following contour in the complex plane: dζ, y (1, ) 1 (1, ) O - x (1, )
25 Unlike the case b > 1 where u (x) blows up only at x = 1 as δ 1 +, when b = 1 the derivative u (x) blows up at every point x in (0, 1] as δ approaches its limiting value! Case b > 1 [M. Stynes & JLG, Lecture Notes in Comput. Sci. Engineering, 2015] We have that lim δ 1 u (x) = + { ±, if x = 1, f b for each x (0, 1). Case b = 1 [M. Stynes & JLG, Lecture Notes in Comput. Sci. Engineering, 2015]. We have that lim δ 1 u (x) = ± for each x (0, 1]. +
26 A sample problem with b > 1 Recall the particular problem D δ u(x) + 2u (x) = 1 for x (0, 1), u(0) 1 δ 1 u (0) = 0.4, u(1) = Exact solution u x Figure: Solution for δ = 1.11
27 A sample problem with b = 1 Consider now the particular problem D δ u(x) + 1u (x) = 1 for x (0, 1), u(0) 1 δ 1 u (0) = 0.4, u(1) = Exact solution u Exact solution u x x Figure: Solution for δ = 1.1 and δ = 1.05
28 A sample problem with b = 1 Consider now the particular problem D δ u(x) + 1u (x) = 1 for x (0, 1), u(0) 1 δ 1 u (0) = 0.4, u(1) = Exact solution u Exact solution u x x Figure: Solution for δ = 1.01 and δ = 1.001
29 Variable-coefficient problem We now consider the following two-point boundary value problem D δ u(x) + b(x)u (x) + c(x)u(x) = f (x) for x (0, 1), u(0) α 0 u (0) = γ 0, u(1) + α 1 u (1) = γ 1, where 1 < δ < 2 and D δ is the Caputo fractional derivative. The constants α 0, α 1, γ 0, γ 1 and functions b, c and f are given. Assume that α 1 0 and α 0 1 δ 1. We will explain when boundary layers at x = 1 occur in the general case where the differential operator has variable coefficients and f is any smooth function.
30 Comparison principle and Barrier functions [M. Stynes & JLG, IMA J. Numer. Anal. 2015] Assume that g C 2 (0, 1], g (x) Cx θ, 0 < x 1 and b, c C[0, 1] with c(x) 0 for all x (0, 1). If D δ g + bg + cg 0 on (0, 1), g(0) α 0 g (0) 0 and g(1) + α 1 g (1) 0, where α 0 1 δ 1 and α 1 0, then g 0 on [0, 1].
31 Comparison principle and Barrier functions [M. Stynes & JLG, IMA J. Numer. Anal. 2015] Assume that g C 2 (0, 1], g (x) Cx θ, 0 < x 1 and b, c C[0, 1] with c(x) 0 for all x (0, 1). If D δ g + bg + cg 0 on (0, 1), g(0) α 0 g (0) 0 and g(1) + α 1 g (1) 0, where α 0 1 δ 1 and α 1 0, then g 0 on [0, 1]. Counterexample Set g(x) = x x x. Then D 1.2 g(x) 0 for 0 < x < 1, g(0) = g(1) = 0, but g(0.8) < 0, so the Dirichlet boundary conditions do not justify a comparison principle. Note that g(0) α 0 g (0) < 0 if α 0 > 0.
32 Comparison principle and Barrier functions [M. Stynes & JLG, IMA J. Numer. Anal. 2015] Assume that g C 2 (0, 1], g (x) Cx θ, 0 < x 1 and b, c C[0, 1] with c(x) 0 for all x (0, 1). If D δ g + bg + cg 0 on (0, 1), g(0) α 0 g (0) 0 and g(1) + α 1 g (1) 0, where α 0 1 δ 1 and α 1 0, then g 0 on [0, 1]. Comparison principle and barrier functions Thus, if ( D δ + b(x)d + c(x)i )(B ± g) 0, x (0, 1) (I α 0 D)(B ± g)(0) 0 (I + α 1 D)(B ± g)(1) 0, = g B.
33 Summary of the variable-coefficient problem In our class of problems define D δ u(x) + b(x)u (x) + c(x)u(x) = f (x) for x (0, 1), u(0) α 0 u (0) = γ 0, u(1) + α 1 u (1) = γ 1, M = max x [0,1] b(x). u u (0) u (1) c > 0 Bounded No Layer No layer (Robin) Layer (Dirichlet) c 0 M < 0 Bounded No Layer No Layer 0 M < 1 Bounded No Layer No Layer M = 1 Unbounded No Layer Layer M > 1 & b > 0 Unbounded No Layer Layer Table: [M. Stynes & JLG, Comput. Methods Appl. Math. 2015]
34 An example: Consider the problem with Then, D δ u(x) + b(x)u (x) + c(x)u(x) = f (x) for x (0, 1), u(0) 1 δ 1 u (0) = γ 0, u(1) + α 1 u (1) = γ 1, M = max b(x) > 1, b(x) > 0, c(x) 0, x [0, 1]. x [0,1] (i) u C δ 1, (ii) u (0) C, (iii) u (1) C δ 1, if α 1 > 0, (iv) (Robin) u (1) CM1/(δ 1), if α 1 = 0 (Dirichlet) δ 1
35 Example: M > 1, b > 0 and α 1 > 0 Example D δ u(x) + (x + 0.2)u (x) = 2(3 + x) for x (0, 1), u(0) 1 δ 1 u (0) = 0.4, u(1) + u (1) = 1.7. The exact solution of this problem is unknown. We compute a numerical solution with a finite difference scheme on a uniform mesh with N = This scheme uses L2 approximation for D δ [K.B Oldham, J. Spanier, 1974] upwinding for u Thus, this scheme is monotone [JLG & M. Stynes, J. Comput. Math. Appl. 2015].
36 Example: M > 1, b > 0 and α 1 > 0 Example Observe that D δ u(x) + (x + 0.2)u (x) = 2(3 + x) for x (0, 1), u(0) 1 δ 1 u (0) = 0.4, u(1) + u (1) = 1.7. M = = 1.2 > 1 and b 0.2 > 0. Table: Verifying sharpness of estimates (i) u u (1) C δ 1 C δ 1 and (iii) δ = 1.01 δ = δ = δ = (δ 1) max 0 j N u j (δ 1) u N u N 1 x N x N
37 Figure: Plot of b(x) = 8x 2 + 6x + 1: b(x) > 1 on part of the interval but b < 0 near x = 1. Example: An interior layer Example D δ u(x) + ( 8x 2 + 6x + 1)u (x) = 1 for x (0, 1), u(0) 1 δ 1 u (0) = 0.4, u(1) = b(x)= 8x 2 +6x x
38 Example: An interior layer 8 90 Computed solution x Computed solution x Figure: Test problem with δ = 1.1 (left) and δ = 1.01 (right)
39 Example: An interior layer Computed solution Computed solution x x Figure: Test problem with δ = (left) and δ = (right)
40 Riemann-Liouville fractional derivative problem Finally, we consider the following BVP D δ RL u(x) + bu (x) = f for x (0, 1), u(0) = 0, u(1) + α 1 u (1) = γ 1, with b and f constants and DRL δ is the Riemann-Liouville fractional derivative of order δ and it is defined by DRL δ g(x) = d 2 [ 1 x ] dx 2 (x t) 1 δ g(t) dt for 0 < x 1. Γ(2 δ) t=0 Existence and uniqueness of solution is discussed in [N. Kopteva, M. Stynes, submitted].
41 Riemann-Liouville fractional derivative problem The Laplace transform L of the Riemann-Liouville fractional derivative is L{D δ RL u} = sδ L{u} C 1 sc 2, with C 1 = [ ] [ ] D δ 1 RL u (0) and C 2 = D δ 2 RL u (0). taking the Laplace transform, one obtains L{u} = Using again that C 1 s(s δ 1 b) + C 2 s δ 1 b f s 2 (s δ 1 b). L { x β 1 E α,β (±ax α ) } = k!sα β s α a,
42 Riemann-Liouville fractional derivative problem we deduce that u(x) = C 1 x δ 1 E δ 1,δ (bx δ 1 ) + C 2 x δ 2 E δ 1,δ 1 (bx δ 1 ) fx δ E δ 1,δ+1 (bx δ 1 ). In this formula the constants C 1 and C 2 must be chosen to satisfy the boundary conditions. The boundary condition u(0) = 0 forces C 2 = 0. Differentiating this expression, we get u (x) = C 1 x δ 2 E δ 1,δ 1 (bx δ 1 ) fx δ 1 E δ 1,δ (bx δ 1 ), and hence, imposing that u(1) + α 1 u (1) = γ 1, C 1 = γ 1 + f [E δ 1,δ+1 (b) + α 1 E δ 1,δ (b)] E δ 1,δ (b) + α 1 E δ 1,δ 1 (b).
43 Exact solution constant-coefficient problem Consequently, E δ 1,1 (bx δ 1 ) 1 u(x) = γ 1 E δ 1,1 (b) 1 + α 1 E δ 1,0 (b) + x f ( ) 1 E δ 1,2 (bx δ 1 ) b + f ( ) Eδ 1,2 (b) 1 + α 1 [E δ 1,1 (b) 1] ( E δ 1,1(bx δ 1 ) 1) b E δ 1,1 (b) 1 + α 1 E δ 1,0 (b) and u (x) = γ 1 x 1 E δ 1,0 (bx δ 1 ) E δ 1,1 (b) 1 + α 1 E δ 1,0 (b) + f ( ) 1 E δ 1,1 (bx δ 1 ) b + f ( ) b x 1 Eδ 1,2 (b) 1 + α 1 [E δ 1,1 (b) 1] E δ 1,0 (bx δ 1 ). E δ 1,1 (b) 1 + α 1 E δ 1,0 (b) Observe that for each δ (1, 2) u (x) as x 0 +
44 Summary of the RL constant-coefficient problem Recall that D δ RL u(x) + bu (x) = f for x (0, 1), u(0) = 0, u(1) + α 1 u (1) = γ 1. u u (1) b < 1 Bounded No layer b = 1 Unbounded Layer: u (1) C δ 1 b > 1 Unbounded Layer: u (1) C (Robin BC) (δ 1) 2 Layer: u (1) Cb1/(δ 1) (δ 1) 2 (Dirichlet BC) Table: [JLG & M. Stynes, proceedings BAIL 2015]
45 An example Consider the test problem: DRL δ u(x) + u (x) = 1 for x (0, 1), u(0) = 0, u(1) = Exact solution u Exact solution u x x Figure: Exact solution for δ = 1.01 (left fig.) and δ = (right fig.).
46 Figure: Exact solution for δ = Another example Consider the test problem: DRL δ u(x) + 1.1u (x) = 1 for x (0, 1), u(0) = 0, u(1) = Exact solution u x
47 Conclusions: We considered a 2pt BVP whose leading term is a Caputo fractional derivative of order δ with 1 < δ < 2. We first considered the special case of a constant-coefficient operator and we showed that when δ is near 1, the solution of the BVP may exhibit a boundary layer at x = 1. Moving on to the general case of a variable-coefficient differential operator, we then determined conditions on the data of the problem under which boundary layers at each endpoint can occur. This analysis showed that a crucial parameter in the presence or absence of a boundary layer at x = 1 is the quantity M := max x [0,1] b(x). We gave an example, where the sign of b is variable, and we observed that the solution can also exhibit an interior layer. We finally considered a Riemann-Liouville fractional derivative 2pt BVP and we determined when a layer at x = 1 can occur.
48 References: [GS15a] [GS15b] [KS] [OS74] J.L. Gracia, M. Stynes, Central difference approximation of convection in Caputo fractional derivative two-point boundary value problems, J. Comput. Math. Appl. 273 (2015) J.L. Gracia, M. Stynes, Boundary layers in a Riemann-Liouville fractional derivative two-point boundary value problem. Accepted in Lecture Notes in Comput. Sci. Engineering (proceedings BAIL 2014). N. Kopteva, M. Stynes, Analysis and numerical solution of a Riemann-Liouville fractional derivative two-point boundary value problem. (Submitted for publication). K.B. Oldham, J. Spanier, The fractional calculus, Academic Press, 1974.
49 References: [P99] I. Podlubny, Fractional Differential Equations. An Introduction to Fractional Derivatives, Fractional Di.erential Equations, to Methods of Their Solution and some of Their Applications, Math. Sci. Eng. 198, Academic Press, San Diego, [SG15a] M. Stynes, J.L. Gracia, Blow-up of solutions and interior layers in a Caputo two-point boundary value problem. Accepted in Lecture Notes in Comput. Sci. Engineering (proceedings BAIL 2014). [SG15b] M. Stynes, J.L. Gracia, A finite difference method for a two-point boundary value problem with a Caputo fractional derivative, IMA J. Numer. Anal. 35 (2015) [SG15c] M. Stynes, J.L. Gracia, Boundary layers in a two-point boundary value problem with a Caputo fractional derivative, Comput. Methods Appl. Math. 15 (2015)
50 Thank you very much for your attention
51 Thank you very much for your attention Boundary layers in a two-point boundary value problem with a Caputo fractional derivative J.L. Gracia and M. Stynes Institute of Mathematics and Applications (IUMA) and Department of Applied Mathematics, University of Zaragoza (Spain) National University of Ireland (Cork) and Beijing Computational Science Research Center (China) Third Fractional Calculus, Probability and Non-Local Operators: Applications and Recent Developments. Bilbao Nov Research partially supported by IUMA, project MTM P and Diputación General de Aragón
Applied 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 informationA PARAMETER ROBUST NUMERICAL METHOD FOR A TWO DIMENSIONAL REACTION-DIFFUSION PROBLEM
A PARAMETER ROBUST NUMERICAL METHOD FOR A TWO DIMENSIONAL REACTION-DIFFUSION PROBLEM C. CLAVERO, J.L. GRACIA, AND E. O RIORDAN Abstract. In this paper a singularly perturbed reaction diffusion partial
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 informationInterior Layers in Singularly Perturbed Problems
Interior Layers in Singularly Perturbed Problems Eugene O Riordan Abstract To construct layer adapted meshes for a class of singularly perturbed problems, whose solutions contain boundary layers, it is
More informationDETERMINATION OF AN UNKNOWN SOURCE TERM IN A SPACE-TIME FRACTIONAL DIFFUSION EQUATION
Journal of Fractional Calculus and Applications, Vol. 6(1) Jan. 2015, pp. 83-90. ISSN: 2090-5858. http://fcag-egypt.com/journals/jfca/ DETERMINATION OF AN UNKNOWN SOURCE TERM IN A SPACE-TIME FRACTIONAL
More informationMaximum principle for the fractional diusion equations and its applications
Maximum principle for the fractional diusion equations and its applications Yuri Luchko Department of Mathematics, Physics, and Chemistry Beuth Technical University of Applied Sciences Berlin Berlin, Germany
More informationThe solutions of time and space conformable fractional heat equations with conformable Fourier transform
Acta Univ. Sapientiae, Mathematica, 7, 2 (25) 3 4 DOI:.55/ausm-25-9 The solutions of time and space conformable fractional heat equations with conformable Fourier transform Yücel Çenesiz Department of
More informationExact Solution of Some Linear Fractional Differential Equations by Laplace Transform. 1 Introduction. 2 Preliminaries and notations
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.16(213) No.1,pp.3-11 Exact Solution of Some Linear Fractional Differential Equations by Laplace Transform Saeed
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 informationTHE NATIONAL UNIVERSITY OF IRELAND, CORK COLÁISTE NA hollscoile, CORCAIGH UNIVERSITY COLLEGE, CORK. Summer Examination 2009.
OLLSCOIL NA héireann, CORCAIGH THE NATIONAL UNIVERSITY OF IRELAND, CORK COLÁISTE NA hollscoile, CORCAIGH UNIVERSITY COLLEGE, CORK Summer Examination 2009 First Engineering MA008 Calculus and Linear Algebra
More informationA COMPUTATIONAL TECHNIQUE FOR SOLVING BOUNDARY VALUE PROBLEM WITH TWO SMALL PARAMETERS
Electronic Journal of Differential Equations, Vol. 2013 (2013), No. 30, pp. 1 10. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu A COMPUTATIONAL
More informationAn Operator Theoretical Approach to Nonlocal Differential Equations
An Operator Theoretical Approach to Nonlocal Differential Equations Joshua Lee Padgett Department of Mathematics and Statistics Texas Tech University Analysis Seminar November 27, 2017 Joshua Lee Padgett
More informationFinite difference method for elliptic problems: I
Finite difference method for elliptic problems: I Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen
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 informationTime fractional Schrödinger equation
Time fractional Schrödinger equation Mark Naber a) Department of Mathematics Monroe County Community College Monroe, Michigan, 48161-9746 The Schrödinger equation is considered with the first order time
More informationA Numerical Scheme for Generalized Fractional Optimal Control Problems
Available at http://pvamuedu/aam Appl Appl Math ISSN: 1932-9466 Vol 11, Issue 2 (December 216), pp 798 814 Applications and Applied Mathematics: An International Journal (AAM) A Numerical Scheme for Generalized
More informationExistence Theory: Green s Functions
Chapter 5 Existence Theory: Green s Functions In this chapter we describe a method for constructing a Green s Function The method outlined is formal (not rigorous) When we find a solution to a PDE by constructing
More informationINTEGRAL TRANSFORMS METHOD TO SOLVE A TIME-SPACE FRACTIONAL DIFFUSION EQUATION. Abstract
INTEGRAL TRANSFORMS METHOD TO SOLVE A TIME-SPACE FRACTIONAL DIFFUSION EQUATION Yanka Nikolova 1, Lyubomir Boyadjiev 2 Abstract The method of integral transforms based on using a fractional generalization
More informationThe Laplace Transform. Background: Improper Integrals
The Laplace Transform Background: Improper Integrals Recall: Definite Integral: a, b real numbers, a b; f continuous on [a, b] b a f(x) dx 1 Improper integrals: Type I Infinite interval of integration
More informationPOSITIVE SOLUTIONS FOR BOUNDARY VALUE PROBLEM OF SINGULAR FRACTIONAL FUNCTIONAL DIFFERENTIAL EQUATION
International Journal of Pure and Applied Mathematics Volume 92 No. 2 24, 69-79 ISSN: 3-88 (printed version); ISSN: 34-3395 (on-line version) url: http://www.ijpam.eu doi: http://dx.doi.org/.2732/ijpam.v92i2.3
More informationNumerical Solution of Space-Time Fractional Convection-Diffusion Equations with Variable Coefficients Using Haar Wavelets
Copyright 22 Tech Science Press CMES, vol.89, no.6, pp.48-495, 22 Numerical Solution of Space-Time Fractional Convection-Diffusion Equations with Variable Coefficients Using Haar Wavelets Jinxia Wei, Yiming
More informationExact Solutions of Fractional-Order Biological Population Model
Commun. Theor. Phys. (Beijing China) 5 (009) pp. 99 996 c Chinese Physical Society and IOP Publishing Ltd Vol. 5 No. 6 December 15 009 Exact Solutions of Fractional-Order Biological Population Model A.M.A.
More informationAP Calculus BC Fall Final Part IIa
AP Calculus BC 18-19 Fall Final Part IIa Calculator Required Name: 1. At time t = 0, there are 120 gallons of oil in a tank. During the time interval 0 t 10 hours, oil flows into the tank at a rate of
More informationThe Laplace Transform. Background: Improper Integrals
The Laplace Transform Background: Improper Integrals Recall: Definite Integral: a, b real numbers, a b; f continuous on [a, b] b a f(x) dx 1 Improper integrals: Type I Infinite interval of integration
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 information[2] (a) Develop and describe the piecewise linear Galerkin finite element approximation of,
269 C, Vese Practice problems [1] Write the differential equation u + u = f(x, y), (x, y) Ω u = 1 (x, y) Ω 1 n + u = x (x, y) Ω 2, Ω = {(x, y) x 2 + y 2 < 1}, Ω 1 = {(x, y) x 2 + y 2 = 1, x 0}, Ω 2 = {(x,
More informationRECONSTRUCTION OF THE THERMAL CONDUCTIVITY COEFFICIENT IN THE SPACE FRACTIONAL HEAT CONDUCTION EQUATION
THERMAL SCIENCE: Year 207, Vol. 2, No. A, pp. 8-88 8 RECONSTRUCTION OF THE THERMAL CONDUCTIVITY COEFFICIENT IN THE SPACE FRACTIONAL HEAT CONDUCTION EQUATION by Rafal BROCIEK * and Damian SLOTA Institute
More informationarxiv: v1 [math.ap] 26 Mar 2013
Analytic solutions of fractional differential equations by operational methods arxiv:134.156v1 [math.ap] 26 Mar 213 Roberto Garra 1 & Federico Polito 2 (1) Dipartimento di Scienze di Base e Applicate per
More informationA coupled system of singularly perturbed semilinear reaction-diffusion equations
A coupled system of singularly perturbed semilinear reaction-diffusion equations J.L. Gracia, F.J. Lisbona, M. Madaune-Tort E. O Riordan September 9, 008 Abstract In this paper singularly perturbed semilinear
More informationFundamental Solutions and Green s functions. Simulation Methods in Acoustics
Fundamental Solutions and Green s functions Simulation Methods in Acoustics Definitions Fundamental solution The solution F (x, x 0 ) of the linear PDE L {F (x, x 0 )} = δ(x x 0 ) x R d Is called the fundamental
More informationUNIVERSITY OF MANITOBA
Question Points Score INSTRUCTIONS TO STUDENTS: This is a 6 hour examination. No extra time will be given. No texts, notes, or other aids are permitted. There are no calculators, cellphones or electronic
More informationAppendix 2 Complex Analysis
Appendix Complex Analsis This section is intended to review several important theorems and definitions from complex analsis. Analtic function Let f be a complex variable function of a complex value, i.e.
More informationMS 2001: Test 1 B Solutions
MS 2001: Test 1 B Solutions Name: Student Number: Answer all questions. Marks may be lost if necessary work is not clearly shown. Remarks by me in italics and would not be required in a test - J.P. Question
More informationResearch Article He s Variational Iteration Method for Solving Fractional Riccati Differential Equation
International Differential Equations Volume 2010, Article ID 764738, 8 pages doi:10.1155/2010/764738 Research Article He s Variational Iteration Method for Solving Fractional Riccati Differential Equation
More informationOn a Two-Point Boundary-Value Problem with Spurious Solutions
On a Two-Point Boundary-Value Problem with Spurious Solutions By C. H. Ou and R. Wong The Carrier Pearson equation εü + u = with boundary conditions u ) = u) = 0 is studied from a rigorous point of view.
More informationSolutions to Math 41 First Exam October 15, 2013
Solutions to Math 41 First Exam October 15, 2013 1. (16 points) Find each of the following its, with justification. If the it does not exist, explain why. If there is an infinite it, then explain whether
More informationTHE FUNDAMENTAL SOLUTION OF THE DISTRIBUTED ORDER FRACTIONAL WAVE EQUATION IN ONE SPACE DIMENSION IS A PROBABILITY DENSITY
THE FUNDAMENTAL SOLUTION OF THE DISTRIBUTED ORDER FRACTIONAL WAVE EQUATION IN ONE SPACE DIMENSION IS A PROBABILITY DENSITY RUDOLF GORENFLO Free University of Berlin Germany Email:gorenflo@mi.fu.berlin.de
More informationParameter-Uniform Numerical Methods for a Class of Singularly Perturbed Problems with a Neumann Boundary Condition
Parameter-Uniform Numerical Methods for a Class of Singularly Perturbed Problems with a Neumann Boundary Condition P. A. Farrell 1,A.F.Hegarty 2, J. J. H. Miller 3, E. O Riordan 4,andG.I. Shishkin 5 1
More informationON THE C-LAGUERRE FUNCTIONS
ON THE C-LAGUERRE FUNCTIONS M. Ishteva, L. Boyadjiev 2 (Submitted by... on... ) MATHEMATIQUES Fonctions Specialles This announcement refers to a fractional extension of the classical Laguerre polynomials.
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 informationCRANK-NICOLSON FINITE DIFFERENCE METHOD FOR SOLVING TIME-FRACTIONAL DIFFUSION EQUATION
Journal of Fractional Calculus and Applications, Vol. 2. Jan. 2012, No. 2, pp. 1-9. ISSN: 2090-5858. http://www.fcaj.webs.com/ CRANK-NICOLSON FINITE DIFFERENCE METHOD FOR SOLVING TIME-FRACTIONAL DIFFUSION
More informationFractional Calculus for Solving Abel s Integral Equations Using Chebyshev Polynomials
Applied Mathematical Sciences, Vol. 5, 211, no. 45, 227-2216 Fractional Calculus for Solving Abel s Integral Equations Using Chebyshev Polynomials Z. Avazzadeh, B. Shafiee and G. B. Loghmani Department
More informationMS 3011 Exercises. December 11, 2013
MS 3011 Exercises December 11, 2013 The exercises are divided into (A) easy (B) medium and (C) hard. If you are particularly interested I also have some projects at the end which will deepen your understanding
More informationIowa State University. Instructor: Alex Roitershtein Summer Homework #5. Solutions
Math 50 Iowa State University Introduction to Real Analysis Department of Mathematics Instructor: Alex Roitershtein Summer 205 Homework #5 Solutions. Let α and c be real numbers, c > 0, and f is defined
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 informationOn Local Asymptotic Stability of q-fractional Nonlinear Dynamical Systems
Available at http://pvamuedu/aam Appl Appl Math ISSN: 1932-9466 Vol 11, Issue 1 (June 2016), pp 174-183 Applications and Applied Mathematics: An International Journal (AAM) On Local Asymptotic Stability
More informationLimits at Infinity. Horizontal Asymptotes. Definition (Limits at Infinity) Horizontal Asymptotes
Limits at Infinity If a function f has a domain that is unbounded, that is, one of the endpoints of its domain is ±, we can determine the long term behavior of the function using a it at infinity. Definition
More informationThe geometric and physical interpretation of fractional order derivatives of polynomial functions
The geometric and physical interpretation of fractional order derivatives of polynomial functions M.H. Tavassoli, A. Tavassoli, M.R. Ostad Rahimi Abstract. In this paper, after a brief mention of the definitions
More informationSection 12.6: Non-homogeneous Problems
Section 12.6: Non-homogeneous Problems 1 Introduction Up to this point all the problems we have considered are we what we call homogeneous problems. This means that for an interval < x < l the problems
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 information1 Review of di erential calculus
Review of di erential calculus This chapter presents the main elements of di erential calculus needed in probability theory. Often, students taking a course on probability theory have problems with concepts
More informationINVERSE PROBLEMS FOR FRACTIONAL DIFFUSION EQUATIONS. A Dissertation LIHUA ZUO
INVERSE PROBLEMS FOR FRACTIONAL DIFFUSION EQUATIONS A Dissertation by LIHUA ZUO Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment of the requirements for the degree
More informationarxiv: v2 [math.ca] 8 Nov 2014
JOURNAL OF THE AMERICAN MATHEMATICAL SOCIETY Volume 00, Number 0, Pages 000 000 S 0894-0347(XX)0000-0 A NEW FRACTIONAL DERIVATIVE WITH CLASSICAL PROPERTIES arxiv:1410.6535v2 [math.ca] 8 Nov 2014 UDITA
More informationLaplace Transform Method Solution of Fractional Ordinary Differential Equations
P P P Faculty P Faculty University of Africa Journal of Sciences (U.A.J.S., Vol.2, 139-160) Laplace Transform Method Solution of Fractional Ordinary Differential Equations Eltayeb. A.M. Yousif P 1 Pand
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 informationIntegro-differential equations: Regularity theory and Pohozaev identities
Integro-differential equations: Regularity theory and Pohozaev identities Xavier Ros Oton Departament Matemàtica Aplicada I, Universitat Politècnica de Catalunya PhD Thesis Advisor: Xavier Cabré Xavier
More informationINTEGRATION WORKSHOP 2004 COMPLEX ANALYSIS EXERCISES
INTEGRATION WORKSHOP 2004 COMPLEX ANALYSIS EXERCISES PHILIP FOTH 1. Cauchy s Formula and Cauchy s Theorem 1. Suppose that γ is a piecewise smooth positively ( counterclockwise ) oriented simple closed
More informationMATH 126 FINAL EXAM. Name:
MATH 126 FINAL EXAM Name: Exam policies: Closed book, closed notes, no external resources, individual work. Please write your name on the exam and on each page you detach. Unless stated otherwise, you
More informationFinal: Solutions Math 118A, Fall 2013
Final: Solutions Math 118A, Fall 2013 1. [20 pts] For each of the following PDEs for u(x, y), give their order and say if they are nonlinear or linear. If they are linear, say if they are homogeneous or
More informationANALYSIS AND APPLICATION OF DIFFUSION EQUATIONS INVOLVING A NEW FRACTIONAL DERIVATIVE WITHOUT SINGULAR KERNEL
Electronic Journal of Differential Equations, Vol. 217 (217), No. 289, pp. 1 6. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ANALYSIS AND APPLICATION OF DIFFUSION EQUATIONS
More informationAbstract. 1. Introduction
Journal of Computational Mathematics Vol.28, No.2, 2010, 273 288. http://www.global-sci.org/jcm doi:10.4208/jcm.2009.10-m2870 UNIFORM SUPERCONVERGENCE OF GALERKIN METHODS FOR SINGULARLY PERTURBED PROBLEMS
More informationThe Chebyshev Collection Method for Solving Fractional Order Klein-Gordon Equation
The Chebyshev Collection Method for Solving Fractional Order Klein-Gordon Equation M. M. KHADER Faculty of Science, Benha University Department of Mathematics Benha EGYPT mohamedmbd@yahoo.com N. H. SWETLAM
More informationPartial Differential Equations
M3M3 Partial Differential Equations Solutions to problem sheet 3/4 1* (i) Show that the second order linear differential operators L and M, defined in some domain Ω R n, and given by Mφ = Lφ = j=1 j=1
More information13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs)
13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs) A prototypical problem we will discuss in detail is the 1D diffusion equation u t = Du xx < x < l, t > finite-length rod u(x,
More informationCollege, Nashik-Road, Dist. - Nashik (MS), India,
Approximate Solution of Space Fractional Partial Differential Equations and Its Applications [1] Kalyanrao Takale, [2] Manisha Datar, [3] Sharvari Kulkarni [1] Department of Mathematics, Gokhale Education
More informationMATH 409 Advanced Calculus I Lecture 10: Continuity. Properties of continuous functions.
MATH 409 Advanced Calculus I Lecture 10: Continuity. Properties of continuous functions. Continuity Definition. Given a set E R, a function f : E R, and a point c E, the function f is continuous at c if
More informationStat 451 Lecture Notes Numerical Integration
Stat 451 Lecture Notes 03 12 Numerical Integration Ryan Martin UIC www.math.uic.edu/~rgmartin 1 Based on Chapter 5 in Givens & Hoeting, and Chapters 4 & 18 of Lange 2 Updated: February 11, 2016 1 / 29
More informationMATH34032 Mid-term Test 10.00am 10.50am, 26th March 2010 Answer all six question [20% of the total mark for this course]
MATH3432: Green s Functions, Integral Equations and the Calculus of Variations 1 MATH3432 Mid-term Test 1.am 1.5am, 26th March 21 Answer all six question [2% of the total mark for this course] Qu.1 (a)
More informationWeak form of Boundary Value Problems. Simulation Methods in Acoustics
Weak form of Boundary Value Problems Simulation Methods in Acoustics Note on finite dimensional description of functions Approximation: N f (x) ˆf (x) = q j φ j (x) j=1 Residual function: r(x) = f (x)
More informationMath 321 Final Examination April 1995 Notation used in this exam: N. (1) S N (f,x) = f(t)e int dt e inx.
Math 321 Final Examination April 1995 Notation used in this exam: N 1 π (1) S N (f,x) = f(t)e int dt e inx. 2π n= N π (2) C(X, R) is the space of bounded real-valued functions on the metric space X, equipped
More informationLinear Fractionally Damped Oscillator. Mark Naber a) Department of Mathematics. Monroe County Community College. Monroe, Michigan,
Linear Fractionally Damped Oscillator Mark Naber a) Department of Mathematics Monroe County Community College Monroe, Michigan, 48161-9746 In this paper the linearly damped oscillator equation is considered
More informationA Singularly Perturbed Convection Diffusion Turning Point Problem with an Interior Layer
A Singularly Perturbed Convection Diffusion Turning Point Problem with an Interior Layer E. O Riordan, J. Quinn School of Mathematical Sciences, Dublin City University, Glasnevin, Dublin 9, Ireland. Abstract
More informationReal Analysis Math 131AH Rudin, Chapter #1. Dominique Abdi
Real Analysis Math 3AH Rudin, Chapter # Dominique Abdi.. If r is rational (r 0) and x is irrational, prove that r + x and rx are irrational. Solution. Assume the contrary, that r+x and rx are rational.
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationError Analyses of Stabilized Pseudo-Derivative MLS Methods
Department of Mathematical Sciences Revision: April 11, 2015 Error Analyses of Stabilized Pseudo-Derivative MLS Methods Don French, M. Osorio and J. Clack Moving least square (MLS) methods are a popular
More informationINTEGRATION: THE FUNDAMENTAL THEOREM OF CALCULUS MR. VELAZQUEZ AP CALCULUS
INTEGRATION: THE FUNDAMENTAL THEOREM OF CALCULUS MR. VELAZQUEZ AP CALCULUS RECALL: ANTIDERIVATIVES When we last spoke of integration, we examined a physics problem where we saw that the area under the
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 informationA Linearised Singularly Perturbed Convection- Diffusion Problem with an Interior Layer
Dublin Institute of Technology ARROW@DIT Articles School of Mathematics 201 A Linearised Singularly Perturbed Convection- Diffusion Problem with an Interior Layer Eugene O'Riordan Dublin City University,
More informationProblem set 3: Solutions Math 207B, Winter Suppose that u(x) is a non-zero solution of the eigenvalue problem. (u ) 2 dx, u 2 dx.
Problem set 3: Solutions Math 27B, Winter 216 1. Suppose that u(x) is a non-zero solution of the eigenvalue problem u = λu < x < 1, u() =, u(1) =. Show that λ = (u ) 2 dx u2 dx. Deduce that every eigenvalue
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 informationMATH 104: INTRODUCTORY ANALYSIS SPRING 2008/09 PROBLEM SET 8 SOLUTIONS
MATH 04: INTRODUCTORY ANALYSIS SPRING 008/09 PROBLEM SET 8 SOLUTIONS. Let f : R R be continuous periodic with period, i.e. f(x + ) = f(x) for all x R. Prove the following: (a) f is bounded above below
More informationDiscontinuous Galerkin methods for fractional diffusion problems
Discontinuous Galerkin methods for fractional diffusion problems Bill McLean Kassem Mustapha School of Maths and Stats, University of NSW KFUPM, Dhahran Leipzig, 7 October, 2010 Outline Sub-diffusion Equation
More informationModule 7: The Laplace Equation
Module 7: The Laplace Equation In this module, we shall study one of the most important partial differential equations in physics known as the Laplace equation 2 u = 0 in Ω R n, (1) where 2 u := n i=1
More informationHigher-order ordinary differential equations
Higher-order ordinary differential equations 1 A linear ODE of general order n has the form a n (x) dn y dx n +a n 1(x) dn 1 y dx n 1 + +a 1(x) dy dx +a 0(x)y = f(x). If f(x) = 0 then the equation is called
More informationEXISTENCE OF SOLUTIONS TO P-LAPLACE EQUATIONS WITH LOGARITHMIC NONLINEARITY
Electronic Journal of Differential Equations, Vol. 2009(2009), No. 87, pp. 1 10. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu EXISTENCE OF SOLUTIONS
More informationRiemann integral and volume are generalized to unbounded functions and sets. is an admissible set, and its volume is a Riemann integral, 1l E,
Tel Aviv University, 26 Analysis-III 9 9 Improper integral 9a Introduction....................... 9 9b Positive integrands................... 9c Special functions gamma and beta......... 4 9d Change of
More informationConformable variational iteration method
NTMSCI 5, No. 1, 172-178 (217) 172 New Trends in Mathematical Sciences http://dx.doi.org/1.2852/ntmsci.217.135 Conformable variational iteration method Omer Acan 1,2 Omer Firat 3 Yildiray Keskin 1 Galip
More information6.1 Matrices. Definition: A Matrix A is a rectangular array of the form. A 11 A 12 A 1n A 21. A 2n. A m1 A m2 A mn A 22.
61 Matrices Definition: A Matrix A is a rectangular array of the form A 11 A 12 A 1n A 21 A 22 A 2n A m1 A m2 A mn The size of A is m n, where m is the number of rows and n is the number of columns The
More informationMAGIC058 & MATH64062: Partial Differential Equations 1
MAGIC58 & MATH6462: Partial Differential Equations Section 6 Fourier transforms 6. The Fourier integral formula We have seen from section 4 that, if a function f(x) satisfies the Dirichlet conditions,
More informationMATH 409 Advanced Calculus I Lecture 16: Mean value theorem. Taylor s formula.
MATH 409 Advanced Calculus I Lecture 16: Mean value theorem. Taylor s formula. Points of local extremum Let f : E R be a function defined on a set E R. Definition. We say that f attains a local maximum
More informationTHE GAMMA FUNCTION AND THE ZETA FUNCTION
THE GAMMA FUNCTION AND THE ZETA FUNCTION PAUL DUNCAN Abstract. The Gamma Function and the Riemann Zeta Function are two special functions that are critical to the study of many different fields of mathematics.
More informationP(E t, Ω)dt, (2) 4t has an advantage with respect. to the compactly supported mollifiers, i.e., the function W (t)f satisfies a semigroup law:
Introduction Functions of bounded variation, usually denoted by BV, have had and have an important role in several problems of calculus of variations. The main features that make BV functions suitable
More informationApplied Mathematics Letters. A reproducing kernel method for solving nonlocal fractional boundary value problems
Applied Mathematics Letters 25 (2012) 818 823 Contents lists available at SciVerse ScienceDirect Applied Mathematics Letters journal homepage: www.elsevier.com/locate/aml A reproducing kernel method for
More informationKey Words: critical point, critical zero point, multiplicity, level sets Mathematics Subject Classification. 35J25; 35B38.
Critical points of solutions to a kind of linear elliptic equations in multiply connected domains Haiyun Deng 1, Hairong Liu 2, Xiaoping Yang 3 1 School of Science, Nanjing University of Science and Technology,
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 information1 Implicit Differentiation
MATH 1010E University Mathematics Lecture Notes (week 5) Martin Li 1 Implicit Differentiation Sometimes a function is defined implicitly by an equation of the form f(x, y) = 0, which we think of as a relationship
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 informationSolutions for Math 225 Assignment #4 1
Solutions for Math 225 Assignment #4 () Let B {(3, 4), (4, 5)} and C {(, ), (0, )} be two ordered bases of R 2 (a) Find the change-of-basis matrices P C B and P B C (b) Find v] B if v] C ] (c) Find v]
More informationζ (s) = s 1 s {u} [u] ζ (s) = s 0 u 1+sdu, {u} Note how the integral runs from 0 and not 1.
Problem Sheet 3. From Theorem 3. we have ζ (s) = + s s {u} u+sdu, (45) valid for Res > 0. i) Deduce that for Res >. [u] ζ (s) = s u +sdu ote the integral contains [u] in place of {u}. ii) Deduce that for
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 information