variance.model = list(model = "fGARCH", garchOrder = c(1, 1),submodel = "GARCH", external.regressors = NULL, variance.targeting = FALSE);
mean.model = list(armaOrder = c(1, 0), include.mean = TRUE, archm = FALSE, archpow = 1, arfima = FALSE, external.regressors = NULL, archex = FALSE);
myspec=ugarchspec(variance.model = variance.model, mean.model = mean.model, distribution.model = "norm");
myfit1=ugarchfit(myspec,data=x1,solver="solnp");
E1=as.numeric(residuals(myfit1,standardize=TRUE));