Again, that is not a well-defined "function". The effect of a time t shock depends upon the GARCH covariance matrix at t, producing a different "M" term for each subsequent entry.
Elder uses the unconditional variance for the size of the shock. There for the IRF as he defines should not be dependent on the point in time. The simulation then depends only the coefficient estimates.