No BSD License
function [h,n] = impz(b,a); n = 0:4/(1-max(abs(roots(a)))); delta = zeros(size(n)); delta(1) = 1; % unit impulse h = filter(b,a,delta);
Contact us at files@mathworks.com