“…For the numerical solution of such intial value problems, many methods have been proposed; cf., e.g., [2,11,12,17,21]. A common feature of a large part of these methods is that, in order to correctly handle the non-locality of the fractional differential operators, they require a relatively large amount of time and/or computer memory [14].…”