Numerical treatment of fractional differential models