A high order numerical method for solving nonlinear fractional differential equation with non-uniform meshes

Abstract

We introduce a high-order numerical method for solving nonlinear fractional differential equation with non-uniform meshes. We first transform the fractional nonlinear differential equation into the equivalent Volterra integral equation. Then we approximate the integral by using the quadratic interpolation polynomials. On the first subinterval [t0, t1, we approximate the integral with the quadratic interpolation polynomials defined on the nodes t0, t1, t2 and in the other subinterval (Formula presented), we approximate the integral with the quadratic interpolation polynomials defined on the nodes (Formula presented). A high-order numerical method is obtained. Then we apply this numerical method with the non-uniform meshes with the step size (Formula presented) where (Formula presented). Numerical results show that this method with the non-uniform meshes has the higher convergence order than the standard numerical methods obtained by using the rectangle and the trapzoid rules with the same non-uniform meshes.