How to fit GARCH(1,1) simulation to ARCH(1)

Question

I was asked to simulate 2000 data points for a standard GARCH(1,1) process. I did this using the code:

spec = garchSpec(model = list(omega = 1e-6, alpha = 0.1