Solve the IRFs for a DSGE model with a Markov Switching process

The short answer is: