A new Jacobi operational matrix: An application for solving fractional differential equations

In this paper, we derived the shifted Jacobi operational matrix (JOM) of fractional derivatives which is applied together with spectral tau method for numerical solution of general linear multi-term fractional differential equations (FDEs). A new approach implementing shifted Jacobi operational matrix in combination with the shifted Jacobi collocation technique is introduced for the numerical solution of nonlinear multi-term FDEs. The main characteristic behind this approach is that it reduces such problems to those of solving a system of algebraic equations which greatly simplifying the problem. The proposed methods are applied for solving linear and nonlinear multi-term FDEs subject to initial or boundary conditions, and the exact solutions are obtained for some tested problems. Special attention is given to the comparison of the numerical results obtained by the new algorithm with those found by other known methods.     

An integral operational matrix based on Jacobi polynomials for solving fractional-order differential equations

This paper aims to construct a general formulation for the Jacobi operational matrix of fractional integral operator. Fractional calculus has been used to model physical and engineering processes that are found to be best described by fractional differential equations. Therefore, a reliable and efficient technique for the solution of them is too important. For the concept of fractional derivative we will adopt Caputo’s definition by using Riemann–Liouville fractional integral operator. Our main aim is to generalize the Jacobi integral operational matrix to the fractional calculus. These matrices together with the Tau method are then utilized to reduce the solution of this problem to the solution of a system of algebraic equations. The method is applied to solve linear and nonlinear fractional differential equations. Illustrative examples are included to demonstrate the validity and applicability of the presented technique.     

The convergence and stability analysis of the Jacobi collocation method for solving nonlinear fractional differential equations with integral boundary conditions

In this paper, we apply the Jacobi collocation method for solving nonlinear fractional differential equations with integral boundary conditions. Due to existence of integral boundary conditions, after reformulation of this equation in the integral form, the method is proposed for solving the obtained integral equation. Also, the convergence and stability analysis of the proposed method are studied in two main theorems. Furthermore, the optimum degree of convergence in the L₂ norm is obtained for this method. Furthermore, some numerical examples are presented in order to illustrate the performance of the presented method. Finally, an application of the model in control theory is introduced. Copyright © 2015 John Wiley & Sons, Ltd.     

The alternative Legendre Tau method for solving nonlinear multi-order fractional differential equations

In this paper, the alternative Legendre polynomials (ALPs) are used to approximate the solution of a class of nonlinear multi-order fractional differential equations (FDEs). First, the operational matrix of fractional integration of an arbitrary order and the product operational matrix are derived for ALPs. These matrices together with the spectral Tau method are then utilized to reduce the solution of the mentioned equations into the one of solving a system of nonlinear algebraic equations with unknown ALP coefficients of the exact solution. The fractional derivatives are considered in the Caputo sense and the fractional integration is described in the Riemann-Liouville sense. Numerical examples illustrate that the present method is very effective for linear and nonlinear multi-order FDEs and high accuracy solutions can be obtained only using a small number of ALPs.     

Jacobi spectral Galerkin method for the integrated forms of second-order differential equations

Efficiently direct solvers based on the Jacobi-Galerkin method for the integrated forms of second-order elliptic equations in one and two space variables are presented. They are based on appropriate base functions for the Galerkin formulation which lead to discrete systems with specially structured matrices that can be efficiently inverted. The homogeneous Dirichlet boundary conditions are satisfied exactly by expanding the unknown variable into a polynomial basis of functions which are built upon the Jacobi polynomials. The direct solution algorithm here developed for the homogeneous Dirichlet problem in two-dimensions relies upon a diagonalization process. Numerical results indicating the high accuracy and effectiveness of these algorithms are presented.     

A Jacobi spectral Galerkin method for the integrated forms of fourth‐order elliptic differential equations

This article analyzes the solution of the integrated forms of fourth‐order elliptic differential equations on a rectilinear domain using a spectral Galerkin method. The spatial approximation is based on Jacobi polynomials P (x), with α, β ∈ (?1, ∞) and n the polynomial degree. For α = β, one recovers the ultraspherical polynomials (symmetric Jacobi polynomials) and for α = β = ?½, α = β = 0, the Chebyshev of the first and second kinds and Legendre polynomials respectively; and for the nonsymmetric Jacobi polynomials, the two important special cases α = ?β = ±½ (Chebyshev polynomials of the third and fourth kinds) are also recovered. The two‐dimensional version of the approximations is obtained by tensor products of the one‐dimensional bases. The various matrix systems resulting from these discretizations are carefully investigated, especially their condition number. An algebraic preconditioning yields a condition number of O(N), N being the polynomial degree of approximation, which is an improvement with respect to the well‐known condition number O(N⁸) of spectral methods for biharmonic elliptic operators. The numerical complexity of the solver is proportional to N^d+1 for a d‐dimensional problem. This operational count is the best one can achieve with a spectral method. The numerical results illustrate the theory and constitute a convincing argument for the feasibility of the method. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

Efficient spectral-Galerkin algorithms for direct solution for second-order differential equations using Jacobi polynomials

It is well known that spectral methods (tau, Galerkin, collocation) have a condition number of ( is the number of retained modes of polynomial approximations). This paper presents some efficient spectral algorithms, which have a condition number of , based on the Jacobi–Galerkin methods of second-order elliptic equations in one and two space variables. The key to the efficiency of these algorithms is to construct appropriate base functions, which lead to systems with specially structured matrices that can be efficiently inverted. The complexities of the algorithms are a small multiple of operations for a -dimensional domain with unknowns, while the convergence rates of the algorithms are exponentials with smooth solutions.     

