An implicit Keller Box numerical scheme for the solution of fractional subdiffusion equations

摘要

In this work, we present a new implicit numerical scheme for fractional subdiffusion equations. In this approach we use the Keller Box method [1] to spatially discretise the fractional subdiffusion equation and we use a modified L1 scheme (ML1), similar to the L1 scheme originally developed by Oldham and Spanier [2], to approximate the fractional derivative. The stability of the proposed method was investigated by using Von-Neumann stability analysis. We have proved the method is unconditionally stable when 0<λq