alpha}}+h^2{\right )}.$$ We have that $\mathcal{E}^m=\mathcal{K}^{-1}\mathcal{R}^m$. Hence $${\left \|} \mathcal{E}^m{\right \|}\leq 
{\left \|} \mathcal{K}^{-1}{\right \|}{\left \|} \mathcal{R}^m{\right \|}\leq {\dfrac{1}{{\sigma}_0^{({\alpha})}}}{\left \|} \mathcal{R}^m{\right \|}\leq 2A m^{\alpha}{\tau}^{\alpha}{\left (}{\tau}^{2-{\alpha}}+h^2{\right )}.$$

From we obtain the estimate for the error on the grid $\mathcal{G}$ $${\left \|} \mathcal{E}^m{\right \|}\leq  2A M^{\alpha}{\tau}^{\alpha}{\left (}{\tau}^{2-{\alpha}}+h^2{\right )}\leq  2A T^{\alpha}{\left (}{\tau}^{2-{\alpha}}+h^2{\right )},$$ for all $m=1,2,\cdots,M$.

[99]{} M. Abramowitz, I. A. Stegun: [*Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables,*]{} Dover, New York, 1964. J. H. Ahlberg, E. N. Nilson: [*Convergence Properties of the Spline Fit,*]{} Journal of the Society for Industrial and Applied Mathematics, [**11(1)**]{}, (1963), 95–104. D. Baleanu, K. Diethelm, E. Scalas and J. J. Trujillo: [*Fractional Calculus: Models and Numerical Methods,*]{} World Scientific, (2012). K. Diethelm: [*The Analysis of Fractional Differential Equations: An Application-Oriented Exposition Using Differential Operators of Caputo Type,*]{} Springer; 2010. Y. Dimitrov: [*A second order approximation for the Caputo fractional derivative,*]{} [**arXiv:1502.00719**]{} (2015). Y. Dimitrov: [*Higher-Order Numerical Solutions of the Fractional Relaxation-Oscillation Equation using Fractional Integration,*]{} [**arXiv:1603.08733**]{} (2016). Y. Dimitrov: [*Approximations for the Caputo derivative (I)*]{}, In preparation (2016). H. Ding, C. Li: [*High
alpha}}+h^2{\right) } .$$ We have that $ \mathcal{E}^m=\mathcal{K}^{-1}\mathcal{R}^m$. Hence $ $ { \left \| } \mathcal{E}^m{\right \|}\leq 
 { \left \| } \mathcal{K}^{-1}{\right \|}{\left \| } \mathcal{R}^m{\right \|}\leq { \dfrac{1}{{\sigma}_0^{({\alpha})}}}{\left \| } \mathcal{R}^m{\right \|}\leq 2A m^{\alpha}{\tau}^{\alpha}{\left (} { \tau}^{2-{\alpha}}+h^2{\right) } .$$ 

 From we obtain the estimate for the mistake on the power system $ \mathcal{G}$ $ $ { \left \| } \mathcal{E}^m{\right \|}\leq   2A M^{\alpha}{\tau}^{\alpha}{\left (} { \tau}^{2-{\alpha}}+h^2{\right) } \leq   2A T^{\alpha}{\left (} { \tau}^{2-{\alpha}}+h^2{\right) }, $ $ for all $ m=1,2,\cdots, M$. 

 [ 99 ] { } M. Abramowitz, I. A. Stegun: [ * Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables, * ] { } Dover, New York, 1964. J. H. Ahlberg, E. N. Nilson: [ * Convergence Properties of the Spline Fit, * ] { } Journal of the Society for Industrial and Applied Mathematics, [ * * 11(1) * * ] { }, (1963), 95–104. D. Baleanu, K. Diethelm, E. Scalas and J. J. Trujillo: [ * Fractional Calculus: Models and Numerical Methods, * ] { } World Scientific, (2012). K. Diethelm: [ * The Analysis of Fractional Differential Equations: An Application - Oriented Exposition use Differential Operators of Caputo Type, * ] { } Springer; 2010. Y. Dimitrov: [ * A second orderliness approximation for the Caputo fractional derivative, * ] { } [ * * arXiv:1502.00719 * * ] { } (2015). Y. Dimitrov: [ * high - Order Numerical Solutions of the Fractional Relaxation - Oscillation Equation use Fractional Integration, * ] { } [ * * arXiv:1603.08733 * * ] { } (2016). Y. Dimitrov: [ * Approximations for the Caputo derivative (iodine) * ] { }, In preparation (2016). H. Ding, C. Li: [ * High
alpja}}+h^2{\right )}.$$ We have that $\mxthcal{E}^m=\mathcal{K}^{-1}\mathcan{R}^m$. Hehce $${\left \|} \mathcal{E}^m{\right \|}\leq 
{\left \|} \methcql{K}^{-1}{\ritht \|}{\left \|} \mathcal{R}^m{\rigft \|}\leq {\dfgac{1}{{\sigma}_0^{({\aopha})}}}{\owft \|} \mathcem{R}^m{\righb \|}\leq 2W m^{\an'ha}{\tau}^{\alpha}{\left (}{\tau}^{2-{\alpha}}+h^2{\sight )}.$$

From we mbgann the estimate for the error on the grid $\mstjcal{G}$ $${\left \|} \majhcal{T}^m{\ryght \|}\leq  2A M^{\alpha}{\tau}^{\alpha}{\left (}{\tau}^{2-{\alphz}}+h^2{\right )}\leq  2A T^{\alphs}{\left (}{\tau}^{2-{\alpha}}+h^2{\right )},$$ for wll $l=1,2,\cdots,M$.

[99]{} M. Abramowihz, I. A. Stegoh: [*Hwbdbook of Maghematical Functions wjth Formulas, Graphs, and Mathemagical Tables,*]{} Docee, Nff York, 1964. J. H. Ahlbvrg, E. N. Nilsok: [*Convesgence Lroperties of bhe S'linw Fit,*]{} Journal of the Vociety for Industtial and A[pmied Mathematics, [**11(1)**]{}, (1963), 95–104. D. Baneang, K. Auetfelj, X. Sdalas wnd J. J. Trujimlo: [*Fractiobal Calculus: Models age Numerical Msthods,*]{} Wjrld Scientific, (2012). K. Diethelm: [*The Analysiv or Fractional Differentiql Equations: An Applifation-Oriqnted Exposition Using Differential Operators of Waputk Tyke,*]{} Sprivteg; 2010. Y. Dimitrov: [*A second order approximation fow tne Caputo fractijnal derivayige,*]{} [**wrXiv:1502.00719**]{} (2015). Y. Dimijrov: [*Hiyged-Order Numerical Sllutionf of rhe Fractyonak Relaxation-Oscillation Equqtion using Yraxtional Integratiou,*]{} [**arXiv:1603.08733**]{} (2016). Y. Dnmitrof: [*Apptoximations for the Capbto dedivative (I)*]{}, Ln preparzgion (2016). H. Ding, C. Lk: [*Hpgh
alpha}}+h^2{\right )}.$$ We have that $\mathcal{E}^m=\mathcal{K}^{-1}\mathcal{R}^m$. Hence \mathcal{E}^m{\right {\left \|} \|}{\left \|} \mathcal{R}^m{\right 2A (}{\tau}^{2-{\alpha}}+h^2{\right )}.$$ From obtain the estimate the error on the grid $\mathcal{G}$ \|} \mathcal{E}^m{\right \|}\leq 2A M^{\alpha}{\tau}^{\alpha}{\left (}{\tau}^{2-{\alpha}}+h^2{\right )}\leq 2A T^{\alpha}{\left (}{\tau}^{2-{\alpha}}+h^2{\right )},$$ for all [99]{} M. Abramowitz, I. A. Stegun: [*Handbook of Mathematical Functions with Formulas, Graphs, Mathematical Dover, York, J. H. Ahlberg, E. N. Nilson: [*Convergence Properties of the Spline Fit,*]{} Journal of the Society Industrial and Applied Mathematics, [**11(1)**]{}, (1963), 95–104. D. K. Diethelm, E. Scalas J. J. Trujillo: [*Fractional Calculus: and Methods,*]{} World (2012). Diethelm: Analysis of Fractional Equations: An Application-Oriented Exposition Using Differential Operators of Caputo Type,*]{} Springer; 2010. Y. Dimitrov: [*A second order for the derivative,*]{} [**arXiv:1502.00719**]{} Y. [*Higher-Order Solutions of the Equation using Fractional Integration,*]{} [**arXiv:1603.08733**]{} (2016). for the Caputo derivative (I)*]{}, In preparation (2016). Ding, C. [*High
alpha}}+h^2{\right )}.$$ We have that $\mathCal{E}^m=\mathcAl{K}^{-1}\maThcAl{R}^M$. HEnce $${\Left \|} \Mathcal{E}^m{\right \|}\LEq 
{\leFt \|} \mathcal{K}^{-1}{\right \|}{\left \|} \matHcal{R}^M{\rIGht \|}\lEQ {\dFrac{1}{{\sIgma}_0^{({\alpHA})}}}{\lEFT \|} \maThCaL{R}^m{\RiGHt \|}\Leq 2A m^{\AlpHa}{\tau}^{\alPha}{\left (}{\tau}^{2-{\AlpHa}}+H^2{\right )}.$$

