On some explicit Adams multistep methods for fractional differential equations

In this paper we present a family of explicit formulas for the numerical solution of differential equations of fractional order. The proposed methods are obtained by modifying, in a suitable way, Fractional-Adams–Moulton methods and they represent a way for extending classical Adams–Bashforth multistep methods to the fractional case. The attention is hence focused on the investigation of stability properties. Intervals of stability for k$k$-step methods, k=1,…,5$k=1,\dots ,5$, are computed and plots of stability regions in the complex plane are presented.     

Numerical approach to differential equations of fractional order

In this paper, the variational iteration method and the Adomian decomposition method are implemented to give approximate solutions for linear and nonlinear systems of differential equations of fractional order. The two methods in applied mathematics can be used as alternative methods for obtaining analytic and approximate solutions for different types of differential equations. In these schemes, the solution takes the form of a convergent series with easily computable components. This paper presents a numerical comparison between the two methods for solving systems of fractional differential equations. Numerical results show that the two approaches are easy to implement and accurate when applied to differential equations of fractional order.     

Weighted finite difference methods for a class of space fractional partial differential equations with variable coefficients

A class of weighted finite difference methods (WFDMs) for solving a class of initial-boundary value problems of space fractional partial differential equations with variable coefficients is presented. Their stability and convergence properties are considered. It is proven that the WFDMs are unconditionally-stable for , and conditionally-stable for , here r is the weighting parameter subjected to 0≤r≤1. Some convergence results are given. These methods, problems and results generalize the corresponding methods, problems and results given in [7], [8] and [10]. Some numerical examples are provided to show the effectiveness of the methods with different weighting parameters.     

On accurate product integration rules for linear fractional differential equations

This paper addresses the numerical solution of linear fractional differential equations with a forcing term. Competitive and highly accurate Product Integration rules are derived by starting from an equivalent formulation in terms of a Volterra integral equation with a generalized Mittag-Leffler function in the kernel. The error analysis is reported and aspects related to the computational complexity are treated. Numerical tests confirming the theoretical findings are presented.     

Solving systems of fractional differential equations using differential transform method

This paper presents approximate analytical solutions for systems of fractional differential equations using the differential transform method. The fractional derivatives are described in the Caputo sense. The application of differential transform method, developed for differential equations of integer order, is extended to derive approximate analytical solutions of systems of fractional differential equations. The solutions of our model equations are calculated in the form of convergent series with easily computable components. Some examples are solved as illustrations, using symbolic computation. The numerical results show that the approach is easy to implement and accurate when applied to systems of fractional differential equations. The method introduces a promising tool for solving many linear and nonlinear fractional differential equations.     

A Fourier method and an extrapolation technique for Stokes’ first problem for a heated generalized second grade fluid with fractional derivative

In this article, we consider Stokes’ first problem for a heated generalized second grade fluid with fractional derivative (SFP-HGSGF). Implicit and explicit numerical approximation schemes for the SFP-HGSGF are presented. The stability and convergence of the numerical schemes are discussed using a Fourier method. In addition, the solvability of the implicit numerical approximation scheme is also analyzed. A Richardson extrapolation technique for improving the order of convergence of the implicit scheme is proposed. Finally, a numerical test is given. The numerical results demonstrate the good performance of our theoretical analysis.     

Application of the fractional differential transform method to fractional-order integro-differential equations with nonlocal boundary conditions

In this paper, the fractional differential transform method is developed to solve fractional integro-differential equations with nonlocal boundary conditions. The method is described and illustrated with numerical examples. The results reveal that the method is accurate and easy to apply.     

The asymptotic stability of one-parameter methods for neutral differential equations

This paper deals with the asymptotic stability of theoretical solutions and numerical methods for systems of neutral differential equationsx=Ax(t–)+Bx(t)+Cx(t–), whereA, B, andC are constant complexN ×N matrices, and >0. A necessary and sufficient condition such that the differential equations are asymptotically stable is derived. We also focus on the numerical stability properties of adaptations of one-parameter methods. Further, we investigate carefully the characterization of the stability region.     

Short memory principle and a predictor–corrector approach for fractional differential equations

Fractional differential equations are increasingly used to model problems in acoustics and thermal systems, rheology and modelling of materials and mechanical systems, signal processing and systems identification, control and robotics, and other areas of application. This paper further analyses the underlying structure of fractional differential equations. From a new point of view, we apprehend the short memory principle of fractional calculus and farther apply a Adams-type predictor–corrector approach for the numerical solution of fractional differential equation. And the detailed error analysis is presented. Combining the short memory principle and the predictor–corrector approach, we gain a good numerical approximation of the true solution of fractional differential equation at reasonable computational cost. A numerical example is provided and compared with the exact analytical solution for illustrating the effectiveness of the short memory principle.     

Analytical solution of the linear fractional differential equation by Adomian decomposition method

In this paper, we consider the n-term linear fractional-order differential equation with constant coefficients and obtain the solution of this kind of fractional differential equations by Adomian decomposition method. With the equivalent transmutation, we show that the solution by Adomian decomposition method is the same as the solution by the Green's function. Finally, we illustrate our result with some examples.     

Runge-Kutta methods for quadratic ordinary differential equations

