GE.f90 produces the IRF for the general equilibrium version of the model.  irfzshock.m produces Figure 14.