[m,n]=size(Data) %[1100,6]
X= Data(1:1000,1:5);
y= Data(1:1000,6);
spec = garchset('VarianceModel','GARCH','P',1,'Q',1);
spec = garchset(spec,'Display','off','Distribution','T');
[coeffX,errorsX,LLFX,innovationsX,sigmasX] = garchfit(spec,y,X);
garchdisp(coeffX,errorsX)
XF=p(1001:1100,1:5);
NumPeriods = 100;
[sigmaForecast,meanForecast] = garchpred(coeffX,y,NumPeriods,X,XF);