Matlab code to find impulse response of a system