From we OBtAin the estiMatE for the error On tHe grid $\MaThcAL{G}$ $${\lefT \|} \maThcal{e}^m{\righT \|}\Leq  2A M^{\aLpha}{\tau}^{\alPhA}{\Left (}{\taU}^{2-{\Alpha}}+h^2{\rIGHt )}\Leq  2A t^{\alpha}{\left (}{\tau}^{2-{\alphA}}+H^2{\rIGht )},$$ for all $m=1,2,\cdotS,M$.

[99]{} M. AbrAmOWiTZ, i. A. STegUn: [*Handbook Of matheMAtical FUNcTIONs wITh Formulas, GraPhs, and MatheMAtiCal TabLeS,*]{} DoVEr, New YOrk, 1964. J. H. ahLBerG, E. N. Nilson: [*CoNverGence PropErties OF the SplINe Fit,*]{} JoUrnal oF thE SoCietY FoR INduStRIal ANd appLIed mathematIcS, [**11(1)**]{}, (1963), 95–104. D. baleaNu, K. DIETHElm, E. scaLas aNd J. J. TRujillo: [*FractiOnaL CalCUluS: ModeLs and numeRiCal MeThods,*]{} WOrld SCiEntific, (2012). K. DiethelM: [*The analysis oF FrAcTioNaL DiffERentiaL EqUatIons: An APplicatIOn-ORiENTEd exposition Using DifFeRENtIal OperaTors of cApUtO type,*]{} SpriNgEr; 2010. Y. dimiTROv: [*A seCond ORdEr approxImatioN FoR tHe CaputO fRactioNaL deRivAtive,*]{} [**ARXiv:1502.00719**]{} (2015). y. DimitRov: [*HigheR-OrdeR numerical SolutIOns of the FractIOnAL reLAxatIon-oscillation equaTIon uSing fRaCtiONal InTegraTiON,*]{} [**aRxiv:1603.08733**]{} (2016). Y. Dimitrov: [*ApproxiMaTions fOr the caputo derivatIve (I)*]{}, In prepARATion (2016). H. DinG, C. Li: [*hIgH
alpha}}+h^2{\right )}.$$ W e have tha t $\m ath cal {E }^m= \mat hcal{K}^{-1}\m a thca l{R}^m$. Hence $${\lef t \|}  \ m athc a l{ E}^m{ \right\ |} \ l eq
{ \l eft  \ | }\math cal {K}^{-1 }{\right \ |}{ \l eft \|} \mat h ca l{R}^m{\ri ght  \|}\leq {\d fra c{1}{{ \s igm a }_0^{ ({\ alpha })}}}{ \ left \ |} \mathc al { R}^m{\ r ight \| } \ le q 2A  m^{\alpha}{\tau} ^ {\ a lpha}{\left (} {\tau} ^{ 2 -{ \ a lph a}} +h^2{\righ t)}.$$ 
From w e  o b t a int he estimate f or the erro r  on  the g ri d $ \ mathca l{G}$  $ $ {\l eft \|} \ma thca l{E}^m{\r ight \ | }\leq 2 A M^{\a lpha}{ \ta u}^ {\al p ha }{ \le ft (}{ \ ta u}^ { 2-{ \alpha}} +h ^2 {\rig ht ) } \ l e q  2 A T ^{\a lpha} {\left (}{\ta u}^ {2-{ \ alp ha}}+ h^2{\ righ t)},$$  for a ll $m =1 ,2,\cdots,M$.
[99] {} M. Abr amo wi tz,  I . A.S tegun:  [* Han dbook o f Mathe m ati ca l F un ctions with Formul as , Gr aphs, an d Math e ma ti c al Table s, *]{ } Do v e r, Ne w Yo r k,  1964. J . H. A h lb er g, E. N .Nilson :[*C onv ergen c e Pr operti es of th e Spl i ne Fit,*]{} Jo u rnal of the S o ci e t yf or I ndu strial andAppl i ed M athe m at ics ,  [**1 1(1)* *] { }, (1963), 95–104. D.Ba leanu,  K. D iethelm, E. S calas andJ . J. Truji llo: [* F ractional Calc ulus:  Models an d  Numeric al Me thods,*] {} WorldS c ientific , ( 201 2).  K. D ie thelm: [*TheA n alys is  of Fra cti onal Di ffe ren tia l E qu ations: A n Applic at io n- Or ien ted E x position  U sin gDif feren t ial Op erato rs o fCa p uto  Type,* ] {} S prin ge r;  201 0.Y.  Dimi trov :  [* A secon d order a ppr o xima ti on  for th e Caputo frac ti onal deriv at ive ,*]{}[ * *arXiv:1 502.00719**]{} (2015).Y . Dimit rov : [*H ighe r-Order N ume ricalSol u tionsof the  Frac ti ona l Relax a t io n-O sc illation E q u ati on us in g Fr actiona l Integration,*]{} [** arXiv:1603.08 733 **]{ } (2 016 ) .Y . D im i tro v :  [*Approximatio ns for the  C a pu to derivat i ve(I )*]{},In prep arati o n (2016 ). H. Din g, C. Li:  [ *Hig h 
alpha}}+h^2{\right )}.$$_We have_that $\mathcal{E}^m=\mathcal{K}^{-1}\mathcal{R}^m$. Hence $${\left_\|} \mathcal{E}^m{\right_\|}\leq_
{\left \|}_\mathcal{K}^{-1}{\right_\|}{\left \|} \mathcal{R}^m{\right_\|}\leq {\dfrac{1}{{\sigma}_0^{({\alpha})}}}{\left \|}_\mathcal{R}^m{\right \|}\leq 2A m^{\alpha}{\tau}^{\alpha}{\left_(}{\tau}^{2-{\alpha}}+h^2{\right )}.$$

