sim1.m was adapted to deal with models with leads and lags of more than 1 period. However, the reshaping of the results was not adjusted