syms n z a=[5 0 5]; b=[1 -2.5 1] [r,p,k]=residuez(a,b); zplane(a,b); w=(z^2+z)/(z^2-.5*z-.5); %y=ztrnas iztrans(w)