A computational matrix method for solving systems of high order fractional differential equations

In this paper, we introduced an accurate computational matrix method for solving systems of high order fractional differential equations. The proposed method is based on the derived relation between the Chebyshev coefficient matrix A of the truncated Chebyshev solution u(t)$\mathbf{u}(\mathbf{t})$ and the Chebyshev coefficient matrix A^(ν)${\mathbf{A}}^{(\nu )}$ of the fractional derivative u^(ν)${\mathbf{u}}^{(\nu )}$. The fractional derivatives are presented in terms of Caputo sense. The matrix method for the approximate solution for the systems of high order fractional differential equations (FDEs) in terms of Chebyshev collocation points is presented. The systems of FDEs and their conditions (initial or boundary) are transformed to matrix equations, which corresponds to system of algebraic equations with unknown Chebyshev coefficients. The remaining set of algebraic equations is solved numerically to yield the Chebyshev coefficients. Several numerical examples for real problems are provided to confirm the accuracy and effectiveness of the present method.     

A procedure for solving some second-order linear ordinary differential equations

The purpose of this work is to introduce the concept of pseudo-exactness for second-order linear ordinary differential equations (ODEs), and then to try to solve some specific ODEs.     

A new modification of the Adomian decomposition method for solving boundary value problems for higher order nonlinear differential equations

In this paper we propose a new modified recursion scheme for the resolution of multi-order and multi-point boundary value problems for nonlinear ordinary and partial differential equations by the Adomian decomposition method (ADM). Our new approach, including Duan’s convergence parameter, provides a significant computational advantage by allowing for the acceleration of convergence and expansion of the interval of convergence during calculations of the solution components for nonlinear boundary value problems, in particular for such cases when one of the boundary points lies outside the interval of convergence of the usual decomposition series. We utilize the boundary conditions to derive an integral equation before establishing the recursion scheme for the solution components. Thus we can derive a modified recursion scheme without any undetermined coefficients when computing successive solution components, whereas several prior recursion schemes have done so. This modification also avoids solving a sequence of nonlinear algebraic equations for the undetermined coefficients fraught with multiple roots, which is required to complete calculation of the solution by several prior modified recursion schemes using the ADM.     

Efficient Chebyshev spectral methods for solving multi-term fractional orders differential equations

In this paper, we state and prove a new formula expressing explicitly the derivatives of shifted Chebyshev polynomials of any degree and for any fractional-order in terms of shifted Chebyshev polynomials themselves. We develop also a direct solution technique for solving the linear multi-order fractional differential equations (FDEs) with constant coefficients using a spectral tau method. The spatial approximation with its fractional-order derivatives (described in the Caputo sense) are based on shifted Chebyshev polynomials T_L,n(x) with x ∈ (0, L), L > 0 and n is the polynomial degree. We presented a shifted Chebyshev collocation method with shifted Chebyshev–Gauss points used as collocation nodes for solving nonlinear multi-order fractional initial value problems. Several numerical examples are considered aiming to demonstrate the validity and applicability of the proposed techniques and to compare with the existing results.     

Generalized fractional‐order Jacobi functions for solving a nonlinear systems of fractional partial differential equations numerically

In this paper, a collocation spectral numerical algorithm is presented for solving nonlinear systems of fractional partial differential equations subject to different types of conditions. A proposed error analysis investigates the convergence of the mentioned algorithm. Some numerical examples confirm the efficiency and accuracy of the method.     

An iterative method for solving nonlinear functional equations

An iterative method for solving nonlinear functional equations, viz. nonlinear Volterra integral equations, algebraic equations and systems of ordinary differential equation, nonlinear algebraic equations and fractional differential equations has been discussed.     

The solution of high-order nonlinear ordinary differential equations by Chebyshev Series

By the use of the Chebyshev series, a direct computational method for solving the higher order nonlinear differential equations has been developed in this paper. This method transforms the nonlinear differential equation into the matrix equation, which corresponds to a system of nonlinear algebraic equations with unknown Chebyshev coefficients, via Chebyshev collocation points. The solution of this system yields the Chebyshev coefficients of the solution function. An algorithm for this nonlinear system is also proposed in this paper. The method is valid for both initial-value and boundary-value problems. Several examples are presented to illustrate the accuracy and effectiveness of the method.     

On elliptic solutions of nonlinear ordinary differential equations

The problem of constructing and classifying exact elliptic solutions of autonomous nonlinear ordinary differential equations is studied. An algorithm for finding elliptic solutions in explicit form is presented.     

Jacobi rational approximation and spectral method for differential equations of degenerate type

We introduce an orthogonal system on the half line, induced by Jacobi polynomials. Some results on the Jacobi rational approximation are established, which play important roles in designing and analyzing the Jacobi rational spectral method for various differential equations, with the coefficients degenerating at certain points and growing up at infinity. The Jacobi rational spectral method is proposed for a model problem appearing frequently in finance. Its convergence is proved. Numerical results demonstrate the efficiency of this new approach.

     

Solution of nonlinear fractional differential equations using homotopy analysis method

In this article, the homotopy analysis method has been applied to solve nonlinear differential equations of fractional order. The validity of this method has successfully been accomplished by applying it to find the solution of two nonlinear fractional equations. The results obtained by homotopy analysis method have been compared with those exact solutions. The results show that the solution of homotopy analysis method is good agreement with the exact solution.