We consider the pricing of European options under a modified Black–Scholes equation having fractional derivatives in the “spatial” (price) variable. To be specific, the underlying price is assumed to follow a geometric Koponen–Boyarchenko–Levendorski process. This pure jump Lévy process could better capture the real behaviour of market data. Despite many difficulties caused by the “globalness” of the fractional derivatives, we derive an explicit closed-form analytical solution by solving the fractional partial differential equation analytically, using the Fourier transform technique. Based on the newly derived formula, we also examine, in theory, many basic properties of the option price under the current model. On the other hand, for practical purposes, we impose a reliable implementation method for the current formula so that it can be easily used in the trading market. With the numerical results, the impact of different parameters on the option price are also investigated.

This article is intended to fill in the blank of the numerical schemes with second-order convergence accuracy in time for nonlinear Stokes’ first problem for a heated generalized second grade fluid with fractional derivative. A linearized difference scheme is proposed. The time fractional-order derivative is discretized by second-order shifted and weighted Gr¨unwald-Letnikov difference operator. The convergence accuracy in space is improved by performing the average operator. The presented numerical method is unconditionally stable with the global convergence order of in maximum norm, where τ and h are the step sizes in time and space, respectively. Finally, numerical examples are carried out to verify the theoretical results, showing that our scheme is efficient indeed.

The computational work and storage of numerically solving the time fractional PDEs are generally huge for the traditional direct methods since they require total memory and work, where NT and NS represent the total number of time steps and grid points in space, respectively. To overcome this difficulty, we present an efficient algorithm for the evaluation of the Caputo fractional derivative of order α∈(0,1). The algorithm is based on an efficient sum-of-exponentials (SOE) approximation for the kernel t–1–α on the interval [Δt, T] with a uniform absolute error ε. We give the theoretical analysis to show that the number of exponentials Nexp needed is of order for T≫1 or for TH1 for fixed accuracy ε. The resulting algorithm requires only storage and work when numerically solving the time fractional PDEs. Furthermore, we also give the stability and error analysis of the new scheme, and present several numerical examples to demonstrate the performance of our scheme.

We obtain the long-time behaviour to the variance of the distribution process associated with the solution of the telegraph equation. To this end, we use a version of the Karamata-Feller Tauberian theorem.

In this article, we have introduced a Taylor collocation method, which is based on collocation method for solving fractional Riccati differential equation. The fractional derivatives are described in the Caputo sense. This method is based on first taking the truncated Taylor expansions of the solution function in the fractional Riccati differential equation and then substituting their matrix forms into the equation. Using collocation points, the systems of nonlinear algebraic equation is derived. We further solve the system of nonlinear algebraic equation using Maple 13 and thus obtain the coefficients of the generalized Taylor expansion. Illustrative examples are presented to demonstrate the effectiveness of the proposed method.

The inverse stable subordinator provides a probability model for time-fractional differential equations, and leads to explicit solution formulae. This paper reviews properties of the inverse stable subordinator, and applications to a variety of problems in mathematics and physics. Several different governing equations for the inverse stable subordinator have been proposed in the literature. This paper also shows how these equations can be reconciled.

As the generalization of the integer order partial differential equations (PDE), the fractional order PDEs are drawing more and more attention for their applications in fluid flow, finance and other areas. This paper presents high-order accurate Runge-Kutta local discontinuous Galerkin (DG) methods for one- and two-dimensional fractional diffusion equations containing derivatives of fractional order in space. The Caputo derivative is chosen as the representation of spatial derivative, because it may represent the fractional derivative by an integral operator. Some numerical examples show that the convergence orders of the proposed local Pk-DG methods are O(hk+1) both in one and two dimensions, where Pk denotes the space of the real-valued polynomials with degree at most k.

A scheme is developed to study numerical solution of the time-fractional shock wave equation and wave equation under initial conditions by the homotopy perturbation method (HPM). The fractional derivatives are taken in the Caputo sense. The solutions are given in the form of series with easily computable terms. Numerical results are illustrated through the graph.

We prove estimates for the partial derivatives of the solution to a time-fractional diffusion equation posed over a bounded spatial domain. Such estimates are needed for the analysis of effective numerical methods, particularly since the solution is typically less regular than in the familiar case of classical diffusion.

Continuous-time random walks incorporate a random waiting time between random jumps. They are used in physics to model particle motion. A physically realistic rescaling uses two different time scales for the mean waiting time and the deviation from the mean. This paper derives the scaling limits for such processes. These limit processes are governed by fractional partial differential equations that may be useful in physics. A transfer theorem for weak convergence of finite-dimensional distributions of stochastic processes is also obtained.