garch_Model = arima('ARLags',[1,2],'MALags',[1,2]);
garch_Model.Variance = garch('GARCHLags',1,'ARCHLags',1);
[garch_Es*****l,garch_EstParamCov,garch_logL,garch_info] = estimate(garch_Model,return1);
summarize(garch_Es*****l)