irf_mc1 — Monte-Carlo Impulse Response function for VAR
[irf_low,irf_upp]=irf_mc1(res,mres,S,N,p,niter,siz)
results = results tlist returned by VAR
mres = decomposition method
N = # of endogenous variables in the VAR
S = # of periods
p = # of lags in the VAR
niter = # of draws in the simulation
siz = size of the confidence band
irf_low = ((S+1) x T) lower range of impulse response confidence band
irf_upp = ((S+1) x T) upper range of impulse response confidence band