Impulse Response Function Matlab Code