function results=GARCHforecast(data,error)
[parameters, likelihood, stderrors, robustSE, ht, scores] = fattailed_garch(data , 1 , 1 , error);
results=parameters(1)+parameters(2)*data(end)*data(end)+parameters(3)*ht(end);