Efficient spectral collocation method for fractional differential equation with Caputo-Hadamard derivative

Hadamard type fractional calculus involves logarithmic function of an arbitrary exponent as its convolutional kernel, which causes challenge in numerical treatment. In this paper we present a spectral collocation method with mapped Jacobi log orthogonal functions (MJLOFs) as basis functions and obtain an efficient algorithm to solve Hadamard type fractional differential equations. We develop basic approximation theory for the MJLOFs and derive a recurrence relation to evaluate the collocation differentiation matrix for implementing the spectral collocation algorithm. We demonstrate the effectiveness of the new method for the nonlinear initial and boundary problems, i.e, the fractional Helmholtz equation, and the fractional Burgers equation.