Why don't the irfs match the constraint condition

test.mod (3.9 KB)
we can run this mod file successfully, but the irf of hs displays a sawtooth pattern, and we check the equations, find that the irfs do not match equation (20). For example, exp(hs steady state value)* irf (hs)+exp(hb steady state value)* irf (hb) should be zero, but the irfs do not show this result.
I wonder what accounts for this and look forward to your perspective. Thank you.

You did not provide the SS1.