Calculation of the impulse response

Program for calculation of impulse response of strictly proper SISO systems:

*/num = numerator polynomial coefficients of transfer function
*/den = denominator polynomial coefficients of transfer function
(Coefficients of 'num' and 'den' are specified as a row vector, in
decreasing powers of 's')
*/tO = time at which unit impulse input is applied
*/dt = time-step (should be smaller than 1/ (largest natural freq.))
*/tf = final time for impulse response calculation
*/y = impulse response; t= vector of time points