  In this paper we consider the numerical solution of fractional terminal value
problems (FDE-TVPs). In particular, the proposed procedure uses a Newton-type
iteration which is particularly efficient when coupled with a
recently-introduced step-by-step procedure for solving fractional initial value
problems (FDE-IVPs), able to produce spectrally accurate solutions of FDE
problems. Some numerical tests are reported to make evidence of its
effectiveness.
