arxiv: v1 [math.na] 29 Sep 2017
|
|
- Arabella Stewart
- 5 years ago
- Views:
Transcription
1 SUPERCONVERGENCE POINTS FOR THE SPECTRAL INTERPOLATION OF RIESZ FRACTIONAL DERIVATIVES BEICHUAN DENG, ZHIMIN ZHANG, AND XUAN ZHAO arxiv:79.3v [math.na] 9 Sep 7 Abstract. In this paper, superconvergence points are located for the approximation of the Riesz derivative of order α using classical Lobatto-type polynomials when α (, ) and generalized Jacobi functions (GJF) for arbitrary α >, respectively. For the former, superconvergence points are zeros of the Riesz fractional derivative of the leading term in the truncated Legendre-Lobatto expansion. It is observed that the convergence rate for different α at the superconvergence points is at least O(N ) better than the optimal global convergence rate. Furthermore, the interpolation is generalized to the Riesz derivative of order α > with the help of GJF, which deal well with the singularities. The well-posedness, convergence and superconvergence properties are theoretically analyzed. The gain of the convergence rate at the superconvergence points is analyzed to be O(N (α+3)/ ) for α (, ) and O(N ) for α >. Finally, we apply our findings in solving model FDEs and observe that the convergence rates are indeed much better at the predicted superconvergence points. Key words. Superconvergence, Riesz fractional derivative, spectral interpolation, generalized Jacobi functions AMS subject classifications. 65N35, 65M5, 6A33, 4A5, 4A. Introduction..Over the last two decades, the theory of fractional differential equations (FDEs) has been extensively studied. Especially, the Riesz fractional derivative, which appears frequently in spatial fractional models (such as various diffusion models), has been widely studied. Some finite-difference based numerical schemes for approximating Riesz fractional derivatives and solving linear or nonlinear Riesz FDEs are presented in [37, 38, 4]. There are also works that apply finite element methods to solve Riesz FDEs; some of them concentrate on theoretical analysis, [36, 39], and some are mainly about improving algorithms, such as fast algorithms [4, 4]. On the other hand, spectral methods are promising candidates for solving FDEs since their global nature fits well with the nonlocal definition of fractional operators. Using integer-order orthogonal polynomials as basis functions, spectral methods [,,, 6, 3] help enormously with the alleviation of the memory cost for the discretization of fractional derivatives. The authors of [4, 8, 4, 5] designed suitable bases to deal with singularities, which usually appear in fractional problems. In particular, Mao, Chen, and Shen [33] proposed a spectral Petrov-Galerkin method, which is based on generalized Jacobi functions, for solving Riesz FDEs, and provided rigorous error analysis. In this work, we study the superconvergence phenomenon for some spectral interpolation of the Riesz fractional derivative. In the literature, superconvergence of the h- version finite element method has been well studied and understood, see, e.g., [3, 35], This work is supported in part by the National Natural Science Foundation of China under grants NSFC 473, NSFC 9436, and NSAF U534; the US National Science Foundation through grant DMS-494, the Natural Science Youth Foundation of Jiangsu Province of China (No.SBK666); and the Fundamental Research Funds for the Central Universities of China (No.46K49). Department of Mathematics, Wayne State University, Detroit, MI 48, USA. (beichuan.deng@wayne.edu). Corresponding author. Beijing Computational Science Research Center, Beijing 93, China (zmzhang@csrc.ac.cn); and Department of Mathematics, Wayne State University, Detroit, MI 48, USA, (zzhang@math.wayne.edu). School of Mathematics, Southeast University, Nanjing 96, China. (xuanzhao@seu.edu.cn).
2 while there have been some relatively recent superconvergence studies of polynomial spectral interpolation and spectral methods in the case of integer-order derivatives, see, e.g., [,7,8,34]. As for fractional-order derivatives, Zhao and Zhang studied the Riemann-Liouville case recently [9] and found systematically some superconvergence points in the spirit of an earlier work on integral-order spectral methods [8]. A major difficulty in the investigation of superconvergence of spectral methods for fractional problems, compared with integer-order derivatives, is the nonlocality of fractional operators and the complicated forms of fractional derivatives. The second challenge is the construction of a good basis for a spectral scheme. Given a suitable basis, one can then begin the analysis of the approximation error in order to locate the superconvergence points. One objective of this work is to consider Lobatto-type polynomial interpolants of a sufficiently smooth function, and identify those points where the values of the Riesz fractional derivative of order α are superconvergent. Note that x = ± are interpolation points for the Legendre-Lobatto interpolation; this fact guarantees that after taking the derivative of order α (, ), the global error doesn t blow up. In this case, superconvergence points are zeros of the Riesz fractional derivative of the corresponding Legendre-Lobatto polynomials. Furthermore, according to a series of numerical experiments, we observe that the convergence rate is at least O(N ) better than the global convergence rate. Comparing with the Riemann-Liouville fractional derivative, the main difficulty in studying the Riesz case is to deal with the left and right fractional derivatives simultaneously, see the definition of (.) (.3). Fractional derivatives of order α > create stronger singularities at x = ±, which restrain the use of Lobatto-type polynomials. To handle the stronger singularity, the fractional interpolation using generalized Jacobi functions (GJF) is introduced (in Section 4), and its well-posedness, convergence and superconvergence properties are theoretically analyzed. When < α <, if the given function is GJF-interpolated at zeros of the Jacobi polynomial P α, α (x), then the superconvergence points for a fractional derivative of order α are exactly the interpolation points. Moreover, the convergence rate at the superconvergence points is O(N α+3 ) and O(N ) higher than the global convergence rate for < α < and < α <, respectively. To demonstrate the usefulness of our discovery of these superconvergence points, we use the GJF as a basis to solve a model fractional differential equation by both Petrov-Galerkin and spectral collocation methods. We observe that convergence rates at the predicted superconvegrence points are indeed much better than the best possible global rates. The organization of this paper is as follows. In Section, definitions and properties of fractional derivatives, Jacobi polynomials, generalized Jacobi functions, and Gegenbauer polynomials are introduced. Section 3 is about the Legendre-Lobatto polynomial interpolation and Section 4 deals with the GJF fractional interpolation, along with numerical examples. Section 5 considers some applications of superconvergence theory. Finally, we draw some conclusions in Section 6.. Preliminaries. We begin with some basic definitions and properties. Throughout the paper, Z + denotes the set of all positive integers, N denotes the set of all nonnegative integers. M n (R) denotes the space of all of n n matrices defined on real number field.
3 Superconvergent points for Riesz Fractional Derivatives 3.. Definitions and Properties of Fractional Derivatives. First, we recall the definitions and properties of Riesz fractional derivatives. Definition.. Let γ (, ), the left and right fractional integral are defined respectively, as follows: I γ x u(x) := Γ(γ) x u(τ) dτ, x (, ], (x τ) γ xi γ u(τ) u(x) := dτ, x [, ). Γ(γ) x (τ x) γ Then for α (k, k), where k Z +, the left and right Riemann-Liouville derivatives are defined respectively by: Dx α u(x) = D k ( Ix k α u(x)), xd α u(x) = ( ) k D k ( x I k α u(x)), where D k := dk is the k-th (weak) derivative. dt k Definition.. Let γ (, ), the one dimensional Riesz potentials are defined as follows: Io γ u(x) := c sign(x τ)u(τ) Γ(γ) x τ γ dτ = c ( Ix γ x I γ )u(x) Ie γ u(x) := c u(τ) Γ(γ) x τ γ dτ = c ( Ix γ + x I γ )u(x) where sign(x) is the sign function, c = sin(πγ/), c = cos(πγ/). Then for α (k, k), we can therefore define the Riesz fractional derivative: { R D α D u(x) := k Io k α u(x) = c ( Dx α + x D α )u(x), k is odd D k Ie k α u(x) = c ( Dx α + x D α )u(x), k is even (.) Definition.3. For any positive real number α (k, k), we define: R D α o u(x) := D k I k α o u(x) = c ( D α x + ( ) k xd α )u(x), (.) and R D α e u(x) := D k I k α e u(x) = c ( D α x + ( ) k xd α )u(x). (.3) Let us recall the Leibniz rule for fractional derivatives. Lemma.4. (see [6], Chap.) Let α R +, n Z +, and α (n, n). If both f(x) and g(x) along with all their derivatives are continuous in [, ], then the Leibniz rule for the left Riemann-Liouville fractional differentiation takes the following form D α x [f(x)g(x)] = k= Γ(α + ) Γ(k + )Γ(α k + ) ( D α k x f(x))g (k) (x). (.4)
4 4 By changing variables, we can derive the Leibniz rule for the right Riemann- Liouville derivative. Lemma.5. (see [6], Chap.5) Under the same conditions as Lemma.4, the Leibniz rule for the right Riemann-Liouville differentiation takes the following form xd α [f(x)g(x)] = ( ) k Γ(α + ) Γ(k + )Γ(α k + ) ( xd α k f(x))g (k) (x). (.5) k=.. Jacobi Polynomials and Generalized Jacobi Functions..We start from the definition of Generalized Jacobi Functions and Gegenbauer polynomials. Definition.6. Let α >, the Generalized Jacobi Functions (GJF) is defined as follows: J α, α n (x) := ( x ) α P α,α (x) where P α,α n (x) is the n-th Jacobi polynomial with respect to the weight function ω(x) = ( x) α ( + x) α. Definition.7. We define the Gegenbauer polynomials by Jacobi polynomials: C n (λ) (x) := c (λ,n) P (λ,λ ) n (x) = (λ) n (λ + ) P (λ,λ ) n (x). n The Gegenbauer polynomials have the following properties: n and max x {C(λ) n (x)} { d dx C(λ) n (x) = λc (λ+) n (x), (.6) = C (λ) n () = Γ(n+λ) Γ(n+)Γ(λ) N λ Γ(λ), when λ > D λ N λ, when < λ < (.7) where D λ is a positive constant independent of n; and n+ (z) nλ ρ n+ ( + ρ ) λ, z E ρ, (.8) Γ(λ) C (λ) where E ρ is the Berstein ellipse defined in (3.3). When we set λ = α+, we have: P α, α () = Γ(N + + α ) Γ(N + )Γ( α + ) N α, and consequently C ( α+ ) () = Γ(N + + α) Γ(N + )Γ(α + ) N α, c ( α+,) N α, (.9)
5 Superconvergent points for Riesz Fractional Derivatives 5 where C ( α+ ) (x) = c ( α+,)p α, α (x). Then, the following lemma shows the connection between the GJF and Riesz fraction derivatives. Lemma.8. (see [33], Theorem ) Let α (k, k), k Z +, then we have Iν k α J α, α Γ(n + α + k) n (x) = C(k) k P α k, α k n+k (x); (.) n! moreover, for m =,,..., k, R Dν α m J α, α Γ(n + α + m) n (x) = C(k) m P α m, α m n+m (x); (.) n! where ν = o, C(k) = ( ) k, if k is odd; ν = e, C(k) = ( ) k, if k is even. Especially: R D α J α, α Γ(n + α + ) n (x) = C(k) P α, α n (x). (.) n! 3. Legendre-Lobatto Interpolation for < α <. According to Definition., the Riesz fractional derivative of order α is equivalent to the two-sided Riemann- Liouville fractional derivatives of the same order. This leads to singularities at x = ±. In order to get rid of the singularities, we consider interpolating u(x) by Lobatto-type polynomials, in particular the Legendre-Lobatto polynomials. By doing so, x = ± are zero points of multiplicity of the error u(x) u N (x) since x = ± are two of the interpolation points. This guarantees that, after taking derivatives of order α (, ), the global error is finite. In this section, we always assume that u(x) is analytic on [, ], and can be analytically extended to a certain Berstein ellipse. 3.. Interpolation of Analytic Functions. Let {x i } N i= be the (N + ) interpolation points, where x < x < < x N. Then we define ω (x) = N (x x i ). (3.) If {x i } N i= is the set of zero points of (N + ) degree Legendre-Lobatto polynomial, then i= ω (x) L N (x) L (x), (3.) where L n (x) represents the Legendre polynomial of degree n, and the right hand side is exactly the Legendre-Lobatto polynomial of degree (N + ). Suppose that u(x) is analytic on [, ], it is well known that u(x) can be analytically extended to a domain enclosed by the so-callded Berstein ellipse, with the foci ±: E ρ := {z : z = (ρeiθ + ρe iθ ), θ π}, ρ > (3.3) where i = is the imaginary unit, ρ is the sum of semimajor and semiminor axes. Then we have the following bounds for L(E ρ ), the perimeter of the ellipse, and D ρ, the shortest distance from E ρ to [, ] respectively: L(E ρ ) π(ρ + ρ ), and Dρ = (ρ + ρ ).
6 6 For convenience, we define: M u = sup z E ρ u(z). To study the superconvergent property, by introducing the Hermite s contour integral, we have the following point-wise error expression: u(x) u N (x) = ω (x) u(z) dz, x [, ]. (3.4) πi z x ω (z) The following analysis is based on this error expression. E ρ 3.. Theoretical Statements. Parallel to the conclusion in [9], we have the following theorem. Theorem 3.. Let < α <. For the interpolation using collocation points as the zeros of Legendre-Lobatto polynomials {x i } N i=, the α-th Riesz fractional derivative superconverges at {ξi α }, which satisfies R D α ω (ξ α i ) =, i =,,..., N where ω (x) is defined by (3.). Proof. The proof starts with (3.4), according to (.), (.4), (.5), x [, ], we have: D α (u(x) u N (x)) = D α ( ω (x) πi E ρ z x ) u(z) ω (z) dz = c ( Dx α + x D α )( ω (x) πi E ρ z x ) u(z) ω (z) dz = c Γ(α + ) ( Dx α m + ( ) k xd α m )ω (x) u(z) πi E ρ Γ(α m + ) (z x) m+ ω m= (z) dz = Γ(α + ) R Do α m ω (x) u(z) πi Γ(α m + ) (z x) m+ dz (3.5) ω (z) E ρ m= According to the analysis in [9], the decay of the error is dominated by the leading term: E ρ R D α ω (x) (z x) u(z) ω (z) dz. When x = ξi α, i =,,..., N, the leading term vanishes, and the remaining terms have higher convergent rates. Next, we describe a method to compute {ξi α}n i=. For < µ <, we start from D µ xl n (x) = Γ(n + ) Γ(n µ + ) ( + x) µ P µ, µ n (x), x [, ],
7 Superconvergent points for Riesz Fractional Derivatives 7 xd µ L n(x) = Γ(n + ) Γ(n µ + ) ( x) µ P µ,µ n (x), x [, ], to have: R D α (L N L )(x) = c ( D α x + x D α )(L N (x) L (x)) = c ( Dx α L N (x) + x D α L N (x) Dx α L (x) x D α L (x)) Γ(N) = sin(πγ/)γ(n µ) [( + x) µ P µ, µ N (x) + ( x) µ P µ,µ N (x)] Γ(N + ) sin(πγ/)γ(n µ + ) [( + x) µ P µ, µ (x) + ( x) µ P µ,µ (x)] Then the roots of R D α (L N L )(x) can be computed numerically Numerical Statements and Validations Numerical Validations for Superconvergent Points. In this subsection, some numerical examples are presented to show the superconvergence. We consider the function u(x) = ( + x) 9 ( x) 9. Here u(x) is interpolated at N + = zero points of ω (x), where ω (x) is the α=.9 α=.9 α=.7 α=.7 α=.5 α=.5 α=.3 α=.3 α=. α= Fig. 3.. Curves R D α (u u )(x) for different α, zeros of R D α ω (x) are highlighted by. Legendre-Lobatto polynomial of degree. We set α =.,.3,.5,.7, and.9,
8 8 respectively. Fig 3. plots R D α (u u )(x) on [, ] with different α, and the asterisks indicate the superconvergent points, i.e. the zeros of R D α ω (x), predicted by Theorem 3.. We see the errors at those points are much smaller than the global maximal error. In addition, we observe that both the global maximal error and the errors at the superconvergent points increase, when α increases Numerical Observation of Superconvergent Rates. In order to quantify the superconvergence rate, we define the following ratio: ratio = max x R D α (u u N )(x) max i N R D α (u u N )(ξ α i ) (3.6) and use u(x) = ( + x) 9 ( x) 9 as an example to plot the ratio in the log-log chart with N = 8,,, 4, 6 in Fig. 3. (for α =.,.,.,.3,.4,.5) and Fig. 3.3 (for α =.5,.6,.7,.8,.9,.99). Two lines N and N 3 are also plotted as reference slopes. 3 α=. N 3 N α=. α=. α=.3 α=.4 α= Fig. 3.. Superconvergent ratios of (3.6) for different α-derivatives. We see that at the superconvergent points, the convergent rate is at least O(N ) faster than the global rate. 4. GJF Fractional Interpolation For Arbitrary Positive α. 4.. Theoretical Statements. When α >, interpolation by the Lobatto-type polynomials, which provides zeros of multiplicity at x = ± does not work anymore, since it is not able to control the two-sided Riemann-Liouville fractional derivatives of order α >. Inspired by Lemma.8, the GJF fractional interpolation is introduced here. On the other hand, due to singularities at x = ±, some more strict conditions are required for u(x). Let α > be the order of Riesz fractional derivatives, we always assume that ( x ) α u(x) is analytic on [, ], and can be analytically extended to a Berstein ellipse with an appropriate ρ. In this section, we mainly concentrate on the analysis of the situation < α <. The conclusion can be generalized to the
9 Superconvergent points for Riesz Fractional Derivatives 9 α=.5 α=.6 α=.7 α=.8 α=.9 α=.99 N N 3 N 3 N α=.99 α=.9 α= Fig Superconvergent ratios of (3.6) for different α-derivatives. cases α >. Let s start with the definition of GJF fractional interpolation: Definition 4.. Let α (k, k) be a given positive real number, where k Z +. Define α, < α <, α = α k +, < α and k is odd α k +, < α and k is even Suppose v(x) := ( x ) α u(x) C[, ], the goal is to find where v N (x) P N [, ], such that (4.) u N (x) = ( x ) α vn (x), (4.) u N (x k ) = u(x k ), x < x < < x N, (4.3) then u N (x) is called the GJF fractional interpolant with order α In fact, (4.3) is equivalent to find v N (x) P N [, ], such that of u(x). v N (x k ) = v(x k ) = ( x k) α u(xk ), x < x < < x N. (4.4) Since v N (x) is a polynomial approximation, by (3.4), (4.), (4.3), x [, ], the error can be expressed by: u(x) u N (x) = ( x ) α (v(x) vn (x)) = ( x α ) ω (x) v(z) πi E ρ z x ω (z) dz,
10 then we have: R D α (u(x) u N (x)) = R D α [ ( x α ) ω (x) v(z) ] πi E ρ z x ω (z) dz = πi E ρ m= Γ(α + ) R Dν α m [c ( α + Γ(α m + ),) J α, α (x)] v(z) (z x) m+ ω (z) dz (4.5) where ν = o when k is odd and ν = e when k is even. If we set {x i } N i= points of P α, α (x), then we can consider to be zero ω (x) = c ( α +,) P α, α Due to the form in (.), for m N, define: φ m (x) = (x) = C α Γ(N + ) R Γ(N + + α Dν α m [c ) ( α + + (x).,) J α, α (x)], (4.6) φ m = φ m L [,]. In the following analysis, we only consider < α <, so that α = α. Then, according to (.7) and Lemma.8, when < α <, we have: φ (x) = C α+ (x), φ Γ(N + α + ) (N + )α = Γ(N + )Γ(α + ) Γ(α + ), (4.7) φ (x) = α C α N+ (x), φ D (α )/ α 3 (N + ) α, (4.8) when m, α m <, we have: where c = φ m = φ m (x) = Ix m [Io α cj α, α (x)] = Ix m φ (x), (4.9) Γ(N+) Γ(N++α) c ( α+ max { x Γ(m ),), and therefore x On the other hand, when < α <, we have: (x t) m φ (t)dt} m Γ(m) φ (x) L [,]. (4.) φ (x) = C α+ (x), φ Γ(N + α + ) (N + )α = Γ(N + )Γ(α + ) Γ(α + ), (4.) φ (x) = α C α N+ (x), φ Γ(N + α + ) (N + )α = Γ(N + 3)Γ(α + ) Γ(α + ), (4.)
11 Superconvergent points for Riesz Fractional Derivatives and φ (x) = x φ (x)dx + C = I xφ (x) + C, where C = Ie α Γ(N + ) [ Γ(N + + α) c ( α+,)j α, α ( )] Γ(N + ) = c Γ(N + + α) c c ( α+,) N α ( + t) α ( t) α C α+ (t) ( + t) α dt ( + t) α ( t) α α P, α (t)dt = c c ( α+,) N α ( α, α ) α c(, α ) with c = cos(πγ/)γ(γ), and ( α, α ) c( α, α ) represents the weighted inner product (see [3], pp 76): ( α, α ) c( α, α ) := (, P α, α (x)) ω α, α = Γ(N + α + ) ( ) m Γ(N + α + m + ) Γ(N + α + )Γ( α + ) m!(n m + )!Γ(m + 3), m= so it is a constant term that can be numerically estimated, seeing Fig : by (.9), we have: on the other hand, we have: therefore, ( α, α ) α c(, α ) N α 3, C N α 3 = o(φ ); x φ (x)dx x φ (t) dt ( + x)φ φ Γ() when m 3, α m <, then we have: so we have φ m = φ φ + o(φ ), (4.3) φ m (x) = Ix m [Ie α cj α, α (x)] = Ix m φ (x) max { x Γ(m ) x (x t) m 3 φ (t)dt} m Γ(m ) φ. (4.4)
12 - c -3 N (.-3) Fig. 4.. α =., ( α, α ) c( α, α ) (blue) decays faster than N.99 (red). - c N N Fig. 4.. α =.45, ( α, α ) c( α, α ) (blue) decays faster than N.55 (red). c N Fig α =.99, ( α, α ) c( α, α ) (blue) decays as fast as N. (red).
13 Superconvergent points for Riesz Fractional Derivatives 3 Then, when < α <, by (.8), we can continue the estimate of (4.5): R D α (u(x) u N (x)) Γ(N + + α) π N + Γ(N + + α) π N + E ρ m= m= ( α m ( ) α Γ(m + ) φm (x) m v(z) z x m+ d z C α+ (z) ) Γ(m + )φm M v L(E ρ ) (D ρ ) m+ Γ(N + + α)γ( α+ M )( + ρ ) α+ (ρ + ρ ) v Γ(N + )(N + ) α ρ We estimate the infinity sum next. When m, by (4.4), () α ρ ( + ρ ) α+ Γ( α+ ) ( ) α Γ(m + )φm m (D ρ ) m+. m= Γ(m + )φ m Γ(m + )m (D ρ ) m+ Γ(m )(D ρ ) m+ φ = (m )m m 4D ρ (D ρ ) m φ if we strictly have D ρ >, i.e. ρ > 3 +, then c ρ, d ρ >, s.t. (4.5) so we can have: d ρ D ρ <, and (m )m c ρ (d ρ ) m, Γ(m + )φ m (D ρ ) m+ φ c ρ 4D ρ ( d ρ D ρ ) m. (4.6) Therefore, by (4.6), the infinity sum in (4.5) can be estimated below: ( ) α Γ(m + )φm m (D m= ρ ) m+ = φ + α φ D ρ (D ρ ) + ( ) α Γ(m + )φm m (D m= ρ ) m+ φ + α φ D ρ (D ρ ) + c ρφ ( ) α ( d ρ ) m 4D ρ m D ρ m= = D ρ φ + α (D ρ ) φ + c ρ 4D ρ ( + d ρ D ρ ) α φ (4.7) Similarly, when < α <, by (4.9), (4.), (4.6), we have: R D α (u(x) u N (x)) D ρ φ + c ρ D ρ ( + d ρ D ρ ) α φ. (4.8) To sum up the analysis above, we have the following theorem. Theorem 4.. Let u(x) be a function such that ( x ) α u(x) is analytic on and within the complex ellipse E ρ, where ρ > 3 +, and u N (x) be the GJF fractional interpolant of u(x) at {ξj α}n j=, which are zero points of P α, α (x). Then, for < α <, we obtain the following global error estimation: max x R D α (u u N )(x) Γ( α+ M )( + ρ ) α+ (ρ + ρ ) v Γ(α + )(ρ + ρ (N + ) 3α+ ρ () (4.9) )
14 4 and the error estimation at the superconvergent points {ξ α j }N j= : max j N R D α (u u N )(ξj α ) Γ( α+ c M )( + ρ ) α+ (ρ + ρ ) v Γ(α + )(ρ + ρ (N + ) α ρ (), (4.) ) where M v = sup z Eρ {( z ) α u(z)}, c = cρd (α )/ (α ) ( + dρ D ρ ) α. For < α <, the global error estimation and the error estimation at the superconvergent points {ξj α}n j= are given respectively in the following: max x R D α (u u N )(x) Γ( α+ M )( + ρ ) α+ (ρ + ρ ) v Γ(α + )(ρ + ρ (N + ) 3α+ ρ () ; (4.) ) where c = α D ρ max j N R D α (u u N )(ξj α ) Γ( α+ c M )( + ρ ) α+ (ρ + ρ ) v Γ(α + )(ρ + ρ (N + ) 3α 3 ρ (), (4.) ) + cρ ( + dρ D ρ ) α. Proof. We can derive (4.9) from (4.8), (4.7), (4.8) and derive (4.) from (4.7), (4.), (4.), (4.3). As for (4.) and (4.), we have, so the first term in (4.7) vanishes, i.e. φ (ξ α j ) = P α, α (ξα j ) =, j =,..., N, max j N R D α (u u N )(ξj α α ) c ρ,n ( (D ρ ) φ + c ρ ( + d ρ ) α φ ), 4D ρ D ρ where c ρ,n = M v Γ(N++α)Γ( α+ )(+ρ ) α+ (ρ +ρ ) Γ(N+)() α ρ, which establishes (4.). Similarly, (4.) follows from (4.8), (4.), and (4.). When α >, since R D α = D α α ( R D α ) where α α is an even integer, we can generalize the results. Corollary 4.3. Let α (k, k), k Z +, k N, and α be defined in (4.). Under the same assumptions in Theorem 4., we have: when k is odd, max x R D α (u u N )(x) (N k + ) 3α +(k ) ρ (N k+), (4.3) +
15 Superconvergent points for Riesz Fractional Derivatives 5 the superconvergent points {ξi α}n k+ i= are the zero points of P α at those points, we have: when k is even, α +k, +k N k+ (x), and max i N R D α (u u N )(ξi α ) (N k + ) 3α +(k ) ρ (N k+) ; (4.4) max x R D α (u u N )(x) (N k + 3) 3α +(k ) ρ (N k+3), (4.5) the superconvergent points {ξi α}n k+ i= are the zero points of P α at those points, we have: + + α +k, +k N k+3 (x), and max i N R D α (u u N )(ξi α ) (N k + ) 3α +(k 3) ρ (N k+3), (4.6) where the constants only depend on α, ρ. Proof. When k is odd, by (.), (4.6), and Lemma.8, φ m (x) = D k Γ(N + ) Γ(N + + α ) From (.6), the leading term is R D α m o φ (x) = (α + ) (α + k)c ( α + and the second term is +k ) N k+ φ (x) = (α + ) (α 3 + k)c ( α + +k ) N k+ + [c (α,)j α, α (x)] (x), φ N α++(k ) (x), φ N α++(k ) and φ m = o(φ ), for m. Therefore, the leading term vanishes at {ξi α}n k+ i=, zero points of P α α +k, +k N k+ (x). Similar to the proof of Theorem 4., the estimates (4.3) and (4.4) are derived from (4.8). When k is even, φ m (x) = D k Γ(N + ) Γ(N + + α ) R D α m e [c (α,)j α, α and the rest of the proof is similar to the case when k is odd. (x)], 4.. Numerical Validations. To make sure v(x) is smooth enough, in this subsection, we consider the function: u(x) = ( x ) α + (x + 3) where v(x) = +(x+3), and we set α =.4,.7 respectively. It s easy to see that v(x) has two simple poles at z = 3 ± i in the complex plane. Hence it is analytic within the Berstein ellipse with ρ > 3 +. In the numerical example, we set N = so v N (x) is interpolated at zero points {ξj α} j= of P α, α (x). The true solution of R D α u(x) is approximated by the sum of 4 terms. Fig. 4.4 and 4.5 depict graphs of R D α (u u )(x), where u N (x) is the GJF fractional interpolation, and α =.4 and.7, respectively. According to Theorem 4., the interpolation points are predicted as superconvergent points. Similar with Fig. 3., the errors at those superconvergent points are significantly less than the global maximal error.
16 Fig Curves R D α (u u )(x) for the GJF interpolation at pints: α = Fig Curves R D α (u u )(x) for the GJF interpolation at pints: α =.7 5. Applications. In this section, we focus on applications of superconvergence. Let < α <, and we consider the following FDE: { R D α u(x) + u(x) = f(x), x (, ) (5.) u( ) = u() = We provide two methods to solve for the equation: Petrov-Galerkin method and spectral collocation method. Our goal is to observe superconvergence phenomenon in
17 Superconvergent points for Riesz Fractional Derivatives 7 numerical solutions. In the following numerical examples, we set f(x) be the function such that u(x) = ( x ) α +.5x is the true solution. Then we demonstrate the error curve R D α (u u N ) and highlight, by, its value at the superconvergent points predicted in Theorem Petrov-Galerkin Method. For any given < α <, we are looking for such that v P N [, ], we have: u N S α = span{j α, α,, J α, α N }, ( R D α u N, v) ω α, α + (u N, v) ω α α = (f, v) ω α, α. (5.) According to (.), by setting v = P α, α i, i =,,..., N, (5.) is equivalent to find (c, c,..., c N ) T R, such that, for i =,,..., N, N j= c j [d j (P α, α j, P α, α i ) ω α, α + (P α, α j, P α, α i ) ω α,α] = (f, P α, α i ) ω α, α, (5.3) where d j = Γ(j++α) Γ(j+). We observe that the stiffness matrix is diagonal and dominates the system α=.84 α=.84 α=.7 α= Fig. 5.. Curves R D α (u u 3 )(x) for the Petrov-Galerkin method: α =.7 and.84 We plotted error curves R D α (u u N ) in Figures for α =.7,.84, N = 3, 7, respectively. According to Theorem 4., the superconvergence points are predicted to be zeros of P α, α (x). We observe that errors at those points are much smaller than the global maximal error, and moreover, both the global maximal error and errors at the superconvergence points increase, when α increases. Fig. 5.3 depicted the reciprocal of (3.6), for α =.,.5,.9, respectively, where O(N ) is plotted as a reference slope (Since they are too close to each other, only three α cases are shown in Fig. 5.3). We see that, at the superconvergence points predicted by Theorem 4., the convergence rate is O(N ) faster than the optimal global rate.
18 8-7 3 α=.84 α=.84 α=.7 α= Fig. 5.. Curves R D α (u u 7 )(x) for the Petrov-Galerkin method: α =.7 and α=.9 α=.5 α=. α=. α=.5 α=.9 N N Fig Superconvergent ratios of the Petrov-Galerkin method for different α-derivatives. 5.. Spectral Collocation Method. For any given < α <, according to the definition of GJF fractional interpolation with order α, we have: N N u N (x) = ˆl j (x)v j := ( x ) α lj (x)v j, (5.4) j= j=
19 Superconvergent points for Riesz Fractional Derivatives 9 where l j P N [, ] is the Lagrange basis function satisfying l j (x i ) = δ ij, i, j =,,..., N. -5 α=.84 α=.84 α=.7 α= Fig Curves R D α (u u 3 )(x) for the collocation method: α =.7 and.84-7 α=.84 α=.84 α=.7 α= Fig Curves R D α (u u 7 )(x) for the collocation method: α =.7 and.84 Therefore, we are looking for V = (v, v,, v N ) T R, such that ( R D α u N )(x i ) + u N (x i ) = f(x i ), i =,,..., N, (5.5)
20 -3 α=. α=.3 α=.5 α=.7 α=.9 N -3 N -4 N -3 N -4-4 α=. α=.3 α=.5 α=.7-5 α= Fig Superconvergent ratios of the collocation method for different α-derivatives. where u N (x i ) = ( x i ) α v i. Then, (5.5) is equivalent to solve the linear system: (D + Λ)V = F where D, Λ M (R), F R, and for i, j =,,..., N, D(i, j) = ( R D α ˆLi )(x j ), Λ(i, i) = ( x i ) α, F (i) = f(xi ), (5.6) and the differential matrix D can be analytically calculated by (.). Here, we set {x i } N i= to be zeros of P α, α (x). Error curves R D α (u u N ) are plotted in Figures for α =.7,.84, N = 3, 7, respectively. As predicted by Theorem 4., the superconvergence points {x i } N i= are zeros of P α, α (x). We can see the errors at those points are significantly smaller than the global maximal error. Furthermore, the performance of superconvergence points of the collocation method is much better than that of the Petrov-Galerkin method. Errors at superconvergence points for the collocation method are more closer to zeros than the Petrov-Galerkin case as demonstrated by Fig. 5.6, the reciprocal of (3.6) ratios with O(N 3 ) and O(N 4 ) as reference slopes. We observe that the convergence rate at superconvergence points for the collocation method is about O(N 3 ) better than the optimal global rate. One possible reason is that the interpolation points and superconvergence points of α th Riesz derivative are identical. 6. Concluding Remarks. In this work, we investigated superconvergence for u u N under Riesz fractional derivatives. We identified superconvergence points and found the improved convergence rate at those points. When < α <, we consider u N (x) as either polynomial interpolation or GJF fractional interpolation, the improvement in convergence rates are O(N ) and O(N α+3 ), respectively. When α >, only the GJF fractional interpolant is discussed due to the singularity, and the improvement in the convergence rate is O(N ). In particular, when < α <, for
21 Superconvergent points for Riesz Fractional Derivatives the case of GJF fractional interpolation, the superconvergence points are the same as the interpolation points. In addition, when we apply our superconvergence knowledge to the numerical solution of model FDEs, our theory predicts accurately the locations of superconvergence points. Moreover, we notice that for the Petrov-Galerkin method, the convergence improvement at the superconvergence points is only O(N ), which is inferior to O(N ), the improvement for the interpolation; and for the spectral collocation method, the convergence improvement at the superconvergence points is O(N 3 ) to O(N 4 ), which is superior to O(N ), the improvement for the interpolation. It seems that polynomial-based interpolation plays only a limited role in solving FDEs. We believe that GJF-type fractional interpolation is to be preferred in fractional calculus. We hope that our findings can be useful in numerically solving FDEs, especially when using data at the predicted superconvergence points. REFERENCES [] R. Askey, Orthogonal Polynomials and Special Functions, SIAM, Philadelphia, 975. [] S. N. Bernstein, Sur l ordre de la meilleure approximation des foncions continues par des polynomes de degré donné, Mém. Publ. Class Sci. Acad. Belgique (), 4 (9), pp. C3. [3] F. Chen, Q. Xu and J.S. Hesthaven, A multi-domain spectral method for time-fractional differential equations, J. Comput. Phys., 93 (5), pp [4] S. Chen, J. Shen and L-L. Wang, Generalized Jacobi functions and their applications to fractional differential equations, to appear in Math. Comp. arxiv: v [5] P. J. Davis, Interpolation and Approximation, Dover, New York, 975. [6] S.G. Samko, A.A. Kilbas, O.I. Marichev, Fractional Integrals and Derivatives, Theory and Applications, Gordon and Breach Science Publishers, 993. [7] L. Fatone and D. Funaro, Optimal Collocation Nodes for Fractional Derivative Operators, SIAM J. Sci. Comput., 37 (5), pp. A54-A54. [8] C. Huang, Q. Zhou and Z. Zhang, Spectral method for substantial fractional differential equations, arxiv: v [9] M. Ishteva, L. Boyadjiev and R. Scherer, On the Caputo operator of fractional calculus and C-Laguerre functions, Math. Sci. Res., 9 (5) pp.6-7. [] C. P. Li, F. H. Zeng, and F. Liu, Spectral approximations to the fractional integral and derivative, Frac. Calc. Appl. Anal. 5 () [] X. Li and C. J. Xu,A space-time spectral method for the time fractional diffusion equation, SIAM J. Numer. Anal. 47 (9) 8-3. [] B. B. Mandelbrot and J. W. Van Ness, Fractional Brownian Motions, Fractional Noises and Applications SIAM Review, (968), pp [3] M. M. Meerschaert, D. Benson, B. Baeumer, Operator Lévy motion and multiscaling anomalous diffusion, Phys. Rev. E 63 () -7. [4] R. Metzler and J. Klafter, The random walk s guide to anomalous diffusion: a fractional dynamics approach, Phys. Reports, 339 () -77. [5] K. Mustapha and W. McLean, Uniform convergence for a discontinuous Galerkin, timestepping method applied to a fractional diffusion equation, IMA J. Numer. Anal., 3 (), pp [6] I. Podlubny, Fractional differential equations, Academic Press, New York, 999. [7] S. Shen, F. Liu, V. Anh and I. Turner, The fundamental solution and numerical solution of the Riesz fractional advection-dispersion equation, IMA J. Appl. Math., 73 (8), pp [8] H. G. Sun, W. Chen and Y. Q. Chen, Variable-order fractional differential operators in anomalous diffusion modeling, Phys. A 388 (9) [9] 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., (5) 35, [] H. Wang and N. Du, Fast alternating-direction finite difference methods for three-dimensional space-fractional diffusion equations, J. Comput. Phys., 58 (3), pp [] L-L. Wang, X. D. Zhao and Z. Zhang, Superconvergence of Jacobi-Gauss-type spectral interpolation, J. Sci. Comput., 59 (4), pp [] Q. Xu, J. S. Hesthaven, Stable multi-domain spectral penalty methods for fractional partial
22 differential equations, J. Comput. Phys. 57 (4) [3] Q. Yang, I. Turner, F. Liu and M. Ilić, Novel numerical methods for solving the time-space fractional diffusion equation in two dimensions, SIAM J. Sci. Comput., 33 (), pp [4] M. Zayernouri and G. E. Karniadakis, Fractional Sturm-Liouville eigen-problems: Theory and numerical approximations, J. Comput. Phys., 47 (3) 8-3. [5] M. Zayernouri and G. E. Karniadakis, Fractional spectral collocation method, SIAM J. Sci. Comput., 36 (4), pp. A4-A6. [6] F. Zeng, F. Liu, C. P. Li, K. Burrage, I. Turner, and V. Anh, Crank- Nicolson ADI spectral method for the -D Riesz space fractional nonlinear reaction-diffusion equation, SIAM J. Numer. Anal., 5 (4), pp [7] Z. Zhang, Superconvergence of a Chebyshev spectral collocation method, J. Sci. Comput. 34 (8) [8] Z. Zhang, Superconvergence points of polynomial spectral interpolation, SIAM J. Numer. Anal. 5 (), [9] X. Zhao, Z. Zhang, Superconvergence points of fractional spectral interpolation, SIAM J. Sci. Comput., 38 (6), pp. A598-A63. [3] M. Zheng, F. Liu, I. Turner and V. Anh, A novel high order space-time spectral method for the time-fractional Fokker-Planck equation, SIAM J. Sci. Comput, 37 (5), pp. A7- A74. [3] J. Shen, T. Tang, L-L. Wang, Spectral Methods: Algorithms, Analysis and Applications, Springer Series in Computational Mathematics, Vol. 4, Springer,. [3] Q. Lin and J. Lin, Finite Element Methods: Accuracy and Improvement, Math. Monogr. Ser., Science Press, Beijing, 6. [33] Z. Mao, S. Chen and J. Shen, Efficient and accurate spectral method using generalized Jacobi functions for solving Riesz fractional differential equations, Applied Numerical Mathematics, 6(6), pp [34] Z. Xie, L. Wang and X. Zhao, On exponential convergence of Gegenbauer interpolation and spectral differentiation, Math. Comp., 8() pp [35] L. B. Wahlbin, Superconvergence in Galerkin Finite Element Methods, Lecture Notes in Math. 65, Springer-Verlag, Berlin, 995. [36] J. P. Roop, Computational aspects of FEM approximation of fractional advection dispersion equations on bounded domains in R, J. Comput. Appl. Math, 93()(6), pp [37] K. Deng, W. Deng, Finite difference/predictor-corrector approximations for the space and time fractional Fokker-Planck equation, Appl. Math. Lett., 5()() pp [38] S. Shen, F. Liu, V. Anh, I. Terner, J. Chen, A novel numerical approximation for the space fractional advection-dispersion equation, IMA J. Appl. Math. 79(3)(4), pp [39] W. Bu, Y. Tang, J. Yang, Galerkin finite element method for two-dimensional Riesz space fractional diffusion equations, J. Comput. Phys. 76(4), pp [4] X. Zhao, Z. Sun, Z. Hao, A fourth-order compact ADI scheme for two-dimensional nonlinear space fractional Schrodinger equation, SIAM J. Sci. Comput. 36(4), pp [4] S. Lei, H. Sun, A circulant preconditioner for fractional diffusion equations, J. Comput. Phys. 4(3), pp [4] H. Pang, H. Sun, Multigrid method for fractional diffusion equations, J. Comput. Phys. 3(), pp
SUPERCONVERGENCE POINTS OF FRACTIONAL SPECTRAL INTERPOLATION
SIAM J. SCI. COMPUT. Vol. 38, No. 1, pp. A598 A613 c 16 Society for Industrial and Applied Mathematics SUPERCONVERGENCE POINTS OF FRACTIONAL SPECTRAL INTERPOLATION XUAN ZHAO AND ZHIMIN ZHANG Abstract.
More informationRational Chebyshev pseudospectral method for long-short wave equations
Journal of Physics: Conference Series PAPER OPE ACCESS Rational Chebyshev pseudospectral method for long-short wave equations To cite this article: Zeting Liu and Shujuan Lv 07 J. Phys.: Conf. Ser. 84
More informationSuperconvergence of discontinuous Galerkin methods for 1-D linear hyperbolic equations with degenerate variable coefficients
Superconvergence of discontinuous Galerkin methods for -D linear hyperbolic equations with degenerate variable coefficients Waixiang Cao Chi-Wang Shu Zhimin Zhang Abstract In this paper, we study the superconvergence
More informationSuperconvergence of Jacobi Gauss-Type Spectral Interpolation
J Sci Comput (204 59:667 687 DOI 0.007/s095-03-9777-x Superconvergence of Jacobi Gauss-Type Spectral Interpolation Li-Lian Wang Xiaodan Zhao Zhimin Zhang Received: 23 April 203 / Accepted: 9 August 203
More informationA Legendre Computational Matrix Method for Solving High-Order Fractional Differential Equations
Mathematics A Legendre Computational Matrix Method for Solving High-Order Fractional Differential Equations Mohamed Meabed KHADER * and Ahmed Saied HENDY Department of Mathematics, Faculty of Science,
More informationEFFICIENT SPECTRAL COLLOCATION METHOD FOR SOLVING MULTI-TERM FRACTIONAL DIFFERENTIAL EQUATIONS BASED ON THE GENERALIZED LAGUERRE POLYNOMIALS
Journal of Fractional Calculus and Applications, Vol. 3. July 212, No.13, pp. 1-14. ISSN: 29-5858. http://www.fcaj.webs.com/ EFFICIENT SPECTRAL COLLOCATION METHOD FOR SOLVING MULTI-TERM FRACTIONAL DIFFERENTIAL
More informationJACOBI SPECTRAL GALERKIN METHODS FOR VOLTERRA INTEGRAL EQUATIONS WITH WEAKLY SINGULAR KERNEL
Bull. Korean Math. Soc. 53 (016), No. 1, pp. 47 6 http://dx.doi.org/10.4134/bkms.016.53.1.47 JACOBI SPECTRAL GALERKIN METHODS FOR VOLTERRA INTEGRAL EQUATIONS WITH WEAKLY SINGULAR KERNEL Yin Yang Abstract.
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 information1462. Jacobi pseudo-spectral Galerkin method for second kind Volterra integro-differential equations with a weakly singular kernel
1462. Jacobi pseudo-spectral Galerkin method for second kind Volterra integro-differential equations with a weakly singular kernel Xiaoyong Zhang 1, Junlin Li 2 1 Shanghai Maritime University, Shanghai,
More informationHandling the fractional Boussinesq-like equation by fractional variational iteration method
6 ¹ 5 Jun., COMMUN. APPL. MATH. COMPUT. Vol.5 No. Å 6-633()-46-7 Handling the fractional Boussinesq-like equation by fractional variational iteration method GU Jia-lei, XIA Tie-cheng (College of Sciences,
More information256 Summary. D n f(x j ) = f j+n f j n 2n x. j n=1. α m n = 2( 1) n (m!) 2 (m n)!(m + n)!. PPW = 2π k x 2 N + 1. i=0?d i,j. N/2} N + 1-dim.
56 Summary High order FD Finite-order finite differences: Points per Wavelength: Number of passes: D n f(x j ) = f j+n f j n n x df xj = m α m dx n D n f j j n= α m n = ( ) n (m!) (m n)!(m + n)!. PPW =
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 informationA SPECTRAL COLLOCATION METHOD FOR EIGENVALUE PROBLEMS OF COMPACT INTEGRAL OPERATORS
JOURNAL OF INTEGRAL EQUATIONS AND APPLICATIONS Volume 5, Number 1, Spring 013 A SPECTRAL COLLOCATION METHOD FOR EIGENVALUE PROBLEMS OF COMPACT INTEGRAL OPERATORS CAN HUANG, HAILONG GUO AND ZHIMIN ZHANG
More informationON SPECTRAL METHODS FOR VOLTERRA INTEGRAL EQUATIONS AND THE CONVERGENCE ANALYSIS * 1. Introduction
Journal of Computational Mathematics, Vol.6, No.6, 008, 85 837. ON SPECTRAL METHODS FOR VOLTERRA INTEGRAL EQUATIONS AND THE CONVERGENCE ANALYSIS * Tao Tang Department of Mathematics, Hong Kong Baptist
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 informationAn Implicit Difference-ADI Method for the Two-dimensional Space-time Fractional Diffusion Equation
Iranian Journal of Mathematical Sciences and Informatics Vol., No. 2 (206), pp 7-86 DOI: 0.7508/ijmsi.206.02.005 An Implicit Difference-ADI Method for the Two-dimensional Space-time Fractional Diffusion
More informationJacobi spectral collocation method for the approximate solution of multidimensional nonlinear Volterra integral equation
Wei et al. SpringerPlus (06) 5:70 DOI 0.86/s40064-06-3358-z RESEARCH Jacobi spectral collocation method for the approximate solution of multidimensional nonlinear Volterra integral equation Open Access
More informationFractional Spectral and Spectral Element Methods
Fractional Calculus, Probability and Non-local Operators: Applications and Recent Developments Nov. 6th - 8th 2013, BCAM, Bilbao, Spain Fractional Spectral and Spectral Element Methods (Based on PhD thesis
More informationPART IV Spectral Methods
PART IV Spectral Methods Additional References: R. Peyret, Spectral methods for incompressible viscous flow, Springer (2002), B. Mercier, An introduction to the numerical analysis of spectral methods,
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 informationFractional Partial Differential Equation:
Fractional Partial Differential Equation: Mathematical and Numerical Analysis Issues Hong Wang Department of Mathematics, University of South Carolina IBM Visiting Fellow, Division of Applied Mathematics,
More informationOn the remainder term of Gauss Radau quadratures for analytic functions
Journal of Computational and Applied Mathematics 218 2008) 281 289 www.elsevier.com/locate/cam On the remainder term of Gauss Radau quadratures for analytic functions Gradimir V. Milovanović a,1, Miodrag
More informationFast solver for fractional differential equations based on Hierarchical Matrices
International Conference on Fractional Differentiation and its Applications, Novi Sad, Seria, July 18-20, 2016 Fast solver for fractional differential equations ased on Hierarchical Matrices Xuan Zhao
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 informationON THE FRACTIONAL CAUCHY PROBLEM ASSOCIATED WITH A FELLER SEMIGROUP
Dedicated to Professor Gheorghe Bucur on the occasion of his 7th birthday ON THE FRACTIONAL CAUCHY PROBLEM ASSOCIATED WITH A FELLER SEMIGROUP EMIL POPESCU Starting from the usual Cauchy problem, we give
More informationIMPROVEMENTS OF COMPOSITION RULE FOR THE CANAVATI FRACTIONAL DERIVATIVES AND APPLICATIONS TO OPIAL-TYPE INEQUALITIES
Dynamic Systems and Applications ( 383-394 IMPROVEMENTS OF COMPOSITION RULE FOR THE CANAVATI FRACTIONAL DERIVATIVES AND APPLICATIONS TO OPIAL-TYPE INEQUALITIES M ANDRIĆ, J PEČARIĆ, AND I PERIĆ Faculty
More informationRecurrence Relations and Fast Algorithms
Recurrence Relations and Fast Algorithms Mark Tygert Research Report YALEU/DCS/RR-343 December 29, 2005 Abstract We construct fast algorithms for decomposing into and reconstructing from linear combinations
More informationarxiv: v1 [math.ca] 28 Feb 2014
Communications in Nonlinear Science and Numerical Simulation. Vol.18. No.11. (213) 2945-2948. arxiv:142.7161v1 [math.ca] 28 Feb 214 No Violation of the Leibniz Rule. No Fractional Derivative. Vasily E.
More informationNUMERICAL SOLUTION FOR CLASS OF ONE DIMENSIONAL PARABOLIC PARTIAL INTEGRO DIFFERENTIAL EQUATIONS VIA LEGENDRE SPECTRAL COLLOCATION METHOD
Journal of Fractional Calculus and Applications, Vol. 5(3S) No., pp. 1-11. (6th. Symp. Frac. Calc. Appl. August, 01). ISSN: 090-5858. http://fcag-egypt.com/journals/jfca/ NUMERICAL SOLUTION FOR CLASS OF
More informationPositive solutions for a class of fractional boundary value problems
Nonlinear Analysis: Modelling and Control, Vol. 21, No. 1, 1 17 ISSN 1392-5113 http://dx.doi.org/1.15388/na.216.1.1 Positive solutions for a class of fractional boundary value problems Jiafa Xu a, Zhongli
More informationExistence of Solutions for Nonlocal Boundary Value Problems of Nonlinear Fractional Differential Equations
Advances in Dynamical Systems and Applications ISSN 973-5321, Volume 7, Number 1, pp. 31 4 (212) http://campus.mst.edu/adsa Existence of Solutions for Nonlocal Boundary Value Problems of Nonlinear Fractional
More informationarxiv: v1 [math.na] 4 Nov 2015
ON WELL-CONDITIONED SPECTRAL COLLOCATION AND SPECTRAL METHODS BY THE INTEGRAL REFORMULATION KUI DU arxiv:529v [mathna 4 Nov 25 Abstract Well-conditioned spectral collocation and spectral methods have recently
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 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 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 informationLECTURE 16 GAUSS QUADRATURE In general for Newton-Cotes (equispaced interpolation points/ data points/ integration points/ nodes).
CE 025 - Lecture 6 LECTURE 6 GAUSS QUADRATURE In general for ewton-cotes (equispaced interpolation points/ data points/ integration points/ nodes). x E x S fx dx hw' o f o + w' f + + w' f + E 84 f 0 f
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 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 informationApplication of high-order spectral method for the time fractional mobile/immobile equation
Computational Methods for Differential Equations http://cmde.tabrizu.ac.ir Vol. 4, No. 4, 2016, pp. 309-322 Application of high-order spectral method for the time fractional mobile/immobile equation Hossein
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 informationFINITE DIFFERENCE APPROXIMATIONS FOR A FRACTIONAL ADVECTION DIFFUSION PROBLEM ERCÍLIA SOUSA
Pré-Publicações do Departamento de Matemática Universidade de Coimbra Preprint Number 8 26 FINITE DIFFERENCE APPROXIMATIONS FOR A FRACTIONAL ADVECTION DIFFUSION PROBLEM ERCÍLIA SOUSA Abstract: The use
More informationc 2017 Society for Industrial and Applied Mathematics
SIAM J. SCI. COMPUT. Vol. 39, No. 5, pp. A198 A1950 c 017 Society for Industrial and Applied Mathematics HERMITE SPECTRAL METHODS FOR FRACTIONAL PDEs IN UNBOUNDED DOMAINS ZHIPING MAO AND JIE SHEN Abstract.
More informationA truncation regularization method for a time fractional diffusion equation with an in-homogeneous source
ITM Web of Conferences, 7 18) ICM 18 https://doi.org/1.151/itmconf/187 A truncation regularization method for a time fractional diffusion equation with an in-homogeneous source Luu Vu Cam Hoan 1,,, Ho
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 informationCONVERGENCE ANALYSIS OF THE JACOBI SPECTRAL-COLLOCATION METHODS FOR VOLTERRA INTEGRAL EQUATIONS WITH A WEAKLY SINGULAR KERNEL
COVERGECE AALYSIS OF THE JACOBI SPECTRAL-COLLOCATIO METHODS FOR VOLTERRA ITEGRAL EQUATIOS WITH A WEAKLY SIGULAR KEREL YAPIG CHE AD TAO TAG Abstract. In this paper, a Jacobi-collocation spectral method
More informationComputational Non-Polynomial Spline Function for Solving Fractional Bagely-Torvik Equation
Math. Sci. Lett. 6, No. 1, 83-87 (2017) 83 Mathematical Sciences Letters An International Journal http://dx.doi.org/10.18576/msl/060113 Computational Non-Polynomial Spline Function for Solving Fractional
More informationAccuracy Enhancement Using Spectral Postprocessing for Differential Equations and Integral Equations
COMMUICATIOS I COMPUTATIOAL PHYSICS Vol. 5, o. -4, pp. 779-79 Commun. Comput. Phys. February 009 Accuracy Enhancement Using Spectral Postprocessing for Differential Equations and Integral Equations Tao
More informationSuperconvergence of the Direct Discontinuous Galerkin Method for Convection-Diffusion Equations
Superconvergence of the Direct Discontinuous Galerkin Method for Convection-Diffusion Equations Waixiang Cao, Hailiang Liu, Zhimin Zhang,3 Division of Applied and Computational Mathematics, Beijing Computational
More informationA Petrov-Galerkin Finite Element Method for Solving the Time-fractional Diffusion Equation with Interface
Journal of Mathematics Research; Vol., No. 4; August 8 ISSN 96-979 E-ISSN 96-989 Published by Canadian Center of Science and Education A Petrov-Galerkin Finite Element Method for Solving the Time-fractional
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 informationOn The Leibniz Rule And Fractional Derivative For Differentiable And Non-Differentiable Functions
On The Leibniz Rule And Fractional Derivative For Differentiable And Non-Differentiable Functions Xiong Wang Center of Chaos and Complex Network, Department of Electronic Engineering, City University of
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
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 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 informationNumerical Analysis Comprehensive Exam Questions
Numerical Analysis Comprehensive Exam Questions 1. Let f(x) = (x α) m g(x) where m is an integer and g(x) C (R), g(α). Write down the Newton s method for finding the root α of f(x), and study the order
More informationA Fast O(N log N) Finite Difference Method for the One-Dimensional Space-Fractional Diffusion Equation
Mathematics 2015, 3, 1032-1044; doi:103390/math3041032 OPEN ACCESS mathematics ISSN 2227-7390 wwwmdpicom/journal/mathematics Article A Fast O(N log N) Finite Difference Method for the One-Dimensional Space-Fractional
More informationA NEW COLLOCATION METHOD FOR SOLVING CERTAIN HADAMARD FINITE-PART INTEGRAL EQUATION
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 6, Number, Pages 4 54 c 9 Institute for Scientific Computing and Information A NEW COLLOCATION METHOD FOR SOLVING CERTAIN HADAMARD FINITE-PART
More informationGeneralized Differential Transform Method to Space- Time Fractional Non-linear Schrodinger Equation
International Journal of Latest Engineering Research and Applications (IJLERA) ISSN: 455-737 Volume, Issue, December 7, PP 7-3 Generalized Differential Transform Method to Space- Time Fractional Non-linear
More informationApplied Numerical Mathematics
Applied Numerical Mathematics 131 (2018) 1 15 Contents lists available at ScienceDirect Applied Numerical Mathematics www.elsevier.com/locate/apnum The spectral-galerkin approximation of nonlinear eigenvalue
More informationFRACTIONAL FOURIER TRANSFORM AND FRACTIONAL DIFFUSION-WAVE EQUATIONS
FRACTIONAL FOURIER TRANSFORM AND FRACTIONAL DIFFUSION-WAVE EQUATIONS L. Boyadjiev*, B. Al-Saqabi** Department of Mathematics, Faculty of Science, Kuwait University *E-mail: boyadjievl@yahoo.com **E-mail:
More informationIsogeometric Analysis for the Fractional Laplacian
Isogeometric Analysis for the Fractional Laplacian Kailai Xu & Eric Darve CME500 May 14, 2018 1 / 26 Outline Fractional Calculus Isogeometric Analysis Numerical Experiment Regularity for the Fractional
More informationSPECTRAL METHODS: ORTHOGONAL POLYNOMIALS
SPECTRAL METHODS: ORTHOGONAL POLYNOMIALS 31 October, 2007 1 INTRODUCTION 2 ORTHOGONAL POLYNOMIALS Properties of Orthogonal Polynomials 3 GAUSS INTEGRATION Gauss- Radau Integration Gauss -Lobatto Integration
More informationA computationally effective predictor-corrector method for simulating fractional order dynamical control system
ANZIAM J. 47 (EMA25) pp.168 184, 26 168 A computationally effective predictor-corrector method for simulating fractional order dynamical control system. Yang F. Liu (Received 14 October 25; revised 24
More informationPreconditioners for ill conditioned (block) Toeplitz systems: facts a
Preconditioners for ill conditioned (block) Toeplitz systems: facts and ideas Department of Informatics, Athens University of Economics and Business, Athens, Greece. Email:pvassal@aueb.gr, pvassal@uoi.gr
More informationEXISTENCE OF SOLUTIONS FOR KIRCHHOFF TYPE EQUATIONS WITH UNBOUNDED POTENTIAL. 1. Introduction In this article, we consider the Kirchhoff type problem
Electronic Journal of Differential Equations, Vol. 207 (207), No. 84, pp. 2. ISSN: 072-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE OF SOLUTIONS FOR KIRCHHOFF TYPE EQUATIONS
More informationTHE differential equations of fractional order arise in
IAENG International Journal of Applied Mathematics 45:4 IJAM_45_4_4 Applications of Fixed Point Theorems for Coupled Systems of Fractional Integro-Differential Equations Involving Convergent Series Mohamed
More informationNontrivial Solutions for Boundary Value Problems of Nonlinear Differential Equation
Advances in Dynamical Systems and Applications ISSN 973-532, Volume 6, Number 2, pp. 24 254 (2 http://campus.mst.edu/adsa Nontrivial Solutions for Boundary Value Problems of Nonlinear Differential Equation
More informationNUMERICAL SOLUTION OF FRACTIONAL DIFFUSION-WAVE EQUATION WITH TWO SPACE VARIABLES BY MATRIX METHOD. Mridula Garg, Pratibha Manohar.
NUMERICAL SOLUTION OF FRACTIONAL DIFFUSION-WAVE EQUATION WITH TWO SPACE VARIABLES BY MATRIX METHOD Mridula Garg, Pratibha Manohar Abstract In the present paper we solve space-time fractional diffusion-wave
More informationCONSTRUCTION OF THE HALF-LINE POTENTIAL FROM THE JOST FUNCTION
CONSTRUCTION OF THE HALF-LINE POTENTIAL FROM THE JOST FUNCTION Tuncay Aktosun Department of Mathematics and Statistics Mississippi State University Mississippi State, MS 39762 Abstract: For the one-dimensional
More informationNUMERICAL SOLUTION OF DELAY DIFFERENTIAL EQUATIONS VIA HAAR WAVELETS
TWMS J Pure Appl Math V5, N2, 24, pp22-228 NUMERICAL SOLUTION OF DELAY DIFFERENTIAL EQUATIONS VIA HAAR WAVELETS S ASADI, AH BORZABADI Abstract In this paper, Haar wavelet benefits are applied to the delay
More informationOn Non-degeneracy of Solutions to SU(3) Toda System
On Non-degeneracy of Solutions to SU3 Toda System Juncheng Wei Chunyi Zhao Feng Zhou March 31 010 Abstract We prove that the solution to the SU3 Toda system u + e u e v = 0 in R v e u + e v = 0 in R e
More informationNumerical study of time-fractional hyperbolic partial differential equations
Available online at wwwisr-publicationscom/jmcs J Math Computer Sci, 7 7, 53 65 Research Article Journal Homepage: wwwtjmcscom - wwwisr-publicationscom/jmcs Numerical study of time-fractional hyperbolic
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 informationEXISTENCE AND UNIQUENESS OF SOLUTIONS FOR A SYSTEM OF FRACTIONAL DIFFERENTIAL EQUATIONS 1. Yong Zhou. Abstract
EXISTENCE AND UNIQUENESS OF SOLUTIONS FOR A SYSTEM OF FRACTIONAL DIFFERENTIAL EQUATIONS 1 Yong Zhou Abstract In this paper, the initial value problem is discussed for a system of fractional differential
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 informationOSGOOD TYPE REGULARITY CRITERION FOR THE 3D NEWTON-BOUSSINESQ EQUATION
Electronic Journal of Differential Equations, Vol. 013 (013), No. 3, pp. 1 6. ISSN: 107-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu OSGOOD TYPE REGULARITY
More informationGauss Hermite interval quadrature rule
Computers and Mathematics with Applications 54 (2007) 544 555 www.elsevier.com/locate/camwa Gauss Hermite interval quadrature rule Gradimir V. Milovanović, Alesandar S. Cvetović Department of Mathematics,
More informationSome lecture notes for Math 6050E: PDEs, Fall 2016
Some lecture notes for Math 65E: PDEs, Fall 216 Tianling Jin December 1, 216 1 Variational methods We discuss an example of the use of variational methods in obtaining existence of solutions. Theorem 1.1.
More informationOn the simplest expression of the perturbed Moore Penrose metric generalized inverse
Annals of the University of Bucharest (mathematical series) 4 (LXII) (2013), 433 446 On the simplest expression of the perturbed Moore Penrose metric generalized inverse Jianbing Cao and Yifeng Xue Communicated
More informationNonlocal problems for the generalized Bagley-Torvik fractional differential equation
Nonlocal problems for the generalized Bagley-Torvik fractional differential equation Svatoslav Staněk Workshop on differential equations Malá Morávka, 28. 5. 212 () s 1 / 32 Overview 1) Introduction 2)
More informationSufficient conditions for functions to form Riesz bases in L 2 and applications to nonlinear boundary-value problems
Electronic Journal of Differential Equations, Vol. 200(200), No. 74, pp. 0. ISSN: 072-669. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu (login: ftp) Sufficient conditions
More informationZERO DISTRIBUTION OF POLYNOMIALS ORTHOGONAL ON THE RADIAL RAYS IN THE COMPLEX PLANE* G. V. Milovanović, P. M. Rajković and Z. M.
FACTA UNIVERSITATIS (NIŠ) Ser. Math. Inform. 12 (1997), 127 142 ZERO DISTRIBUTION OF POLYNOMIALS ORTHOGONAL ON THE RADIAL RAYS IN THE COMPLEX PLANE* G. V. Milovanović, P. M. Rajković and Z. M. Marjanović
More informationA SPACE-TIME SPECTRAL PETROV-GALERKIN SPECTRAL METHOD FOR TIME FRACTIONAL DIFFUSION EQUATION
A SPACE-IE SPECRAL PEROV-GALERKI SPECRAL EHOD FOR IE FRACIOAL DIFFUSIO EQUAIO Chang-ao Sheng and Jie Shen,, Abstract. We develop in this paper a space-time Petrov-Galerkin spectral method for linear and
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 informationMATH 543: FUCHSIAN DIFFERENTIAL EQUATIONS HYPERGEOMETRIC FUNCTION
SET 7 MATH 543: FUCHSIAN DIFFERENTIAL EQUATIONS HYERGEOMETRIC FUNCTION References: DK and Sadri Hassan. Historical Notes: lease read the book Linear Differential Equations and the Group Theory by Jeremy
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 informationarxiv: v2 [math.ap] 30 Jul 2012
Blow up for some semilinear wave equations in multi-space dimensions Yi Zhou Wei Han. arxiv:17.536v [math.ap] 3 Jul 1 Abstract In this paper, we discuss a new nonlinear phenomenon. We find that in n space
More informationNonlocal Fractional Semilinear Delay Differential Equations in Separable Banach spaces
IOSR Journal of Mathematics (IOSR-JM) e-issn: 2278-5728, p-issn:2319-765x. Volume 10, Issue 1 Ver. IV. (Feb. 2014), PP 49-55 Nonlocal Fractional Semilinear Delay Differential Equations in Separable Banach
More informationCrank-Nicolson Scheme for Space Fractional Heat Conduction Equation with Mixed Boundary Condition
Crank-Nicolson Scheme for Space Fractional Heat Conduction Equation with Mied Boundary Condition Rafał Brociek Institute of Mathematics Silesian University of Technology Kaszubska 23, 44-00 Gliwice, Poland
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 informationDownloaded 01/21/14 to Redistribution subject to SIAM license or copyright; see
SIAM J. SCI. COMPUT. Vol. 6, o. 1, pp. A40 A62 c 2014 Society for Industrial and Applied Mathematics FRACTIOAL SPECTRAL COLLOCATIO METHOD MOHSE ZAYEROURI AD GEORGE EM KARIADAKIS Abstract. We develop an
More informationc 2017 Society for Industrial and Applied Mathematics
SIAM J. SCI. COMPUT. Vol. 39, No., pp. A24 A228 c 27 Society for Industrial and Applied Mathematics RATIONAL APPROXIMATION TO THE FRACTIONAL LAPLACIAN OPERATOR IN REACTION-DIFFUSION PROBLEMS LIDIA ACETO
More information1 Assignment 1: Nonlinear dynamics (due September
Assignment : Nonlinear dynamics (due September 4, 28). Consider the ordinary differential equation du/dt = cos(u). Sketch the equilibria and indicate by arrows the increase or decrease of the solutions.
More informationOn the positivity of linear weights in WENO approximations. Abstract
On the positivity of linear weights in WENO approximations Yuanyuan Liu, Chi-Wang Shu and Mengping Zhang 3 Abstract High order accurate weighted essentially non-oscillatory (WENO) schemes have been used
More informationThe local fractional Hilbert transform in fractal space
The local fractional ilbert transform in fractal space Guang-Sheng Chen Department of Computer Engineering, Guangxi Modern Vocational Technology College, echi,guangxi, 547000, P.. China E-mail address:
More informationNumerical differentiation by means of Legendre polynomials in the presence of square summable noise
www.oeaw.ac.at Numerical differentiation by means of Legendre polynomials in the presence of square summable noise S. Lu, V. Naumova, S. Pereverzyev RICAM-Report 2012-15 www.ricam.oeaw.ac.at Numerical
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 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 informationOrthonormal series expansion and finite spherical Hankel transform of generalized functions
Malaya Journal of Matematik 2(1)(2013) 77 82 Orthonormal series expansion and finite spherical Hankel transform of generalized functions S.K. Panchal a, a Department of Mathematics, Dr. Babasaheb Ambedkar
More informationA derivative-free nonmonotone line search and its application to the spectral residual method
IMA Journal of Numerical Analysis (2009) 29, 814 825 doi:10.1093/imanum/drn019 Advance Access publication on November 14, 2008 A derivative-free nonmonotone line search and its application to the spectral
More information