function perexpr;
num=[0 0 1];
den=[1 3 3 1];
w=tf(num,den);
step(w);