From we_obtain_the estimate for the error on the grid $\mathcal{G}$ $${\left \|} \mathcal{E}^m{\right \|}\leq _2A_M^{\alpha}{\tau}^{\alpha}{\left (}{\tau}^{2-{\alpha}}+h^2{\right_)}\leq__2A T^{\alpha}{\left (}{\tau}^{2-{\alpha}}+h^2{\right )},$$ for_all $m=1,2,\cdots,M$.

[99]{} M. Abramowitz, I._A. Stegun:_[*Handbook of Mathematical Functions with Formulas, Graphs, and_Mathematical_Tables,*]{} Dover, New_York, 1964. J. H. Ahlberg, E. N. Nilson: [*Convergence_Properties of the Spline Fit,*]{} Journal_of the Society_for_Industrial_and Applied Mathematics, [**11(1)**]{},_(1963), 95–104. D. Baleanu, K. Diethelm,_E. Scalas and J. J. Trujillo:_[*Fractional Calculus: Models and Numerical Methods,*]{} World_Scientific, (2012). K. Diethelm: [*The Analysis_of Fractional Differential Equations: An_Application-Oriented Exposition_Using Differential Operators of Caputo_Type,*]{} Springer; 2010._Y. Dimitrov:_[*A second order_approximation for the Caputo fractional derivative,*]{}_[**arXiv:1502.00719**]{} (2015). Y._Dimitrov: [*Higher-Order Numerical Solutions of the_Fractional_Relaxation-Oscillation Equation using_Fractional_Integration,*]{}_[**arXiv:1603.08733**]{} (2016)._Y. Dimitrov: [*Approximations_for_the Caputo_derivative_(I)*]{}, In preparation (2016). H. Ding,_C._Li: [*High