Many systems of ordinary differential equations are quadratic: the derivative can be expressed as a quadratic function of the dependent variable. We demonstrate that this feature can be exploited in the numerical solution by Runge-Kutta methods, since the quadratic structure serves to decrease the number of order conditions. We discuss issues related to construction design and implementation and present a number of new methods of Runge-Kutta and Runge-Kutta-Nyström type that display superior behaviour when applied to quadratic ordinary differential equations.     

High-order adaptive finite element-singly implicit Runge-Kutta methods for parabolic differential equations

We describe an adaptive mesh refinement finite element method-of-lines procedure for solving one-dimensional parabolic partial differential equations. Solutions are calculated using Galerkin's method with a piecewise hierarchical polynomial basis in space and singly implicit Runge-Kutta (SIRK) methods in time. A modified SIRK formulation eliminates a linear systems solution that is required by the traditional SIRK formulation and leads to a new reduced-order interpolation formula. Stability and temporal error estimation techniques allow acceptance of approximate solutions at intermediate stages, yielding increased efficiency when solving partial differential equations. A priori energy estimates of the local discretization error are obtained for a nonlinear scalar problem. A posteriori estimates of local spatial discretization errors, obtained by order variation, are used with the a priori error estimates to control the adaptive mesh refinement strategy. Computational results suggest convergence of the a posteriori error estimate to the exact discretization error and verify the utility of the adaptive technique.This research was partially supported by the U.S. Air Force Office of Scientific Research, Air Force Systems Command, USAF, under Grant Number AFOSR-90-0194; the U.S. Army Research Office under Contract Number DAAL 03-91-G-0215; by the National Science Foundation under Grant Number CDA-8805910; and by a grant from the Committee on Research, Tulane University.     

The composite Milstein methods for the numerical solution of Ito stochastic differential equations

In this paper, we present the composite Milstein methods for the strong solution of Ito stochastic differential equations. These methods are a combination of semi-implicit and implicit Milstein methods. We give a criterion for choosing either the implicit or the semi-implicit scheme at each step of our numerical solution. The stability and convergence properties are investigated and discussed for the linear test equation. The convergence properties for the nonlinear case are shown numerically to be the same as the linear case. The stability properties of the composite Milstein methods are found to be more superior compared to those of the Milstein, the Euler and even better than the composite Euler method. This superiority in stability makes the methods a better candidate for the solution of stiff SDEs.     

Numerical solution of fractional order differential equations by extrapolation

We present an extrapolation type algorithm for the numerical solution of fractional order differential equations. It is based on the new result that the sequence of approximate solutions of these equations, computed by means of a recently published algorithm by Diethelm [6], possesses an asymptotic expansion with respect to the stepsize. From this we conclude that the application of extrapolation is justified, and we obtain a very efficient differential equation solver with practically no additional numerical costs. This is also illustrated by a number of numerical examples. This revised version was published online in August 2006 with corrections to the Cover Date.     

The monotonic property and stability of solutions of fractional differential equations

In this paper we improve on the monotone property of Lemma 1.7.3 in Lakshmikantham et al. (2009) [5] for the case g(t,u)=λu with a nonnegative real number λ. We also investigate the Mittag-Leffler stability of solutions of fractional differential equations by using the fractional comparison principle.     

A stability property ofA-stable natural Runge-Kutta methods for systems of delay differential equations

A natural Runge-Kutta method is a special type of Runge-Kutta method for delay differential equations (DDEs); it is known that any one-step collocation method is equivalent to one of such methods. In this paper, we consider a linear constant-coefficient system of DDEs with a constant delay, and discuss the application of natural Runge-Kutta methods to the system. We show that anA-stable method preserves the asymptotic stability property of the analytical solutions of the system.     

A method for solving differential equations of fractional order

In this paper, we consider Caputo type fractional differential equations of order 0<α<1$0<\alpha <1$ with initial condition x(0)=_x0$x\left(0\right)=x_0$. We introduce a technique to find the exact solutions of fractional differential equations by using the solutions of integer order differential equations. Generalization of the technique to finite systems is also given. Finally, we give some examples to illustrate the applications of our results.     

On the asymptotic stability properties of Runge-Kutta methods for delay differential equations

Summary. In this paper asymptotic stability properties of Runge-Kutta (R-K) methods for delay differential equations (DDEs) are considered with respect to the following test equation: where and is a continuous real-valued function. In the last few years, stability properties of R-K methods applied to DDEs have been studied by numerous authors who have considered regions of asymptotic stability for “any positive delay” (and thus independent of the specific value of ). In this work we direct attention at the dependence of stability regions on a fixed delay . In particular, natural Runge-Kutta methods for DDEs are extensively examined. Received April 15, 1996 / Revised version received August 8, 1996     

Generalized Cauchy type problems for nonlinear fractional differential equations with composite fractional derivative operator

This paper is devoted to proving the existence and uniqueness of solutions to Cauchy type problems for fractional differential equations with composite fractional derivative operator on a finite interval of the real axis in spaces of summable functions. An approach based on the equivalence of the nonlinear Cauchy type problem to a nonlinear Volterra integral equation of the second kind and applying a variant of the Banach’s fixed point theorem to prove uniqueness and existence of the solution is presented. The Cauchy type problems for integro-differential equations of Volterra type with composite fractional derivative operator, which contain the generalized Mittag-Leffler function in the kernel, are considered. Using the method of successive approximation, and the Laplace transform method, explicit solutions of the open problem proposed by Srivastava and Tomovski (2009) [11] are established in terms of the multinomial Mittag-Leffler function.