%% stepresp.m Plots the step response of a system 
% numerator and denominator of transfer function
num=[1];
den=[1 3 9];
sys=tf(num,den);
% time vector
tfinal=10;
step(sys,tfinal);
