Abstract
Multi-order fractional differential equations are motivated by their flexibility to describe complex multi-rate physical processes. This paper is concerned with the convergence behavior of a spectral collocation method when used to approximate solutions of nonlinear multi-order fractional initial value problems. The collocation scheme and its convergence analysis are developed based on the novel spectral collocation method which has been recently presented by Wang et al. (J Sci Comput 76(1):166-188, 2018) for single fractional-order boundary value problems. The proposed method is an indirect approach since we act on the equivalent Volterra integral equation of the second kind. More precisely, the spectral rate of convergence for the proposed method is established in the L2- and L infinity-norms. The method enjoys high accuracy for problems with smooth solutions. Exponentially rapid convergence is observed with a small number of degree of freedoms and for all samples of fractional orders. Numerical examples are presented to support the theoretical finding.